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

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

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


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2297.8088	4635.6177	0.2192	0.2192
GTRI	20	-2297.8373	4635.6746	0.2130	0.4322
HKYG	16	-2301.8668	4635.7337	0.2068	0.6390
HKYI	16	-2301.8748	4635.7496	0.2052	0.8442
GTRIG	21	-2297.8100	4637.6200	0.0805	0.9247
HKYIG	17	-2301.8778	4637.7556	0.0753	0.9999
GTR	19	-2307.6877	4653.3755	0.0000	1.0000
HKY	15	-2311.6942	4653.3884	0.0000	1.0000
F81I	15	-2337.8751	4705.7501	0.0000	1.0000
F81G	15	-2337.9001	4705.8003	0.0000	1.0000
F81IG	16	-2337.9202	4707.8404	0.0000	1.0000
F81	14	-2347.7198	4723.4397	0.0000	1.0000
K2PI	13	-2388.3508	4802.7015	0.0000	1.0000
K2PG	13	-2388.3735	4802.7470	0.0000	1.0000
K2PIG	14	-2388.4018	4804.8036	0.0000	1.0000
SYMI	17	-2387.3866	4808.7733	0.0000	1.0000
SYMG	17	-2387.4007	4808.8015	0.0000	1.0000
SYMIG	18	-2387.4233	4810.8465	0.0000	1.0000
K2P	12	-2398.4819	4820.9638	0.0000	1.0000
SYM	16	-2397.3970	4826.7939	0.0000	1.0000
JC69I	12	-2418.2049	4860.4097	0.0000	1.0000
JC69G	12	-2418.2532	4860.5063	0.0000	1.0000
JC69IG	13	-2418.2919	4862.5837	0.0000	1.0000
JC69	11	-2427.6913	4877.3827	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2301.8668	4716.2556	0.4942	0.4942
HKYI	16	-2301.8748	4716.2716	0.4903	0.9845
HKYIG	17	-2301.8778	4723.3102	0.0145	0.9991
HKY	15	-2311.6942	4728.8777	0.0009	1.0000
GTRG	20	-2297.8088	4736.2701	0.0000	1.0000
GTRI	20	-2297.8373	4736.3271	0.0000	1.0000
GTRIG	21	-2297.8100	4743.3051	0.0000	1.0000
GTR	19	-2307.6877	4748.9953	0.0000	1.0000
F81I	15	-2337.8751	4781.2395	0.0000	1.0000
F81G	15	-2337.9001	4781.2896	0.0000	1.0000
F81IG	16	-2337.9202	4788.3624	0.0000	1.0000
F81	14	-2347.7198	4793.8964	0.0000	1.0000
K2PI	13	-2388.3508	4868.1257	0.0000	1.0000
K2PG	13	-2388.3735	4868.1711	0.0000	1.0000
K2PIG	14	-2388.4018	4875.2604	0.0000	1.0000
K2P	12	-2398.4819	4881.3553	0.0000	1.0000
SYMI	17	-2387.3866	4894.3279	0.0000	1.0000
SYMG	17	-2387.4007	4894.3561	0.0000	1.0000
SYMIG	18	-2387.4233	4901.4338	0.0000	1.0000
SYM	16	-2397.3970	4907.3159	0.0000	1.0000
JC69I	12	-2418.2049	4920.8012	0.0000	1.0000
JC69G	12	-2418.2532	4920.8978	0.0000	1.0000
JC69IG	13	-2418.2919	4928.0078	0.0000	1.0000
JC69	11	-2427.6913	4932.7415	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 683-tomodel-phylip/chrZ_12021.phylip = 1 - 1133;
 Partition Dummy = 1:683-tomodel-phylip/chrZ_12021.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 (HKYG)]

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


[Mrbayes block for the best BIC model (HKYG)]

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

