
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
16:13 04/04/2011

Input data from file "683-tomodel-phylip/chr2_24766.phylip" (7 taxa, 539 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	-1291.9888	2617.0197	0.4633	0.4633
HKYI	16	-1292.2809	2617.6040	0.3460	0.8093
HKYIG	17	-1291.9268	2619.0282	0.1697	0.9790
GTRG	20	-1291.5832	2624.7880	0.0095	0.9886
GTRI	20	-1291.8389	2625.2995	0.0074	0.9959
GTRIG	21	-1291.5296	2626.8465	0.0034	0.9994
HKY	15	-1299.7253	2630.3683	0.0006	0.9999
K2PG	13	-1305.0577	2636.8087	0.0000	1.0000
K2PI	13	-1305.4176	2637.5286	0.0000	1.0000
GTR	19	-1299.1944	2637.8531	0.0000	1.0000
K2PIG	14	-1304.9894	2638.7803	0.0000	1.0000
SYMG	17	-1303.8252	2642.8251	0.0000	1.0000
SYMI	17	-1304.1416	2643.4579	0.0000	1.0000
SYMIG	18	-1303.7753	2644.8660	0.0000	1.0000
K2P	12	-1312.4250	2649.4431	0.0000	1.0000
F81G	15	-1311.2717	2653.4612	0.0000	1.0000
F81I	15	-1311.5526	2654.0229	0.0000	1.0000
SYM	16	-1310.5433	2654.1287	0.0000	1.0000
F81IG	16	-1311.2008	2655.4437	0.0000	1.0000
F81	14	-1318.7130	2666.2276	0.0000	1.0000
JC69G	12	-1325.9768	2676.5468	0.0000	1.0000
JC69I	12	-1326.2593	2677.1118	0.0000	1.0000
JC69IG	13	-1325.9049	2678.5032	0.0000	1.0000
JC69	11	-1332.7581	2688.0171	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-1291.9888	2615.9776	0.4548	0.4548
HKYI	16	-1292.2809	2616.5618	0.3396	0.7944
HKYIG	17	-1291.9268	2617.8535	0.1780	0.9724
GTRG	20	-1291.5832	2623.1664	0.0125	0.9849
GTRI	20	-1291.8389	2623.6779	0.0097	0.9946
GTRIG	21	-1291.5296	2625.0593	0.0048	0.9994
HKY	15	-1299.7253	2629.4506	0.0005	0.9999
K2PG	13	-1305.0577	2636.1154	0.0000	1.0000
GTR	19	-1299.1944	2636.3887	0.0000	1.0000
K2PI	13	-1305.4176	2636.8353	0.0000	1.0000
K2PIG	14	-1304.9894	2637.9787	0.0000	1.0000
SYMG	17	-1303.8252	2641.6505	0.0000	1.0000
SYMI	17	-1304.1416	2642.2833	0.0000	1.0000
SYMIG	18	-1303.7753	2643.5506	0.0000	1.0000
K2P	12	-1312.4250	2648.8499	0.0000	1.0000
F81G	15	-1311.2717	2652.5434	0.0000	1.0000
SYM	16	-1310.5433	2653.0866	0.0000	1.0000
F81I	15	-1311.5526	2653.1051	0.0000	1.0000
F81IG	16	-1311.2008	2654.4016	0.0000	1.0000
F81	14	-1318.7130	2665.4260	0.0000	1.0000
JC69G	12	-1325.9768	2675.9537	0.0000	1.0000
JC69I	12	-1326.2593	2676.5187	0.0000	1.0000
JC69IG	13	-1325.9049	2677.8098	0.0000	1.0000
JC69	11	-1332.7581	2687.5162	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1291.9888	2684.6130	0.5408	0.5408
HKYI	16	-1292.2809	2685.1973	0.4038	0.9446
HKYIG	17	-1291.9268	2690.7787	0.0248	0.9694
K2PG	13	-1305.0577	2691.8817	0.0143	0.9837
K2PI	13	-1305.4176	2692.6016	0.0100	0.9936
HKY	15	-1299.7253	2693.7963	0.0055	0.9991
K2PIG	14	-1304.9894	2698.0348	0.0007	0.9998
K2P	12	-1312.4250	2700.3265	0.0002	1.0000
GTRG	20	-1291.5832	2708.9607	0.0000	1.0000
GTRI	20	-1291.8389	2709.4722	0.0000	1.0000
SYMG	17	-1303.8252	2714.5756	0.0000	1.0000
GTRIG	21	-1291.5296	2715.1433	0.0000	1.0000
SYMI	17	-1304.1416	2715.2084	0.0000	1.0000
F81G	15	-1311.2717	2716.8891	0.0000	1.0000
F81I	15	-1311.5526	2717.4509	0.0000	1.0000
GTR	19	-1299.1944	2717.8933	0.0000	1.0000
SYMIG	18	-1303.7753	2720.7655	0.0000	1.0000
SYM	16	-1310.5433	2721.7220	0.0000	1.0000
F81IG	16	-1311.2008	2723.0370	0.0000	1.0000
F81	14	-1318.7130	2725.4821	0.0000	1.0000
JC69G	12	-1325.9768	2727.4303	0.0000	1.0000
JC69I	12	-1326.2593	2727.9952	0.0000	1.0000
JC69IG	13	-1325.9049	2733.5761	0.0000	1.0000
JC69	11	-1332.7581	2734.7030	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 683-tomodel-phylip/chr2_24766.phylip = 1 - 539;
 Partition Dummy = 1:683-tomodel-phylip/chr2_24766.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 683-tomodel-phylip/chr2_24766.phylip = 1 - 539;
 Partition Dummy = 1:683-tomodel-phylip/chr2_24766.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 683-tomodel-phylip/chr2_24766.phylip = 1 - 539;
 Partition Dummy = 1:683-tomodel-phylip/chr2_24766.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

