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

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

Minimum AIC  model: HKYI
Minimum AICc model: HKYI
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
HKYI	16	-1306.9030	2646.5044	0.2697	0.2697
HKYG	16	-1306.9054	2646.5092	0.2691	0.5388
HKY	15	-1308.0767	2646.7688	0.2363	0.7751
HKYIG	17	-1306.9142	2648.6150	0.0939	0.8689
GTRI	20	-1304.6205	2650.3248	0.0399	0.9089
GTRG	20	-1304.6228	2650.3296	0.0398	0.9487
GTR	19	-1305.7366	2650.4525	0.0375	0.9862
GTRIG	21	-1304.6287	2652.4512	0.0138	1.0000
F81I	15	-1317.8624	2666.3401	0.0000	1.0000
F81G	15	-1317.8636	2666.3425	0.0000	1.0000
F81	14	-1319.0016	2666.5411	0.0000	1.0000
F81IG	16	-1317.8698	2668.4380	0.0000	1.0000
K2PI	13	-1339.5070	2705.4795	0.0000	1.0000
K2PG	13	-1339.5078	2705.4811	0.0000	1.0000
K2P	12	-1340.7361	2705.8707	0.0000	1.0000
K2PIG	14	-1339.5181	2707.5739	0.0000	1.0000
SYMI	17	-1338.1472	2711.0810	0.0000	1.0000
SYMG	17	-1338.1496	2711.0859	0.0000	1.0000
SYM	16	-1339.2847	2711.2678	0.0000	1.0000
SYMIG	18	-1338.1576	2713.1955	0.0000	1.0000
JC69I	12	-1349.3973	2723.1932	0.0000	1.0000
JC69G	12	-1349.3985	2723.1956	0.0000	1.0000
JC69	11	-1350.5775	2723.4918	0.0000	1.0000
JC69IG	13	-1349.4054	2725.2763	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-1306.9030	2645.8060	0.2641	0.2641
HKYG	16	-1306.9054	2645.8109	0.2635	0.5276
HKY	15	-1308.0767	2646.1534	0.2220	0.7496
HKYIG	17	-1306.9142	2647.8284	0.0961	0.8457
GTRI	20	-1304.6205	2649.2409	0.0474	0.8931
GTRG	20	-1304.6228	2649.2457	0.0473	0.9404
GTR	19	-1305.7366	2649.4731	0.0422	0.9827
GTRIG	21	-1304.6287	2651.2574	0.0173	1.0000
F81I	15	-1317.8624	2665.7247	0.0000	1.0000
F81G	15	-1317.8636	2665.7272	0.0000	1.0000
F81	14	-1319.0016	2666.0033	0.0000	1.0000
F81IG	16	-1317.8698	2667.7396	0.0000	1.0000
K2PI	13	-1339.5070	2705.0140	0.0000	1.0000
K2PG	13	-1339.5078	2705.0156	0.0000	1.0000
K2P	12	-1340.7361	2705.4722	0.0000	1.0000
K2PIG	14	-1339.5181	2707.0362	0.0000	1.0000
SYMI	17	-1338.1472	2710.2944	0.0000	1.0000
SYMG	17	-1338.1496	2710.2993	0.0000	1.0000
SYM	16	-1339.2847	2710.5695	0.0000	1.0000
SYMIG	18	-1338.1576	2712.3152	0.0000	1.0000
JC69I	12	-1349.3973	2722.7947	0.0000	1.0000
JC69G	12	-1349.3985	2722.7971	0.0000	1.0000
JC69	11	-1350.5775	2723.1550	0.0000	1.0000
JC69IG	13	-1349.4054	2724.8108	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKY	15	-1308.0767	2716.3474	0.8106	0.8106
HKYI	16	-1306.9030	2720.6796	0.0929	0.9035
HKYG	16	-1306.9054	2720.6845	0.0927	0.9962
HKYIG	17	-1306.9142	2727.3815	0.0033	0.9995
F81	14	-1319.0016	2731.5177	0.0004	0.9999
F81I	15	-1317.8624	2735.9187	0.0000	0.9999
F81G	15	-1317.8636	2735.9211	0.0000	1.0000
GTR	19	-1305.7366	2738.3855	0.0000	1.0000
F81IG	16	-1317.8698	2742.6132	0.0000	1.0000
GTRI	20	-1304.6205	2742.8329	0.0000	1.0000
GTRG	20	-1304.6228	2742.8377	0.0000	1.0000
GTRIG	21	-1304.6287	2749.5290	0.0000	1.0000
K2P	12	-1340.7361	2761.6274	0.0000	1.0000
K2PI	13	-1339.5070	2765.8488	0.0000	1.0000
K2PG	13	-1339.5078	2765.8504	0.0000	1.0000
K2PIG	14	-1339.5181	2772.5505	0.0000	1.0000
JC69	11	-1350.5775	2774.6306	0.0000	1.0000
JC69I	12	-1349.3973	2778.9499	0.0000	1.0000
JC69G	12	-1349.3985	2778.9523	0.0000	1.0000
SYM	16	-1339.2847	2785.4430	0.0000	1.0000
JC69IG	13	-1349.4054	2785.6456	0.0000	1.0000
SYMI	17	-1338.1472	2789.8476	0.0000	1.0000
SYMG	17	-1338.1496	2789.8524	0.0000	1.0000
SYMIG	18	-1338.1576	2796.5480	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/chrZ_11372.phylip = 1 - 796;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_11372.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/chrZ_11372.phylip = 1 - 796;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_11372.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 (HKY)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chrZ_11372.phylip = 1 - 796;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_11372.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

