
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:40 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/chr5_3273.phylip" (19 taxa, 816 characters)

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


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
GTRG	44	-3977.8827	8048.9015	0.8423	0.8423
GTRIG	45	-3978.4377	8052.2519	0.1577	1.0000
HKYG	40	-3997.1173	8078.4668	0.0000	1.0000
HKYIG	41	-3998.3888	8083.2273	0.0000	1.0000
SYMG	41	-4018.3891	8123.2278	0.0000	1.0000
SYMIG	42	-4017.5372	8123.7471	0.0000	1.0000
K2PG	37	-4024.7321	8127.0785	0.0000	1.0000
GTRI	44	-4017.1116	8127.3595	0.0000	1.0000
K2PIG	38	-4024.7321	8129.2789	0.0000	1.0000
HKYI	40	-4031.8199	8147.8720	0.0000	1.0000
SYMI	41	-4048.7656	8183.9808	0.0000	1.0000
K2PI	37	-4056.9495	8191.5134	0.0000	1.0000
GTR	43	-4053.6015	8198.1046	0.0000	1.0000
HKY	39	-4075.0352	8232.0909	0.0000	1.0000
SYM	40	-4105.0247	8294.2816	0.0000	1.0000
K2P	36	-4113.0270	8301.4738	0.0000	1.0000
F81G	39	-4120.8575	8323.7356	0.0000	1.0000
F81IG	40	-4120.8410	8325.9143	0.0000	1.0000
JC69G	36	-4142.4685	8360.3568	0.0000	1.0000
JC69IG	37	-4145.3041	8368.2226	0.0000	1.0000
F81I	39	-4147.9387	8377.8980	0.0000	1.0000
JC69I	36	-4167.8990	8411.2177	0.0000	1.0000
F81	38	-4193.7918	8467.3983	0.0000	1.0000
JC69	35	-4217.6981	8508.6270	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-3977.8827	8043.7653	0.8256	0.8256
GTRIG	45	-3978.4377	8046.8753	0.1744	1.0000
HKYG	40	-3997.1173	8074.2346	0.0000	1.0000
HKYIG	41	-3998.3888	8078.7776	0.0000	1.0000
SYMG	41	-4018.3891	8118.7782	0.0000	1.0000
SYMIG	42	-4017.5372	8119.0744	0.0000	1.0000
GTRI	44	-4017.1116	8122.2233	0.0000	1.0000
K2PG	37	-4024.7321	8123.4641	0.0000	1.0000
K2PIG	38	-4024.7321	8125.4642	0.0000	1.0000
HKYI	40	-4031.8199	8143.6398	0.0000	1.0000
SYMI	41	-4048.7656	8179.5312	0.0000	1.0000
K2PI	37	-4056.9495	8187.8990	0.0000	1.0000
GTR	43	-4053.6015	8193.2031	0.0000	1.0000
HKY	39	-4075.0352	8228.0703	0.0000	1.0000
SYM	40	-4105.0247	8290.0494	0.0000	1.0000
K2P	36	-4113.0270	8298.0540	0.0000	1.0000
F81G	39	-4120.8575	8319.7150	0.0000	1.0000
F81IG	40	-4120.8410	8321.6820	0.0000	1.0000
JC69G	36	-4142.4685	8356.9370	0.0000	1.0000
JC69IG	37	-4145.3041	8364.6082	0.0000	1.0000
F81I	39	-4147.9387	8373.8774	0.0000	1.0000
JC69I	36	-4167.8990	8407.7979	0.0000	1.0000
F81	38	-4193.7918	8463.5836	0.0000	1.0000
JC69	35	-4217.6981	8505.3962	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	44	-3977.8827	8250.7596	0.9774	0.9774
GTRIG	45	-3978.4377	8258.5740	0.0196	0.9971
HKYG	40	-3997.1173	8262.4111	0.0029	1.0000
HKYIG	41	-3998.3888	8271.6586	0.0000	1.0000
K2PG	37	-4024.7321	8297.5274	0.0000	1.0000
K2PIG	38	-4024.7321	8304.2320	0.0000	1.0000
SYMG	41	-4018.3891	8311.6592	0.0000	1.0000
SYMIG	42	-4017.5372	8316.6598	0.0000	1.0000
GTRI	44	-4017.1116	8329.2175	0.0000	1.0000
HKYI	40	-4031.8199	8331.8163	0.0000	1.0000
K2PI	37	-4056.9495	8361.9623	0.0000	1.0000
SYMI	41	-4048.7656	8372.4122	0.0000	1.0000
GTR	43	-4053.6015	8395.4929	0.0000	1.0000
HKY	39	-4075.0352	8411.5425	0.0000	1.0000
K2P	36	-4113.0270	8467.4130	0.0000	1.0000
SYM	40	-4105.0247	8478.2260	0.0000	1.0000
F81G	39	-4120.8575	8503.1871	0.0000	1.0000
F81IG	40	-4120.8410	8509.8586	0.0000	1.0000
JC69G	36	-4142.4685	8526.2959	0.0000	1.0000
JC69IG	37	-4145.3041	8538.6716	0.0000	1.0000
F81I	39	-4147.9387	8557.3496	0.0000	1.0000
JC69I	36	-4167.8990	8577.1568	0.0000	1.0000
F81	38	-4193.7918	8642.3514	0.0000	1.0000
JC69	35	-4217.6981	8670.0507	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/chr5_3273.phylip = 1 - 816;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_3273.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 (GTRG)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_3273.phylip = 1 - 816;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_3273.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 BIC model (GTRG)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_3273.phylip = 1 - 816;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_3273.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;


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

