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

Input data from file "591-tomodel-phylip/chrZ_9895.phylip" (7 taxa, 1084 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	-3862.8417	7765.6834	0.5465	0.5465
GTRIG	21	-3862.8583	7767.7167	0.1977	0.7442
HKYI	16	-3867.9991	7767.9981	0.1718	0.9160
HKYIG	17	-3868.0185	7770.0371	0.0620	0.9779
GTRG	20	-3866.3019	7772.6038	0.0172	0.9951
HKYG	16	-3871.5587	7775.1174	0.0049	1.0000
GTR	19	-3904.0257	7846.0514	0.0000	1.0000
HKY	15	-3909.4947	7848.9894	0.0000	1.0000
F81I	15	-3955.5298	7941.0596	0.0000	1.0000
F81IG	16	-3955.5331	7943.0662	0.0000	1.0000
F81G	15	-3957.8197	7945.6394	0.0000	1.0000
SYMI	17	-3961.7225	7957.4449	0.0000	1.0000
SYMIG	18	-3961.5695	7959.1390	0.0000	1.0000
SYMG	17	-3962.6546	7959.3091	0.0000	1.0000
K2PI	13	-3968.5154	7963.0307	0.0000	1.0000
K2PIG	14	-3968.5106	7965.0211	0.0000	1.0000
K2PG	13	-3970.5279	7967.0557	0.0000	1.0000
F81	14	-3990.8356	8009.6712	0.0000	1.0000
SYM	16	-3999.0053	8030.0106	0.0000	1.0000
K2P	12	-4006.1063	8036.2126	0.0000	1.0000
JC69I	12	-4036.5391	8097.0782	0.0000	1.0000
JC69IG	13	-4036.5347	8099.0694	0.0000	1.0000
JC69G	12	-4038.3637	8100.7273	0.0000	1.0000
JC69	11	-4071.4699	8164.9398	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-3867.9991	7847.8128	0.9448	0.9448
HKYIG	17	-3868.0185	7854.8401	0.0281	0.9730
HKYG	16	-3871.5587	7854.9320	0.0269	0.9999
GTRI	20	-3862.8417	7865.4517	0.0001	1.0000
GTRG	20	-3866.3019	7872.3720	0.0000	1.0000
GTRIG	21	-3862.8583	7872.4734	0.0000	1.0000
HKY	15	-3909.4947	7923.8156	0.0000	1.0000
GTR	19	-3904.0257	7940.8313	0.0000	1.0000
F81I	15	-3955.5298	8015.8858	0.0000	1.0000
F81G	15	-3957.8197	8020.4656	0.0000	1.0000
F81IG	16	-3955.5331	8022.8808	0.0000	1.0000
K2PI	13	-3968.5154	8027.8801	0.0000	1.0000
K2PG	13	-3970.5279	8031.9051	0.0000	1.0000
K2PIG	14	-3968.5106	8034.8589	0.0000	1.0000
SYMI	17	-3961.7225	8042.2479	0.0000	1.0000
SYMG	17	-3962.6546	8044.1121	0.0000	1.0000
SYMIG	18	-3961.5695	8048.9304	0.0000	1.0000
F81	14	-3990.8356	8079.5090	0.0000	1.0000
K2P	12	-4006.1063	8096.0736	0.0000	1.0000
SYM	16	-3999.0053	8109.8253	0.0000	1.0000
JC69I	12	-4036.5391	8156.9392	0.0000	1.0000
JC69G	12	-4038.3637	8160.5883	0.0000	1.0000
JC69IG	13	-4036.5347	8163.9188	0.0000	1.0000
JC69	11	-4071.4699	8219.8123	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/chrZ_9895.phylip = 1 - 1084;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_9895.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 591-tomodel-phylip/chrZ_9895.phylip = 1 - 1084;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_9895.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 591-tomodel-phylip/chrZ_9895.phylip = 1 - 1084;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_9895.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

