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

Input data from file "591-tomodel-phylip/chr5_11247.phylip" (7 taxa, 1203 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	-3942.6183	7925.2367	0.7213	0.7213
GTRIG	21	-3942.6082	7927.2164	0.2681	0.9894
GTRI	20	-3946.8383	7933.6767	0.0106	1.0000
HKYG	16	-3965.1372	7962.2743	0.0000	1.0000
HKYIG	17	-3965.0797	7964.1593	0.0000	1.0000
HKYI	16	-3968.2827	7968.5653	0.0000	1.0000
GTR	19	-3983.0971	8004.1942	0.0000	1.0000
SYMG	17	-4002.5543	8039.1086	0.0000	1.0000
SYMIG	18	-4002.5020	8041.0041	0.0000	1.0000
HKY	15	-4005.5208	8041.0417	0.0000	1.0000
SYMI	17	-4005.1078	8044.2156	0.0000	1.0000
K2PG	13	-4009.8410	8045.6820	0.0000	1.0000
K2PIG	14	-4009.7400	8047.4801	0.0000	1.0000
K2PI	13	-4012.4069	8050.8138	0.0000	1.0000
SYM	16	-4042.8179	8117.6358	0.0000	1.0000
K2P	12	-4052.3784	8128.7568	0.0000	1.0000
F81G	15	-4062.1770	8154.3540	0.0000	1.0000
F81IG	16	-4062.1039	8156.2079	0.0000	1.0000
F81I	15	-4064.3915	8158.7830	0.0000	1.0000
JC69G	12	-4098.0204	8220.0408	0.0000	1.0000
JC69IG	13	-4097.9300	8221.8600	0.0000	1.0000
JC69I	12	-4100.0219	8224.0439	0.0000	1.0000
F81	14	-4098.9700	8225.9400	0.0000	1.0000
JC69	11	-4135.9714	8293.9428	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	20	-3942.6183	8027.0882	0.9578	0.9578
GTRIG	21	-3942.6082	8034.1604	0.0279	0.9857
GTRI	20	-3946.8383	8035.5282	0.0141	0.9998
HKYG	16	-3965.1372	8043.7555	0.0002	1.0000
HKYI	16	-3968.2827	8050.0465	0.0000	1.0000
HKYIG	17	-3965.0797	8050.7331	0.0000	1.0000
GTR	19	-3983.0971	8100.9531	0.0000	1.0000
K2PG	13	-4009.8410	8111.8854	0.0000	1.0000
K2PI	13	-4012.4069	8117.0173	0.0000	1.0000
HKY	15	-4005.5208	8117.4303	0.0000	1.0000
K2PIG	14	-4009.7400	8118.7761	0.0000	1.0000
SYMG	17	-4002.5543	8125.6824	0.0000	1.0000
SYMI	17	-4005.1078	8130.7893	0.0000	1.0000
SYMIG	18	-4002.5020	8132.6704	0.0000	1.0000
K2P	12	-4052.3784	8189.8677	0.0000	1.0000
SYM	16	-4042.8179	8199.1170	0.0000	1.0000
F81G	15	-4062.1770	8230.7426	0.0000	1.0000
F81I	15	-4064.3915	8235.1716	0.0000	1.0000
F81IG	16	-4062.1039	8237.6891	0.0000	1.0000
JC69G	12	-4098.0204	8281.1517	0.0000	1.0000
JC69I	12	-4100.0219	8285.1548	0.0000	1.0000
JC69IG	13	-4097.9300	8288.0634	0.0000	1.0000
F81	14	-4098.9700	8297.2360	0.0000	1.0000
JC69	11	-4135.9714	8349.9611	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/chr5_11247.phylip = 1 - 1203;
 Partition Dummy = 1:591-tomodel-phylip/chr5_11247.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/chr5_11247.phylip = 1 - 1203;
 Partition Dummy = 1:591-tomodel-phylip/chr5_11247.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/chr5_11247.phylip = 1 - 1203;
 Partition Dummy = 1:591-tomodel-phylip/chr5_11247.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

