
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:56 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_10849.phylip" (19 taxa, 445 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	-1039.8430	2167.8049	0.3875	0.3875
HKYI	40	-1039.8434	2167.8057	0.3874	0.7749
HKYIG	41	-1039.8586	2170.2632	0.1134	0.8883
GTRG	44	-1037.0462	2171.9924	0.0478	0.9360
GTRI	44	-1037.1734	2172.2467	0.0421	0.9781
GTRIG	45	-1036.9918	2174.3595	0.0146	0.9927
F81I	39	-1046.1588	2178.0213	0.0023	0.9951
F81G	39	-1046.2031	2178.1099	0.0022	0.9973
HKY	39	-1046.3874	2178.4786	0.0019	0.9992
F81IG	40	-1046.2408	2180.6004	0.0006	0.9998
GTR	43	-1043.8513	2183.1390	0.0002	1.0000
F81	38	-1052.7215	2188.7435	0.0000	1.0000
K2PI	37	-1064.5937	2210.0964	0.0000	1.0000
K2PG	37	-1064.6513	2210.2117	0.0000	1.0000
K2PIG	38	-1064.6512	2212.6030	0.0000	1.0000
JC69I	36	-1069.4745	2217.4784	0.0000	1.0000
JC69G	36	-1069.5440	2217.6174	0.0000	1.0000
SYMI	41	-1063.9293	2218.4045	0.0000	1.0000
SYMG	41	-1063.9456	2218.4371	0.0000	1.0000
JC69IG	37	-1069.5992	2220.1075	0.0000	1.0000
SYMIG	42	-1063.9705	2220.9261	0.0000	1.0000
K2P	36	-1071.3700	2221.2694	0.0000	1.0000
JC69	35	-1076.1940	2228.5493	0.0000	1.0000
SYM	40	-1070.5862	2229.2911	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-1039.8430	2159.6861	0.3269	0.3269
HKYI	40	-1039.8434	2159.6869	0.3268	0.6536
HKYIG	41	-1039.8586	2161.7173	0.1184	0.7720
GTRG	44	-1037.0462	2162.0924	0.0981	0.8702
GTRI	44	-1037.1734	2162.3467	0.0864	0.9566
GTRIG	45	-1036.9918	2163.9836	0.0381	0.9947
F81I	39	-1046.1588	2170.3176	0.0016	0.9963
F81G	39	-1046.2031	2170.4062	0.0015	0.9979
HKY	39	-1046.3874	2170.7749	0.0013	0.9992
F81IG	40	-1046.2408	2172.4816	0.0005	0.9997
GTR	43	-1043.8513	2173.7026	0.0003	1.0000
F81	38	-1052.7215	2181.4430	0.0000	1.0000
K2PI	37	-1064.5937	2203.1873	0.0000	1.0000
K2PG	37	-1064.6513	2203.3026	0.0000	1.0000
K2PIG	38	-1064.6512	2205.3025	0.0000	1.0000
SYMI	41	-1063.9293	2209.8586	0.0000	1.0000
SYMG	41	-1063.9456	2209.8912	0.0000	1.0000
JC69I	36	-1069.4745	2210.9490	0.0000	1.0000
JC69G	36	-1069.5440	2211.0879	0.0000	1.0000
SYMIG	42	-1063.9705	2211.9411	0.0000	1.0000
JC69IG	37	-1069.5992	2213.1984	0.0000	1.0000
K2P	36	-1071.3700	2214.7400	0.0000	1.0000
SYM	40	-1070.5862	2221.1723	0.0000	1.0000
JC69	35	-1076.1940	2222.3880	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-1039.8430	2323.6091	0.4642	0.4642
HKYI	40	-1039.8434	2323.6099	0.4640	0.9282
HKYIG	41	-1039.8586	2329.7383	0.0217	0.9499
F81I	39	-1046.1588	2330.1425	0.0177	0.9676
F81G	39	-1046.2031	2330.2311	0.0169	0.9845
HKY	39	-1046.3874	2330.5998	0.0141	0.9986
F81IG	40	-1046.2408	2336.4045	0.0008	0.9994
F81	38	-1052.7215	2337.1699	0.0005	0.9999
GTRG	44	-1037.0462	2342.4077	0.0000	1.0000
GTRI	44	-1037.1734	2342.6620	0.0000	1.0000
GTRIG	45	-1036.9918	2348.3969	0.0000	1.0000
GTR	43	-1043.8513	2349.9198	0.0000	1.0000
K2PI	37	-1064.5937	2354.8161	0.0000	1.0000
K2PG	37	-1064.6513	2354.9314	0.0000	1.0000
JC69I	36	-1069.4745	2358.4797	0.0000	1.0000
JC69G	36	-1069.5440	2358.6186	0.0000	1.0000
K2PIG	38	-1064.6512	2361.0293	0.0000	1.0000
K2P	36	-1071.3700	2362.2707	0.0000	1.0000
JC69IG	37	-1069.5992	2364.8271	0.0000	1.0000
JC69	35	-1076.1940	2365.8206	0.0000	1.0000
SYMI	41	-1063.9293	2377.8797	0.0000	1.0000
SYMG	41	-1063.9456	2377.9122	0.0000	1.0000
SYMIG	42	-1063.9705	2384.0602	0.0000	1.0000
SYM	40	-1070.5862	2385.0953	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_10849.phylip = 1 - 445;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_10849.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_10849.phylip = 1 - 445;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_10849.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_10849.phylip = 1 - 445;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_10849.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

