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

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

Minimum AIC  model: SYMG
Minimum AICc model: SYMG
Minimum BIC  model: K2PG


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
SYMG	17	-1256.6201	2548.7626	0.2841	0.2841
K2PG	13	-1261.5470	2549.9906	0.1538	0.4379
SYMIG	18	-1256.4475	2550.6007	0.1133	0.5512
SYMI	17	-1257.7615	2551.0454	0.0907	0.6419
K2PI	13	-1262.2440	2551.3846	0.0766	0.7185
HKYG	16	-1259.1080	2551.5659	0.0699	0.7885
K2PIG	14	-1261.4079	2551.8528	0.0606	0.8491
GTRG	20	-1255.0204	2552.1460	0.0523	0.9014
HKYI	16	-1259.8377	2553.0253	0.0337	0.9351
HKYIG	17	-1258.9687	2553.4597	0.0271	0.9623
GTRIG	21	-1254.8595	2554.0406	0.0203	0.9826
GTRI	20	-1256.1190	2554.3433	0.0174	1.0000
SYM	16	-1272.4998	2578.3496	0.0000	1.0000
JC69G	12	-1276.8543	2578.4751	0.0000	1.0000
F81G	15	-1273.7497	2578.6875	0.0000	1.0000
JC69I	12	-1277.3546	2579.4757	0.0000	1.0000
F81I	15	-1274.2258	2579.6396	0.0000	1.0000
JC69IG	13	-1276.7899	2580.4763	0.0000	1.0000
F81IG	16	-1273.6784	2580.7067	0.0000	1.0000
K2P	12	-1278.1544	2581.0753	0.0000	1.0000
GTR	19	-1270.5892	2581.0784	0.0000	1.0000
HKY	15	-1275.4268	2582.0418	0.0000	1.0000
JC69	11	-1292.7393	2608.1257	0.0000	1.0000
F81	14	-1289.8423	2608.7216	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
SYMG	17	-1256.6201	2547.2403	0.2957	0.2957
SYMIG	18	-1256.4475	2548.8949	0.1293	0.4249
K2PG	13	-1261.5470	2549.0940	0.1170	0.5419
SYMI	17	-1257.7615	2549.5230	0.0944	0.6364
GTRG	20	-1255.0204	2550.0408	0.0729	0.7093
HKYG	16	-1259.1080	2550.2161	0.0668	0.7760
K2PI	13	-1262.2440	2550.4880	0.0583	0.8343
K2PIG	14	-1261.4079	2550.8158	0.0495	0.8838
HKYI	16	-1259.8377	2551.6754	0.0322	0.9160
GTRIG	21	-1254.8595	2551.7190	0.0315	0.9475
HKYIG	17	-1258.9687	2551.9373	0.0282	0.9757
GTRI	20	-1256.1190	2552.2380	0.0243	1.0000
SYM	16	-1272.4998	2576.9997	0.0000	1.0000
F81G	15	-1273.7497	2577.4994	0.0000	1.0000
JC69G	12	-1276.8543	2577.7085	0.0000	1.0000
F81I	15	-1274.2258	2578.4515	0.0000	1.0000
JC69I	12	-1277.3546	2578.7091	0.0000	1.0000
GTR	19	-1270.5892	2579.1784	0.0000	1.0000
F81IG	16	-1273.6784	2579.3569	0.0000	1.0000
JC69IG	13	-1276.7899	2579.5798	0.0000	1.0000
K2P	12	-1278.1544	2580.3087	0.0000	1.0000
HKY	15	-1275.4268	2580.8537	0.0000	1.0000
JC69	11	-1292.7393	2607.4786	0.0000	1.0000
F81	14	-1289.8423	2607.6845	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-1261.5470	2601.6173	0.6421	0.6421
K2PI	13	-1262.2440	2603.0113	0.3198	0.9620
K2PIG	14	-1261.4079	2607.3793	0.0360	0.9980
HKYG	16	-1259.1080	2614.8601	0.0009	0.9988
SYMG	17	-1256.6201	2615.9246	0.0005	0.9993
HKYI	16	-1259.8377	2616.3195	0.0004	0.9998
SYMI	17	-1257.7615	2618.2073	0.0002	0.9999
HKYIG	17	-1258.9687	2620.6217	0.0000	1.0000
SYMIG	18	-1256.4475	2621.6195	0.0000	1.0000
JC69G	12	-1276.8543	2626.1916	0.0000	1.0000
JC69I	12	-1277.3546	2627.1922	0.0000	1.0000
K2P	12	-1278.1544	2628.7918	0.0000	1.0000
GTRG	20	-1255.0204	2630.8459	0.0000	1.0000
JC69IG	13	-1276.7899	2632.1031	0.0000	1.0000
GTRI	20	-1256.1190	2633.0431	0.0000	1.0000
GTRIG	21	-1254.8595	2636.5643	0.0000	1.0000
F81G	15	-1273.7497	2638.1032	0.0000	1.0000
F81I	15	-1274.2258	2639.0553	0.0000	1.0000
HKY	15	-1275.4268	2641.4575	0.0000	1.0000
SYM	16	-1272.4998	2641.6438	0.0000	1.0000
F81IG	16	-1273.6784	2644.0009	0.0000	1.0000
JC69	11	-1292.7393	2651.9214	0.0000	1.0000
GTR	19	-1270.5892	2655.9433	0.0000	1.0000
F81	14	-1289.8423	2664.2481	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr11_2599.phylip = 1 - 420;
 Partition Dummy = 1:591-tomodel-phylip/chr11_2599.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=Fixed(Equal);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr11_2599.phylip = 1 - 420;
 Partition Dummy = 1:591-tomodel-phylip/chr11_2599.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=Fixed(Equal);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr11_2599.phylip = 1 - 420;
 Partition Dummy = 1:591-tomodel-phylip/chr11_2599.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

