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

Input data from file "683-tomodel-phylip/chr1_24644.phylip" (7 taxa, 1150 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	-2546.8584	5125.7169	0.2478	0.2478
GTRG	20	-2543.0177	5126.0354	0.2114	0.4592
HKYI	16	-2547.1373	5126.2745	0.1875	0.6467
GTRI	20	-2543.3296	5126.6591	0.1547	0.8014
HKYIG	17	-2546.6921	5127.3843	0.1077	0.9091
GTRIG	21	-2542.8616	5127.7231	0.0909	1.0000
GTR	19	-2560.7920	5159.5839	0.0000	1.0000
HKY	15	-2565.6161	5161.2322	0.0000	1.0000
F81G	15	-2598.2676	5226.5353	0.0000	1.0000
F81I	15	-2598.5003	5227.0007	0.0000	1.0000
F81IG	16	-2598.1260	5228.2519	0.0000	1.0000
F81	14	-2615.6702	5259.3404	0.0000	1.0000
K2PG	13	-2618.0492	5262.0985	0.0000	1.0000
K2PI	13	-2618.1703	5262.3405	0.0000	1.0000
SYMG	17	-2614.7399	5263.4799	0.0000	1.0000
SYMI	17	-2614.9240	5263.8480	0.0000	1.0000
K2PIG	14	-2617.9903	5263.9806	0.0000	1.0000
SYMIG	18	-2614.6168	5265.2336	0.0000	1.0000
K2P	12	-2638.2178	5300.4356	0.0000	1.0000
SYM	16	-2634.3915	5300.7831	0.0000	1.0000
JC69G	12	-2662.8115	5349.6231	0.0000	1.0000
JC69I	12	-2662.9842	5349.9684	0.0000	1.0000
JC69IG	13	-2662.7157	5351.4315	0.0000	1.0000
JC69	11	-2681.1277	5384.2553	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2546.8584	5206.4771	0.5582	0.5582
HKYI	16	-2547.1373	5207.0348	0.4224	0.9805
HKYIG	17	-2546.6921	5213.1921	0.0194	1.0000
GTRG	20	-2543.0177	5226.9857	0.0000	1.0000
GTRI	20	-2543.3296	5227.6095	0.0000	1.0000
GTRIG	21	-2542.8616	5233.7210	0.0000	1.0000
HKY	15	-2565.6161	5236.9449	0.0000	1.0000
GTR	19	-2560.7920	5255.4868	0.0000	1.0000
F81G	15	-2598.2676	5302.2481	0.0000	1.0000
F81I	15	-2598.5003	5302.7134	0.0000	1.0000
F81IG	16	-2598.1260	5309.0122	0.0000	1.0000
K2PG	13	-2618.0492	5327.7162	0.0000	1.0000
K2PI	13	-2618.1703	5327.9583	0.0000	1.0000
F81	14	-2615.6702	5330.0057	0.0000	1.0000
K2PIG	14	-2617.9903	5334.6459	0.0000	1.0000
SYMG	17	-2614.7399	5349.2877	0.0000	1.0000
SYMI	17	-2614.9240	5349.6558	0.0000	1.0000
SYMIG	18	-2614.6168	5356.0889	0.0000	1.0000
K2P	12	-2638.2178	5361.0058	0.0000	1.0000
SYM	16	-2634.3915	5381.5433	0.0000	1.0000
JC69G	12	-2662.8115	5410.1933	0.0000	1.0000
JC69I	12	-2662.9842	5410.5386	0.0000	1.0000
JC69IG	13	-2662.7157	5417.0492	0.0000	1.0000
JC69	11	-2681.1277	5439.7780	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/chr1_24644.phylip = 1 - 1150;
 Partition Dummy = 1:683-tomodel-phylip/chr1_24644.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/chr1_24644.phylip = 1 - 1150;
 Partition Dummy = 1:683-tomodel-phylip/chr1_24644.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/chr1_24644.phylip = 1 - 1150;
 Partition Dummy = 1:683-tomodel-phylip/chr1_24644.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

