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

Minimum AIC  model: HKYIG
Minimum AICc model: HKYIG
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
HKYIG	41	-1678.7227	3445.4039	0.3177	0.3177
HKYI	40	-1679.9737	3445.6123	0.2863	0.6040
GTRI	44	-1676.0169	3446.9208	0.1488	0.7528
HKYG	40	-1680.8615	3447.3879	0.1178	0.8706
GTRG	44	-1676.6140	3448.1150	0.0819	0.9525
GTRIG	45	-1675.9964	3449.2053	0.0475	1.0000
K2PI	37	-1706.1391	3491.1098	0.0000	1.0000
K2PIG	38	-1705.3706	3491.8426	0.0000	1.0000
SYMI	41	-1701.9523	3491.8630	0.0000	1.0000
SYMIG	42	-1701.4945	3493.2489	0.0000	1.0000
K2PG	37	-1709.0347	3496.9009	0.0000	1.0000
SYMG	41	-1704.9908	3497.9401	0.0000	1.0000
F81I	39	-1727.3449	3538.0691	0.0000	1.0000
F81G	39	-1727.8214	3539.0221	0.0000	1.0000
F81IG	40	-1726.9368	3539.5386	0.0000	1.0000
JC69I	36	-1750.8956	3578.3607	0.0000	1.0000
JC69IG	37	-1750.1932	3579.2180	0.0000	1.0000
JC69G	36	-1753.4357	3583.4409	0.0000	1.0000
GTR	43	-1750.6001	3593.7696	0.0000	1.0000
HKY	39	-1756.6178	3596.6149	0.0000	1.0000
SYM	40	-1781.8217	3649.3083	0.0000	1.0000
K2P	36	-1787.3807	3651.3308	0.0000	1.0000
F81	38	-1801.6302	3684.3620	0.0000	1.0000
JC69	35	-1828.7384	3731.7919	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYIG	41	-1678.7227	3439.4454	0.2830	0.2830
HKYI	40	-1679.9737	3439.9473	0.2202	0.5032
GTRI	44	-1676.0169	3440.0339	0.2109	0.7141
GTRG	44	-1676.6140	3441.2280	0.1161	0.8302
HKYG	40	-1680.8615	3441.7230	0.0906	0.9208
GTRIG	45	-1675.9964	3441.9928	0.0792	1.0000
SYMI	41	-1701.9523	3485.9046	0.0000	1.0000
K2PI	37	-1706.1391	3486.2781	0.0000	1.0000
K2PIG	38	-1705.3706	3486.7411	0.0000	1.0000
SYMIG	42	-1701.4945	3486.9889	0.0000	1.0000
SYMG	41	-1704.9908	3491.9816	0.0000	1.0000
K2PG	37	-1709.0347	3492.0693	0.0000	1.0000
F81I	39	-1727.3449	3532.6897	0.0000	1.0000
F81G	39	-1727.8214	3533.6428	0.0000	1.0000
F81IG	40	-1726.9368	3533.8736	0.0000	1.0000
JC69I	36	-1750.8956	3573.7912	0.0000	1.0000
JC69IG	37	-1750.1932	3574.3864	0.0000	1.0000
JC69G	36	-1753.4357	3578.8715	0.0000	1.0000
GTR	43	-1750.6001	3587.2002	0.0000	1.0000
HKY	39	-1756.6178	3591.2356	0.0000	1.0000
SYM	40	-1781.8217	3643.6434	0.0000	1.0000
K2P	36	-1787.3807	3646.7613	0.0000	1.0000
F81	38	-1801.6302	3679.2605	0.0000	1.0000
JC69	35	-1828.7384	3727.4768	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	40	-1679.9737	3617.1361	0.6443	0.6443
HKYG	40	-1680.8615	3618.9118	0.2652	0.9095
HKYIG	41	-1678.7227	3621.0639	0.0904	0.9999
GTRI	44	-1676.0169	3634.9415	0.0001	0.9999
GTRG	44	-1676.6140	3636.1357	0.0000	1.0000
GTRIG	45	-1675.9964	3641.3302	0.0000	1.0000
K2PI	37	-1706.1391	3650.1778	0.0000	1.0000
K2PIG	38	-1705.3706	3655.0704	0.0000	1.0000
K2PG	37	-1709.0347	3655.9689	0.0000	1.0000
SYMI	41	-1701.9523	3667.5231	0.0000	1.0000
SYMIG	42	-1701.4945	3673.0371	0.0000	1.0000
SYMG	41	-1704.9908	3673.6001	0.0000	1.0000
F81I	39	-1727.3449	3705.4488	0.0000	1.0000
F81G	39	-1727.8214	3706.4018	0.0000	1.0000
F81IG	40	-1726.9368	3711.0624	0.0000	1.0000
JC69I	36	-1750.8956	3733.2611	0.0000	1.0000
JC69IG	37	-1750.1932	3738.2860	0.0000	1.0000
JC69G	36	-1753.4357	3738.3414	0.0000	1.0000
HKY	39	-1756.6178	3763.9946	0.0000	1.0000
GTR	43	-1750.6001	3777.6781	0.0000	1.0000
K2P	36	-1787.3807	3806.2312	0.0000	1.0000
SYM	40	-1781.8217	3820.8322	0.0000	1.0000
F81	38	-1801.6302	3847.5898	0.0000	1.0000
JC69	35	-1828.7384	3882.5170	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_25705.phylip = 1 - 620;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_25705.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_25705.phylip = 1 - 620;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_25705.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=invgamma;
 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) 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/chr1_25705.phylip = 1 - 620;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_25705.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

