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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
Minimum BIC  model: HKYG


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
HKYG	40	-3378.9136	6842.4147	0.4627	0.4627
GTRG	44	-3375.0865	6843.7426	0.2382	0.7009
HKYIG	41	-3378.6275	6844.0785	0.2014	0.9022
GTRIG	45	-3374.8463	6845.5235	0.0978	1.0000
HKYI	40	-3394.1840	6872.9554	0.0000	1.0000
GTRI	44	-3390.8041	6875.1777	0.0000	1.0000
K2PG	37	-3411.3881	6900.6927	0.0000	1.0000
SYMG	41	-3407.2305	6901.2846	0.0000	1.0000
K2PIG	38	-3411.1620	6902.4579	0.0000	1.0000
SYMIG	42	-3407.0503	6903.1666	0.0000	1.0000
SYMI	41	-3420.0637	6926.9510	0.0000	1.0000
K2PI	37	-3424.5253	6926.9671	0.0000	1.0000
HKY	39	-3463.1665	7008.6906	0.0000	1.0000
GTR	43	-3460.2756	7011.8658	0.0000	1.0000
F81G	39	-3467.0544	7016.4664	0.0000	1.0000
F81IG	40	-3468.8461	7022.2797	0.0000	1.0000
F81I	39	-3477.7259	7037.8093	0.0000	1.0000
JC69G	36	-3492.3089	7060.3230	0.0000	1.0000
JC69IG	37	-3492.1767	7062.2698	0.0000	1.0000
K2P	36	-3496.6159	7068.9369	0.0000	1.0000
SYM	40	-3493.1219	7070.8312	0.0000	1.0000
JC69I	36	-3503.1562	7082.0176	0.0000	1.0000
F81	38	-3543.1490	7166.4319	0.0000	1.0000
JC69	35	-3569.9913	7213.4825	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-3378.9136	6837.8273	0.3670	0.3670
GTRG	44	-3375.0865	6838.1730	0.3088	0.6758
HKYIG	41	-3378.6275	6839.2550	0.1798	0.8556
GTRIG	45	-3374.8463	6839.6925	0.1444	1.0000
HKYI	40	-3394.1840	6868.3679	0.0000	1.0000
GTRI	44	-3390.8041	6869.6081	0.0000	1.0000
SYMG	41	-3407.2305	6896.4610	0.0000	1.0000
K2PG	37	-3411.3881	6896.7762	0.0000	1.0000
SYMIG	42	-3407.0503	6898.1007	0.0000	1.0000
K2PIG	38	-3411.1620	6898.3240	0.0000	1.0000
SYMI	41	-3420.0637	6922.1274	0.0000	1.0000
K2PI	37	-3424.5253	6923.0507	0.0000	1.0000
HKY	39	-3463.1665	7004.3330	0.0000	1.0000
GTR	43	-3460.2756	7006.5512	0.0000	1.0000
F81G	39	-3467.0544	7012.1088	0.0000	1.0000
F81IG	40	-3468.8461	7017.6923	0.0000	1.0000
F81I	39	-3477.7259	7033.4518	0.0000	1.0000
JC69G	36	-3492.3089	7056.6178	0.0000	1.0000
JC69IG	37	-3492.1767	7058.3534	0.0000	1.0000
K2P	36	-3496.6159	7065.2317	0.0000	1.0000
SYM	40	-3493.1219	7066.2438	0.0000	1.0000
JC69I	36	-3503.1562	7078.3125	0.0000	1.0000
F81	38	-3543.1490	7162.2980	0.0000	1.0000
JC69	35	-3569.9913	7209.9825	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-3378.9136	7022.9490	0.9537	0.9537
HKYIG	41	-3378.6275	7029.0047	0.0462	0.9999
GTRG	44	-3375.0865	7041.8068	0.0001	1.0000
GTRIG	45	-3374.8463	7047.9544	0.0000	1.0000
HKYI	40	-3394.1840	7053.4896	0.0000	1.0000
K2PG	37	-3411.3881	7068.0138	0.0000	1.0000
GTRI	44	-3390.8041	7073.2419	0.0000	1.0000
K2PIG	38	-3411.1620	7074.1896	0.0000	1.0000
SYMG	41	-3407.2305	7086.2107	0.0000	1.0000
SYMIG	42	-3407.0503	7092.4784	0.0000	1.0000
K2PI	37	-3424.5253	7094.2882	0.0000	1.0000
SYMI	41	-3420.0637	7111.8771	0.0000	1.0000
HKY	39	-3463.1665	7184.8266	0.0000	1.0000
F81G	39	-3467.0544	7192.6025	0.0000	1.0000
F81IG	40	-3468.8461	7202.8139	0.0000	1.0000
GTR	43	-3460.2756	7205.5569	0.0000	1.0000
F81I	39	-3477.7259	7213.9454	0.0000	1.0000
JC69G	36	-3492.3089	7223.2273	0.0000	1.0000
JC69IG	37	-3492.1767	7229.5909	0.0000	1.0000
K2P	36	-3496.6159	7231.8412	0.0000	1.0000
JC69I	36	-3503.1562	7244.9219	0.0000	1.0000
SYM	40	-3493.1219	7251.3655	0.0000	1.0000
F81	38	-3543.1490	7338.1635	0.0000	1.0000
JC69	35	-3569.9913	7371.9640	0.0000	1.0000



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

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

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


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

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


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

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


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

