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

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

Minimum AIC  model: GTRIG
Minimum AICc model: GTRIG
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
GTRIG	65	-3648.5184	7440.7649	0.5130	0.5130
GTRG	64	-3650.5064	7442.3036	0.2377	0.7507
SYMIG	62	-3653.8763	7444.1921	0.0925	0.8432
SYMG	61	-3655.7236	7445.4727	0.0487	0.8919
K2PIG	58	-3659.5752	7445.9795	0.0378	0.9297
HKYIG	61	-3656.1502	7446.3259	0.0318	0.9615
K2PG	57	-3661.3563	7447.1580	0.0210	0.9825
HKYG	60	-3657.9537	7447.5264	0.0175	1.0000
GTRI	64	-3660.3877	7462.0661	0.0000	1.0000
SYMI	61	-3667.8559	7469.7372	0.0000	1.0000
HKYI	60	-3670.2273	7472.0737	0.0000	1.0000
K2PI	57	-3674.6810	7473.8075	0.0000	1.0000
JC69IG	57	-3748.4034	7621.2523	0.0000	1.0000
F81IG	60	-3745.3938	7622.4067	0.0000	1.0000
JC69G	56	-3750.2144	7622.4982	0.0000	1.0000
F81G	59	-3747.2259	7623.6721	0.0000	1.0000
F81I	59	-3759.5391	7648.2984	0.0000	1.0000
JC69I	56	-3763.1261	7648.3216	0.0000	1.0000
GTR	63	-3777.9608	7694.7829	0.0000	1.0000
SYM	60	-3786.5214	7704.6618	0.0000	1.0000
HKY	59	-3788.6250	7706.4702	0.0000	1.0000
K2P	56	-3794.2003	7710.4699	0.0000	1.0000
F81	58	-3875.5736	7877.9764	0.0000	1.0000
JC69	55	-3879.8722	7879.4452	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-3648.5184	7427.0369	0.6365	0.6365
GTRG	64	-3650.5064	7429.0129	0.2370	0.8735
SYMIG	62	-3653.8763	7431.7526	0.0602	0.9337
SYMG	61	-3655.7236	7433.4472	0.0258	0.9595
HKYIG	61	-3656.1502	7434.3005	0.0168	0.9764
K2PIG	58	-3659.5752	7435.1504	0.0110	0.9874
HKYG	60	-3657.9537	7435.9073	0.0075	0.9949
K2PG	57	-3661.3563	7436.7125	0.0050	1.0000
GTRI	64	-3660.3877	7448.7754	0.0000	1.0000
SYMI	61	-3667.8559	7457.7117	0.0000	1.0000
HKYI	60	-3670.2273	7460.4546	0.0000	1.0000
K2PI	57	-3674.6810	7463.3620	0.0000	1.0000
F81IG	60	-3745.3938	7610.7877	0.0000	1.0000
JC69IG	57	-3748.4034	7610.8068	0.0000	1.0000
JC69G	56	-3750.2144	7612.4288	0.0000	1.0000
F81G	59	-3747.2259	7612.4518	0.0000	1.0000
F81I	59	-3759.5391	7637.0781	0.0000	1.0000
JC69I	56	-3763.1261	7638.2522	0.0000	1.0000
GTR	63	-3777.9608	7681.9217	0.0000	1.0000
SYM	60	-3786.5214	7693.0428	0.0000	1.0000
HKY	59	-3788.6250	7695.2500	0.0000	1.0000
K2P	56	-3794.2003	7700.4005	0.0000	1.0000
F81	58	-3875.5736	7867.1473	0.0000	1.0000
JC69	55	-3879.8722	7869.7444	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-3661.3563	7695.3865	0.8139	0.8139
K2PIG	58	-3659.5752	7698.3625	0.1838	0.9977
HKYG	60	-3657.9537	7708.1957	0.0013	0.9991
SYMG	61	-3655.7236	7710.2738	0.0005	0.9996
HKYIG	61	-3656.1502	7711.1270	0.0003	0.9999
SYMIG	62	-3653.8763	7713.1173	0.0001	1.0000
GTRG	64	-3650.5064	7719.4538	0.0000	1.0000
GTRIG	65	-3648.5184	7722.0159	0.0000	1.0000
K2PI	57	-3674.6810	7722.0359	0.0000	1.0000
HKYI	60	-3670.2273	7732.7430	0.0000	1.0000
SYMI	61	-3667.8559	7734.5383	0.0000	1.0000
GTRI	64	-3660.3877	7739.2163	0.0000	1.0000
JC69G	56	-3750.2144	7866.5646	0.0000	1.0000
JC69IG	57	-3748.4034	7869.4807	0.0000	1.0000
F81G	59	-3747.2259	7880.2020	0.0000	1.0000
F81IG	60	-3745.3938	7883.0760	0.0000	1.0000
JC69I	56	-3763.1261	7892.3880	0.0000	1.0000
F81I	59	-3759.5391	7904.8284	0.0000	1.0000
K2P	56	-3794.2003	7954.5364	0.0000	1.0000
HKY	59	-3788.6250	7963.0002	0.0000	1.0000
SYM	60	-3786.5214	7965.3312	0.0000	1.0000
GTR	63	-3777.9608	7967.8245	0.0000	1.0000
JC69	55	-3879.8722	8119.3421	0.0000	1.0000
F81	58	-3875.5736	8130.3594	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/6949_hg18.phylip = 1 - 691;
 Partition Dummy = 1:261-aggressive-phylip/6949_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/6949_hg18.phylip = 1 - 691;
 Partition Dummy = 1:261-aggressive-phylip/6949_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

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

