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

Input data from file "591-tomodel-phylip/chr5_11146.phylip" (7 taxa, 711 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	-2185.7240	4404.2320	0.3353	0.3353
HKYI	16	-2185.7320	4404.2479	0.3326	0.6679
GTRI	20	-2182.4859	4406.1892	0.1260	0.7939
HKYIG	17	-2185.8269	4406.5369	0.1059	0.8998
GTRG	20	-2183.3438	4407.9050	0.0534	0.9532
GTRIG	21	-2182.4145	4408.1700	0.0468	1.0000
HKY	15	-2216.1878	4463.0662	0.0000	1.0000
GTR	19	-2214.7293	4468.5584	0.0000	1.0000
SYMI	17	-2226.4105	4487.7042	0.0000	1.0000
SYMG	17	-2226.6721	4488.2273	0.0000	1.0000
K2PI	13	-2231.8490	4490.2201	0.0000	1.0000
K2PG	13	-2232.0626	4490.6473	0.0000	1.0000
SYMIG	18	-2227.2819	4491.5522	0.0000	1.0000
K2PIG	14	-2232.1515	4492.9065	0.0000	1.0000
F81I	15	-2245.9683	4522.6273	0.0000	1.0000
F81IG	16	-2245.3813	4523.5464	0.0000	1.0000
F81G	15	-2246.5600	4523.8107	0.0000	1.0000
SYM	16	-2258.0563	4548.8964	0.0000	1.0000
K2P	12	-2262.2577	4548.9623	0.0000	1.0000
F81	14	-2273.5480	4575.6995	0.0000	1.0000
JC69I	12	-2284.0253	4592.4976	0.0000	1.0000
JC69G	12	-2284.5787	4593.6045	0.0000	1.0000
JC69IG	13	-2284.0274	4594.5771	0.0000	1.0000
JC69	11	-2312.1092	4646.5960	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2185.7240	4403.4481	0.3151	0.3151
HKYI	16	-2185.7320	4403.4640	0.3126	0.6278
GTRI	20	-2182.4859	4404.9718	0.1471	0.7749
HKYIG	17	-2185.8269	4405.6538	0.1046	0.8795
GTRG	20	-2183.3438	4406.6876	0.0624	0.9419
GTRIG	21	-2182.4145	4406.8289	0.0581	1.0000
HKY	15	-2216.1878	4462.3755	0.0000	1.0000
GTR	19	-2214.7293	4467.4586	0.0000	1.0000
SYMI	17	-2226.4105	4486.8211	0.0000	1.0000
SYMG	17	-2226.6721	4487.3442	0.0000	1.0000
K2PI	13	-2231.8490	4489.6979	0.0000	1.0000
K2PG	13	-2232.0626	4490.1251	0.0000	1.0000
SYMIG	18	-2227.2819	4490.5637	0.0000	1.0000
K2PIG	14	-2232.1515	4492.3031	0.0000	1.0000
F81I	15	-2245.9683	4521.9367	0.0000	1.0000
F81IG	16	-2245.3813	4522.7626	0.0000	1.0000
F81G	15	-2246.5600	4523.1200	0.0000	1.0000
SYM	16	-2258.0563	4548.1126	0.0000	1.0000
K2P	12	-2262.2577	4548.5153	0.0000	1.0000
F81	14	-2273.5480	4575.0961	0.0000	1.0000
JC69I	12	-2284.0253	4592.0506	0.0000	1.0000
JC69G	12	-2284.5787	4593.1575	0.0000	1.0000
JC69IG	13	-2284.0274	4594.0548	0.0000	1.0000
JC69	11	-2312.1092	4646.2183	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2185.7240	4476.5149	0.4936	0.4936
HKYI	16	-2185.7320	4476.5308	0.4897	0.9833
HKYIG	17	-2185.8269	4483.2872	0.0167	1.0000
GTRI	20	-2182.4859	4496.3052	0.0000	1.0000
GTRG	20	-2183.3438	4498.0211	0.0000	1.0000
GTRIG	21	-2182.4145	4502.7291	0.0000	1.0000
HKY	15	-2216.1878	4530.8756	0.0000	1.0000
K2PI	13	-2231.8490	4549.0646	0.0000	1.0000
K2PG	13	-2232.0626	4549.4918	0.0000	1.0000
GTR	19	-2214.7293	4554.2253	0.0000	1.0000
K2PIG	14	-2232.1515	4556.2365	0.0000	1.0000
SYMI	17	-2226.4105	4564.4545	0.0000	1.0000
SYMG	17	-2226.6721	4564.9776	0.0000	1.0000
SYMIG	18	-2227.2819	4572.7638	0.0000	1.0000
F81I	15	-2245.9683	4590.4367	0.0000	1.0000
F81G	15	-2246.5600	4591.6201	0.0000	1.0000
F81IG	16	-2245.3813	4595.8293	0.0000	1.0000
K2P	12	-2262.2577	4603.3154	0.0000	1.0000
SYM	16	-2258.0563	4621.1793	0.0000	1.0000
F81	14	-2273.5480	4639.0295	0.0000	1.0000
JC69I	12	-2284.0253	4646.8507	0.0000	1.0000
JC69G	12	-2284.5787	4647.9575	0.0000	1.0000
JC69IG	13	-2284.0274	4653.4216	0.0000	1.0000
JC69	11	-2312.1092	4696.4517	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/chr5_11146.phylip = 1 - 711;
 Partition Dummy = 1:591-tomodel-phylip/chr5_11146.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/chr5_11146.phylip = 1 - 711;
 Partition Dummy = 1:591-tomodel-phylip/chr5_11146.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/chr5_11146.phylip = 1 - 711;
 Partition Dummy = 1:591-tomodel-phylip/chr5_11146.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

