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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	20	-1467.9503	2977.6011	0.3542	0.3542
GTRI	20	-1468.5511	2978.8026	0.1942	0.5484
HKYG	16	-1472.9774	2979.0472	0.1719	0.7203
GTRIG	21	-1467.9485	2979.7712	0.1197	0.8399
HKYI	16	-1473.5075	2980.1074	0.1012	0.9411
HKYIG	17	-1472.9803	2981.1920	0.0588	0.9999
GTR	19	-1477.4451	2994.4256	0.0001	1.0000
HKY	15	-1482.5988	2996.1595	0.0000	1.0000
K2PG	13	-1497.0802	3020.8870	0.0000	1.0000
K2PI	13	-1497.8336	3022.3938	0.0000	1.0000
K2PIG	14	-1497.0697	3022.9794	0.0000	1.0000
SYMG	17	-1495.8953	3027.0220	0.0000	1.0000
F81G	15	-1498.5399	3028.0417	0.0000	1.0000
SYMI	17	-1496.7343	3028.7000	0.0000	1.0000
SYMIG	18	-1495.8823	3029.1436	0.0000	1.0000
F81I	15	-1499.1787	3029.3192	0.0000	1.0000
F81IG	16	-1498.5387	3030.1697	0.0000	1.0000
K2P	12	-1507.4432	3039.5079	0.0000	1.0000
F81	14	-1507.4404	3043.7207	0.0000	1.0000
SYM	16	-1505.5501	3044.1926	0.0000	1.0000
JC69G	12	-1519.4153	3063.4521	0.0000	1.0000
JC69I	12	-1520.1511	3064.9237	0.0000	1.0000
JC69IG	13	-1519.4072	3065.5410	0.0000	1.0000
JC69	11	-1528.8343	3080.1934	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1467.9503	2975.9007	0.3821	0.3821
GTRI	20	-1468.5511	2977.1022	0.2095	0.5916
GTRIG	21	-1467.9485	2977.8970	0.1408	0.7324
HKYG	16	-1472.9774	2977.9549	0.1368	0.8692
HKYI	16	-1473.5075	2979.0150	0.0805	0.9497
HKYIG	17	-1472.9803	2979.9607	0.0502	0.9999
GTR	19	-1477.4451	2992.8902	0.0001	1.0000
HKY	15	-1482.5988	2995.1976	0.0000	1.0000
K2PG	13	-1497.0802	3020.1605	0.0000	1.0000
K2PI	13	-1497.8336	3021.6673	0.0000	1.0000
K2PIG	14	-1497.0697	3022.1394	0.0000	1.0000
SYMG	17	-1495.8953	3025.7906	0.0000	1.0000
F81G	15	-1498.5399	3027.0797	0.0000	1.0000
SYMI	17	-1496.7343	3027.4686	0.0000	1.0000
SYMIG	18	-1495.8823	3027.7646	0.0000	1.0000
F81I	15	-1499.1787	3028.3573	0.0000	1.0000
F81IG	16	-1498.5387	3029.0774	0.0000	1.0000
K2P	12	-1507.4432	3038.8863	0.0000	1.0000
F81	14	-1507.4404	3042.8807	0.0000	1.0000
SYM	16	-1505.5501	3043.1002	0.0000	1.0000
JC69G	12	-1519.4153	3062.8306	0.0000	1.0000
JC69I	12	-1520.1511	3064.3022	0.0000	1.0000
JC69IG	13	-1519.4072	3064.8145	0.0000	1.0000
JC69	11	-1528.8343	3079.6685	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1472.9774	3045.8616	0.6116	0.6116
HKYI	16	-1473.5075	3046.9217	0.3600	0.9716
HKYIG	17	-1472.9803	3052.1115	0.0269	0.9985
HKY	15	-1482.5988	3058.8601	0.0009	0.9994
GTRG	20	-1467.9503	3060.7840	0.0004	0.9998
GTRI	20	-1468.5511	3061.9855	0.0002	1.0000
GTRIG	21	-1467.9485	3067.0245	0.0000	1.0000
GTR	19	-1477.4451	3073.5294	0.0000	1.0000
K2PG	13	-1497.0802	3075.3347	0.0000	1.0000
K2PI	13	-1497.8336	3076.8414	0.0000	1.0000
K2PIG	14	-1497.0697	3081.5578	0.0000	1.0000
K2P	12	-1507.4432	3089.8163	0.0000	1.0000
F81G	15	-1498.5399	3090.7422	0.0000	1.0000
F81I	15	-1499.1787	3092.0198	0.0000	1.0000
F81IG	16	-1498.5387	3096.9841	0.0000	1.0000
SYMG	17	-1495.8953	3097.9415	0.0000	1.0000
SYMI	17	-1496.7343	3099.6195	0.0000	1.0000
F81	14	-1507.4404	3102.2991	0.0000	1.0000
SYMIG	18	-1495.8823	3104.1596	0.0000	1.0000
SYM	16	-1505.5501	3111.0069	0.0000	1.0000
JC69G	12	-1519.4153	3113.7606	0.0000	1.0000
JC69I	12	-1520.1511	3115.2322	0.0000	1.0000
JC69IG	13	-1519.4072	3119.9887	0.0000	1.0000
JC69	11	-1528.8343	3126.3544	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr1_32299.phylip = 1 - 515;
 Partition Dummy = 1:591-tomodel-phylip/chr1_32299.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr1_32299.phylip = 1 - 515;
 Partition Dummy = 1:591-tomodel-phylip/chr1_32299.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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);
END;


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

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

