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

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

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
Minimum BIC  model: HKY


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
GTRI	20	-1492.0315	3025.6600	0.3020	0.3020
GTRG	20	-1492.4179	3026.4328	0.2052	0.5072
HKYI	16	-1497.1187	3027.2637	0.1354	0.6427
HKYG	16	-1497.1588	3027.3440	0.1301	0.7728
GTRIG	21	-1492.0411	3027.8422	0.1014	0.8742
GTR	19	-1494.9086	3029.2594	0.0499	0.9242
HKYIG	17	-1497.0865	3029.3299	0.0482	0.9724
HKY	15	-1499.9055	3030.7150	0.0241	0.9965
F81I	15	-1502.7989	3036.5017	0.0013	0.9978
F81G	15	-1502.8362	3036.5763	0.0013	0.9991
F81IG	16	-1502.8053	3038.6371	0.0005	0.9996
F81	14	-1505.5903	3039.9701	0.0002	0.9998
K2PI	13	-1508.0213	3042.7256	0.0001	0.9999
K2PG	13	-1508.2652	3043.2133	0.0000	0.9999
SYMI	17	-1504.4459	3044.0487	0.0000	0.9999
SYMG	17	-1504.7998	3044.7564	0.0000	1.0000
K2PIG	14	-1508.2275	3045.2444	0.0000	1.0000
SYMIG	18	-1504.4552	3046.2058	0.0000	1.0000
K2P	12	-1510.8764	3046.3371	0.0000	1.0000
SYM	16	-1507.1834	3047.3932	0.0000	1.0000
JC69I	12	-1512.8349	3050.2541	0.0000	1.0000
JC69G	12	-1513.1339	3050.8521	0.0000	1.0000
JC69IG	13	-1512.8432	3052.3693	0.0000	1.0000
JC69	11	-1515.6581	3053.8096	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-1492.0315	3024.0630	0.3278	0.3278
GTRG	20	-1492.4179	3024.8358	0.2227	0.5505
GTRIG	21	-1492.0411	3026.0822	0.1194	0.6699
HKYI	16	-1497.1187	3026.2373	0.1105	0.7805
HKYG	16	-1497.1588	3026.3176	0.1062	0.8866
GTR	19	-1494.9086	3027.8173	0.0502	0.9368
HKYIG	17	-1497.0865	3028.1730	0.0420	0.9788
HKY	15	-1499.9055	3029.8111	0.0185	0.9973
F81I	15	-1502.7989	3035.5977	0.0010	0.9983
F81G	15	-1502.8362	3035.6724	0.0010	0.9993
F81IG	16	-1502.8053	3037.6106	0.0004	0.9997
F81	14	-1505.5903	3039.1806	0.0002	0.9998
K2PI	13	-1508.0213	3042.0426	0.0000	0.9999
K2PG	13	-1508.2652	3042.5304	0.0000	0.9999
SYMI	17	-1504.4459	3042.8918	0.0000	0.9999
SYMG	17	-1504.7998	3043.5995	0.0000	1.0000
K2PIG	14	-1508.2275	3044.4549	0.0000	1.0000
SYMIG	18	-1504.4552	3044.9104	0.0000	1.0000
K2P	12	-1510.8764	3045.7529	0.0000	1.0000
SYM	16	-1507.1834	3046.3667	0.0000	1.0000
JC69I	12	-1512.8349	3049.6699	0.0000	1.0000
JC69G	12	-1513.1339	3050.2678	0.0000	1.0000
JC69IG	13	-1512.8432	3051.6864	0.0000	1.0000
JC69	11	-1515.6581	3053.3161	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKY	15	-1499.9055	3094.3778	0.3126	0.3126
HKYI	16	-1497.1187	3095.1085	0.2170	0.5296
HKYG	16	-1497.1588	3095.1888	0.2084	0.7380
K2P	12	-1510.8764	3097.4063	0.0688	0.8068
K2PI	13	-1508.0213	3098.0005	0.0511	0.8579
K2PG	13	-1508.2652	3098.4882	0.0400	0.8979
F81	14	-1505.5903	3099.4429	0.0248	0.9227
F81I	15	-1502.7989	3100.1645	0.0173	0.9401
F81G	15	-1502.8362	3100.2391	0.0167	0.9567
JC69	11	-1515.6581	3100.6650	0.0135	0.9702
JC69I	12	-1512.8349	3101.3232	0.0097	0.9799
HKYIG	17	-1497.0865	3101.3487	0.0096	0.9895
JC69G	12	-1513.1339	3101.9212	0.0072	0.9967
K2PIG	14	-1508.2275	3104.7172	0.0018	0.9985
F81IG	16	-1502.8053	3106.4818	0.0007	0.9992
JC69IG	13	-1512.8432	3107.6442	0.0004	0.9996
GTR	19	-1494.9086	3109.6018	0.0002	0.9998
GTRI	20	-1492.0315	3110.1520	0.0001	0.9999
GTRG	20	-1492.4179	3110.9248	0.0001	1.0000
SYM	16	-1507.1834	3115.2379	0.0000	1.0000
SYMI	17	-1504.4459	3116.0674	0.0000	1.0000
GTRIG	21	-1492.0411	3116.4756	0.0000	1.0000
SYMG	17	-1504.7998	3116.7752	0.0000	1.0000
SYMIG	18	-1504.4552	3122.3905	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr5_1649.phylip = 1 - 547;
 Partition Dummy = 1:591-tomodel-phylip/chr5_1649.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr5_1649.phylip = 1 - 547;
 Partition Dummy = 1:591-tomodel-phylip/chr5_1649.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr5_1649.phylip = 1 - 547;
 Partition Dummy = 1:591-tomodel-phylip/chr5_1649.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0);
END;


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

