
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_8655.phylip" (19 taxa, 507 characters)

Minimum AIC  model: HKYI
Minimum AICc model: HKYI
Minimum BIC  model: HKYI


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
HKYI	40	-1062.1210	2211.2806	0.4592	0.4592
HKYIG	41	-1061.3603	2212.1270	0.3008	0.7600
HKYG	40	-1063.5446	2214.1277	0.1106	0.8706
GTRG	44	-1058.8569	2214.2852	0.1022	0.9728
GTRI	44	-1060.6684	2217.9082	0.0167	0.9895
GTRIG	45	-1059.9293	2218.8391	0.0105	1.0000
F81G	39	-1074.0949	2232.8707	0.0000	1.0000
K2PI	37	-1077.6538	2235.3034	0.0000	1.0000
K2PIG	38	-1076.9221	2236.1776	0.0000	1.0000
F81I	39	-1076.4945	2237.6699	0.0000	1.0000
F81IG	40	-1075.5951	2238.2288	0.0000	1.0000
K2PG	37	-1079.1398	2238.2753	0.0000	1.0000
SYMG	41	-1076.2074	2241.8213	0.0000	1.0000
HKY	39	-1078.7811	2242.2431	0.0000	1.0000
GTR	43	-1074.4168	2243.0064	0.0000	1.0000
SYMI	41	-1077.7666	2244.9396	0.0000	1.0000
SYMIG	42	-1076.9892	2245.7629	0.0000	1.0000
JC69G	36	-1089.0407	2255.7495	0.0000	1.0000
JC69I	36	-1091.3630	2260.3941	0.0000	1.0000
F81	38	-1089.0739	2260.4812	0.0000	1.0000
JC69IG	37	-1090.4453	2260.8864	0.0000	1.0000
K2P	36	-1094.3516	2266.3712	0.0000	1.0000
SYM	40	-1091.4849	2270.0085	0.0000	1.0000
JC69	35	-1104.0633	2283.4769	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	40	-1062.1210	2204.2420	0.3780	0.3780
HKYIG	41	-1061.3603	2204.7206	0.2975	0.6755
GTRG	44	-1058.8569	2205.7138	0.1811	0.8566
HKYG	40	-1063.5446	2207.0891	0.0910	0.9476
GTRI	44	-1060.6684	2209.3368	0.0296	0.9772
GTRIG	45	-1059.9293	2209.8587	0.0228	1.0000
F81G	39	-1074.0949	2226.1897	0.0000	1.0000
K2PI	37	-1077.6538	2229.3077	0.0000	1.0000
K2PIG	38	-1076.9221	2229.8443	0.0000	1.0000
F81I	39	-1076.4945	2230.9890	0.0000	1.0000
F81IG	40	-1075.5951	2231.1902	0.0000	1.0000
K2PG	37	-1079.1398	2232.2796	0.0000	1.0000
SYMG	41	-1076.2074	2234.4148	0.0000	1.0000
GTR	43	-1074.4168	2234.8336	0.0000	1.0000
HKY	39	-1078.7811	2235.5621	0.0000	1.0000
SYMI	41	-1077.7666	2237.5331	0.0000	1.0000
SYMIG	42	-1076.9892	2237.9784	0.0000	1.0000
JC69G	36	-1089.0407	2250.0814	0.0000	1.0000
F81	38	-1089.0739	2254.1478	0.0000	1.0000
JC69I	36	-1091.3630	2254.7260	0.0000	1.0000
JC69IG	37	-1090.4453	2254.8906	0.0000	1.0000
K2P	36	-1094.3516	2260.7031	0.0000	1.0000
SYM	40	-1091.4849	2262.9698	0.0000	1.0000
JC69	35	-1104.0633	2278.1265	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	40	-1062.1210	2373.3824	0.7469	0.7469
HKYG	40	-1063.5446	2376.2296	0.1799	0.9268
HKYIG	41	-1061.3603	2378.0895	0.0710	0.9978
K2PI	37	-1077.6538	2385.7626	0.0015	0.9993
K2PG	37	-1079.1398	2388.7345	0.0003	0.9997
K2PIG	38	-1076.9221	2390.5277	0.0001	0.9998
F81G	39	-1074.0949	2391.1017	0.0001	0.9999
GTRG	44	-1058.8569	2391.7683	0.0001	1.0000
GTRI	44	-1060.6684	2395.3913	0.0000	1.0000
F81I	39	-1076.4945	2395.9009	0.0000	1.0000
GTRIG	45	-1059.9293	2400.1417	0.0000	1.0000
F81IG	40	-1075.5951	2400.3306	0.0000	1.0000
HKY	39	-1078.7811	2400.4740	0.0000	1.0000
JC69G	36	-1089.0407	2402.3078	0.0000	1.0000
JC69I	36	-1091.3630	2406.9524	0.0000	1.0000
SYMG	41	-1076.2074	2407.7838	0.0000	1.0000
SYMI	41	-1077.7666	2410.9021	0.0000	1.0000
JC69IG	37	-1090.4453	2411.3455	0.0000	1.0000
K2P	36	-1094.3516	2412.9295	0.0000	1.0000
F81	38	-1089.0739	2414.8313	0.0000	1.0000
SYMIG	42	-1076.9892	2415.5759	0.0000	1.0000
GTR	43	-1074.4168	2416.6596	0.0000	1.0000
JC69	35	-1104.0633	2426.1244	0.0000	1.0000
SYM	40	-1091.4849	2432.1103	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (HKYI)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_8655.phylip = 1 - 507;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_8655.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best AICc model (HKYI)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_8655.phylip = 1 - 507;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_8655.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_8655.phylip = 1 - 507;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_8655.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

