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

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

Minimum AIC  model: GTRIG
Minimum AICc model: GTRIG
Minimum BIC  model: HKYG


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	21	-2456.7942	4955.5883	0.3946	0.3946
GTRG	20	-2458.2470	4956.4940	0.2509	0.6456
HKYIG	17	-2461.7463	4957.4926	0.1523	0.7979
HKYG	16	-2463.1157	4958.2314	0.1053	0.9031
GTRI	20	-2459.6143	4959.2285	0.0639	0.9671
HKYI	16	-2464.2778	4960.5555	0.0329	1.0000
K2PIG	14	-2491.4172	5010.8345	0.0000	1.0000
K2PG	13	-2492.4630	5010.9261	0.0000	1.0000
K2PI	13	-2493.5364	5013.0728	0.0000	1.0000
SYMG	17	-2489.5910	5013.1820	0.0000	1.0000
SYMIG	18	-2488.6643	5013.3285	0.0000	1.0000
SYMI	17	-2490.7286	5015.4573	0.0000	1.0000
F81G	15	-2494.8701	5019.7402	0.0000	1.0000
F81IG	16	-2494.0027	5020.0053	0.0000	1.0000
F81I	15	-2495.8965	5021.7929	0.0000	1.0000
GTR	19	-2497.6198	5033.2396	0.0000	1.0000
HKY	15	-2502.7980	5035.5960	0.0000	1.0000
JC69G	12	-2519.7010	5063.4021	0.0000	1.0000
JC69IG	13	-2518.9362	5063.8725	0.0000	1.0000
JC69I	12	-2520.6529	5065.3058	0.0000	1.0000
K2P	12	-2529.1231	5082.2463	0.0000	1.0000
SYM	16	-2525.1646	5082.3292	0.0000	1.0000
F81	14	-2531.7171	5091.4343	0.0000	1.0000
JC69	11	-2555.4675	5132.9350	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2463.1157	5035.6630	0.6936	0.6936
HKYI	16	-2464.2778	5037.9871	0.2170	0.9106
HKYIG	17	-2461.7463	5039.7637	0.0893	0.9999
GTRG	20	-2458.2470	5053.2835	0.0001	1.0000
GTRI	20	-2459.6143	5056.0180	0.0000	1.0000
GTRIG	21	-2456.7942	5057.2173	0.0000	1.0000
K2PG	13	-2492.4630	5073.8393	0.0000	1.0000
K2PI	13	-2493.5364	5075.9860	0.0000	1.0000
K2PIG	14	-2491.4172	5078.5872	0.0000	1.0000
F81G	15	-2494.8701	5092.3324	0.0000	1.0000
F81I	15	-2495.8965	5094.3851	0.0000	1.0000
SYMG	17	-2489.5910	5095.4531	0.0000	1.0000
F81IG	16	-2494.0027	5097.4369	0.0000	1.0000
SYMI	17	-2490.7286	5097.7284	0.0000	1.0000
SYMIG	18	-2488.6643	5100.4391	0.0000	1.0000
HKY	15	-2502.7980	5108.1881	0.0000	1.0000
JC69G	12	-2519.7010	5121.4758	0.0000	1.0000
JC69I	12	-2520.6529	5123.3796	0.0000	1.0000
GTR	19	-2497.6198	5125.1897	0.0000	1.0000
JC69IG	13	-2518.9362	5126.7857	0.0000	1.0000
K2P	12	-2529.1231	5140.3200	0.0000	1.0000
F81	14	-2531.7171	5159.1870	0.0000	1.0000
SYM	16	-2525.1646	5159.7608	0.0000	1.0000
JC69	11	-2555.4675	5186.1692	0.0000	1.0000



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

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

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


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

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


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chrZ_6570.phylip = 1 - 934;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_6570.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

