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

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

Minimum AIC  model: GTRI
Minimum AICc model: HKYI
Minimum BIC  model: K2PI


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-4187.7237	8415.4474	0.3175	0.3175
HKYI	16	-4191.8020	8415.6041	0.2936	0.6111
GTRIG	21	-4187.7016	8417.4033	0.1194	0.7306
HKYIG	17	-4191.8026	8417.6052	0.1080	0.8385
GTRG	20	-4189.3158	8418.6317	0.0646	0.9031
HKYG	16	-4193.6599	8419.3197	0.0458	0.9489
K2PI	13	-4197.2596	8420.5192	0.0251	0.9741
K2PIG	14	-4197.2585	8422.5169	0.0093	0.9833
SYMI	17	-4194.4041	8422.8082	0.0080	0.9913
K2PG	13	-4199.0346	8424.0693	0.0043	0.9956
SYMIG	18	-4194.3954	8424.7909	0.0030	0.9986
SYMG	17	-4196.1336	8426.2672	0.0014	1.0000
GTR	19	-4221.1492	8480.2984	0.0000	1.0000
HKY	15	-4225.3998	8480.7996	0.0000	1.0000
K2P	12	-4230.3529	8484.7057	0.0000	1.0000
SYM	16	-4228.0560	8488.1120	0.0000	1.0000
F81I	15	-4275.8610	8581.7220	0.0000	1.0000
F81IG	16	-4275.8914	8583.7827	0.0000	1.0000
F81G	15	-4278.3835	8586.7669	0.0000	1.0000
JC69I	12	-4283.6251	8591.2502	0.0000	1.0000
JC69IG	13	-4283.6569	8593.3138	0.0000	1.0000
JC69G	12	-4288.2905	8600.5810	0.0000	1.0000
F81	14	-4304.6986	8637.3972	0.0000	1.0000
JC69	11	-4311.7979	8645.5959	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PI	13	-4197.2596	8485.5708	0.8283	0.8283
K2PG	13	-4199.0346	8489.1209	0.1404	0.9687
K2PIG	14	-4197.2585	8492.5726	0.0250	0.9937
HKYI	16	-4191.8020	8495.6677	0.0053	0.9990
HKYG	16	-4193.6599	8499.3833	0.0008	0.9998
HKYIG	17	-4191.8026	8502.6728	0.0002	1.0000
SYMI	17	-4194.4041	8507.8757	0.0000	1.0000
SYMG	17	-4196.1336	8511.3348	0.0000	1.0000
SYMIG	18	-4194.3954	8514.8624	0.0000	1.0000
GTRI	20	-4187.7237	8515.5269	0.0000	1.0000
GTRG	20	-4189.3158	8518.7111	0.0000	1.0000
GTRIG	21	-4187.7016	8522.4867	0.0000	1.0000
K2P	12	-4230.3529	8544.7534	0.0000	1.0000
HKY	15	-4225.3998	8555.8593	0.0000	1.0000
SYM	16	-4228.0560	8568.1756	0.0000	1.0000
GTR	19	-4221.1492	8575.3739	0.0000	1.0000
JC69I	12	-4283.6251	8651.2978	0.0000	1.0000
F81I	15	-4275.8610	8656.7816	0.0000	1.0000
JC69IG	13	-4283.6569	8658.3655	0.0000	1.0000
JC69G	12	-4288.2905	8660.6286	0.0000	1.0000
F81G	15	-4278.3835	8661.8265	0.0000	1.0000
F81IG	16	-4275.8914	8663.8463	0.0000	1.0000
JC69	11	-4311.7979	8700.6396	0.0000	1.0000
F81	14	-4304.6986	8707.4528	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr9_6291.phylip = 1 - 1101;
 Partition Dummy = 1:591-tomodel-phylip/chr9_6291.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr9_6291.phylip = 1 - 1101;
 Partition Dummy = 1:591-tomodel-phylip/chr9_6291.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 (K2PI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr9_6291.phylip = 1 - 1101;
 Partition Dummy = 1:591-tomodel-phylip/chr9_6291.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

