
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:47 08/01/2011

Input data from file "261-aggressive-phylip/7474_hg18.phylip" (29 taxa, 951 characters)

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	64	-2092.8821	4323.1547	0.7642	0.7642
GTRIG	65	-2092.9129	4325.5207	0.2341	0.9984
HKYG	60	-2103.8811	4335.9870	0.0012	0.9996
HKYIG	61	-2103.9120	4338.3324	0.0004	1.0000
GTR	63	-2130.1400	4395.3714	0.0000	1.0000
K2PG	57	-2137.1902	4395.7846	0.0000	1.0000
GTRI	64	-2130.1654	4397.7214	0.0000	1.0000
K2PIG	58	-2137.2211	4398.1148	0.0000	1.0000
SYMG	61	-2134.7508	4400.0101	0.0000	1.0000
SYMIG	62	-2134.7817	4402.3607	0.0000	1.0000
HKY	59	-2139.4106	4404.7672	0.0000	1.0000
HKYI	60	-2139.4363	4407.0973	0.0000	1.0000
F81G	59	-2147.9833	4421.9128	0.0000	1.0000
F81IG	60	-2148.0144	4424.2536	0.0000	1.0000
K2P	56	-2173.8221	4466.7850	0.0000	1.0000
K2PI	57	-2173.8478	4469.0999	0.0000	1.0000
SYM	60	-2171.2578	4470.7402	0.0000	1.0000
JC69G	56	-2175.9080	4470.9569	0.0000	1.0000
SYMI	61	-2171.2835	4473.0754	0.0000	1.0000
JC69IG	57	-2175.9390	4473.2823	0.0000	1.0000
F81	58	-2182.8480	4489.3687	0.0000	1.0000
F81I	59	-2182.8738	4491.6938	0.0000	1.0000
JC69	55	-2211.5190	4539.9206	0.0000	1.0000
JC69I	56	-2211.5448	4542.2306	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-2092.8821	4313.7642	0.7364	0.7364
GTRIG	65	-2092.9129	4315.8258	0.2627	0.9991
HKYG	60	-2103.8811	4327.7622	0.0007	0.9998
HKYIG	61	-2103.9120	4329.8240	0.0002	1.0000
GTR	63	-2130.1400	4386.2801	0.0000	1.0000
GTRI	64	-2130.1654	4388.3309	0.0000	1.0000
K2PG	57	-2137.1902	4388.3803	0.0000	1.0000
K2PIG	58	-2137.2211	4390.4421	0.0000	1.0000
SYMG	61	-2134.7508	4391.5017	0.0000	1.0000
SYMIG	62	-2134.7817	4393.5634	0.0000	1.0000
HKY	59	-2139.4106	4396.8211	0.0000	1.0000
HKYI	60	-2139.4363	4398.8726	0.0000	1.0000
F81G	59	-2147.9833	4413.9667	0.0000	1.0000
F81IG	60	-2148.0144	4416.0289	0.0000	1.0000
K2P	56	-2173.8221	4459.6441	0.0000	1.0000
K2PI	57	-2173.8478	4461.6957	0.0000	1.0000
SYM	60	-2171.2578	4462.5155	0.0000	1.0000
JC69G	56	-2175.9080	4463.8160	0.0000	1.0000
SYMI	61	-2171.2835	4464.5670	0.0000	1.0000
JC69IG	57	-2175.9390	4465.8780	0.0000	1.0000
F81	58	-2182.8480	4481.6960	0.0000	1.0000
F81I	59	-2182.8738	4483.7476	0.0000	1.0000
JC69	55	-2211.5190	4533.0380	0.0000	1.0000
JC69I	56	-2211.5448	4535.0896	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-2103.8811	4619.2131	0.9094	0.9094
GTRG	64	-2092.8821	4624.6451	0.0601	0.9695
HKYIG	61	-2103.9120	4626.1323	0.0286	0.9981
GTRIG	65	-2092.9129	4631.5642	0.0019	1.0000
K2PG	57	-2137.1902	4665.2586	0.0000	1.0000
K2PIG	58	-2137.2211	4672.1780	0.0000	1.0000
HKY	59	-2139.4106	4683.4144	0.0000	1.0000
SYMG	61	-2134.7508	4687.8100	0.0000	1.0000
HKYI	60	-2139.4363	4690.3234	0.0000	1.0000
GTR	63	-2130.1400	4692.3034	0.0000	1.0000
SYMIG	62	-2134.7817	4694.7293	0.0000	1.0000
GTRI	64	-2130.1654	4699.2118	0.0000	1.0000
F81G	59	-2147.9833	4700.5600	0.0000	1.0000
F81IG	60	-2148.0144	4707.4797	0.0000	1.0000
K2P	56	-2173.8221	4731.6649	0.0000	1.0000
JC69G	56	-2175.9080	4735.8368	0.0000	1.0000
K2PI	57	-2173.8478	4738.5740	0.0000	1.0000
JC69IG	57	-2175.9390	4742.7563	0.0000	1.0000
SYM	60	-2171.2578	4753.9664	0.0000	1.0000
SYMI	61	-2171.2835	4760.8753	0.0000	1.0000
F81	58	-2182.8480	4763.4319	0.0000	1.0000
F81I	59	-2182.8738	4770.3410	0.0000	1.0000
JC69	55	-2211.5190	4800.2012	0.0000	1.0000
JC69I	56	-2211.5448	4807.1104	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (GTRG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/7474_hg18.phylip = 1 - 951;
 Partition Dummy = 1:261-aggressive-phylip/7474_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


[Mrbayes block for the best AICc model (GTRG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/7474_hg18.phylip = 1 - 951;
 Partition Dummy = 1:261-aggressive-phylip/7474_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/7474_hg18.phylip = 1 - 951;
 Partition Dummy = 1:261-aggressive-phylip/7474_hg18.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

