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

Minimum AIC  model: HKYG
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	-907.8899	1902.8588	0.3058	0.3058
K2PI	37	-908.1122	1903.3035	0.2449	0.5507
HKYG	40	-904.2728	1904.0173	0.1714	0.7220
HKYI	40	-904.4063	1904.2843	0.1499	0.8720
K2PIG	38	-907.8878	1905.6260	0.0767	0.9486
HKYIG	41	-904.1978	1906.7180	0.0444	0.9931
GTRG	44	-903.1168	1913.2722	0.0017	0.9947
GTRI	44	-903.1825	1913.4034	0.0016	0.9963
SYMG	41	-907.5481	1913.4185	0.0016	0.9979
SYMI	41	-907.7118	1913.7458	0.0013	0.9992
GTRIG	45	-902.9887	1915.9775	0.0004	0.9996
SYMIG	42	-907.5250	1916.2500	0.0004	1.0000
K2P	36	-920.1942	1924.7218	0.0000	1.0000
HKY	39	-916.0481	1924.7441	0.0000	1.0000
GTR	43	-914.6247	1933.3547	0.0000	1.0000
SYM	40	-919.5174	1934.5065	0.0000	1.0000
JC69G	36	-945.7081	1975.7496	0.0000	1.0000
JC69I	36	-945.7714	1975.8762	0.0000	1.0000
JC69IG	37	-945.6284	1978.3359	0.0000	1.0000
F81G	39	-942.9911	1978.6300	0.0000	1.0000
F81I	39	-943.0030	1978.6538	0.0000	1.0000
F81IG	40	-942.8877	1981.2472	0.0000	1.0000
JC69	35	-956.4265	1994.4659	0.0000	1.0000
F81	38	-953.4398	1996.7300	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-904.2728	1888.5456	0.2770	0.2770
HKYI	40	-904.4063	1888.8126	0.2424	0.5194
K2PG	37	-907.8899	1889.7797	0.1495	0.6689
K2PI	37	-908.1122	1890.2244	0.1197	0.7885
HKYIG	41	-904.1978	1890.3957	0.1098	0.8984
K2PIG	38	-907.8878	1891.7756	0.0551	0.9535
GTRG	44	-903.1168	1894.2337	0.0161	0.9696
GTRI	44	-903.1825	1894.3649	0.0151	0.9847
GTRIG	45	-902.9887	1895.9775	0.0067	0.9914
SYMG	41	-907.5481	1897.0962	0.0039	0.9953
SYMI	41	-907.7118	1897.4236	0.0033	0.9985
SYMIG	42	-907.5250	1899.0500	0.0015	1.0000
HKY	39	-916.0481	1910.0962	0.0000	1.0000
K2P	36	-920.1942	1912.3885	0.0000	1.0000
GTR	43	-914.6247	1915.2495	0.0000	1.0000
SYM	40	-919.5174	1919.0348	0.0000	1.0000
JC69G	36	-945.7081	1963.4163	0.0000	1.0000
JC69I	36	-945.7714	1963.5429	0.0000	1.0000
F81G	39	-942.9911	1963.9821	0.0000	1.0000
F81I	39	-943.0030	1964.0059	0.0000	1.0000
JC69IG	37	-945.6284	1965.2568	0.0000	1.0000
F81IG	40	-942.8877	1965.7755	0.0000	1.0000
JC69	35	-956.4265	1982.8530	0.0000	1.0000
F81	38	-953.4398	1982.8795	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	37	-907.8899	2020.5152	0.5314	0.5314
K2PI	37	-908.1122	2020.9598	0.4255	0.9569
K2PIG	38	-907.8878	2026.0444	0.0335	0.9904
HKYG	40	-904.2728	2029.8812	0.0049	0.9953
HKYI	40	-904.4063	2030.1482	0.0043	0.9996
HKYIG	41	-904.1978	2035.2646	0.0003	0.9999
K2P	36	-920.1942	2039.5905	0.0000	1.0000
SYMG	41	-907.5481	2041.9652	0.0000	1.0000
SYMI	41	-907.7118	2042.2925	0.0000	1.0000
SYMIG	42	-907.5250	2047.4524	0.0000	1.0000
HKY	39	-916.0481	2047.8984	0.0000	1.0000
GTRG	44	-903.1168	2049.7028	0.0000	1.0000
GTRI	44	-903.1825	2049.8341	0.0000	1.0000
GTRIG	45	-902.9887	2054.9800	0.0000	1.0000
SYM	40	-919.5174	2060.3703	0.0000	1.0000
GTR	43	-914.6247	2067.1852	0.0000	1.0000
JC69G	36	-945.7081	2090.6183	0.0000	1.0000
JC69I	36	-945.7714	2090.7449	0.0000	1.0000
JC69IG	37	-945.6284	2095.9923	0.0000	1.0000
F81G	39	-942.9911	2101.7843	0.0000	1.0000
F81I	39	-943.0030	2101.8081	0.0000	1.0000
JC69	35	-956.4265	2106.5216	0.0000	1.0000
F81IG	40	-942.8877	2107.1111	0.0000	1.0000
F81	38	-953.4398	2117.1483	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/chr3_18148.phylip = 1 - 253;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_18148.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 (K2PG)]

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

