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

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

Minimum AIC  model: HKYI
Minimum AICc model: HKYI
Minimum BIC  model: HKYI


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
HKYI	16	-1481.8138	2996.6295	0.4398	0.4398
HKYG	16	-1482.1102	2997.2222	0.3270	0.7668
HKYIG	17	-1481.8217	2998.7726	0.1506	0.9175
HKY	15	-1485.4923	3001.8670	0.0321	0.9495
GTRI	20	-1480.4877	3002.5338	0.0230	0.9725
GTRG	20	-1480.7689	3003.0962	0.0173	0.9898
GTRIG	21	-1480.4955	3004.7086	0.0077	0.9976
GTR	19	-1484.2211	3007.8496	0.0016	0.9992
K2PI	13	-1492.0988	3010.8643	0.0004	0.9995
K2PG	13	-1492.3556	3011.3779	0.0003	0.9998
K2PIG	14	-1492.1056	3012.9819	0.0001	0.9999
K2P	12	-1495.7715	3016.1135	0.0000	1.0000
SYMI	17	-1490.9600	3017.0491	0.0000	1.0000
SYMG	17	-1491.2151	3017.5593	0.0000	1.0000
SYMIG	18	-1490.9670	3019.1983	0.0000	1.0000
SYM	16	-1494.7095	3022.4209	0.0000	1.0000
F81I	15	-1512.7078	3056.2980	0.0000	1.0000
F81G	15	-1512.9291	3056.7406	0.0000	1.0000
F81IG	16	-1512.7138	3058.4295	0.0000	1.0000
F81	14	-1515.6642	3060.0991	0.0000	1.0000
JC69I	12	-1521.7992	3068.1688	0.0000	1.0000
JC69G	12	-1522.0020	3068.5744	0.0000	1.0000
JC69IG	13	-1521.8041	3070.2750	0.0000	1.0000
JC69	11	-1524.9287	3072.3391	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-1481.8138	2995.6276	0.4293	0.4293
HKYG	16	-1482.1102	2996.2204	0.3192	0.7485
HKYIG	17	-1481.8217	2997.6434	0.1567	0.9051
GTRI	20	-1480.4877	3000.9754	0.0296	0.9347
HKY	15	-1485.4923	3000.9846	0.0295	0.9642
GTRG	20	-1480.7689	3001.5378	0.0224	0.9866
GTRIG	21	-1480.4955	3002.9911	0.0108	0.9974
GTR	19	-1484.2211	3006.4422	0.0019	0.9993
K2PI	13	-1492.0988	3010.1977	0.0003	0.9996
K2PG	13	-1492.3556	3010.7112	0.0002	0.9998
K2PIG	14	-1492.1056	3012.2113	0.0001	0.9999
K2P	12	-1495.7715	3015.5431	0.0000	1.0000
SYMI	17	-1490.9600	3015.9200	0.0000	1.0000
SYMG	17	-1491.2151	3016.4301	0.0000	1.0000
SYMIG	18	-1490.9670	3017.9340	0.0000	1.0000
SYM	16	-1494.7095	3021.4191	0.0000	1.0000
F81I	15	-1512.7078	3055.4156	0.0000	1.0000
F81G	15	-1512.9291	3055.8582	0.0000	1.0000
F81IG	16	-1512.7138	3057.4277	0.0000	1.0000
F81	14	-1515.6642	3059.3284	0.0000	1.0000
JC69I	12	-1521.7992	3067.5985	0.0000	1.0000
JC69G	12	-1522.0020	3068.0041	0.0000	1.0000
JC69IG	13	-1521.8041	3069.6083	0.0000	1.0000
JC69	11	-1524.9287	3071.8573	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1481.8138	3064.8746	0.2876	0.2876
HKYG	16	-1482.1102	3065.4673	0.2139	0.5015
HKY	15	-1485.4923	3065.9037	0.1719	0.6734
K2PI	13	-1492.0988	3066.4609	0.1301	0.8036
K2PG	13	-1492.3556	3066.9744	0.1007	0.9042
K2P	12	-1495.7715	3067.4783	0.0782	0.9825
HKYIG	17	-1481.8217	3071.2183	0.0121	0.9945
K2PIG	14	-1492.1056	3072.8024	0.0055	1.0000
GTRI	20	-1480.4877	3087.5341	0.0000	1.0000
GTRG	20	-1480.7689	3088.0965	0.0000	1.0000
GTR	19	-1484.2211	3088.6730	0.0000	1.0000
SYMI	17	-1490.9600	3089.4949	0.0000	1.0000
SYMG	17	-1491.2151	3090.0050	0.0000	1.0000
SYM	16	-1494.7095	3090.6660	0.0000	1.0000
GTRIG	21	-1480.4955	3093.8778	0.0000	1.0000
SYMIG	18	-1490.9670	3095.8368	0.0000	1.0000
JC69	11	-1524.9287	3119.4646	0.0000	1.0000
JC69I	12	-1521.7992	3119.5337	0.0000	1.0000
F81	14	-1515.6642	3119.9196	0.0000	1.0000
JC69G	12	-1522.0020	3119.9393	0.0000	1.0000
F81I	15	-1512.7078	3120.3347	0.0000	1.0000
F81G	15	-1512.9291	3120.7773	0.0000	1.0000
JC69IG	13	-1521.8041	3125.8715	0.0000	1.0000
F81IG	16	-1512.7138	3126.6746	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr4_13330.phylip = 1 - 560;
 Partition Dummy = 1:591-tomodel-phylip/chr4_13330.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr4_13330.phylip = 1 - 560;
 Partition Dummy = 1:591-tomodel-phylip/chr4_13330.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr4_13330.phylip = 1 - 560;
 Partition Dummy = 1:591-tomodel-phylip/chr4_13330.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


-------------------------------------------------------------
End of Output

