
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12: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/chr6_5852.phylip" (19 taxa, 1229 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	-6263.8367	12610.4344	0.7589	0.7589
GTRG	44	-6261.3019	12613.9484	0.1310	0.8899
GTRIG	45	-6260.4064	12614.3124	0.1092	0.9990
HKYIG	41	-6269.4274	12623.7563	0.0010	1.0000
HKYI	40	-6300.5634	12683.8878	0.0000	1.0000
GTRI	44	-6296.3832	12684.1111	0.0000	1.0000
SYMG	41	-6312.7970	12710.4954	0.0000	1.0000
SYMIG	42	-6312.2768	12711.5992	0.0000	1.0000
K2PIG	38	-6319.2335	12716.9577	0.0000	1.0000
K2PG	37	-6326.8229	12730.0068	0.0000	1.0000
SYMI	41	-6352.3914	12789.6843	0.0000	1.0000
K2PI	37	-6365.4293	12807.2197	0.0000	1.0000
F81G	39	-6419.4498	12919.5236	0.0000	1.0000
F81IG	40	-6419.2127	12921.1863	0.0000	1.0000
F81I	39	-6450.5733	12981.7706	0.0000	1.0000
JC69G	36	-6458.5082	12991.2513	0.0000	1.0000
JC69IG	37	-6458.2995	12992.9601	0.0000	1.0000
JC69I	36	-6493.9780	13062.1909	0.0000	1.0000
GTR	43	-6525.7704	13140.7340	0.0000	1.0000
HKY	39	-6531.7661	13144.1562	0.0000	1.0000
SYM	40	-6562.0920	13206.9450	0.0000	1.0000
K2P	36	-6577.7007	13229.6364	0.0000	1.0000
F81	38	-6659.6751	13397.8409	0.0000	1.0000
JC69	35	-6695.1884	13462.4892	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-6263.8367	12607.6735	0.6942	0.6942
GTRG	44	-6261.3019	12610.6038	0.1604	0.8546
GTRIG	45	-6260.4064	12610.8128	0.1445	0.9990
HKYIG	41	-6269.4274	12620.8548	0.0010	1.0000
GTRI	44	-6296.3832	12680.7665	0.0000	1.0000
HKYI	40	-6300.5634	12681.1269	0.0000	1.0000
SYMG	41	-6312.7970	12707.5939	0.0000	1.0000
SYMIG	42	-6312.2768	12708.5536	0.0000	1.0000
K2PIG	38	-6319.2335	12714.4670	0.0000	1.0000
K2PG	37	-6326.8229	12727.6458	0.0000	1.0000
SYMI	41	-6352.3914	12786.7829	0.0000	1.0000
K2PI	37	-6365.4293	12804.8586	0.0000	1.0000
F81G	39	-6419.4498	12916.8995	0.0000	1.0000
F81IG	40	-6419.2127	12918.4254	0.0000	1.0000
F81I	39	-6450.5733	12979.1466	0.0000	1.0000
JC69G	36	-6458.5082	12989.0164	0.0000	1.0000
JC69IG	37	-6458.2995	12990.5991	0.0000	1.0000
JC69I	36	-6493.9780	13059.9560	0.0000	1.0000
GTR	43	-6525.7704	13137.5408	0.0000	1.0000
HKY	39	-6531.7661	13141.5321	0.0000	1.0000
SYM	40	-6562.0920	13204.1840	0.0000	1.0000
K2P	36	-6577.7007	13227.4015	0.0000	1.0000
F81	38	-6659.6751	13395.3501	0.0000	1.0000
JC69	35	-6695.1884	13460.3769	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-6263.8367	12812.2317	0.9999	0.9999
HKYIG	41	-6269.4274	12830.5270	0.0001	1.0000
GTRG	44	-6261.3019	12835.6178	0.0000	1.0000
GTRIG	45	-6260.4064	12840.9408	0.0000	1.0000
HKYI	40	-6300.5634	12885.6851	0.0000	1.0000
GTRI	44	-6296.3832	12905.7806	0.0000	1.0000
K2PIG	38	-6319.2335	12908.7973	0.0000	1.0000
K2PG	37	-6326.8229	12916.8622	0.0000	1.0000
SYMG	41	-6312.7970	12917.2661	0.0000	1.0000
SYMIG	42	-6312.2768	12923.3398	0.0000	1.0000
K2PI	37	-6365.4293	12994.0750	0.0000	1.0000
SYMI	41	-6352.3914	12996.4551	0.0000	1.0000
F81G	39	-6419.4498	13116.3438	0.0000	1.0000
F81IG	40	-6419.2127	13122.9836	0.0000	1.0000
JC69G	36	-6458.5082	13173.1188	0.0000	1.0000
F81I	39	-6450.5733	13178.5909	0.0000	1.0000
JC69IG	37	-6458.2995	13179.8155	0.0000	1.0000
JC69I	36	-6493.9780	13244.0584	0.0000	1.0000
HKY	39	-6531.7661	13340.9764	0.0000	1.0000
GTR	43	-6525.7704	13357.4409	0.0000	1.0000
SYM	40	-6562.0920	13408.7423	0.0000	1.0000
K2P	36	-6577.7007	13411.5039	0.0000	1.0000
F81	38	-6659.6751	13589.6805	0.0000	1.0000
JC69	35	-6695.1884	13639.3653	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/chr6_5852.phylip = 1 - 1229;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_5852.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/chr6_5852.phylip = 1 - 1229;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_5852.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/chr6_5852.phylip = 1 - 1229;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_5852.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

