
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:52 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/chr7_6333.phylip" (19 taxa, 450 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	-1885.9960	3860.0115	0.6783	0.6783
HKYIG	41	-1885.9960	3862.4331	0.2021	0.8804
GTRG	44	-1883.4390	3864.6557	0.0665	0.9469
HKYI	40	-1889.0955	3866.2105	0.0306	0.9774
GTRIG	45	-1883.4390	3867.1254	0.0193	0.9968
GTRI	44	-1886.4839	3870.7456	0.0032	1.0000
HKY	39	-1896.8677	3879.3451	0.0000	1.0000
GTR	43	-1894.6675	3884.6552	0.0000	1.0000
K2PG	37	-1908.0586	3896.9424	0.0000	1.0000
K2PIG	38	-1908.0587	3899.3290	0.0000	1.0000
K2PI	37	-1911.1066	3903.0385	0.0000	1.0000
SYMG	41	-1906.6845	3903.8101	0.0000	1.0000
SYMIG	42	-1906.6845	3906.2438	0.0000	1.0000
SYMI	41	-1909.6036	3909.6484	0.0000	1.0000
K2P	36	-1918.8378	3916.1259	0.0000	1.0000
SYM	40	-1917.4075	3922.8345	0.0000	1.0000
F81G	39	-1941.4529	3968.5156	0.0000	1.0000
F81IG	40	-1941.4529	3970.9254	0.0000	1.0000
F81I	39	-1944.0303	3973.6703	0.0000	1.0000
F81	38	-1951.2985	3985.8086	0.0000	1.0000
JC69G	36	-1958.1564	3994.7632	0.0000	1.0000
JC69IG	37	-1958.1565	3997.1383	0.0000	1.0000
JC69I	36	-1960.8050	4000.0604	0.0000	1.0000
JC69	35	-1968.1897	4012.4664	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-1885.9960	3851.9919	0.5723	0.5723
HKYIG	41	-1885.9960	3853.9919	0.2105	0.7828
GTRG	44	-1883.4390	3854.8779	0.1352	0.9180
GTRIG	45	-1883.4390	3856.8779	0.0497	0.9677
HKYI	40	-1889.0955	3858.1909	0.0258	0.9935
GTRI	44	-1886.4839	3860.9678	0.0064	1.0000
HKY	39	-1896.8677	3871.7353	0.0000	1.0000
GTR	43	-1894.6675	3875.3350	0.0000	1.0000
K2PG	37	-1908.0586	3890.1172	0.0000	1.0000
K2PIG	38	-1908.0587	3892.1173	0.0000	1.0000
SYMG	41	-1906.6845	3895.3689	0.0000	1.0000
K2PI	37	-1911.1066	3896.2132	0.0000	1.0000
SYMIG	42	-1906.6845	3897.3691	0.0000	1.0000
SYMI	41	-1909.6036	3901.2073	0.0000	1.0000
K2P	36	-1918.8378	3909.6756	0.0000	1.0000
SYM	40	-1917.4075	3914.8150	0.0000	1.0000
F81G	39	-1941.4529	3960.9059	0.0000	1.0000
F81IG	40	-1941.4529	3962.9059	0.0000	1.0000
F81I	39	-1944.0303	3966.0605	0.0000	1.0000
F81	38	-1951.2985	3978.5970	0.0000	1.0000
JC69G	36	-1958.1564	3988.3129	0.0000	1.0000
JC69IG	37	-1958.1565	3990.3130	0.0000	1.0000
JC69I	36	-1960.8050	3993.6100	0.0000	1.0000
JC69	35	-1968.1897	4006.3795	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-1885.9960	4016.3618	0.9152	0.9152
HKYIG	41	-1885.9960	4022.4711	0.0431	0.9583
HKYI	40	-1889.0955	4022.5608	0.0412	0.9996
HKY	39	-1896.8677	4031.9960	0.0004	0.9999
GTRG	44	-1883.4390	4035.6848	0.0001	1.0000
GTRI	44	-1886.4839	4041.7747	0.0000	1.0000
GTRIG	45	-1883.4390	4041.7941	0.0000	1.0000
K2PG	37	-1908.0586	4042.1594	0.0000	1.0000
K2PI	37	-1911.1066	4048.2554	0.0000	1.0000
K2PIG	38	-1908.0587	4048.2687	0.0000	1.0000
GTR	43	-1894.6675	4052.0326	0.0000	1.0000
K2P	36	-1918.8378	4057.6085	0.0000	1.0000
SYMG	41	-1906.6845	4063.8481	0.0000	1.0000
SYMI	41	-1909.6036	4069.6864	0.0000	1.0000
SYMIG	42	-1906.6845	4069.9575	0.0000	1.0000
SYM	40	-1917.4075	4079.1849	0.0000	1.0000
F81G	39	-1941.4529	4121.1665	0.0000	1.0000
F81I	39	-1944.0303	4126.3212	0.0000	1.0000
F81IG	40	-1941.4529	4127.2758	0.0000	1.0000
F81	38	-1951.2985	4134.7484	0.0000	1.0000
JC69G	36	-1958.1564	4136.2458	0.0000	1.0000
JC69I	36	-1960.8050	4141.5429	0.0000	1.0000
JC69IG	37	-1958.1565	4142.3552	0.0000	1.0000
JC69	35	-1968.1897	4150.2031	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/chr7_6333.phylip = 1 - 450;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_6333.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/chr7_6333.phylip = 1 - 450;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_6333.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/chr7_6333.phylip = 1 - 450;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_6333.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

