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

Minimum AIC  model: GTRIG
Minimum AICc model: GTRIG
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
GTRIG	45	-1987.7377	4072.1635	0.6289	0.6289
GTRG	44	-1989.5150	4073.4172	0.3360	0.9649
GTRI	44	-1991.7743	4077.9357	0.0351	1.0000
SYMIG	42	-2003.3806	4096.5683	0.0000	1.0000
SYMG	41	-2005.0349	4097.5980	0.0000	1.0000
SYMI	41	-2006.4607	4100.4495	0.0000	1.0000
HKYIG	41	-2006.6249	4100.7779	0.0000	1.0000
HKYG	40	-2008.4201	4102.0966	0.0000	1.0000
HKYI	40	-2009.7481	4104.7525	0.0000	1.0000
K2PIG	38	-2014.7234	4110.1817	0.0000	1.0000
K2PG	37	-2017.4381	4113.3611	0.0000	1.0000
K2PI	37	-2018.4238	4115.3324	0.0000	1.0000
F81IG	40	-2034.5867	4154.4298	0.0000	1.0000
F81G	39	-2036.1129	4155.2178	0.0000	1.0000
F81I	39	-2037.3360	4157.6640	0.0000	1.0000
JC69IG	37	-2041.4974	4161.4796	0.0000	1.0000
JC69G	36	-2043.0474	4162.3368	0.0000	1.0000
JC69I	36	-2044.1949	4164.6319	0.0000	1.0000
GTR	43	-2042.1315	4176.3564	0.0000	1.0000
SYM	40	-2058.8978	4203.0521	0.0000	1.0000
HKY	39	-2063.2551	4209.5023	0.0000	1.0000
K2P	36	-2072.4828	4221.2077	0.0000	1.0000
F81	38	-2090.8607	4262.4561	0.0000	1.0000
JC69	35	-2098.3224	4270.6512	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	45	-1987.7377	4065.4753	0.6633	0.6633
GTRG	44	-1989.5150	4067.0301	0.3049	0.9682
GTRI	44	-1991.7743	4071.5486	0.0318	1.0000
SYMIG	42	-2003.3806	4090.7612	0.0000	1.0000
SYMG	41	-2005.0349	4092.0699	0.0000	1.0000
SYMI	41	-2006.4607	4094.9214	0.0000	1.0000
HKYIG	41	-2006.6249	4095.2498	0.0000	1.0000
HKYG	40	-2008.4201	4096.8402	0.0000	1.0000
HKYI	40	-2009.7481	4099.4961	0.0000	1.0000
K2PIG	38	-2014.7234	4105.4468	0.0000	1.0000
K2PG	37	-2017.4381	4108.8763	0.0000	1.0000
K2PI	37	-2018.4238	4110.8475	0.0000	1.0000
F81IG	40	-2034.5867	4149.1733	0.0000	1.0000
F81G	39	-2036.1129	4150.2258	0.0000	1.0000
F81I	39	-2037.3360	4152.6720	0.0000	1.0000
JC69IG	37	-2041.4974	4156.9948	0.0000	1.0000
JC69G	36	-2043.0474	4158.0948	0.0000	1.0000
JC69I	36	-2044.1949	4160.3898	0.0000	1.0000
GTR	43	-2042.1315	4170.2630	0.0000	1.0000
SYM	40	-2058.8978	4197.7957	0.0000	1.0000
HKY	39	-2063.2551	4204.5103	0.0000	1.0000
K2P	36	-2072.4828	4216.9656	0.0000	1.0000
F81	38	-2090.8607	4257.7213	0.0000	1.0000
JC69	35	-2098.3224	4266.6449	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	44	-1989.5150	4265.0207	0.7389	0.7389
GTRIG	45	-1987.7377	4267.9658	0.1695	0.9083
GTRI	44	-1991.7743	4269.5392	0.0772	0.9855
K2PG	37	-2017.4381	4275.3684	0.0042	0.9896
K2PIG	38	-2014.7234	4276.4387	0.0024	0.9921
SYMG	41	-2005.0349	4276.5611	0.0023	0.9944
HKYG	40	-2008.4201	4276.8317	0.0020	0.9964
K2PI	37	-2018.4238	4277.3396	0.0016	0.9980
SYMI	41	-2006.4607	4279.4127	0.0006	0.9985
HKYI	40	-2009.7481	4279.4876	0.0005	0.9991
HKYIG	41	-2006.6249	4279.7411	0.0005	0.9995
SYMIG	42	-2003.3806	4279.7522	0.0005	1.0000
JC69G	36	-2043.0474	4320.0871	0.0000	1.0000
JC69I	36	-2044.1949	4322.3822	0.0000	1.0000
JC69IG	37	-2041.4974	4323.4869	0.0000	1.0000
F81G	39	-2036.1129	4325.7175	0.0000	1.0000
F81I	39	-2037.3360	4328.1637	0.0000	1.0000
F81IG	40	-2034.5867	4329.1648	0.0000	1.0000
GTR	43	-2042.1315	4363.7539	0.0000	1.0000
SYM	40	-2058.8978	4377.7871	0.0000	1.0000
K2P	36	-2072.4828	4378.9580	0.0000	1.0000
HKY	39	-2063.2551	4380.0020	0.0000	1.0000
JC69	35	-2098.3224	4424.1374	0.0000	1.0000
F81	38	-2090.8607	4428.7132	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/chr6_9559.phylip = 1 - 665;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_9559.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/chr6_9559.phylip = 1 - 665;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_9559.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 (GTRG)]

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

