
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
21:54 04/04/2011

Input data from file "591-tomodel-phylip/chrZ_7827.phylip" (7 taxa, 1137 characters)

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: GTRG


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-3182.1696	6404.3392	0.5426	0.5426
GTRIG	21	-3182.0031	6406.0062	0.2358	0.7784
GTRI	20	-3183.0654	6406.1308	0.2216	1.0000
HKYG	16	-3200.8751	6433.7502	0.0000	1.0000
HKYI	16	-3201.2831	6434.5662	0.0000	1.0000
HKYIG	17	-3200.7200	6435.4399	0.0000	1.0000
GTR	19	-3211.3758	6460.7516	0.0000	1.0000
SYMG	17	-3214.1353	6462.2706	0.0000	1.0000
SYMI	17	-3214.3789	6462.7578	0.0000	1.0000
SYMIG	18	-3213.9909	6463.9819	0.0000	1.0000
K2PG	13	-3226.5494	6479.0988	0.0000	1.0000
K2PI	13	-3226.6654	6479.3308	0.0000	1.0000
K2PIG	14	-3226.3721	6480.7442	0.0000	1.0000
HKY	15	-3229.9295	6489.8591	0.0000	1.0000
SYM	16	-3242.8419	6517.6838	0.0000	1.0000
K2P	12	-3256.4822	6536.9643	0.0000	1.0000
F81G	15	-3261.7149	6553.4297	0.0000	1.0000
F81I	15	-3261.8052	6553.6103	0.0000	1.0000
F81IG	16	-3261.5435	6555.0870	0.0000	1.0000
JC69G	12	-3284.8426	6593.6852	0.0000	1.0000
JC69I	12	-3284.9177	6593.8355	0.0000	1.0000
JC69IG	13	-3284.8261	6595.6522	0.0000	1.0000
F81	14	-3290.0399	6608.0799	0.0000	1.0000
JC69	11	-3313.4989	6648.9978	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	20	-3182.1696	6505.0622	0.6850	0.6850
GTRI	20	-3183.0654	6506.8537	0.2797	0.9647
GTRIG	21	-3182.0031	6511.7653	0.0240	0.9887
HKYG	16	-3200.8751	6514.3286	0.0067	0.9953
HKYI	16	-3201.2831	6515.1446	0.0044	0.9998
HKYIG	17	-3200.7200	6521.0544	0.0002	1.0000
K2PG	13	-3226.5494	6544.5688	0.0000	1.0000
K2PI	13	-3226.6654	6544.8007	0.0000	1.0000
SYMG	17	-3214.1353	6547.8851	0.0000	1.0000
SYMI	17	-3214.3789	6548.3724	0.0000	1.0000
K2PIG	14	-3226.3721	6551.2502	0.0000	1.0000
SYMIG	18	-3213.9909	6554.6326	0.0000	1.0000
GTR	19	-3211.3758	6556.4385	0.0000	1.0000
HKY	15	-3229.9295	6565.4013	0.0000	1.0000
K2P	12	-3256.4822	6597.3981	0.0000	1.0000
SYM	16	-3242.8419	6598.2622	0.0000	1.0000
F81G	15	-3261.7149	6628.9719	0.0000	1.0000
F81I	15	-3261.8052	6629.1525	0.0000	1.0000
F81IG	16	-3261.5435	6635.6654	0.0000	1.0000
JC69G	12	-3284.8426	6654.1190	0.0000	1.0000
JC69I	12	-3284.9177	6654.2693	0.0000	1.0000
JC69IG	13	-3284.8261	6661.1221	0.0000	1.0000
F81	14	-3290.0399	6678.5860	0.0000	1.0000
JC69	11	-3313.4989	6704.3954	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chrZ_7827.phylip = 1 - 1137;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_7827.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chrZ_7827.phylip = 1 - 1137;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_7827.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chrZ_7827.phylip = 1 - 1137;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_7827.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

