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

Input data from file "683-tomodel-phylip/chr5_3213.phylip" (7 taxa, 1020 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	-2510.2874	5052.5749	0.3422	0.3422
GTRG	20	-2506.3057	5052.6115	0.3360	0.6781
GTRI	20	-2506.7442	5053.4884	0.2167	0.8949
HKYG	16	-2511.9972	5055.9944	0.0619	0.9568
HKYIG	17	-2511.6194	5057.2388	0.0332	0.9900
GTRIG	21	-2508.8182	5059.6363	0.0100	1.0000
K2PG	13	-2531.6733	5089.3465	0.0000	1.0000
K2PI	13	-2532.0000	5090.0000	0.0000	1.0000
SYMG	17	-2528.0069	5090.0137	0.0000	1.0000
SYMI	17	-2528.3615	5090.7230	0.0000	1.0000
K2PIG	14	-2531.4092	5090.8185	0.0000	1.0000
SYMIG	18	-2530.4586	5096.9173	0.0000	1.0000
HKY	15	-2542.7518	5115.5037	0.0000	1.0000
GTR	19	-2539.5580	5117.1160	0.0000	1.0000
F81G	15	-2547.9552	5125.9104	0.0000	1.0000
F81I	15	-2548.3218	5126.6437	0.0000	1.0000
F81IG	16	-2547.8020	5127.6040	0.0000	1.0000
K2P	12	-2563.0095	5150.0190	0.0000	1.0000
SYM	16	-2560.6876	5153.3751	0.0000	1.0000
JC69G	12	-2569.5894	5163.1788	0.0000	1.0000
JC69I	12	-2569.9275	5163.8551	0.0000	1.0000
JC69IG	13	-2569.4602	5164.9205	0.0000	1.0000
F81	14	-2576.3493	5180.6985	0.0000	1.0000
JC69	11	-2597.3342	5216.6684	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-2510.2874	5131.4158	0.8408	0.8408
HKYG	16	-2511.9972	5134.8353	0.1521	0.9930
HKYIG	17	-2511.6194	5141.0073	0.0069	0.9999
GTRG	20	-2506.3057	5151.1626	0.0000	0.9999
GTRI	20	-2506.7442	5152.0396	0.0000	1.0000
K2PG	13	-2531.6733	5153.4048	0.0000	1.0000
K2PI	13	-2532.0000	5154.0583	0.0000	1.0000
K2PIG	14	-2531.4092	5159.8043	0.0000	1.0000
GTRIG	21	-2508.8182	5163.1150	0.0000	1.0000
SYMG	17	-2528.0069	5173.7822	0.0000	1.0000
SYMI	17	-2528.3615	5174.4915	0.0000	1.0000
SYMIG	18	-2530.4586	5185.6133	0.0000	1.0000
HKY	15	-2542.7518	5189.4170	0.0000	1.0000
F81G	15	-2547.9552	5199.8238	0.0000	1.0000
F81I	15	-2548.3218	5200.5570	0.0000	1.0000
F81IG	16	-2547.8020	5206.4449	0.0000	1.0000
K2P	12	-2563.0095	5209.1497	0.0000	1.0000
GTR	19	-2539.5580	5210.7396	0.0000	1.0000
JC69G	12	-2569.5894	5222.3095	0.0000	1.0000
JC69I	12	-2569.9275	5222.9858	0.0000	1.0000
JC69IG	13	-2569.4602	5228.9787	0.0000	1.0000
SYM	16	-2560.6876	5232.2160	0.0000	1.0000
F81	14	-2576.3493	5249.6843	0.0000	1.0000
JC69	11	-2597.3342	5270.8715	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/chr5_3213.phylip = 1 - 1020;
 Partition Dummy = 1:683-tomodel-phylip/chr5_3213.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/chr5_3213.phylip = 1 - 1020;
 Partition Dummy = 1:683-tomodel-phylip/chr5_3213.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/chr5_3213.phylip = 1 - 1020;
 Partition Dummy = 1:683-tomodel-phylip/chr5_3213.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

