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

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

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


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-4206.0038	8452.0076	0.4372	0.4372
GTRG	20	-4206.1580	8452.3159	0.3747	0.8119
GTRIG	21	-4205.8485	8453.6970	0.1878	0.9997
HKYI	16	-4217.9969	8467.9937	0.0001	0.9998
HKYG	16	-4218.4646	8468.9292	0.0001	0.9999
HKYIG	17	-4217.9195	8469.8390	0.0001	1.0000
GTR	19	-4227.4491	8492.8982	0.0000	1.0000
HKY	15	-4240.3576	8510.7152	0.0000	1.0000
SYMI	17	-4240.4209	8514.8417	0.0000	1.0000
SYMG	17	-4240.7711	8515.5422	0.0000	1.0000
SYMIG	18	-4240.4298	8516.8595	0.0000	1.0000
K2PI	13	-4250.4620	8526.9240	0.0000	1.0000
K2PG	13	-4250.8433	8527.6866	0.0000	1.0000
K2PIG	14	-4250.3552	8528.7104	0.0000	1.0000
SYM	16	-4259.8050	8551.6100	0.0000	1.0000
K2P	12	-4270.3465	8564.6929	0.0000	1.0000
F81I	15	-4369.7539	8769.5079	0.0000	1.0000
F81G	15	-4370.3066	8770.6131	0.0000	1.0000
F81IG	16	-4369.7293	8771.4586	0.0000	1.0000
F81	14	-4386.2852	8800.5704	0.0000	1.0000
JC69I	12	-4393.2100	8810.4201	0.0000	1.0000
JC69G	12	-4393.7805	8811.5610	0.0000	1.0000
JC69IG	13	-4393.1862	8812.3723	0.0000	1.0000
JC69	11	-4409.4932	8840.9864	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-4217.9969	8548.8925	0.5304	0.5304
HKYG	16	-4218.4646	8549.8280	0.3323	0.8627
GTRI	20	-4206.0038	8553.1311	0.0637	0.9264
GTRG	20	-4206.1580	8553.4394	0.0546	0.9810
HKYIG	17	-4217.9195	8555.7940	0.0168	0.9978
GTRIG	21	-4205.8485	8559.8767	0.0022	1.0000
HKY	15	-4240.3576	8586.5578	0.0000	1.0000
GTR	19	-4227.4491	8588.9656	0.0000	1.0000
K2PI	13	-4250.4620	8592.6542	0.0000	1.0000
K2PG	13	-4250.8433	8593.4168	0.0000	1.0000
K2PIG	14	-4250.3552	8599.4969	0.0000	1.0000
SYMI	17	-4240.4209	8600.7967	0.0000	1.0000
SYMG	17	-4240.7711	8601.4972	0.0000	1.0000
SYMIG	18	-4240.4298	8607.8707	0.0000	1.0000
K2P	12	-4270.3465	8625.3670	0.0000	1.0000
SYM	16	-4259.8050	8632.5088	0.0000	1.0000
F81I	15	-4369.7539	8845.3505	0.0000	1.0000
F81G	15	-4370.3066	8846.4557	0.0000	1.0000
F81IG	16	-4369.7293	8852.3574	0.0000	1.0000
JC69I	12	-4393.2100	8871.0942	0.0000	1.0000
F81	14	-4386.2852	8871.3568	0.0000	1.0000
JC69G	12	-4393.7805	8872.2351	0.0000	1.0000
JC69IG	13	-4393.1862	8878.1026	0.0000	1.0000
JC69	11	-4409.4932	8896.6043	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 683-tomodel-phylip/chr1_22454.phylip = 1 - 1160;
 Partition Dummy = 1:683-tomodel-phylip/chr1_22454.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 (GTRI)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr1_22454.phylip = 1 - 1160;
 Partition Dummy = 1:683-tomodel-phylip/chr1_22454.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 BIC model (HKYI)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr1_22454.phylip = 1 - 1160;
 Partition Dummy = 1:683-tomodel-phylip/chr1_22454.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

