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

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

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


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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-3071.7527	6175.5055	0.3887	0.3887
HKYI	16	-3072.0227	6176.0455	0.2967	0.6854
HKYIG	17	-3071.7569	6177.5138	0.1424	0.8277
GTRG	20	-3069.2922	6178.5844	0.0834	0.9111
GTRI	20	-3069.6567	6179.3134	0.0579	0.9690
GTRIG	21	-3069.2821	6180.5642	0.0310	1.0000
K2PG	13	-3089.6600	6205.3201	0.0000	1.0000
K2PI	13	-3089.8484	6205.6967	0.0000	1.0000
K2PIG	14	-3089.6739	6207.3477	0.0000	1.0000
SYMG	17	-3086.9237	6207.8473	0.0000	1.0000
SYMI	17	-3087.1545	6208.3089	0.0000	1.0000
HKY	15	-3089.1632	6208.3265	0.0000	1.0000
SYMIG	18	-3086.9342	6209.8685	0.0000	1.0000
GTR	19	-3087.2903	6212.5807	0.0000	1.0000
K2P	12	-3106.5204	6237.0408	0.0000	1.0000
SYM	16	-3104.1025	6240.2051	0.0000	1.0000
F81G	15	-3144.8533	6319.7066	0.0000	1.0000
F81I	15	-3145.0134	6320.0269	0.0000	1.0000
F81IG	16	-3144.8610	6321.7220	0.0000	1.0000
JC69G	12	-3159.8504	6343.7009	0.0000	1.0000
JC69I	12	-3159.9847	6343.9693	0.0000	1.0000
JC69IG	13	-3159.8605	6345.7210	0.0000	1.0000
F81	14	-3159.7418	6347.4837	0.0000	1.0000
JC69	11	-3174.5618	6371.1236	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-3071.7527	6256.7049	0.5575	0.5575
HKYI	16	-3072.0227	6257.2449	0.4256	0.9831
HKYIG	17	-3071.7569	6263.7882	0.0161	0.9993
K2PG	13	-3089.6600	6271.2946	0.0004	0.9997
K2PI	13	-3089.8484	6271.6712	0.0003	1.0000
K2PIG	14	-3089.6739	6278.3972	0.0000	1.0000
GTRG	20	-3069.2922	6280.0836	0.0000	1.0000
GTRI	20	-3069.6567	6280.8126	0.0000	1.0000
HKY	15	-3089.1632	6284.4509	0.0000	1.0000
GTRIG	21	-3069.2821	6287.1385	0.0000	1.0000
SYMG	17	-3086.9237	6294.1217	0.0000	1.0000
SYMI	17	-3087.1545	6294.5833	0.0000	1.0000
K2P	12	-3106.5204	6297.9403	0.0000	1.0000
SYMIG	18	-3086.9342	6301.2178	0.0000	1.0000
GTR	19	-3087.2903	6309.0050	0.0000	1.0000
SYM	16	-3104.1025	6321.4045	0.0000	1.0000
F81G	15	-3144.8533	6395.8310	0.0000	1.0000
F81I	15	-3145.0134	6396.1513	0.0000	1.0000
F81IG	16	-3144.8610	6402.9214	0.0000	1.0000
JC69G	12	-3159.8504	6404.6005	0.0000	1.0000
JC69I	12	-3159.9847	6404.8689	0.0000	1.0000
JC69IG	13	-3159.8605	6411.6955	0.0000	1.0000
F81	14	-3159.7418	6418.5332	0.0000	1.0000
JC69	11	-3174.5618	6426.9482	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr2_8663.phylip = 1 - 1182;
 Partition Dummy = 1:683-tomodel-phylip/chr2_8663.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 AICc model (HKYG)]

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

