
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
13:07 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/chrZ_3271.phylip" (19 taxa, 571 characters)

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	44	-1882.5097	3860.5480	0.3893	0.3893
HKYG	40	-1887.4836	3861.1559	0.2873	0.6765
GTRIG	45	-1882.1294	3862.1446	0.1752	0.8518
HKYIG	41	-1886.9844	3862.4793	0.1482	1.0000
GTRI	44	-1893.4267	3882.3820	0.0000	1.0000
HKYI	40	-1898.5262	3883.2411	0.0000	1.0000
K2PG	37	-1903.0231	3885.3220	0.0000	1.0000
K2PIG	38	-1902.3472	3886.2658	0.0000	1.0000
SYMG	41	-1898.9270	3886.3645	0.0000	1.0000
SYMIG	42	-1898.2863	3887.4135	0.0000	1.0000
K2PI	37	-1905.2914	3889.8586	0.0000	1.0000
SYMI	41	-1901.0885	3890.6874	0.0000	1.0000
HKY	39	-1932.5365	3948.9487	0.0000	1.0000
GTR	43	-1928.8743	3950.9289	0.0000	1.0000
F81G	39	-1942.9550	3969.7857	0.0000	1.0000
SYM	40	-1942.0782	3970.3451	0.0000	1.0000
K2P	36	-1946.7459	3970.4805	0.0000	1.0000
F81I	39	-1945.7725	3975.4207	0.0000	1.0000
F81IG	40	-1949.2785	3984.7456	0.0000	1.0000
JC69G	36	-1964.2108	4005.4103	0.0000	1.0000
JC69IG	37	-1963.7535	4006.7829	0.0000	1.0000
JC69I	36	-1966.2504	4009.4896	0.0000	1.0000
F81	38	-1979.6274	4040.8262	0.0000	1.0000
JC69	35	-2000.2646	4075.2394	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-1882.5097	3853.0195	0.4663	0.4663
GTRIG	45	-1882.1294	3854.2588	0.2509	0.7172
HKYG	40	-1887.4836	3854.9672	0.1761	0.8933
HKYIG	41	-1886.9844	3855.9689	0.1067	1.0000
GTRI	44	-1893.4267	3874.8534	0.0000	1.0000
HKYI	40	-1898.5262	3877.0524	0.0000	1.0000
SYMG	41	-1898.9270	3879.8541	0.0000	1.0000
K2PG	37	-1903.0231	3880.0462	0.0000	1.0000
SYMIG	42	-1898.2863	3880.5726	0.0000	1.0000
K2PIG	38	-1902.3472	3880.6944	0.0000	1.0000
SYMI	41	-1901.0885	3884.1770	0.0000	1.0000
K2PI	37	-1905.2914	3884.5828	0.0000	1.0000
HKY	39	-1932.5365	3943.0730	0.0000	1.0000
GTR	43	-1928.8743	3943.7487	0.0000	1.0000
F81G	39	-1942.9550	3963.9100	0.0000	1.0000
SYM	40	-1942.0782	3964.1565	0.0000	1.0000
K2P	36	-1946.7459	3965.4918	0.0000	1.0000
F81I	39	-1945.7725	3969.5450	0.0000	1.0000
F81IG	40	-1949.2785	3978.5569	0.0000	1.0000
JC69G	36	-1964.2108	4000.4215	0.0000	1.0000
JC69IG	37	-1963.7535	4001.5071	0.0000	1.0000
JC69I	36	-1966.2504	4004.5008	0.0000	1.0000
F81	38	-1979.6274	4035.2548	0.0000	1.0000
JC69	35	-2000.2646	4070.5291	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-1887.4836	4028.8627	0.9326	0.9326
HKYIG	41	-1886.9844	4034.2119	0.0643	0.9969
K2PG	37	-1903.0231	4040.8996	0.0023	0.9991
GTRG	44	-1882.5097	4044.3046	0.0004	0.9995
K2PI	37	-1905.2914	4045.4362	0.0002	0.9998
K2PIG	38	-1902.3472	4045.8952	0.0002	1.0000
GTRIG	45	-1882.1294	4049.8914	0.0000	1.0000
HKYI	40	-1898.5262	4050.9480	0.0000	1.0000
SYMG	41	-1898.9270	4058.0970	0.0000	1.0000
SYMI	41	-1901.0885	4062.4199	0.0000	1.0000
SYMIG	42	-1898.2863	4063.1629	0.0000	1.0000
GTRI	44	-1893.4267	4066.1386	0.0000	1.0000
HKY	39	-1932.5365	4112.6212	0.0000	1.0000
K2P	36	-1946.7459	4121.9978	0.0000	1.0000
GTR	43	-1928.8743	4130.6864	0.0000	1.0000
F81G	39	-1942.9550	4133.4582	0.0000	1.0000
SYM	40	-1942.0782	4138.0520	0.0000	1.0000
F81I	39	-1945.7725	4139.0932	0.0000	1.0000
F81IG	40	-1949.2785	4152.4525	0.0000	1.0000
JC69G	36	-1964.2108	4156.9275	0.0000	1.0000
JC69I	36	-1966.2504	4161.0069	0.0000	1.0000
JC69IG	37	-1963.7535	4162.3605	0.0000	1.0000
F81	38	-1979.6274	4200.4556	0.0000	1.0000
JC69	35	-2000.2646	4222.6877	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/chrZ_3271.phylip = 1 - 571;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_3271.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/chrZ_3271.phylip = 1 - 571;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_3271.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 (HKYG)]

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

