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

Minimum AIC  model: GTRIG
Minimum AICc model: GTRIG
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
GTRIG	45	-4767.3369	9628.5001	0.8658	0.8658
GTRG	44	-4770.3158	9632.2881	0.1303	0.9961
HKYIG	41	-4777.1851	9639.5415	0.0035	0.9995
HKYG	40	-4780.2707	9643.5588	0.0005	1.0000
K2PIG	38	-4794.7372	9668.1962	0.0000	1.0000
SYMIG	42	-4790.9172	9669.1634	0.0000	1.0000
K2PG	37	-4797.6653	9671.9103	0.0000	1.0000
SYMG	41	-4793.5607	9672.2926	0.0000	1.0000
GTRI	44	-4792.9975	9677.6515	0.0000	1.0000
HKYI	40	-4800.6971	9684.4117	0.0000	1.0000
K2PI	37	-4816.0355	9708.6508	0.0000	1.0000
SYMI	41	-4812.0950	9709.3613	0.0000	1.0000
F81IG	40	-4882.0173	9847.0521	0.0000	1.0000
F81G	39	-4884.5013	9849.8702	0.0000	1.0000
JC69IG	37	-4893.1528	9862.8855	0.0000	1.0000
JC69G	36	-4895.5611	9865.5639	0.0000	1.0000
F81I	39	-4901.9522	9884.7720	0.0000	1.0000
JC69I	36	-4912.2345	9898.9108	0.0000	1.0000
GTR	43	-4954.7941	9999.0789	0.0000	1.0000
HKY	39	-4967.4425	10015.7527	0.0000	1.0000
SYM	40	-4975.7922	10034.6019	0.0000	1.0000
K2P	36	-4983.8662	10042.1742	0.0000	1.0000
F81	38	-5066.3694	10211.4605	0.0000	1.0000
JC69	35	-5077.6910	10227.6897	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	45	-4767.3369	9624.6739	0.8761	0.8761
GTRG	44	-4770.3158	9628.6316	0.1211	0.9972
HKYIG	41	-4777.1851	9636.3702	0.0025	0.9997
HKYG	40	-4780.2707	9640.5413	0.0003	1.0000
K2PIG	38	-4794.7372	9665.4744	0.0000	1.0000
SYMIG	42	-4790.9172	9665.8343	0.0000	1.0000
SYMG	41	-4793.5607	9669.1213	0.0000	1.0000
K2PG	37	-4797.6653	9669.3305	0.0000	1.0000
GTRI	44	-4792.9975	9673.9950	0.0000	1.0000
HKYI	40	-4800.6971	9681.3943	0.0000	1.0000
K2PI	37	-4816.0355	9706.0710	0.0000	1.0000
SYMI	41	-4812.0950	9706.1900	0.0000	1.0000
F81IG	40	-4882.0173	9844.0346	0.0000	1.0000
F81G	39	-4884.5013	9847.0026	0.0000	1.0000
JC69IG	37	-4893.1528	9860.3056	0.0000	1.0000
JC69G	36	-4895.5611	9863.1221	0.0000	1.0000
F81I	39	-4901.9522	9881.9044	0.0000	1.0000
JC69I	36	-4912.2345	9896.4690	0.0000	1.0000
GTR	43	-4954.7941	9995.5882	0.0000	1.0000
HKY	39	-4967.4425	10012.8851	0.0000	1.0000
SYM	40	-4975.7922	10031.5845	0.0000	1.0000
K2P	36	-4983.8662	10039.7324	0.0000	1.0000
F81	38	-5066.3694	10208.7388	0.0000	1.0000
JC69	35	-5077.6910	10225.3820	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-4780.2707	9841.6694	0.5955	0.5955
HKYIG	41	-4777.1851	9842.5265	0.3879	0.9834
GTRG	44	-4770.3158	9849.8724	0.0099	0.9933
GTRIG	45	-4767.3369	9850.9429	0.0058	0.9990
K2PG	37	-4797.6653	9855.3740	0.0006	0.9996
K2PIG	38	-4794.7372	9856.5461	0.0004	1.0000
SYMG	41	-4793.5607	9875.2776	0.0000	1.0000
SYMIG	42	-4790.9172	9877.0188	0.0000	1.0000
HKYI	40	-4800.6971	9882.5223	0.0000	1.0000
K2PI	37	-4816.0355	9892.1145	0.0000	1.0000
GTRI	44	-4792.9975	9895.2359	0.0000	1.0000
SYMI	41	-4812.0950	9912.3463	0.0000	1.0000
F81G	39	-4884.5013	10043.1025	0.0000	1.0000
JC69G	36	-4895.5611	10044.1374	0.0000	1.0000
F81IG	40	-4882.0173	10045.1627	0.0000	1.0000
JC69IG	37	-4893.1528	10046.3491	0.0000	1.0000
JC69I	36	-4912.2345	10077.4843	0.0000	1.0000
F81I	39	-4901.9522	10078.0042	0.0000	1.0000
HKY	39	-4967.4425	10208.9849	0.0000	1.0000
GTR	43	-4954.7941	10211.8008	0.0000	1.0000
K2P	36	-4983.8662	10220.7477	0.0000	1.0000
SYM	40	-4975.7922	10232.7125	0.0000	1.0000
F81	38	-5066.3694	10399.8104	0.0000	1.0000
JC69	35	-5077.6910	10401.3690	0.0000	1.0000



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

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

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


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

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

