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

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


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
K2PG	37	-1913.0955	3904.4776	0.2374	0.2374
K2PI	37	-1913.6118	3905.5102	0.1417	0.3791
GTRG	44	-1905.8603	3905.8223	0.1212	0.5003
SYMG	41	-1909.3433	3905.9687	0.1127	0.6130
K2PIG	38	-1912.8915	3906.3082	0.0951	0.7081
SYMI	41	-1909.9857	3907.2537	0.0593	0.7673
GTRI	44	-1906.5944	3907.2905	0.0582	0.8255
GTRIG	45	-1905.6262	3907.6413	0.0488	0.8743
SYMIG	42	-1909.1312	3907.8108	0.0448	0.9192
HKYG	40	-1911.4864	3907.9959	0.0409	0.9601
HKYI	40	-1912.0361	3909.0952	0.0236	0.9836
HKYIG	41	-1911.2737	3909.8296	0.0163	1.0000
JC69G	36	-1925.3366	3926.7281	0.0000	1.0000
JC69I	36	-1926.0837	3928.2222	0.0000	1.0000
JC69IG	37	-1925.0287	3928.3440	0.0000	1.0000
F81G	39	-1923.9303	3930.6312	0.0000	1.0000
F81I	39	-1924.6863	3932.1433	0.0000	1.0000
F81IG	40	-1923.6327	3932.2883	0.0000	1.0000
GTR	43	-1933.1564	3958.1344	0.0000	1.0000
K2P	36	-1941.1770	3958.4087	0.0000	1.0000
SYM	40	-1936.8282	3958.6793	0.0000	1.0000
HKY	39	-1939.1907	3961.1521	0.0000	1.0000
JC69	35	-1953.4070	3980.6437	0.0000	1.0000
F81	38	-1951.6877	3983.9005	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-1905.8603	3899.7206	0.1929	0.1929
K2PG	37	-1913.0955	3900.1910	0.1524	0.3453
SYMG	41	-1909.3433	3900.6865	0.1190	0.4643
GTRI	44	-1906.5944	3901.1888	0.0926	0.5568
K2PI	37	-1913.6118	3901.2236	0.0910	0.6478
GTRIG	45	-1905.6262	3901.2525	0.0897	0.7374
K2PIG	38	-1912.8915	3901.7830	0.0688	0.8062
SYMI	41	-1909.9857	3901.9715	0.0626	0.8688
SYMIG	42	-1909.1312	3902.2624	0.0541	0.9229
HKYG	40	-1911.4864	3902.9729	0.0379	0.9608
HKYI	40	-1912.0361	3904.0722	0.0219	0.9827
HKYIG	41	-1911.2737	3904.5474	0.0173	1.0000
JC69G	36	-1925.3366	3922.6733	0.0000	1.0000
JC69IG	37	-1925.0287	3924.0574	0.0000	1.0000
JC69I	36	-1926.0837	3924.1674	0.0000	1.0000
F81G	39	-1923.9303	3925.8606	0.0000	1.0000
F81IG	40	-1923.6327	3927.2653	0.0000	1.0000
F81I	39	-1924.6863	3927.3726	0.0000	1.0000
GTR	43	-1933.1564	3952.3128	0.0000	1.0000
SYM	40	-1936.8282	3953.6563	0.0000	1.0000
K2P	36	-1941.1770	3954.3539	0.0000	1.0000
HKY	39	-1939.1907	3956.3814	0.0000	1.0000
JC69	35	-1953.4070	3976.8139	0.0000	1.0000
F81	38	-1951.6877	3979.3753	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	37	-1913.0955	4068.2625	0.6083	0.6083
K2PI	37	-1913.6118	4069.2951	0.3630	0.9712
K2PIG	38	-1912.8915	4074.3970	0.0283	0.9995
HKYG	40	-1911.4864	4084.6718	0.0002	0.9997
HKYI	40	-1912.0361	4085.7711	0.0001	0.9998
JC69G	36	-1925.3366	4086.2023	0.0001	0.9999
SYMG	41	-1909.3433	4086.9279	0.0001	0.9999
JC69I	36	-1926.0837	4087.6964	0.0000	1.0000
SYMI	41	-1909.9857	4088.2128	0.0000	1.0000
HKYIG	41	-1911.2737	4090.7887	0.0000	1.0000
JC69IG	37	-1925.0287	4092.1289	0.0000	1.0000
SYMIG	42	-1909.1312	4093.0462	0.0000	1.0000
GTRG	44	-1905.8603	4099.5894	0.0000	1.0000
GTRI	44	-1906.5944	4101.0576	0.0000	1.0000
F81G	39	-1923.9303	4103.0170	0.0000	1.0000
F81I	39	-1924.6863	4104.5290	0.0000	1.0000
GTRIG	45	-1905.6262	4105.6637	0.0000	1.0000
F81IG	40	-1923.6327	4108.9642	0.0000	1.0000
K2P	36	-1941.1770	4117.8829	0.0000	1.0000
HKY	39	-1939.1907	4133.5378	0.0000	1.0000
SYM	40	-1936.8282	4135.3552	0.0000	1.0000
JC69	35	-1953.4070	4135.8005	0.0000	1.0000
GTR	43	-1933.1564	4147.6391	0.0000	1.0000
F81	38	-1951.6877	4151.9892	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/chr8_4014.phylip = 1 - 694;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr8_4014.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 (K2PG)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr8_4014.phylip = 1 - 694;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr8_4014.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


[Mrbayes block for the best BIC model (K2PG)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr8_4014.phylip = 1 - 694;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr8_4014.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

