
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:57 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/chr2_21308.phylip" (19 taxa, 425 characters)

Minimum AIC  model: GTRG
Minimum AICc model: K2PG
Minimum BIC  model: K2PG


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
K2PG	37	-1200.3168	2481.8998	0.2134	0.2134
K2PI	37	-1200.7878	2482.8417	0.1332	0.3466
K2PIG	38	-1199.6880	2483.0548	0.1198	0.4663
GTRG	44	-1192.4120	2483.2450	0.1089	0.5752
SYMG	41	-1196.3655	2483.7232	0.0857	0.6610
GTRIG	45	-1191.4834	2483.8902	0.0789	0.7398
GTRI	44	-1192.7992	2484.0195	0.0739	0.8138
SYMI	41	-1196.6793	2484.3507	0.0626	0.8764
SYMIG	42	-1195.7011	2484.8576	0.0486	0.9250
HKYG	40	-1198.4893	2485.5203	0.0349	0.9599
HKYI	40	-1198.9780	2486.4976	0.0214	0.9813
HKYIG	41	-1197.8910	2486.7741	0.0186	1.0000
JC69G	36	-1213.6075	2506.0809	0.0000	1.0000
JC69I	36	-1215.5689	2510.0037	0.0000	1.0000
JC69IG	37	-1214.5036	2510.2734	0.0000	1.0000
F81G	39	-1212.2524	2510.6087	0.0000	1.0000
F81I	39	-1214.2021	2514.5082	0.0000	1.0000
F81IG	40	-1213.1544	2514.8504	0.0000	1.0000
K2P	36	-1229.7625	2538.3910	0.0000	1.0000
SYM	40	-1226.2194	2540.9805	0.0000	1.0000
GTR	43	-1223.3256	2542.5829	0.0000	1.0000
HKY	39	-1228.5284	2543.1606	0.0000	1.0000
JC69	35	-1242.8613	2562.2008	0.0000	1.0000
F81	38	-1241.5639	2566.8066	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-1192.4120	2472.8239	0.2094	0.2094
GTRIG	45	-1191.4834	2472.9667	0.1950	0.4044
GTRI	44	-1192.7992	2473.5985	0.1422	0.5466
K2PG	37	-1200.3168	2474.6336	0.0847	0.6314
SYMG	41	-1196.3655	2474.7311	0.0807	0.7121
SYMI	41	-1196.6793	2475.3586	0.0590	0.7710
K2PIG	38	-1199.6880	2475.3760	0.0585	0.8295
SYMIG	42	-1195.7011	2475.4021	0.0577	0.8872
K2PI	37	-1200.7878	2475.5756	0.0529	0.9401
HKYG	40	-1198.4893	2476.9787	0.0262	0.9664
HKYIG	41	-1197.8910	2477.7819	0.0176	0.9839
HKYI	40	-1198.9780	2477.9560	0.0161	1.0000
JC69G	36	-1213.6075	2499.2150	0.0000	1.0000
F81G	39	-1212.2524	2502.5048	0.0000	1.0000
JC69IG	37	-1214.5036	2503.0073	0.0000	1.0000
JC69I	36	-1215.5689	2503.1378	0.0000	1.0000
F81IG	40	-1213.1544	2506.3087	0.0000	1.0000
F81I	39	-1214.2021	2506.4043	0.0000	1.0000
K2P	36	-1229.7625	2531.5251	0.0000	1.0000
SYM	40	-1226.2194	2532.4389	0.0000	1.0000
GTR	43	-1223.3256	2532.6511	0.0000	1.0000
HKY	39	-1228.5284	2535.0567	0.0000	1.0000
JC69	35	-1242.8613	2555.7226	0.0000	1.0000
F81	38	-1241.5639	2559.1279	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	37	-1200.3168	2624.5609	0.5824	0.5824
K2PI	37	-1200.7878	2625.5029	0.3636	0.9460
K2PIG	38	-1199.6880	2629.3554	0.0530	0.9990
HKYG	40	-1198.4893	2639.0622	0.0004	0.9994
HKYI	40	-1198.9780	2640.0395	0.0003	0.9996
SYMG	41	-1196.3655	2640.8667	0.0002	0.9998
SYMI	41	-1196.6793	2641.4942	0.0001	0.9999
HKYIG	41	-1197.8910	2643.9176	0.0000	1.0000
JC69G	36	-1213.6075	2645.0902	0.0000	1.0000
SYMIG	42	-1195.7011	2645.5899	0.0000	1.0000
JC69I	36	-1215.5689	2649.0130	0.0000	1.0000
GTRG	44	-1192.4120	2651.1158	0.0000	1.0000
GTRI	44	-1192.7992	2651.8904	0.0000	1.0000
JC69IG	37	-1214.5036	2652.9346	0.0000	1.0000
GTRIG	45	-1191.4834	2655.3107	0.0000	1.0000
F81G	39	-1212.2524	2660.5363	0.0000	1.0000
F81I	39	-1214.2021	2664.4358	0.0000	1.0000
F81IG	40	-1213.1544	2668.3923	0.0000	1.0000
K2P	36	-1229.7625	2677.4003	0.0000	1.0000
HKY	39	-1228.5284	2693.0882	0.0000	1.0000
SYM	40	-1226.2194	2694.5224	0.0000	1.0000
JC69	35	-1242.8613	2697.5458	0.0000	1.0000
GTR	43	-1223.3256	2706.8910	0.0000	1.0000
F81	38	-1241.5639	2713.1072	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_21308.phylip = 1 - 425;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_21308.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_21308.phylip = 1 - 425;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_21308.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_21308.phylip = 1 - 425;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_21308.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

