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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	40	-1118.0351	2324.4162	0.3782	0.3782
HKYI	40	-1118.2206	2324.7873	0.3141	0.6923
GTRG	44	-1114.2279	2326.6357	0.1247	0.8170
HKYIG	41	-1118.0012	2326.7882	0.1155	0.9325
GTRI	44	-1115.2980	2328.7760	0.0428	0.9753
GTRIG	45	-1115.0307	2330.7314	0.0161	0.9913
K2PG	37	-1126.2528	2333.6065	0.0038	0.9952
K2PI	37	-1126.4508	2334.0027	0.0031	0.9983
K2PIG	38	-1126.2021	2335.9079	0.0012	0.9995
SYMG	41	-1124.4641	2339.7138	0.0002	0.9997
SYMI	41	-1124.6679	2340.1215	0.0001	0.9998
HKY	39	-1127.6129	2341.1446	0.0001	0.9999
SYMIG	42	-1124.4182	2342.0742	0.0001	1.0000
GTR	43	-1123.7060	2343.1146	0.0000	1.0000
K2P	36	-1136.3421	2351.3946	0.0000	1.0000
SYM	40	-1133.7077	2355.7615	0.0000	1.0000
F81G	39	-1137.5382	2360.9952	0.0000	1.0000
F81I	39	-1138.1179	2362.1547	0.0000	1.0000
F81IG	40	-1138.1084	2364.5629	0.0000	1.0000
JC69G	36	-1145.4693	2369.6490	0.0000	1.0000
JC69I	36	-1145.5146	2369.7395	0.0000	1.0000
JC69IG	37	-1145.4825	2372.0661	0.0000	1.0000
F81	38	-1146.9721	2377.4480	0.0000	1.0000
JC69	35	-1154.5675	2385.4666	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-1118.0351	2316.0702	0.2883	0.2883
HKYI	40	-1118.2206	2316.4413	0.2395	0.5278
GTRG	44	-1114.2279	2316.4557	0.2378	0.7656
HKYIG	41	-1118.0012	2318.0025	0.1097	0.8753
GTRI	44	-1115.2980	2318.5961	0.0815	0.9569
GTRIG	45	-1115.0307	2320.0613	0.0392	0.9961
K2PG	37	-1126.2528	2326.5055	0.0016	0.9976
K2PI	37	-1126.4508	2326.9017	0.0013	0.9989
K2PIG	38	-1126.2021	2328.4041	0.0006	0.9995
SYMG	41	-1124.4641	2330.9281	0.0002	0.9997
SYMI	41	-1124.6679	2331.3358	0.0001	0.9998
SYMIG	42	-1124.4182	2332.8364	0.0001	0.9999
HKY	39	-1127.6129	2333.2259	0.0001	1.0000
GTR	43	-1123.7060	2333.4120	0.0000	1.0000
K2P	36	-1136.3421	2344.6843	0.0000	1.0000
SYM	40	-1133.7077	2347.4154	0.0000	1.0000
F81G	39	-1137.5382	2353.0764	0.0000	1.0000
F81I	39	-1138.1179	2354.2359	0.0000	1.0000
F81IG	40	-1138.1084	2356.2168	0.0000	1.0000
JC69G	36	-1145.4693	2362.9387	0.0000	1.0000
JC69I	36	-1145.5146	2363.0292	0.0000	1.0000
JC69IG	37	-1145.4825	2364.9651	0.0000	1.0000
F81	38	-1146.9721	2369.9442	0.0000	1.0000
JC69	35	-1154.5675	2379.1349	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	37	-1126.2528	2477.2082	0.3783	0.3783
K2PI	37	-1126.4508	2477.6043	0.3104	0.6887
HKYG	40	-1118.0351	2478.9919	0.1551	0.8438
HKYI	40	-1118.2206	2479.3630	0.1288	0.9726
K2PIG	38	-1126.2021	2483.1798	0.0191	0.9917
HKYIG	41	-1118.0012	2484.9973	0.0077	0.9994
K2P	36	-1136.3421	2491.3139	0.0003	0.9997
HKY	39	-1127.6129	2492.0746	0.0002	0.9999
GTRG	44	-1114.2279	2495.6697	0.0000	1.0000
GTRI	44	-1115.2980	2497.8100	0.0000	1.0000
SYMG	41	-1124.4641	2497.9229	0.0000	1.0000
SYMI	41	-1124.6679	2498.3306	0.0000	1.0000
GTRIG	45	-1115.0307	2503.3483	0.0000	1.0000
SYMIG	42	-1124.4182	2503.9042	0.0000	1.0000
GTR	43	-1123.7060	2508.5530	0.0000	1.0000
JC69G	36	-1145.4693	2509.5683	0.0000	1.0000
JC69I	36	-1145.5146	2509.6588	0.0000	1.0000
SYM	40	-1133.7077	2510.3372	0.0000	1.0000
F81G	39	-1137.5382	2511.9252	0.0000	1.0000
F81I	39	-1138.1179	2513.0846	0.0000	1.0000
JC69IG	37	-1145.4825	2515.6677	0.0000	1.0000
F81IG	40	-1138.1084	2519.1386	0.0000	1.0000
JC69	35	-1154.5675	2521.6915	0.0000	1.0000
F81	38	-1146.9721	2524.7199	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/chr2_22041.phylip = 1 - 434;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_22041.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/chr2_22041.phylip = 1 - 434;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_22041.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 (K2PG)]

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

