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

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

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


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-3492.1531	7024.3062	0.5304	0.5304
GTRIG	21	-3492.1663	7026.3326	0.1925	0.7229
HKYG	16	-3497.4261	7026.8523	0.1485	0.8714
HKYIG	17	-3497.4620	7028.9240	0.0527	0.9241
GTRI	20	-3494.4942	7028.9885	0.0510	0.9751
HKYI	16	-3499.2124	7030.4247	0.0249	1.0000
K2PG	13	-3534.3155	7094.6311	0.0000	1.0000
K2PIG	14	-3534.3185	7096.6369	0.0000	1.0000
SYMG	17	-3531.5495	7097.0991	0.0000	1.0000
SYMIG	18	-3531.5561	7099.1121	0.0000	1.0000
K2PI	13	-3537.1191	7100.2382	0.0000	1.0000
SYMI	17	-3534.5545	7103.1089	0.0000	1.0000
GTR	19	-3540.6187	7119.2374	0.0000	1.0000
HKY	15	-3544.9773	7119.9546	0.0000	1.0000
F81G	15	-3551.9281	7133.8561	0.0000	1.0000
F81IG	16	-3551.9280	7135.8561	0.0000	1.0000
F81I	15	-3555.0139	7140.0278	0.0000	1.0000
K2P	12	-3576.0031	7176.0062	0.0000	1.0000
SYM	16	-3572.5766	7177.1532	0.0000	1.0000
JC69IG	13	-3580.4282	7186.8564	0.0000	1.0000
JC69G	12	-3582.2843	7188.5686	0.0000	1.0000
JC69I	12	-3582.6558	7189.3116	0.0000	1.0000
F81	14	-3594.4504	7216.9007	0.0000	1.0000
JC69	11	-3619.8805	7261.7611	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-3497.4261	7104.5557	0.8339	0.8339
HKYI	16	-3499.2124	7108.1281	0.1398	0.9737
HKYIG	17	-3497.4620	7111.4839	0.0261	0.9998
GTRG	20	-3492.1531	7121.4354	0.0002	1.0000
GTRI	20	-3494.4942	7126.1177	0.0000	1.0000
GTRIG	21	-3492.1663	7128.3183	0.0000	1.0000
K2PG	13	-3534.3155	7157.7651	0.0000	1.0000
K2PI	13	-3537.1191	7163.3722	0.0000	1.0000
K2PIG	14	-3534.3185	7164.6274	0.0000	1.0000
SYMG	17	-3531.5495	7179.6589	0.0000	1.0000
SYMI	17	-3534.5545	7185.6688	0.0000	1.0000
SYMIG	18	-3531.5561	7186.5284	0.0000	1.0000
HKY	15	-3544.9773	7192.8015	0.0000	1.0000
F81G	15	-3551.9281	7206.7031	0.0000	1.0000
GTR	19	-3540.6187	7211.5102	0.0000	1.0000
F81I	15	-3555.0139	7212.8747	0.0000	1.0000
F81IG	16	-3551.9280	7213.5595	0.0000	1.0000
K2P	12	-3576.0031	7234.2838	0.0000	1.0000
JC69G	12	-3582.2843	7246.8461	0.0000	1.0000
JC69I	12	-3582.6558	7247.5891	0.0000	1.0000
JC69IG	13	-3580.4282	7249.9904	0.0000	1.0000
SYM	16	-3572.5766	7254.8566	0.0000	1.0000
F81	14	-3594.4504	7284.8912	0.0000	1.0000
JC69	11	-3619.8805	7315.1822	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/chr4_16780.phylip = 1 - 950;
 Partition Dummy = 1:591-tomodel-phylip/chr4_16780.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/chr4_16780.phylip = 1 - 950;
 Partition Dummy = 1:591-tomodel-phylip/chr4_16780.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 (HKYG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr4_16780.phylip = 1 - 950;
 Partition Dummy = 1:591-tomodel-phylip/chr4_16780.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

