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

Input data from file "683-tomodel-phylip/chr2_23261.phylip" (7 taxa, 584 characters)

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


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	-1296.8416	2626.6426	0.2853	0.2853
HKYG	16	-1296.9535	2626.8664	0.2551	0.5404
GTRI	20	-1293.4584	2628.4087	0.1180	0.6584
GTRG	20	-1293.5520	2628.5959	0.1074	0.7658
HKYIG	17	-1296.9535	2628.9883	0.0883	0.8541
GTRIG	21	-1293.5519	2630.7480	0.0366	0.8908
K2PI	13	-1302.3843	2631.4071	0.0263	0.9171
K2PG	13	-1302.4929	2631.6243	0.0236	0.9408
SYMI	17	-1298.3041	2631.6894	0.0229	0.9636
SYMG	17	-1298.4161	2631.9135	0.0205	0.9841
K2PIG	14	-1302.4929	2633.7239	0.0083	0.9924
SYMIG	18	-1298.4162	2634.0429	0.0071	0.9994
HKY	15	-1304.4580	2639.7611	0.0004	0.9998
GTR	19	-1301.6768	2642.7010	0.0001	0.9999
K2P	12	-1309.5539	2643.6541	0.0001	1.0000
SYM	16	-1306.0236	2645.0066	0.0000	1.0000
F81I	15	-1317.3740	2665.5932	0.0000	1.0000
F81G	15	-1317.4611	2665.7672	0.0000	1.0000
F81IG	16	-1317.4611	2667.8816	0.0000	1.0000
JC69I	12	-1322.6332	2669.8129	0.0000	1.0000
JC69G	12	-1322.7247	2669.9957	0.0000	1.0000
JC69IG	13	-1322.7247	2672.0880	0.0000	1.0000
F81	14	-1324.6001	2677.9383	0.0000	1.0000
JC69	11	-1329.6996	2681.8607	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-1296.8416	2625.6831	0.2631	0.2631
HKYG	16	-1296.9535	2625.9070	0.2352	0.4982
GTRI	20	-1293.4584	2626.9167	0.1420	0.6402
GTRG	20	-1293.5520	2627.1039	0.1293	0.7695
HKYIG	17	-1296.9535	2627.9070	0.0865	0.8560
GTRIG	21	-1293.5519	2629.1039	0.0476	0.9036
SYMI	17	-1298.3041	2630.6082	0.0224	0.9260
K2PI	13	-1302.3843	2630.7685	0.0207	0.9467
SYMG	17	-1298.4161	2630.8323	0.0200	0.9667
K2PG	13	-1302.4929	2630.9857	0.0186	0.9853
SYMIG	18	-1298.4162	2632.8323	0.0074	0.9926
K2PIG	14	-1302.4929	2632.9858	0.0068	0.9995
HKY	15	-1304.4580	2638.9160	0.0004	0.9998
GTR	19	-1301.6768	2641.3535	0.0001	0.9999
K2P	12	-1309.5539	2643.1077	0.0000	1.0000
SYM	16	-1306.0236	2644.0472	0.0000	1.0000
F81I	15	-1317.3740	2664.7481	0.0000	1.0000
F81G	15	-1317.4611	2664.9222	0.0000	1.0000
F81IG	16	-1317.4611	2666.9221	0.0000	1.0000
JC69I	12	-1322.6332	2669.2665	0.0000	1.0000
JC69G	12	-1322.7247	2669.4493	0.0000	1.0000
JC69IG	13	-1322.7247	2671.4494	0.0000	1.0000
F81	14	-1324.6001	2677.2001	0.0000	1.0000
JC69	11	-1329.6996	2681.3992	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PI	13	-1302.3843	2687.5773	0.5029	0.5029
K2PG	13	-1302.4929	2687.7945	0.4512	0.9541
K2PIG	14	-1302.4929	2694.1644	0.0187	0.9728
K2P	12	-1309.5539	2695.5466	0.0094	0.9821
HKYI	16	-1296.8416	2695.6016	0.0091	0.9912
HKYG	16	-1296.9535	2695.8254	0.0081	0.9994
HKYIG	17	-1296.9535	2702.1953	0.0003	0.9997
HKY	15	-1304.4580	2704.4645	0.0001	0.9998
SYMI	17	-1298.3041	2704.8965	0.0001	0.9999
SYMG	17	-1298.4161	2705.1206	0.0001	1.0000
SYMIG	18	-1298.4162	2711.4905	0.0000	1.0000
SYM	16	-1306.0236	2713.9656	0.0000	1.0000
GTRI	20	-1293.4584	2714.3147	0.0000	1.0000
GTRG	20	-1293.5520	2714.5020	0.0000	1.0000
GTRIG	21	-1293.5519	2720.8718	0.0000	1.0000
JC69I	12	-1322.6332	2721.7053	0.0000	1.0000
JC69G	12	-1322.7247	2721.8881	0.0000	1.0000
GTR	19	-1301.6768	2724.3816	0.0000	1.0000
JC69IG	13	-1322.7247	2728.2581	0.0000	1.0000
JC69	11	-1329.6996	2729.4681	0.0000	1.0000
F81I	15	-1317.3740	2730.2966	0.0000	1.0000
F81G	15	-1317.4611	2730.4707	0.0000	1.0000
F81IG	16	-1317.4611	2736.8406	0.0000	1.0000
F81	14	-1324.6001	2738.3787	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 683-tomodel-phylip/chr2_23261.phylip = 1 - 584;
 Partition Dummy = 1:683-tomodel-phylip/chr2_23261.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 683-tomodel-phylip/chr2_23261.phylip = 1 - 584;
 Partition Dummy = 1:683-tomodel-phylip/chr2_23261.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 (K2PI)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr2_23261.phylip = 1 - 584;
 Partition Dummy = 1:683-tomodel-phylip/chr2_23261.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

