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

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

Minimum AIC  model: HKYG
Minimum AICc model: K2PG
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
K2PG	13	-2217.8935	4462.3769	0.2266	0.2266
HKYG	16	-2214.8123	4462.5106	0.2120	0.4386
GTRG	20	-2210.8752	4463.1274	0.1557	0.5943
SYMG	17	-2214.3437	4463.6858	0.1178	0.7121
K2PIG	14	-2217.8661	4464.4139	0.0818	0.7939
HKYIG	17	-2214.7889	4464.5762	0.0755	0.8693
GTRIG	21	-2210.8489	4465.2150	0.0548	0.9242
SYMIG	18	-2214.3108	4465.7393	0.0422	0.9664
K2PI	13	-2220.7910	4468.1720	0.0125	0.9789
HKYI	16	-2217.7477	4468.3814	0.0113	0.9901
GTRI	20	-2214.1621	4469.7013	0.0058	0.9959
SYMI	17	-2217.7098	4470.4180	0.0041	1.0000
K2P	12	-2244.6949	4513.8946	0.0000	1.0000
HKY	15	-2241.7455	4514.2715	0.0000	1.0000
GTR	19	-2237.6340	4514.5119	0.0000	1.0000
SYM	16	-2240.9333	4514.7526	0.0000	1.0000
F81G	15	-2262.4547	4555.6899	0.0000	1.0000
JC69G	12	-2266.2881	4557.0810	0.0000	1.0000
F81IG	16	-2262.4558	4557.7976	0.0000	1.0000
JC69IG	13	-2266.2874	4559.1647	0.0000	1.0000
F81I	15	-2264.6677	4560.1160	0.0000	1.0000
JC69I	12	-2268.5509	4561.6066	0.0000	1.0000
F81	14	-2286.2370	4601.1558	0.0000	1.0000
JC69	11	-2289.9895	4602.4055	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2214.8123	4461.6246	0.2036	0.2036
GTRG	20	-2210.8752	4461.7504	0.1912	0.3948
K2PG	13	-2217.8935	4461.7869	0.1877	0.5825
SYMG	17	-2214.3437	4462.6874	0.1197	0.7022
HKYIG	17	-2214.7889	4463.5779	0.0767	0.7789
GTRIG	21	-2210.8489	4463.6977	0.0722	0.8511
K2PIG	14	-2217.8661	4463.7321	0.0710	0.9221
SYMIG	18	-2214.3108	4464.6216	0.0455	0.9676
HKYI	16	-2217.7477	4467.4954	0.0108	0.9784
K2PI	13	-2220.7910	4467.5820	0.0104	0.9887
GTRI	20	-2214.1621	4468.3243	0.0071	0.9959
SYMI	17	-2217.7098	4469.4196	0.0041	1.0000
GTR	19	-2237.6340	4513.2681	0.0000	1.0000
K2P	12	-2244.6949	4513.3898	0.0000	1.0000
HKY	15	-2241.7455	4513.4910	0.0000	1.0000
SYM	16	-2240.9333	4513.8666	0.0000	1.0000
F81G	15	-2262.4547	4554.9095	0.0000	1.0000
JC69G	12	-2266.2881	4556.5762	0.0000	1.0000
F81IG	16	-2262.4558	4556.9116	0.0000	1.0000
JC69IG	13	-2266.2874	4558.5748	0.0000	1.0000
F81I	15	-2264.6677	4559.3355	0.0000	1.0000
JC69I	12	-2268.5509	4561.1018	0.0000	1.0000
F81	14	-2286.2370	4600.4740	0.0000	1.0000
JC69	11	-2289.9895	4601.9790	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-2217.8935	4519.6019	0.9110	0.9110
K2PI	13	-2220.7910	4525.3970	0.0503	0.9613
K2PIG	14	-2217.8661	4525.9944	0.0373	0.9985
HKYG	16	-2214.8123	4532.7815	0.0013	0.9998
SYMG	17	-2214.3437	4538.2916	0.0001	0.9999
HKYI	16	-2217.7477	4538.6523	0.0001	0.9999
HKYIG	17	-2214.7889	4539.1821	0.0001	1.0000
SYMIG	18	-2214.3108	4544.6731	0.0000	1.0000
SYMI	17	-2217.7098	4545.0238	0.0000	1.0000
GTRG	20	-2210.8752	4550.6965	0.0000	1.0000
GTRIG	21	-2210.8489	4557.0911	0.0000	1.0000
GTRI	20	-2214.1621	4557.2704	0.0000	1.0000
K2P	12	-2244.6949	4566.7575	0.0000	1.0000
HKY	15	-2241.7455	4580.2006	0.0000	1.0000
SYM	16	-2240.9333	4585.0235	0.0000	1.0000
GTR	19	-2237.6340	4597.7669	0.0000	1.0000
JC69G	12	-2266.2881	4609.9438	0.0000	1.0000
JC69I	12	-2268.5509	4614.4695	0.0000	1.0000
JC69IG	13	-2266.2874	4616.3897	0.0000	1.0000
F81G	15	-2262.4547	4621.6190	0.0000	1.0000
F81I	15	-2264.6677	4626.0451	0.0000	1.0000
F81IG	16	-2262.4558	4628.0685	0.0000	1.0000
JC69	11	-2289.9895	4650.8993	0.0000	1.0000
F81	14	-2286.2370	4662.7363	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 591-tomodel-phylip/chrZ_6615.phylip = 1 - 631;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_6615.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 (K2PG)]

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


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

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

