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

Input data from file "261-aggressive-phylip/3385_hg18.phylip" (29 taxa, 698 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	60	-1670.5329	3472.5572	0.7385	0.7385
HKYIG	61	-1670.5340	3474.9611	0.2220	0.9605
GTRG	64	-1668.8938	3478.9314	0.0305	0.9910
GTRIG	65	-1668.8949	3481.3658	0.0090	1.0000
F81G	59	-1688.0637	3505.2245	0.0000	1.0000
K2PG	57	-1690.9679	3506.2671	0.0000	1.0000
F81IG	60	-1688.0649	3507.6211	0.0000	1.0000
K2PIG	58	-1690.9692	3508.6488	0.0000	1.0000
SYMG	61	-1688.0469	3509.9869	0.0000	1.0000
SYMIG	62	-1688.0480	3512.3984	0.0000	1.0000
HKYI	60	-1692.0095	3515.5104	0.0000	1.0000
HKY	59	-1695.1040	3519.3051	0.0000	1.0000
GTRI	64	-1690.5877	3522.3191	0.0000	1.0000
GTR	63	-1693.8992	3526.5176	0.0000	1.0000
JC69G	56	-1707.1974	3536.3542	0.0000	1.0000
JC69IG	57	-1707.1986	3538.7284	0.0000	1.0000
F81I	59	-1709.1389	3547.3750	0.0000	1.0000
K2PI	57	-1712.6510	3549.6332	0.0000	1.0000
F81	58	-1711.9418	3550.5942	0.0000	1.0000
K2P	56	-1715.7872	3553.5338	0.0000	1.0000
SYMI	61	-1710.1746	3554.2422	0.0000	1.0000
SYM	60	-1713.6152	3558.7217	0.0000	1.0000
JC69I	56	-1728.4545	3578.8684	0.0000	1.0000
JC69	55	-1731.3425	3582.2800	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-1670.5329	3461.0659	0.6682	0.6682
HKYIG	61	-1670.5340	3463.0680	0.2456	0.9138
GTRG	64	-1668.8938	3465.7877	0.0630	0.9768
GTRIG	65	-1668.8949	3467.7899	0.0232	1.0000
F81G	59	-1688.0637	3494.1273	0.0000	1.0000
K2PG	57	-1690.9679	3495.9359	0.0000	1.0000
F81IG	60	-1688.0649	3496.1297	0.0000	1.0000
K2PIG	58	-1690.9692	3497.9383	0.0000	1.0000
SYMG	61	-1688.0469	3498.0938	0.0000	1.0000
SYMIG	62	-1688.0480	3500.0960	0.0000	1.0000
HKYI	60	-1692.0095	3504.0191	0.0000	1.0000
HKY	59	-1695.1040	3508.2079	0.0000	1.0000
GTRI	64	-1690.5877	3509.1754	0.0000	1.0000
GTR	63	-1693.8992	3513.7984	0.0000	1.0000
JC69G	56	-1707.1974	3526.3948	0.0000	1.0000
JC69IG	57	-1707.1986	3528.3972	0.0000	1.0000
F81I	59	-1709.1389	3536.2778	0.0000	1.0000
K2PI	57	-1712.6510	3539.3019	0.0000	1.0000
F81	58	-1711.9418	3539.8837	0.0000	1.0000
SYMI	61	-1710.1746	3542.3491	0.0000	1.0000
K2P	56	-1715.7872	3543.5744	0.0000	1.0000
SYM	60	-1713.6152	3547.2303	0.0000	1.0000
JC69I	56	-1728.4545	3568.9089	0.0000	1.0000
JC69	55	-1731.3425	3572.6850	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1670.5329	3733.9590	0.9635	0.9635
HKYIG	61	-1670.5340	3740.5094	0.0364	1.0000
K2PG	57	-1690.9679	3755.1844	0.0000	1.0000
GTRG	64	-1668.8938	3756.8737	0.0000	1.0000
K2PIG	58	-1690.9692	3761.7350	0.0000	1.0000
F81G	59	-1688.0637	3762.4722	0.0000	1.0000
GTRIG	65	-1668.8949	3763.4241	0.0000	1.0000
F81IG	60	-1688.0649	3769.0229	0.0000	1.0000
SYMG	61	-1688.0469	3775.5352	0.0000	1.0000
HKY	59	-1695.1040	3776.5528	0.0000	1.0000
HKYI	60	-1692.0095	3776.9122	0.0000	1.0000
JC69G	56	-1707.1974	3781.0950	0.0000	1.0000
SYMIG	62	-1688.0480	3782.0856	0.0000	1.0000
JC69IG	57	-1707.1986	3787.6456	0.0000	1.0000
K2P	56	-1715.7872	3798.2746	0.0000	1.0000
K2PI	57	-1712.6510	3798.5504	0.0000	1.0000
GTRI	64	-1690.5877	3800.2614	0.0000	1.0000
GTR	63	-1693.8992	3800.3362	0.0000	1.0000
F81	58	-1711.9418	3803.6804	0.0000	1.0000
F81I	59	-1709.1389	3804.6227	0.0000	1.0000
SYMI	61	-1710.1746	3819.7905	0.0000	1.0000
SYM	60	-1713.6152	3820.1234	0.0000	1.0000
JC69	55	-1731.3425	3822.8371	0.0000	1.0000
JC69I	56	-1728.4545	3823.6092	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/3385_hg18.phylip = 1 - 698;
 Partition Dummy = 1:261-aggressive-phylip/3385_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/3385_hg18.phylip = 1 - 698;
 Partition Dummy = 1:261-aggressive-phylip/3385_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 (HKYG)]

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

