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

Input data from file "591-tomodel-phylip/chrZ_6695.phylip" (7 taxa, 699 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	-2233.1294	4499.0565	0.5884	0.5884
HKYIG	17	-2233.1279	4501.1545	0.2061	0.7945
HKYI	16	-2234.4201	4501.6378	0.1619	0.9563
GTRG	20	-2231.9769	4505.1927	0.0274	0.9837
GTRIG	21	-2231.9772	4507.3191	0.0094	0.9931
GTRI	20	-2233.5210	4508.2810	0.0058	0.9990
K2PG	13	-2243.1461	4512.8236	0.0006	0.9996
K2PIG	14	-2243.1466	4514.9072	0.0002	0.9998
K2PI	13	-2244.5598	4515.6509	0.0001	0.9999
SYMG	17	-2241.6130	4518.1246	0.0000	1.0000
SYMIG	18	-2241.6108	4520.2274	0.0000	1.0000
SYMI	17	-2242.9111	4520.7210	0.0000	1.0000
HKY	15	-2247.0933	4524.8894	0.0000	1.0000
GTR	19	-2245.6218	4530.3628	0.0000	1.0000
K2P	12	-2257.4618	4539.3785	0.0000	1.0000
SYM	16	-2256.1680	4545.1336	0.0000	1.0000
F81G	15	-2283.6269	4597.9567	0.0000	1.0000
F81IG	16	-2283.6313	4600.0603	0.0000	1.0000
F81I	15	-2284.8958	4600.4943	0.0000	1.0000
JC69G	12	-2293.0228	4610.5004	0.0000	1.0000
JC69IG	13	-2293.0236	4612.5786	0.0000	1.0000
JC69I	12	-2294.2458	4612.9464	0.0000	1.0000
F81	14	-2295.8515	4620.3169	0.0000	1.0000
JC69	11	-2305.6573	4633.6989	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2233.1294	4498.2588	0.5758	0.5758
HKYIG	17	-2233.1279	4500.2558	0.2121	0.7879
HKYI	16	-2234.4201	4500.8401	0.1584	0.9463
GTRG	20	-2231.9769	4503.9537	0.0334	0.9797
GTRIG	21	-2231.9772	4505.9543	0.0123	0.9920
GTRI	20	-2233.5210	4507.0420	0.0071	0.9991
K2PG	13	-2243.1461	4512.2922	0.0005	0.9996
K2PIG	14	-2243.1466	4514.2931	0.0002	0.9998
K2PI	13	-2244.5598	4515.1195	0.0001	0.9999
SYMG	17	-2241.6130	4517.2259	0.0000	1.0000
SYMIG	18	-2241.6108	4519.2215	0.0000	1.0000
SYMI	17	-2242.9111	4519.8223	0.0000	1.0000
HKY	15	-2247.0933	4524.1866	0.0000	1.0000
GTR	19	-2245.6218	4529.2435	0.0000	1.0000
K2P	12	-2257.4618	4538.9236	0.0000	1.0000
SYM	16	-2256.1680	4544.3359	0.0000	1.0000
F81G	15	-2283.6269	4597.2539	0.0000	1.0000
F81IG	16	-2283.6313	4599.2626	0.0000	1.0000
F81I	15	-2284.8958	4599.7915	0.0000	1.0000
JC69G	12	-2293.0228	4610.0456	0.0000	1.0000
JC69IG	13	-2293.0236	4612.0473	0.0000	1.0000
JC69I	12	-2294.2458	4612.4916	0.0000	1.0000
F81	14	-2295.8515	4619.7029	0.0000	1.0000
JC69	11	-2305.6573	4633.3146	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2233.1294	4571.0533	0.4219	0.4219
K2PG	13	-2243.1461	4571.4376	0.3481	0.7701
HKYI	16	-2234.4201	4573.6345	0.1161	0.8861
K2PI	13	-2244.5598	4574.2650	0.0847	0.9708
HKYIG	17	-2233.1279	4577.5999	0.0160	0.9868
K2PIG	14	-2243.1466	4577.9882	0.0132	1.0000
HKY	15	-2247.0933	4592.4313	0.0000	1.0000
K2P	12	-2257.4618	4593.5194	0.0000	1.0000
SYMG	17	-2241.6130	4594.5700	0.0000	1.0000
GTRG	20	-2231.9769	4594.9467	0.0000	1.0000
SYMI	17	-2242.9111	4597.1664	0.0000	1.0000
GTRI	20	-2233.5210	4598.0351	0.0000	1.0000
SYMIG	18	-2241.6108	4601.1153	0.0000	1.0000
GTRIG	21	-2231.9772	4601.4970	0.0000	1.0000
GTR	19	-2245.6218	4615.6869	0.0000	1.0000
SYM	16	-2256.1680	4617.1304	0.0000	1.0000
JC69G	12	-2293.0228	4664.6414	0.0000	1.0000
F81G	15	-2283.6269	4665.4986	0.0000	1.0000
JC69I	12	-2294.2458	4667.0874	0.0000	1.0000
F81I	15	-2284.8958	4668.0363	0.0000	1.0000
JC69IG	13	-2293.0236	4671.1927	0.0000	1.0000
F81IG	16	-2283.6313	4672.0570	0.0000	1.0000
JC69	11	-2305.6573	4683.3607	0.0000	1.0000
F81	14	-2295.8515	4683.3980	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_6695.phylip = 1 - 699;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_6695.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 591-tomodel-phylip/chrZ_6695.phylip = 1 - 699;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_6695.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 591-tomodel-phylip/chrZ_6695.phylip = 1 - 699;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_6695.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

