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

Input data from file "261-aggressive-phylip/749_hg18.phylip" (29 taxa, 918 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	-3678.8694	7495.4926	0.9988	0.9988
GTRIG	65	-3684.6786	7509.4277	0.0009	0.9998
HKYG	60	-3691.9824	7512.5062	0.0002	1.0000
HKYIG	61	-3693.2169	7517.2703	0.0000	1.0000
GTRI	64	-3699.0332	7535.8203	0.0000	1.0000
HKYI	60	-3703.6422	7535.8259	0.0000	1.0000
SYMG	61	-3716.1435	7563.1235	0.0000	1.0000
K2PG	57	-3723.0931	7567.8745	0.0000	1.0000
SYMIG	62	-3721.2769	7575.6906	0.0000	1.0000
K2PIG	58	-3728.6485	7581.2644	0.0000	1.0000
SYMI	61	-3731.4286	7593.6937	0.0000	1.0000
K2PI	57	-3738.1381	7597.9647	0.0000	1.0000
GTR	63	-3770.8148	7677.0722	0.0000	1.0000
HKY	59	-3782.1306	7690.5130	0.0000	1.0000
SYM	60	-3799.1781	7726.8977	0.0000	1.0000
K2P	56	-3807.5635	7734.5416	0.0000	1.0000
F81G	59	-3815.6010	7757.4537	0.0000	1.0000
F81IG	60	-3819.9329	7768.4073	0.0000	1.0000
F81I	59	-3828.7214	7783.6946	0.0000	1.0000
JC69G	56	-3837.8821	7795.1788	0.0000	1.0000
JC69IG	57	-3842.7759	7807.2402	0.0000	1.0000
JC69I	56	-3850.7314	7820.8774	0.0000	1.0000
F81	58	-3897.3895	7918.7465	0.0000	1.0000
JC69	55	-3918.5165	7954.1791	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-3678.8694	7485.7388	0.9988	0.9988
GTRIG	65	-3684.6786	7499.3573	0.0011	0.9999
HKYG	60	-3691.9824	7503.9647	0.0001	1.0000
HKYIG	61	-3693.2169	7508.4338	0.0000	1.0000
GTRI	64	-3699.0332	7526.0665	0.0000	1.0000
HKYI	60	-3703.6422	7527.2845	0.0000	1.0000
SYMG	61	-3716.1435	7554.2871	0.0000	1.0000
K2PG	57	-3723.0931	7560.1861	0.0000	1.0000
SYMIG	62	-3721.2769	7566.5537	0.0000	1.0000
K2PIG	58	-3728.6485	7573.2970	0.0000	1.0000
SYMI	61	-3731.4286	7584.8573	0.0000	1.0000
K2PI	57	-3738.1381	7590.2763	0.0000	1.0000
GTR	63	-3770.8148	7667.6296	0.0000	1.0000
HKY	59	-3782.1306	7682.2612	0.0000	1.0000
SYM	60	-3799.1781	7718.3563	0.0000	1.0000
K2P	56	-3807.5635	7727.1270	0.0000	1.0000
F81G	59	-3815.6010	7749.2020	0.0000	1.0000
F81IG	60	-3819.9329	7759.8659	0.0000	1.0000
F81I	59	-3828.7214	7775.4428	0.0000	1.0000
JC69G	56	-3837.8821	7787.7642	0.0000	1.0000
JC69IG	57	-3842.7759	7799.5518	0.0000	1.0000
JC69I	56	-3850.7314	7813.4628	0.0000	1.0000
F81	58	-3897.3895	7910.7791	0.0000	1.0000
JC69	55	-3918.5165	7947.0330	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-3691.9824	7793.2966	0.6260	0.6260
GTRG	64	-3678.8694	7794.3595	0.3679	0.9939
HKYIG	61	-3693.2169	7802.5879	0.0060	1.0000
GTRIG	65	-3684.6786	7812.8001	0.0000	1.0000
HKYI	60	-3703.6422	7816.6163	0.0000	1.0000
GTRI	64	-3699.0332	7834.6871	0.0000	1.0000
K2PG	57	-3723.0931	7835.0514	0.0000	1.0000
SYMG	61	-3716.1435	7848.4411	0.0000	1.0000
K2PIG	58	-3728.6485	7852.9844	0.0000	1.0000
K2PI	57	-3738.1381	7865.1415	0.0000	1.0000
SYMIG	62	-3721.2769	7865.5300	0.0000	1.0000
SYMI	61	-3731.4286	7879.0113	0.0000	1.0000
HKY	59	-3782.1306	7966.7709	0.0000	1.0000
GTR	63	-3770.8148	7971.4281	0.0000	1.0000
K2P	56	-3807.5635	7997.1700	0.0000	1.0000
SYM	60	-3799.1781	8007.6881	0.0000	1.0000
F81G	59	-3815.6010	8033.7116	0.0000	1.0000
F81IG	60	-3819.9329	8049.1977	0.0000	1.0000
JC69G	56	-3837.8821	8057.8072	0.0000	1.0000
F81I	59	-3828.7214	8059.9525	0.0000	1.0000
JC69IG	57	-3842.7759	8074.4171	0.0000	1.0000
JC69I	56	-3850.7314	8083.5059	0.0000	1.0000
F81	58	-3897.3895	8190.4665	0.0000	1.0000
JC69	55	-3918.5165	8212.2538	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/749_hg18.phylip = 1 - 918;
 Partition Dummy = 1:261-aggressive-phylip/749_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/749_hg18.phylip = 1 - 918;
 Partition Dummy = 1:261-aggressive-phylip/749_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/749_hg18.phylip = 1 - 918;
 Partition Dummy = 1:261-aggressive-phylip/749_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

