
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:19 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/chr11_6465.phylip" (19 taxa, 364 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	-1090.2031	2270.5610	0.5645	0.5645
HKYI	40	-1091.0268	2272.2084	0.2477	0.8122
HKYIG	41	-1090.1415	2272.9787	0.1685	0.9808
GTRG	44	-1089.0419	2278.4976	0.0107	0.9914
GTRI	44	-1089.7205	2279.8547	0.0054	0.9968
GTRIG	45	-1088.9767	2280.9722	0.0031	0.9999
HKY	39	-1100.7040	2289.0376	0.0001	1.0000
GTR	43	-1099.1547	2296.1343	0.0000	1.0000
K2PG	37	-1108.8958	2300.4173	0.0000	1.0000
K2PI	37	-1109.9829	2302.5916	0.0000	1.0000
K2PIG	38	-1108.7766	2302.6732	0.0000	1.0000
SYMG	41	-1108.0725	2308.8407	0.0000	1.0000
SYMI	41	-1108.9725	2310.6407	0.0000	1.0000
SYMIG	42	-1107.9667	2311.1857	0.0000	1.0000
K2P	36	-1118.7786	2317.7040	0.0000	1.0000
SYM	40	-1117.4620	2325.0788	0.0000	1.0000
F81G	39	-1123.5852	2334.8001	0.0000	1.0000
F81I	39	-1124.1395	2335.9087	0.0000	1.0000
F81IG	40	-1123.5762	2337.3073	0.0000	1.0000
F81	38	-1132.1411	2349.4021	0.0000	1.0000
JC69G	36	-1138.9738	2358.0944	0.0000	1.0000
JC69I	36	-1139.5284	2359.2036	0.0000	1.0000
JC69IG	37	-1138.9678	2360.5613	0.0000	1.0000
JC69	35	-1147.4446	2372.5721	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-1090.2031	2260.4062	0.5151	0.5151
HKYI	40	-1091.0268	2262.0536	0.2260	0.7412
HKYIG	41	-1090.1415	2262.2831	0.2015	0.9427
GTRG	44	-1089.0419	2266.0838	0.0301	0.9728
GTRI	44	-1089.7205	2267.4409	0.0153	0.9881
GTRIG	45	-1088.9767	2267.9534	0.0118	1.0000
HKY	39	-1100.7040	2279.4080	0.0000	1.0000
GTR	43	-1099.1547	2284.3093	0.0000	1.0000
K2PG	37	-1108.8958	2291.7916	0.0000	1.0000
K2PIG	38	-1108.7766	2293.5532	0.0000	1.0000
K2PI	37	-1109.9829	2293.9659	0.0000	1.0000
SYMG	41	-1108.0725	2298.1451	0.0000	1.0000
SYMIG	42	-1107.9667	2299.9334	0.0000	1.0000
SYMI	41	-1108.9725	2299.9451	0.0000	1.0000
K2P	36	-1118.7786	2309.5572	0.0000	1.0000
SYM	40	-1117.4620	2314.9240	0.0000	1.0000
F81G	39	-1123.5852	2325.1704	0.0000	1.0000
F81I	39	-1124.1395	2326.2791	0.0000	1.0000
F81IG	40	-1123.5762	2327.1525	0.0000	1.0000
F81	38	-1132.1411	2340.2821	0.0000	1.0000
JC69G	36	-1138.9738	2349.9477	0.0000	1.0000
JC69I	36	-1139.5284	2351.0568	0.0000	1.0000
JC69IG	37	-1138.9678	2351.9355	0.0000	1.0000
JC69	35	-1147.4446	2364.8892	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-1090.2031	2416.2924	0.6688	0.6688
HKYI	40	-1091.0268	2417.9397	0.2935	0.9623
HKYIG	41	-1090.1415	2422.0664	0.0373	0.9996
HKY	39	-1100.7040	2431.3970	0.0004	0.9999
K2PG	37	-1108.8958	2435.9863	0.0000	1.0000
GTRG	44	-1089.0419	2437.5586	0.0000	1.0000
K2PI	37	-1109.9829	2438.1606	0.0000	1.0000
GTRI	44	-1089.7205	2438.9157	0.0000	1.0000
K2PIG	38	-1108.7766	2441.6451	0.0000	1.0000
GTRIG	45	-1088.9767	2443.3253	0.0000	1.0000
K2P	36	-1118.7786	2449.8548	0.0000	1.0000
GTR	43	-1099.1547	2451.8869	0.0000	1.0000
SYMG	41	-1108.0725	2457.9284	0.0000	1.0000
SYMI	41	-1108.9725	2459.7284	0.0000	1.0000
SYMIG	42	-1107.9667	2463.6138	0.0000	1.0000
SYM	40	-1117.4620	2470.8102	0.0000	1.0000
F81G	39	-1123.5852	2477.1594	0.0000	1.0000
F81I	39	-1124.1395	2478.2681	0.0000	1.0000
F81IG	40	-1123.5762	2483.0386	0.0000	1.0000
F81	38	-1132.1411	2488.3740	0.0000	1.0000
JC69G	36	-1138.9738	2490.2452	0.0000	1.0000
JC69I	36	-1139.5284	2491.3543	0.0000	1.0000
JC69IG	37	-1138.9678	2496.1302	0.0000	1.0000
JC69	35	-1147.4446	2501.2896	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/chr11_6465.phylip = 1 - 364;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr11_6465.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/chr11_6465.phylip = 1 - 364;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr11_6465.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/chr11_6465.phylip = 1 - 364;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr11_6465.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

