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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: K2PG


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	-3531.1580	7203.3771	0.6314	0.6314
GTRIG	65	-3530.4821	7204.4548	0.3684	0.9998
SYMG	61	-3543.2587	7220.3362	0.0001	0.9999
SYMIG	62	-3542.6207	7221.4668	0.0001	1.0000
HKYG	60	-3546.7036	7224.8269	0.0000	1.0000
HKYIG	61	-3545.9128	7225.6444	0.0000	1.0000
K2PG	57	-3553.4543	7231.1756	0.0000	1.0000
K2PIG	58	-3552.7657	7232.1753	0.0000	1.0000
GTRI	64	-3546.6428	7234.3469	0.0000	1.0000
SYMI	61	-3557.8370	7249.4928	0.0000	1.0000
HKYI	60	-3564.3244	7260.0684	0.0000	1.0000
K2PI	57	-3568.6213	7261.5097	0.0000	1.0000
F81G	59	-3612.0402	7353.1084	0.0000	1.0000
GTR	63	-3607.4713	7353.5821	0.0000	1.0000
F81IG	60	-3611.3481	7354.1158	0.0000	1.0000
JC69G	56	-3617.2886	7356.4749	0.0000	1.0000
JC69IG	57	-3616.6193	7357.5057	0.0000	1.0000
SYM	60	-3618.6194	7368.6585	0.0000	1.0000
F81I	59	-3626.5929	7382.2138	0.0000	1.0000
JC69I	56	-3631.9065	7385.7107	0.0000	1.0000
K2P	56	-3632.0323	7385.9624	0.0000	1.0000
HKY	59	-3628.5884	7386.2048	0.0000	1.0000
F81	58	-3689.2991	7505.2420	0.0000	1.0000
JC69	55	-3694.9768	7509.4892	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-3531.1580	7190.3159	0.5803	0.5803
GTRIG	65	-3530.4821	7190.9643	0.4196	0.9999
SYMG	61	-3543.2587	7208.5175	0.0001	0.9999
SYMIG	62	-3542.6207	7209.2414	0.0000	1.0000
HKYG	60	-3546.7036	7213.4072	0.0000	1.0000
HKYIG	61	-3545.9128	7213.8256	0.0000	1.0000
K2PG	57	-3553.4543	7220.9085	0.0000	1.0000
GTRI	64	-3546.6428	7221.2857	0.0000	1.0000
K2PIG	58	-3552.7657	7221.5314	0.0000	1.0000
SYMI	61	-3557.8370	7237.6740	0.0000	1.0000
HKYI	60	-3564.3244	7248.6487	0.0000	1.0000
K2PI	57	-3568.6213	7251.2426	0.0000	1.0000
GTR	63	-3607.4713	7340.9426	0.0000	1.0000
F81G	59	-3612.0402	7342.0803	0.0000	1.0000
F81IG	60	-3611.3481	7342.6961	0.0000	1.0000
JC69G	56	-3617.2886	7346.5772	0.0000	1.0000
JC69IG	57	-3616.6193	7347.2386	0.0000	1.0000
SYM	60	-3618.6194	7357.2389	0.0000	1.0000
F81I	59	-3626.5929	7371.1858	0.0000	1.0000
HKY	59	-3628.5884	7375.1768	0.0000	1.0000
JC69I	56	-3631.9065	7375.8130	0.0000	1.0000
K2P	56	-3632.0323	7376.0647	0.0000	1.0000
F81	58	-3689.2991	7494.5982	0.0000	1.0000
JC69	55	-3694.9768	7499.9536	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-3553.4543	7480.4827	0.5720	0.5720
GTRG	64	-3531.1580	7481.7677	0.3009	0.8729
K2PIG	58	-3552.7657	7485.6596	0.0430	0.9159
SYMG	61	-3543.2587	7486.3074	0.0311	0.9470
HKYG	60	-3546.7036	7486.6432	0.0263	0.9733
GTRIG	65	-3530.4821	7486.9699	0.0223	0.9956
SYMIG	62	-3542.6207	7491.5853	0.0022	0.9978
HKYIG	61	-3545.9128	7491.6156	0.0022	1.0000
K2PI	57	-3568.6213	7510.8168	0.0000	1.0000
GTRI	64	-3546.6428	7512.7374	0.0000	1.0000
SYMI	61	-3557.8370	7515.4640	0.0000	1.0000
HKYI	60	-3564.3244	7521.8847	0.0000	1.0000
JC69G	56	-3617.2886	7601.5975	0.0000	1.0000
JC69IG	57	-3616.6193	7606.8128	0.0000	1.0000
F81G	59	-3612.0402	7610.7624	0.0000	1.0000
F81IG	60	-3611.3481	7615.9321	0.0000	1.0000
GTR	63	-3607.4713	7627.8404	0.0000	1.0000
SYM	60	-3618.6194	7630.4749	0.0000	1.0000
JC69I	56	-3631.9065	7630.8333	0.0000	1.0000
K2P	56	-3632.0323	7631.0850	0.0000	1.0000
F81I	59	-3626.5929	7639.8678	0.0000	1.0000
HKY	59	-3628.5884	7643.8589	0.0000	1.0000
JC69	55	-3694.9768	7750.4199	0.0000	1.0000
F81	58	-3689.2991	7758.7263	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/6339_hg18.phylip = 1 - 702;
 Partition Dummy = 1:261-aggressive-phylip/6339_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/6339_hg18.phylip = 1 - 702;
 Partition Dummy = 1:261-aggressive-phylip/6339_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 (K2PG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/6339_hg18.phylip = 1 - 702;
 Partition Dummy = 1:261-aggressive-phylip/6339_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


-------------------------------------------------------------
End of Output

