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

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

Minimum AIC  model: GTRIG
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	-2273.6417	4694.6323	0.5413	0.5413
GTRIG	65	-2272.4881	4694.9762	0.4558	0.9971
HKYG	60	-2284.9502	4706.7668	0.0013	0.9984
HKYIG	61	-2283.6921	4706.8530	0.0012	0.9996
GTRI	64	-2280.8163	4708.9814	0.0004	1.0000
HKYI	60	-2292.0077	4720.8817	0.0000	1.0000
K2PIG	58	-2301.0214	4733.7400	0.0000	1.0000
K2PG	57	-2302.4567	4734.0439	0.0000	1.0000
SYMIG	62	-2297.6735	4737.4303	0.0000	1.0000
SYMG	61	-2299.0087	4737.4862	0.0000	1.0000
K2PI	57	-2309.8712	4748.8729	0.0000	1.0000
SYMI	61	-2306.2091	4751.8870	0.0000	1.0000
GTR	63	-2329.8796	4804.4693	0.0000	1.0000
HKY	59	-2341.7618	4817.7994	0.0000	1.0000
SYM	60	-2356.2819	4849.4302	0.0000	1.0000
K2P	56	-2361.8479	4850.2712	0.0000	1.0000
F81G	59	-2358.0253	4850.3264	0.0000	1.0000
F81IG	60	-2356.8985	4850.6633	0.0000	1.0000
F81I	59	-2364.5563	4863.3884	0.0000	1.0000
JC69G	56	-2372.8313	4872.2380	0.0000	1.0000
JC69IG	57	-2371.6421	4872.4147	0.0000	1.0000
JC69I	56	-2379.4843	4885.5439	0.0000	1.0000
F81	58	-2413.3425	4958.3823	0.0000	1.0000
JC69	55	-2429.0813	4982.1944	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-2272.4881	4674.9762	0.5378	0.5378
GTRG	64	-2273.6417	4675.2835	0.4612	0.9989
HKYIG	61	-2283.6921	4689.3842	0.0004	0.9993
GTRI	64	-2280.8163	4689.6326	0.0004	0.9997
HKYG	60	-2284.9502	4689.9004	0.0003	1.0000
HKYI	60	-2292.0077	4704.0154	0.0000	1.0000
K2PIG	58	-2301.0214	4718.0428	0.0000	1.0000
K2PG	57	-2302.4567	4718.9134	0.0000	1.0000
SYMIG	62	-2297.6735	4719.3470	0.0000	1.0000
SYMG	61	-2299.0087	4720.0174	0.0000	1.0000
K2PI	57	-2309.8712	4733.7425	0.0000	1.0000
SYMI	61	-2306.2091	4734.4182	0.0000	1.0000
GTR	63	-2329.8796	4785.7593	0.0000	1.0000
HKY	59	-2341.7618	4801.5236	0.0000	1.0000
SYM	60	-2356.2819	4832.5638	0.0000	1.0000
F81IG	60	-2356.8985	4833.7970	0.0000	1.0000
F81G	59	-2358.0253	4834.0505	0.0000	1.0000
K2P	56	-2361.8479	4835.6958	0.0000	1.0000
F81I	59	-2364.5563	4847.1125	0.0000	1.0000
JC69IG	57	-2371.6421	4857.2843	0.0000	1.0000
JC69G	56	-2372.8313	4857.6627	0.0000	1.0000
JC69I	56	-2379.4843	4870.9686	0.0000	1.0000
F81	58	-2413.3425	4942.6850	0.0000	1.0000
JC69	55	-2429.0813	4968.1625	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-2284.9502	4942.1739	0.6495	0.6495
GTRG	64	-2273.6417	4944.3752	0.2161	0.8656
HKYIG	61	-2283.6921	4945.8622	0.1027	0.9683
GTRIG	65	-2272.4881	4948.2724	0.0308	0.9991
HKYI	60	-2292.0077	4956.2888	0.0006	0.9996
K2PG	57	-2302.4567	4958.5732	0.0002	0.9998
GTRI	64	-2280.8163	4958.7243	0.0002	1.0000
K2PIG	58	-2301.0214	4961.9072	0.0000	1.0000
K2PI	57	-2309.8712	4973.4023	0.0000	1.0000
SYMG	61	-2299.0087	4976.4954	0.0000	1.0000
SYMIG	62	-2297.6735	4980.0296	0.0000	1.0000
SYMI	61	-2306.2091	4990.8962	0.0000	1.0000
HKY	59	-2341.7618	5049.5925	0.0000	1.0000
GTR	63	-2329.8796	5050.6464	0.0000	1.0000
K2P	56	-2361.8479	5071.1511	0.0000	1.0000
F81G	59	-2358.0253	5082.1194	0.0000	1.0000
SYM	60	-2356.2819	5084.8373	0.0000	1.0000
F81IG	60	-2356.8985	5086.0704	0.0000	1.0000
JC69G	56	-2372.8313	5093.1179	0.0000	1.0000
F81I	59	-2364.5563	5095.1814	0.0000	1.0000
JC69IG	57	-2371.6421	5096.9441	0.0000	1.0000
JC69I	56	-2379.4843	5106.4238	0.0000	1.0000
F81	58	-2413.3425	5186.5494	0.0000	1.0000
JC69	55	-2429.0813	5199.4132	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/12070_hg18.phylip = 1 - 495;
 Partition Dummy = 1:261-aggressive-phylip/12070_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 (GTRG)]

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

