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

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

Minimum AIC  model: GTRG
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	-1843.6768	3833.0377	0.4136	0.4136
GTRG	64	-1838.3410	3834.2906	0.2210	0.6346
HKYI	60	-1844.8108	3835.3057	0.1331	0.7677
HKYIG	61	-1843.4854	3835.6046	0.1146	0.8823
GTRIG	65	-1838.1475	3836.9378	0.0588	0.9411
GTRI	64	-1839.6663	3836.9411	0.0587	0.9998
HKY	59	-1853.0582	3848.8716	0.0002	1.0000
K2PG	57	-1859.0673	3855.0930	0.0000	1.0000
K2PIG	58	-1858.9109	3857.6685	0.0000	1.0000
K2PI	57	-1860.4215	3857.8013	0.0000	1.0000
GTR	63	-1852.0480	3858.6918	0.0000	1.0000
SYMG	61	-1857.1182	3862.8701	0.0000	1.0000
SYMIG	62	-1857.0313	3865.6668	0.0000	1.0000
SYMI	61	-1858.7991	3866.2321	0.0000	1.0000
K2P	56	-1873.1903	3880.4706	0.0000	1.0000
SYM	60	-1871.1332	3887.9507	0.0000	1.0000
F81I	59	-1894.1263	3931.0078	0.0000	1.0000
F81G	59	-1896.8840	3936.5232	0.0000	1.0000
F81IG	60	-1896.7381	3939.1603	0.0000	1.0000
JC69I	56	-1904.2633	3942.6165	0.0000	1.0000
JC69G	56	-1906.9377	3947.9654	0.0000	1.0000
JC69IG	57	-1906.8007	3950.5597	0.0000	1.0000
F81	58	-1909.1724	3958.1915	0.0000	1.0000
JC69	55	-1919.8331	3970.9075	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1838.3410	3804.6821	0.4594	0.4594
GTRIG	65	-1838.1475	3806.2950	0.2051	0.6644
GTRI	64	-1839.6663	3807.3326	0.1221	0.7865
HKYG	60	-1843.6768	3807.3535	0.1208	0.9073
HKYIG	61	-1843.4854	3808.9708	0.0538	0.9611
HKYI	60	-1844.8108	3809.6215	0.0389	1.0000
HKY	59	-1853.0582	3824.1164	0.0000	1.0000
GTR	63	-1852.0480	3830.0960	0.0000	1.0000
K2PG	57	-1859.0673	3832.1347	0.0000	1.0000
K2PIG	58	-1858.9109	3833.8218	0.0000	1.0000
K2PI	57	-1860.4215	3834.8429	0.0000	1.0000
SYMG	61	-1857.1182	3836.2363	0.0000	1.0000
SYMIG	62	-1857.0313	3838.0625	0.0000	1.0000
SYMI	61	-1858.7991	3839.5983	0.0000	1.0000
K2P	56	-1873.1903	3858.3806	0.0000	1.0000
SYM	60	-1871.1332	3862.2664	0.0000	1.0000
F81I	59	-1894.1263	3906.2525	0.0000	1.0000
F81G	59	-1896.8840	3911.7679	0.0000	1.0000
F81IG	60	-1896.7381	3913.4761	0.0000	1.0000
JC69I	56	-1904.2633	3920.5266	0.0000	1.0000
JC69G	56	-1906.9377	3925.8754	0.0000	1.0000
JC69IG	57	-1906.8007	3927.6013	0.0000	1.0000
F81	58	-1909.1724	3934.3448	0.0000	1.0000
JC69	55	-1919.8331	3949.6661	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1843.6768	4038.1398	0.7181	0.7181
HKYI	60	-1844.8108	4040.4078	0.2311	0.9492
HKYIG	61	-1843.4854	4043.6036	0.0467	0.9960
GTRG	64	-1838.3410	4050.8542	0.0012	0.9972
HKY	59	-1853.0582	4051.0562	0.0011	0.9983
K2PG	57	-1859.0673	4051.3817	0.0010	0.9993
GTRI	64	-1839.6663	4053.5047	0.0003	0.9996
K2PI	57	-1860.4215	4054.0900	0.0002	0.9999
GTRIG	65	-1838.1475	4056.3135	0.0001	0.9999
K2PIG	58	-1858.9109	4056.9152	0.0001	1.0000
SYMG	61	-1857.1182	4070.8691	0.0000	1.0000
GTR	63	-1852.0480	4072.4217	0.0000	1.0000
K2P	56	-1873.1903	4073.7812	0.0000	1.0000
SYMI	61	-1858.7991	4074.2310	0.0000	1.0000
SYMIG	62	-1857.0313	4076.5417	0.0000	1.0000
SYM	60	-1871.1332	4093.0528	0.0000	1.0000
F81I	59	-1894.1263	4133.1924	0.0000	1.0000
JC69I	56	-1904.2633	4135.9272	0.0000	1.0000
F81G	59	-1896.8840	4138.7078	0.0000	1.0000
JC69G	56	-1906.9377	4141.2760	0.0000	1.0000
F81IG	60	-1896.7381	4144.2624	0.0000	1.0000
JC69IG	57	-1906.8007	4146.8483	0.0000	1.0000
F81	58	-1909.1724	4157.4383	0.0000	1.0000
JC69	55	-1919.8331	4161.2203	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/1589_hg18.phylip = 1 - 346;
 Partition Dummy = 1:261-aggressive-phylip/1589_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 (HKYG)]

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

