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

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

Minimum AIC  model: GTRIG
Minimum AICc model: GTRI
Minimum BIC  model: GTRI


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
GTRI	64	-1917.7954	3975.0666	0.4976	0.4976
GTRIG	65	-1916.6837	3975.2183	0.4612	0.9588
GTRG	64	-1920.2870	3980.0499	0.0412	0.9999
HKYIG	61	-1930.7167	3993.8235	0.0000	1.0000
HKYI	60	-1932.9020	3995.8452	0.0000	1.0000
HKYG	60	-1933.6312	3997.3036	0.0000	1.0000
SYMIG	62	-1945.9216	4026.5888	0.0000	1.0000
SYMI	61	-1947.9651	4028.3202	0.0000	1.0000
SYMG	61	-1948.6130	4029.6160	0.0000	1.0000
K2PIG	58	-1954.1707	4033.7039	0.0000	1.0000
K2PI	57	-1956.2435	4035.5197	0.0000	1.0000
K2PG	57	-1956.7368	4036.5063	0.0000	1.0000
F81IG	60	-1953.9931	4038.0273	0.0000	1.0000
F81I	59	-1955.9894	4039.6775	0.0000	1.0000
F81G	59	-1956.3772	4040.4531	0.0000	1.0000
GTR	63	-1966.3184	4069.7443	0.0000	1.0000
JC69IG	57	-1974.7186	4072.4701	0.0000	1.0000
JC69I	56	-1976.6980	4074.1054	0.0000	1.0000
JC69G	56	-1976.9466	4074.6026	0.0000	1.0000
HKY	59	-1977.3025	4082.3037	0.0000	1.0000
SYM	60	-1991.6313	4113.3037	0.0000	1.0000
K2P	56	-1999.2561	4119.2217	0.0000	1.0000
F81	58	-1999.1745	4123.7116	0.0000	1.0000
JC69	55	-2019.1423	4156.6769	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-1916.6837	3963.3675	0.5080	0.5080
GTRI	64	-1917.7954	3963.5907	0.4544	0.9624
GTRG	64	-1920.2870	3968.5740	0.0376	1.0000
HKYIG	61	-1930.7167	3983.4334	0.0000	1.0000
HKYI	60	-1932.9020	3985.8041	0.0000	1.0000
HKYG	60	-1933.6312	3987.2624	0.0000	1.0000
SYMIG	62	-1945.9216	4015.8433	0.0000	1.0000
SYMI	61	-1947.9651	4017.9301	0.0000	1.0000
SYMG	61	-1948.6130	4019.2259	0.0000	1.0000
K2PIG	58	-1954.1707	4024.3414	0.0000	1.0000
K2PI	57	-1956.2435	4026.4869	0.0000	1.0000
K2PG	57	-1956.7368	4027.4736	0.0000	1.0000
F81IG	60	-1953.9931	4027.9861	0.0000	1.0000
F81I	59	-1955.9894	4029.9789	0.0000	1.0000
F81G	59	-1956.3772	4030.7545	0.0000	1.0000
GTR	63	-1966.3184	4058.6369	0.0000	1.0000
JC69IG	57	-1974.7186	4063.4373	0.0000	1.0000
JC69I	56	-1976.6980	4065.3960	0.0000	1.0000
JC69G	56	-1976.9466	4065.8932	0.0000	1.0000
HKY	59	-1977.3025	4072.6050	0.0000	1.0000
SYM	60	-1991.6313	4103.2625	0.0000	1.0000
K2P	56	-1999.2561	4110.5123	0.0000	1.0000
F81	58	-1999.1745	4114.3491	0.0000	1.0000
JC69	55	-2019.1423	4148.2845	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRI	64	-1917.7954	4262.6008	0.6669	0.6669
HKYI	60	-1932.9020	4266.1260	0.1144	0.7813
GTRIG	65	-1916.6837	4267.0496	0.0721	0.8534
GTRG	64	-1920.2870	4267.5841	0.0552	0.9086
HKYG	60	-1933.6312	4267.5844	0.0552	0.9638
HKYIG	61	-1930.7167	4268.4274	0.0362	1.0000
K2PI	57	-1956.2435	4292.7928	0.0000	1.0000
K2PG	57	-1956.7368	4293.7794	0.0000	1.0000
K2PIG	58	-1954.1707	4295.3193	0.0000	1.0000
SYMI	61	-1947.9651	4302.9241	0.0000	1.0000
SYMG	61	-1948.6130	4304.2199	0.0000	1.0000
SYMIG	62	-1945.9216	4305.5093	0.0000	1.0000
F81I	59	-1955.9894	4305.6288	0.0000	1.0000
F81G	59	-1956.3772	4306.4044	0.0000	1.0000
F81IG	60	-1953.9931	4308.3081	0.0000	1.0000
JC69I	56	-1976.6980	4327.0298	0.0000	1.0000
JC69G	56	-1976.9466	4327.5270	0.0000	1.0000
JC69IG	57	-1974.7186	4329.7432	0.0000	1.0000
HKY	59	-1977.3025	4348.2550	0.0000	1.0000
GTR	63	-1966.3184	4352.9749	0.0000	1.0000
K2P	56	-1999.2561	4372.1461	0.0000	1.0000
SYM	60	-1991.6313	4383.5845	0.0000	1.0000
F81	58	-1999.1745	4385.3270	0.0000	1.0000
JC69	55	-2019.1423	4405.2464	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/506_hg18.phylip = 1 - 790;
 Partition Dummy = 1:261-aggressive-phylip/506_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 (GTRI)]

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


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

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


-------------------------------------------------------------
End of Output

