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

Input data from file "591-tomodel-phylip/chr13_720.phylip" (7 taxa, 729 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	-1730.6283	3494.0207	0.3659	0.3659
HKYI	16	-1730.7542	3494.2724	0.3227	0.6886
HKYIG	17	-1730.6046	3496.0699	0.1313	0.8199
GTRG	20	-1727.9336	3497.0537	0.0803	0.9002
GTRI	20	-1728.0480	3497.2825	0.0716	0.9719
GTRIG	21	-1727.9222	3499.1513	0.0281	1.0000
HKY	15	-1748.7740	3528.2212	0.0000	1.0000
GTR	19	-1746.0474	3531.1668	0.0000	1.0000
F81G	15	-1763.8708	3558.4149	0.0000	1.0000
F81I	15	-1763.8803	3558.4337	0.0000	1.0000
F81IG	16	-1763.8819	3560.5279	0.0000	1.0000
F81	14	-1781.0648	3590.7178	0.0000	1.0000
K2PG	13	-1783.8460	3594.2011	0.0000	1.0000
K2PI	13	-1783.9478	3594.4047	0.0000	1.0000
K2PIG	14	-1783.8283	3596.2449	0.0000	1.0000
SYMG	17	-1781.0857	3597.0322	0.0000	1.0000
SYMI	17	-1781.1538	3597.1684	0.0000	1.0000
SYMIG	18	-1781.0841	3599.1317	0.0000	1.0000
K2P	12	-1802.4971	3629.4300	0.0000	1.0000
SYM	16	-1800.1030	3632.9700	0.0000	1.0000
JC69G	12	-1811.1335	3646.7028	0.0000	1.0000
JC69I	12	-1811.1544	3646.7446	0.0000	1.0000
JC69IG	13	-1811.1383	3648.7856	0.0000	1.0000
JC69	11	-1828.8561	3680.0803	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-1730.6283	3493.2567	0.3482	0.3482
HKYI	16	-1730.7542	3493.5084	0.3070	0.6551
HKYIG	17	-1730.6046	3495.2091	0.1312	0.7863
GTRG	20	-1727.9336	3495.8672	0.0944	0.8807
GTRI	20	-1728.0480	3496.0961	0.0842	0.9649
GTRIG	21	-1727.9222	3497.8444	0.0351	1.0000
HKY	15	-1748.7740	3527.5480	0.0000	1.0000
GTR	19	-1746.0474	3530.0948	0.0000	1.0000
F81G	15	-1763.8708	3557.7417	0.0000	1.0000
F81I	15	-1763.8803	3557.7605	0.0000	1.0000
F81IG	16	-1763.8819	3559.7639	0.0000	1.0000
F81	14	-1781.0648	3590.1295	0.0000	1.0000
K2PG	13	-1783.8460	3593.6920	0.0000	1.0000
K2PI	13	-1783.9478	3593.8956	0.0000	1.0000
K2PIG	14	-1783.8283	3595.6567	0.0000	1.0000
SYMG	17	-1781.0857	3596.1715	0.0000	1.0000
SYMI	17	-1781.1538	3596.3076	0.0000	1.0000
SYMIG	18	-1781.0841	3598.1683	0.0000	1.0000
K2P	12	-1802.4971	3628.9943	0.0000	1.0000
SYM	16	-1800.1030	3632.2059	0.0000	1.0000
JC69G	12	-1811.1335	3646.2670	0.0000	1.0000
JC69I	12	-1811.1544	3646.3089	0.0000	1.0000
JC69IG	13	-1811.1383	3648.2765	0.0000	1.0000
JC69	11	-1828.8561	3679.7121	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1730.6283	3566.7235	0.5209	0.5209
HKYI	16	-1730.7542	3566.9752	0.4593	0.9802
HKYIG	17	-1730.6046	3573.2676	0.0198	1.0000
GTRG	20	-1727.9336	3587.7007	0.0000	1.0000
GTRI	20	-1728.0480	3587.9295	0.0000	1.0000
GTRIG	21	-1727.9222	3594.2695	0.0000	1.0000
HKY	15	-1748.7740	3596.4231	0.0000	1.0000
GTR	19	-1746.0474	3617.3366	0.0000	1.0000
F81G	15	-1763.8708	3626.6168	0.0000	1.0000
F81I	15	-1763.8803	3626.6356	0.0000	1.0000
F81IG	16	-1763.8819	3633.2306	0.0000	1.0000
K2PG	13	-1783.8460	3653.3837	0.0000	1.0000
K2PI	13	-1783.9478	3653.5874	0.0000	1.0000
F81	14	-1781.0648	3654.4130	0.0000	1.0000
K2PIG	14	-1783.8283	3659.9401	0.0000	1.0000
SYMG	17	-1781.0857	3674.2299	0.0000	1.0000
SYMI	17	-1781.1538	3674.3661	0.0000	1.0000
SYMIG	18	-1781.0841	3680.8184	0.0000	1.0000
K2P	12	-1802.4971	3684.0944	0.0000	1.0000
JC69G	12	-1811.1335	3701.3671	0.0000	1.0000
JC69I	12	-1811.1544	3701.4090	0.0000	1.0000
SYM	16	-1800.1030	3705.6727	0.0000	1.0000
JC69IG	13	-1811.1383	3707.9683	0.0000	1.0000
JC69	11	-1828.8561	3730.2205	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/chr13_720.phylip = 1 - 729;
 Partition Dummy = 1:591-tomodel-phylip/chr13_720.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/chr13_720.phylip = 1 - 729;
 Partition Dummy = 1:591-tomodel-phylip/chr13_720.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/chr13_720.phylip = 1 - 729;
 Partition Dummy = 1:591-tomodel-phylip/chr13_720.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

