
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
16:17 04/04/2011

Input data from file "683-tomodel-phylip/chr7_9558.phylip" (7 taxa, 1116 characters)

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


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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-2515.9055	5063.8110	0.4186	0.4186
HKYG	16	-2516.1019	5064.2039	0.3439	0.7625
HKYIG	17	-2516.1020	5066.2040	0.1265	0.8890
GTRI	20	-2513.9571	5067.9141	0.0538	0.9428
GTRG	20	-2514.2088	5068.4176	0.0418	0.9846
GTRIG	21	-2514.2088	5070.4175	0.0154	1.0000
HKY	15	-2532.1445	5094.2890	0.0000	1.0000
GTR	19	-2529.7973	5097.5946	0.0000	1.0000
F81I	15	-2536.7950	5103.5900	0.0000	1.0000
F81G	15	-2537.0256	5104.0513	0.0000	1.0000
F81IG	16	-2537.0257	5106.0513	0.0000	1.0000
F81	14	-2553.0580	5134.1160	0.0000	1.0000
K2PI	13	-2595.8473	5217.6946	0.0000	1.0000
K2PG	13	-2595.9987	5217.9973	0.0000	1.0000
SYMI	17	-2592.7571	5219.5141	0.0000	1.0000
SYMG	17	-2592.9007	5219.8014	0.0000	1.0000
K2PIG	14	-2595.9987	5219.9974	0.0000	1.0000
SYMIG	18	-2592.9007	5221.8014	0.0000	1.0000
JC69I	12	-2611.8225	5247.6451	0.0000	1.0000
K2P	12	-2611.9135	5247.8270	0.0000	1.0000
JC69G	12	-2612.0093	5248.0187	0.0000	1.0000
SYM	16	-2608.6814	5249.3629	0.0000	1.0000
JC69IG	13	-2612.0094	5250.0187	0.0000	1.0000
JC69	11	-2628.0802	5278.1603	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-2515.9055	5144.0911	0.5416	0.5416
HKYG	16	-2516.1019	5144.4840	0.4450	0.9867
HKYIG	17	-2516.1020	5151.5016	0.0133	1.0000
GTRI	20	-2513.9571	5168.2643	0.0000	1.0000
GTRG	20	-2514.2088	5168.7677	0.0000	1.0000
HKY	15	-2532.1445	5169.5516	0.0000	1.0000
GTRIG	21	-2514.2088	5175.7851	0.0000	1.0000
F81I	15	-2536.7950	5178.8526	0.0000	1.0000
F81G	15	-2537.0256	5179.3139	0.0000	1.0000
F81IG	16	-2537.0257	5186.3314	0.0000	1.0000
GTR	19	-2529.7973	5192.9272	0.0000	1.0000
F81	14	-2553.0580	5204.3610	0.0000	1.0000
K2PI	13	-2595.8473	5282.9222	0.0000	1.0000
K2PG	13	-2595.9987	5283.2249	0.0000	1.0000
K2PIG	14	-2595.9987	5290.2425	0.0000	1.0000
SYMI	17	-2592.7571	5304.8117	0.0000	1.0000
SYMG	17	-2592.9007	5305.0990	0.0000	1.0000
JC69I	12	-2611.8225	5307.8551	0.0000	1.0000
K2P	12	-2611.9135	5308.0370	0.0000	1.0000
JC69G	12	-2612.0093	5308.2288	0.0000	1.0000
SYMIG	18	-2592.9007	5312.1165	0.0000	1.0000
JC69IG	13	-2612.0094	5315.2463	0.0000	1.0000
SYM	16	-2608.6814	5329.6430	0.0000	1.0000
JC69	11	-2628.0802	5333.3529	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 683-tomodel-phylip/chr7_9558.phylip = 1 - 1116;
 Partition Dummy = 1:683-tomodel-phylip/chr7_9558.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 683-tomodel-phylip/chr7_9558.phylip = 1 - 1116;
 Partition Dummy = 1:683-tomodel-phylip/chr7_9558.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 683-tomodel-phylip/chr7_9558.phylip = 1 - 1116;
 Partition Dummy = 1:683-tomodel-phylip/chr7_9558.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

