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

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

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


Note: number of parameters compared to sample size
is low (Nchar/Nparams < 40). Consider the use of AICc.


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

Model	df	lnL		AICc		wAICc	cwAICc
HKYG	16	-2711.2439	5455.2045	0.5609	0.5609
HKYIG	17	-2711.2497	5457.3069	0.1960	0.7569
HKYI	16	-2712.5999	5457.9165	0.1445	0.9015
GTRG	20	-2709.2164	5459.5454	0.0640	0.9655
GTRIG	21	-2709.2167	5461.6589	0.0222	0.9877
GTRI	20	-2710.8699	5462.8524	0.0123	1.0000
HKY	15	-2729.7479	5490.1273	0.0000	1.0000
GTR	19	-2727.9843	5494.9739	0.0000	1.0000
K2PG	13	-2749.1731	5524.8239	0.0000	1.0000
K2PI	13	-2749.9010	5526.2796	0.0000	1.0000
K2PIG	14	-2749.1532	5526.8584	0.0000	1.0000
SYMG	17	-2748.7036	5532.2146	0.0000	1.0000
SYMI	17	-2749.3641	5533.5355	0.0000	1.0000
SYMIG	18	-2748.6792	5534.2620	0.0000	1.0000
K2P	12	-2767.9200	5560.2488	0.0000	1.0000
SYM	16	-2767.7016	5568.1199	0.0000	1.0000
F81G	15	-2780.1774	5590.9864	0.0000	1.0000
F81I	15	-2780.5154	5591.6625	0.0000	1.0000
F81IG	16	-2780.0723	5592.8614	0.0000	1.0000
F81	14	-2796.5012	5621.5544	0.0000	1.0000
JC69G	12	-2811.9413	5648.2915	0.0000	1.0000
JC69I	12	-2812.0996	5648.6082	0.0000	1.0000
JC69IG	13	-2811.7728	5650.0233	0.0000	1.0000
JC69	11	-2828.5228	5679.3911	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2711.2439	5454.4877	0.5434	0.5434
HKYIG	17	-2711.2497	5456.4995	0.1987	0.7421
HKYI	16	-2712.5999	5457.1997	0.1400	0.8822
GTRG	20	-2709.2164	5458.4329	0.0756	0.9577
GTRIG	21	-2709.2167	5460.4334	0.0278	0.9855
GTRI	20	-2710.8699	5461.7398	0.0145	1.0000
HKY	15	-2729.7479	5489.4957	0.0000	1.0000
GTR	19	-2727.9843	5493.9686	0.0000	1.0000
K2PG	13	-2749.1731	5524.3462	0.0000	1.0000
K2PI	13	-2749.9010	5525.8019	0.0000	1.0000
K2PIG	14	-2749.1532	5526.3065	0.0000	1.0000
SYMG	17	-2748.7036	5531.4072	0.0000	1.0000
SYMI	17	-2749.3641	5532.7281	0.0000	1.0000
SYMIG	18	-2748.6792	5533.3584	0.0000	1.0000
K2P	12	-2767.9200	5559.8399	0.0000	1.0000
SYM	16	-2767.7016	5567.4032	0.0000	1.0000
F81G	15	-2780.1774	5590.3548	0.0000	1.0000
F81I	15	-2780.5154	5591.0309	0.0000	1.0000
F81IG	16	-2780.0723	5592.1446	0.0000	1.0000
F81	14	-2796.5012	5621.0025	0.0000	1.0000
JC69G	12	-2811.9413	5647.8826	0.0000	1.0000
JC69I	12	-2812.0996	5648.1992	0.0000	1.0000
JC69IG	13	-2811.7728	5649.5456	0.0000	1.0000
JC69	11	-2828.5228	5679.0455	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2711.2439	5528.9542	0.7732	0.7732
HKYI	16	-2712.5999	5531.6662	0.1992	0.9724
HKYIG	17	-2711.2497	5535.6201	0.0276	1.0000
GTRG	20	-2709.2164	5551.5159	0.0000	1.0000
GTRI	20	-2710.8699	5554.8229	0.0000	1.0000
GTRIG	21	-2709.2167	5558.1706	0.0000	1.0000
HKY	15	-2729.7479	5559.3080	0.0000	1.0000
GTR	19	-2727.9843	5582.3975	0.0000	1.0000
K2PG	13	-2749.1731	5584.8501	0.0000	1.0000
K2PI	13	-2749.9010	5586.3059	0.0000	1.0000
K2PIG	14	-2749.1532	5591.4646	0.0000	1.0000
SYMG	17	-2748.7036	5610.5278	0.0000	1.0000
SYMI	17	-2749.3641	5611.8487	0.0000	1.0000
K2P	12	-2767.9200	5615.6898	0.0000	1.0000
SYMIG	18	-2748.6792	5617.1332	0.0000	1.0000
SYM	16	-2767.7016	5641.8696	0.0000	1.0000
F81G	15	-2780.1774	5660.1671	0.0000	1.0000
F81I	15	-2780.5154	5660.8432	0.0000	1.0000
F81IG	16	-2780.0723	5666.6111	0.0000	1.0000
F81	14	-2796.5012	5686.1606	0.0000	1.0000
JC69G	12	-2811.9413	5703.7324	0.0000	1.0000
JC69I	12	-2812.0996	5704.0491	0.0000	1.0000
JC69IG	13	-2811.7728	5710.0496	0.0000	1.0000
JC69	11	-2828.5228	5730.2412	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 591-tomodel-phylip/chrZ_3148.phylip = 1 - 776;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_3148.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 591-tomodel-phylip/chrZ_3148.phylip = 1 - 776;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_3148.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 591-tomodel-phylip/chrZ_3148.phylip = 1 - 776;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_3148.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

