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

Input data from file "261-aggressive-phylip/7695_hg18.phylip" (29 taxa, 461 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	-1971.3260	4091.6620	0.7285	0.7285
GTRIG	65	-1971.0718	4093.8652	0.2421	0.9706
GTRI	64	-1974.9426	4098.8954	0.0196	0.9901
HKYG	60	-1981.8927	4102.0853	0.0040	0.9941
SYMG	61	-1980.8878	4102.7330	0.0029	0.9970
HKYIG	61	-1981.5347	4104.0268	0.0015	0.9985
SYMIG	62	-1980.5240	4104.6761	0.0011	0.9996
K2PG	57	-1988.9613	4108.3295	0.0002	0.9997
SYMI	61	-1984.3680	4109.6934	0.0001	0.9998
HKYI	60	-1985.7013	4109.7025	0.0001	0.9999
K2PIG	58	-1988.4386	4109.9020	0.0001	1.0000
K2PI	57	-1992.7870	4115.9810	0.0000	1.0000
GTR	63	-2022.0132	4190.3388	0.0000	1.0000
SYM	60	-2030.2348	4198.7696	0.0000	1.0000
HKY	59	-2031.9739	4199.6036	0.0000	1.0000
K2P	56	-2038.5535	4204.9089	0.0000	1.0000
F81G	59	-2043.1844	4222.0246	0.0000	1.0000
JC69G	56	-2047.7972	4223.3964	0.0000	1.0000
F81IG	60	-2042.8044	4223.9087	0.0000	1.0000
JC69IG	57	-2047.3851	4225.1772	0.0000	1.0000
F81I	59	-2047.0368	4229.7295	0.0000	1.0000
JC69I	56	-2051.7029	4231.2078	0.0000	1.0000
F81	58	-2091.9842	4316.9933	0.0000	1.0000
JC69	55	-2096.5520	4318.3138	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1971.3260	4070.6519	0.6642	0.6642
GTRIG	65	-1971.0718	4072.1437	0.3150	0.9792
GTRI	64	-1974.9426	4077.8853	0.0178	0.9971
SYMG	61	-1980.8878	4083.7756	0.0009	0.9980
HKYG	60	-1981.8927	4083.7853	0.0009	0.9989
SYMIG	62	-1980.5240	4085.0480	0.0005	0.9994
HKYIG	61	-1981.5347	4085.0694	0.0005	0.9999
SYMI	61	-1984.3680	4090.7360	0.0000	1.0000
HKYI	60	-1985.7013	4091.4025	0.0000	1.0000
K2PG	57	-1988.9613	4091.9226	0.0000	1.0000
K2PIG	58	-1988.4386	4092.8771	0.0000	1.0000
K2PI	57	-1992.7870	4099.5741	0.0000	1.0000
GTR	63	-2022.0132	4170.0265	0.0000	1.0000
SYM	60	-2030.2348	4180.4696	0.0000	1.0000
HKY	59	-2031.9739	4181.9478	0.0000	1.0000
K2P	56	-2038.5535	4189.1069	0.0000	1.0000
F81G	59	-2043.1844	4204.3687	0.0000	1.0000
F81IG	60	-2042.8044	4205.6087	0.0000	1.0000
JC69G	56	-2047.7972	4207.5944	0.0000	1.0000
JC69IG	57	-2047.3851	4208.7703	0.0000	1.0000
F81I	59	-2047.0368	4212.0737	0.0000	1.0000
JC69I	56	-2051.7029	4215.4058	0.0000	1.0000
F81	58	-2091.9842	4299.9684	0.0000	1.0000
JC69	55	-2096.5520	4303.1039	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-1988.9613	4327.5262	0.7876	0.7876
HKYG	60	-1981.8927	4331.7892	0.0935	0.8811
K2PIG	58	-1988.4386	4332.6142	0.0619	0.9429
K2PI	57	-1992.7870	4335.1777	0.0172	0.9601
GTRG	64	-1971.3260	4335.1894	0.0171	0.9772
SYMG	61	-1980.8878	4335.9129	0.0119	0.9891
HKYIG	61	-1981.5347	4337.2067	0.0062	0.9953
HKYI	60	-1985.7013	4339.4064	0.0021	0.9974
GTRIG	65	-1971.0718	4340.8145	0.0010	0.9984
SYMIG	62	-1980.5240	4341.3186	0.0008	0.9992
GTRI	64	-1974.9426	4342.4228	0.0005	0.9996
SYMI	61	-1984.3680	4342.8733	0.0004	1.0000
K2P	56	-2038.5535	4420.5772	0.0000	1.0000
HKY	59	-2031.9739	4425.8183	0.0000	1.0000
SYM	60	-2030.2348	4428.4735	0.0000	1.0000
GTR	63	-2022.0132	4430.4306	0.0000	1.0000
JC69G	56	-2047.7972	4439.0647	0.0000	1.0000
JC69IG	57	-2047.3851	4444.3739	0.0000	1.0000
JC69I	56	-2051.7029	4446.8761	0.0000	1.0000
F81G	59	-2043.1844	4448.2392	0.0000	1.0000
F81IG	60	-2042.8044	4453.6126	0.0000	1.0000
F81I	59	-2047.0368	4455.9442	0.0000	1.0000
JC69	55	-2096.5520	4530.4408	0.0000	1.0000
F81	58	-2091.9842	4539.7055	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/7695_hg18.phylip = 1 - 461;
 Partition Dummy = 1:261-aggressive-phylip/7695_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/7695_hg18.phylip = 1 - 461;
 Partition Dummy = 1:261-aggressive-phylip/7695_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/7695_hg18.phylip = 1 - 461;
 Partition Dummy = 1:261-aggressive-phylip/7695_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

