
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:02 02/24/2011

Input data from file "/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr13_2812.phylip" (19 taxa, 496 characters)

Minimum AIC  model: HKYI
Minimum AICc model: HKYI
Minimum BIC  model: HKYI


Note: number of parameters compared to sample size
is low (Nchar/Nparams < 40). Consider the use of AICc.


(Output sorted by AICc. w: Akaike weight, cw: cumulative w.)

Model	df	lnL		AICc		wAICc	cwAICc
HKYI	40	-1668.4550	3424.1187	0.5668	0.5668
HKYG	40	-1669.1025	3425.4137	0.2966	0.8635
HKYIG	41	-1669.1023	3427.7905	0.0904	0.9539
GTRI	44	-1666.7550	3430.2904	0.0259	0.9798
GTRG	44	-1667.4857	3431.7519	0.0125	0.9922
GTRIG	45	-1666.7513	3432.7027	0.0078	1.0000
K2PI	37	-1683.0510	3446.2417	0.0000	1.0000
K2PG	37	-1683.8278	3447.7953	0.0000	1.0000
K2PIG	38	-1683.0511	3448.5879	0.0000	1.0000
SYMI	41	-1680.7591	3451.1041	0.0000	1.0000
SYMG	41	-1681.6920	3452.9699	0.0000	1.0000
SYMIG	42	-1680.7629	3453.4993	0.0000	1.0000
HKY	39	-1695.0223	3474.8868	0.0000	1.0000
GTR	43	-1693.5211	3481.4140	0.0000	1.0000
K2P	36	-1709.4240	3496.6518	0.0000	1.0000
SYM	40	-1707.4046	3502.0180	0.0000	1.0000
F81I	39	-1721.8221	3528.4862	0.0000	1.0000
F81G	39	-1722.4911	3529.8242	0.0000	1.0000
F81IG	40	-1721.8226	3530.8541	0.0000	1.0000
JC69I	36	-1733.1117	3544.0272	0.0000	1.0000
JC69G	36	-1733.8148	3545.4336	0.0000	1.0000
JC69IG	37	-1733.1133	3546.3663	0.0000	1.0000
F81	38	-1746.9465	3576.3787	0.0000	1.0000
JC69	35	-1758.3574	3592.1930	0.0000	1.0000



(Output sorted by AIC. w: Akaike weight, cw: cumulative w.)

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	40	-1668.4550	3416.9099	0.5259	0.5259
HKYG	40	-1669.1025	3418.2049	0.2752	0.8011
HKYIG	41	-1669.1023	3420.2046	0.1013	0.9024
GTRI	44	-1666.7550	3421.5099	0.0527	0.9551
GTRG	44	-1667.4857	3422.9715	0.0254	0.9805
GTRIG	45	-1666.7513	3423.5027	0.0195	1.0000
K2PI	37	-1683.0510	3440.1020	0.0000	1.0000
K2PG	37	-1683.8278	3441.6555	0.0000	1.0000
K2PIG	38	-1683.0511	3442.1022	0.0000	1.0000
SYMI	41	-1680.7591	3443.5182	0.0000	1.0000
SYMG	41	-1681.6920	3445.3840	0.0000	1.0000
SYMIG	42	-1680.7629	3445.5258	0.0000	1.0000
HKY	39	-1695.0223	3468.0447	0.0000	1.0000
GTR	43	-1693.5211	3473.0423	0.0000	1.0000
K2P	36	-1709.4240	3490.8479	0.0000	1.0000
SYM	40	-1707.4046	3494.8092	0.0000	1.0000
F81I	39	-1721.8221	3521.6441	0.0000	1.0000
F81G	39	-1722.4911	3522.9821	0.0000	1.0000
F81IG	40	-1721.8226	3523.6453	0.0000	1.0000
JC69I	36	-1733.1117	3538.2233	0.0000	1.0000
JC69G	36	-1733.8148	3539.6297	0.0000	1.0000
JC69IG	37	-1733.1133	3540.2265	0.0000	1.0000
F81	38	-1746.9465	3569.8930	0.0000	1.0000
JC69	35	-1758.3574	3586.7147	0.0000	1.0000



(Output sorted by BIC. w: Akaike weight, cw: cumulative w.)

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	40	-1668.4550	3585.1730	0.6433	0.6433
HKYG	40	-1669.1025	3586.4679	0.3367	0.9800
HKYIG	41	-1669.1023	3592.6742	0.0151	0.9951
K2PI	37	-1683.0510	3595.7453	0.0033	0.9983
K2PG	37	-1683.8278	3597.2988	0.0015	0.9998
K2PIG	38	-1683.0511	3601.9520	0.0001	1.0000
GTRI	44	-1666.7550	3606.5993	0.0000	1.0000
GTRG	44	-1667.4857	3608.0608	0.0000	1.0000
GTRIG	45	-1666.7513	3612.7986	0.0000	1.0000
SYMI	41	-1680.7591	3615.9878	0.0000	1.0000
SYMG	41	-1681.6920	3617.8536	0.0000	1.0000
SYMIG	42	-1680.7629	3622.2019	0.0000	1.0000
HKY	39	-1695.0223	3632.1012	0.0000	1.0000
K2P	36	-1709.4240	3642.2846	0.0000	1.0000
GTR	43	-1693.5211	3653.9250	0.0000	1.0000
SYM	40	-1707.4046	3663.0722	0.0000	1.0000
F81I	39	-1721.8221	3685.7006	0.0000	1.0000
F81G	39	-1722.4911	3687.0386	0.0000	1.0000
JC69I	36	-1733.1117	3689.6600	0.0000	1.0000
JC69G	36	-1733.8148	3691.0664	0.0000	1.0000
F81IG	40	-1721.8226	3691.9083	0.0000	1.0000
JC69IG	37	-1733.1133	3695.8698	0.0000	1.0000
F81	38	-1746.9465	3729.7428	0.0000	1.0000
JC69	35	-1758.3574	3733.9449	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (HKYI)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr13_2812.phylip = 1 - 496;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr13_2812.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best AICc model (HKYI)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr13_2812.phylip = 1 - 496;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr13_2812.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best BIC model (HKYI)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr13_2812.phylip = 1 - 496;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr13_2812.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


-------------------------------------------------------------
End of Output

