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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
Minimum BIC  model: HKYG


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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2978.7608	5989.5215	0.4230	0.4230
HKYI	16	-2979.4253	5990.8506	0.2177	0.6407
HKYIG	17	-2978.7597	5991.5194	0.1558	0.7965
GTRG	20	-2976.0922	5992.1844	0.1117	0.9082
GTRI	20	-2976.8774	5993.7549	0.0509	0.9592
GTRIG	21	-2976.0986	5994.1973	0.0408	1.0000
HKY	15	-3003.0855	6036.1710	0.0000	1.0000
GTR	19	-3000.0181	6038.0363	0.0000	1.0000
F81G	15	-3020.2275	6070.4550	0.0000	1.0000
F81I	15	-3020.7093	6071.4187	0.0000	1.0000
F81IG	16	-3020.2418	6072.4837	0.0000	1.0000
SYMG	17	-3027.6771	6089.3543	0.0000	1.0000
K2PG	13	-3032.2366	6090.4732	0.0000	1.0000
SYMIG	18	-3027.6507	6091.3015	0.0000	1.0000
SYMI	17	-3028.8533	6091.7067	0.0000	1.0000
K2PI	13	-3033.1259	6092.2518	0.0000	1.0000
K2PIG	14	-3032.2247	6092.4495	0.0000	1.0000
F81	14	-3043.4336	6114.8673	0.0000	1.0000
SYM	16	-3050.9377	6133.8754	0.0000	1.0000
K2P	12	-3056.0797	6136.1595	0.0000	1.0000
JC69G	12	-3066.9171	6157.8342	0.0000	1.0000
JC69I	12	-3067.5547	6159.1094	0.0000	1.0000
JC69IG	13	-3066.9224	6159.8448	0.0000	1.0000
JC69	11	-3090.1051	6202.2102	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2978.7608	6068.4563	0.6469	0.6469
HKYI	16	-2979.4253	6069.7854	0.3328	0.9798
HKYIG	17	-2978.7597	6075.3876	0.0202	1.0000
GTRG	20	-2976.0922	6090.8528	0.0000	1.0000
GTRI	20	-2976.8774	6092.4233	0.0000	1.0000
GTRIG	21	-2976.0986	6097.7991	0.0000	1.0000
HKY	15	-3003.0855	6110.1724	0.0000	1.0000
GTR	19	-3000.0181	6131.7713	0.0000	1.0000
F81G	15	-3020.2275	6144.4564	0.0000	1.0000
F81I	15	-3020.7093	6145.4200	0.0000	1.0000
F81IG	16	-3020.2418	6151.4184	0.0000	1.0000
K2PG	13	-3032.2366	6154.6077	0.0000	1.0000
K2PI	13	-3033.1259	6156.3863	0.0000	1.0000
K2PIG	14	-3032.2247	6161.5174	0.0000	1.0000
SYMG	17	-3027.6771	6173.2225	0.0000	1.0000
SYMI	17	-3028.8533	6175.5749	0.0000	1.0000
SYMIG	18	-3027.6507	6180.1031	0.0000	1.0000
F81	14	-3043.4336	6183.9352	0.0000	1.0000
K2P	12	-3056.0797	6195.3606	0.0000	1.0000
SYM	16	-3050.9377	6212.8101	0.0000	1.0000
JC69G	12	-3066.9171	6217.0353	0.0000	1.0000
JC69I	12	-3067.5547	6218.3105	0.0000	1.0000
JC69IG	13	-3066.9224	6223.9793	0.0000	1.0000
JC69	11	-3090.1051	6256.4779	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 683-tomodel-phylip/chr2_6754.phylip = 1 - 1026;
 Partition Dummy = 1:683-tomodel-phylip/chr2_6754.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 683-tomodel-phylip/chr2_6754.phylip = 1 - 1026;
 Partition Dummy = 1:683-tomodel-phylip/chr2_6754.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 683-tomodel-phylip/chr2_6754.phylip = 1 - 1026;
 Partition Dummy = 1:683-tomodel-phylip/chr2_6754.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

