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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	60	-659.7255	1452.7118	0.4825	0.4825
HKY	59	-661.5277	1453.8582	0.2720	0.7545
HKYIG	61	-659.7607	1455.2491	0.1357	0.8902
HKYI	60	-661.5624	1456.3857	0.0769	0.9670
GTRG	64	-658.1237	1459.4299	0.0168	0.9838
GTR	63	-659.9867	1460.6620	0.0091	0.9929
GTRIG	65	-658.1590	1462.0036	0.0046	0.9975
GTRI	64	-660.0217	1463.2258	0.0025	1.0000
K2PG	57	-676.2148	1478.3431	0.0000	1.0000
K2P	56	-677.6804	1478.8428	0.0000	1.0000
K2PIG	58	-676.2500	1480.8538	0.0000	1.0000
K2PI	57	-677.7152	1481.3439	0.0000	1.0000
F81G	59	-676.2112	1483.2252	0.0000	1.0000
F81	58	-677.7059	1483.7656	0.0000	1.0000
SYMG	61	-674.5430	1484.8138	0.0000	1.0000
SYM	60	-675.9980	1485.2569	0.0000	1.0000
F81IG	60	-676.2466	1485.7541	0.0000	1.0000
F81I	59	-677.7409	1486.2848	0.0000	1.0000
SYMIG	62	-674.5782	1487.3600	0.0000	1.0000
SYMI	61	-676.0328	1487.7934	0.0000	1.0000
JC69G	56	-690.1911	1503.8641	0.0000	1.0000
JC69	55	-691.6036	1504.2664	0.0000	1.0000
JC69IG	57	-690.2264	1506.3664	0.0000	1.0000
JC69I	56	-691.6386	1506.7592	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-659.7255	1439.4509	0.4677	0.4677
HKY	59	-661.5277	1441.0553	0.2097	0.6775
HKYIG	61	-659.7607	1441.5213	0.1661	0.8436
HKYI	60	-661.5624	1443.1249	0.0745	0.9181
GTRG	64	-658.1237	1444.2474	0.0425	0.9606
GTR	63	-659.9867	1445.9734	0.0179	0.9785
GTRIG	65	-658.1590	1446.3181	0.0151	0.9936
GTRI	64	-660.0217	1448.0433	0.0064	1.0000
K2PG	57	-676.2148	1466.4296	0.0000	1.0000
K2P	56	-677.6804	1467.3608	0.0000	1.0000
K2PIG	58	-676.2500	1468.5000	0.0000	1.0000
K2PI	57	-677.7152	1469.4304	0.0000	1.0000
F81G	59	-676.2112	1470.4223	0.0000	1.0000
SYMG	61	-674.5430	1471.0860	0.0000	1.0000
F81	58	-677.7059	1471.4118	0.0000	1.0000
SYM	60	-675.9980	1471.9960	0.0000	1.0000
F81IG	60	-676.2466	1472.4932	0.0000	1.0000
SYMIG	62	-674.5782	1473.1564	0.0000	1.0000
F81I	59	-677.7409	1473.4819	0.0000	1.0000
SYMI	61	-676.0328	1474.0656	0.0000	1.0000
JC69G	56	-690.1911	1492.3821	0.0000	1.0000
JC69	55	-691.6036	1493.2072	0.0000	1.0000
JC69IG	57	-690.2264	1494.4529	0.0000	1.0000
JC69I	56	-691.6386	1495.2772	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKY	59	-661.5277	1701.7389	0.7721	0.7721
HKYG	60	-659.7255	1704.5528	0.1891	0.9611
HKYI	60	-661.5624	1708.2268	0.0301	0.9912
HKYIG	61	-659.7607	1711.0416	0.0074	0.9986
K2P	56	-677.6804	1714.7892	0.0011	0.9997
K2PG	57	-676.2148	1718.2764	0.0002	0.9999
K2PI	57	-677.7152	1721.2772	0.0000	1.0000
GTR	63	-659.9867	1724.3304	0.0000	1.0000
K2PIG	58	-676.2500	1724.7652	0.0000	1.0000
GTRG	64	-658.1237	1727.0228	0.0000	1.0000
F81	58	-677.7059	1727.6770	0.0000	1.0000
GTRI	64	-660.0217	1730.8187	0.0000	1.0000
F81G	59	-676.2112	1731.1059	0.0000	1.0000
GTRIG	65	-658.1590	1733.5118	0.0000	1.0000
F81I	59	-677.7409	1734.1654	0.0000	1.0000
JC69	55	-691.6036	1736.2173	0.0000	1.0000
SYM	60	-675.9980	1737.0979	0.0000	1.0000
F81IG	60	-676.2466	1737.5951	0.0000	1.0000
JC69G	56	-690.1911	1739.8105	0.0000	1.0000
SYMG	61	-674.5430	1740.6063	0.0000	1.0000
JC69I	56	-691.6386	1742.7056	0.0000	1.0000
SYMI	61	-676.0328	1743.5859	0.0000	1.0000
JC69IG	57	-690.2264	1746.2997	0.0000	1.0000
SYMIG	62	-674.5782	1747.0950	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 261-aggressive-phylip/2337_hg18.phylip = 1 - 613;
 Partition Dummy = 1:261-aggressive-phylip/2337_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;


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

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


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

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

