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

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

Minimum AIC  model: HKYI
Minimum AICc model: F81I
Minimum BIC  model: F81


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
F81I	59	-320.6330	806.7828	0.2557	0.2557
F81	58	-322.7495	807.1256	0.2154	0.4710
HKY	59	-321.1819	807.8806	0.1477	0.6187
HKYI	60	-319.2556	807.9706	0.1412	0.7599
F81G	59	-321.4604	808.4376	0.1118	0.8716
HKYG	60	-319.9162	809.2918	0.0729	0.9446
F81IG	60	-320.6330	810.7255	0.0356	0.9802
HKYIG	61	-319.2554	811.9665	0.0191	0.9993
GTR	63	-319.2715	820.1569	0.0003	0.9996
GTRI	64	-317.8133	821.4044	0.0002	0.9998
GTRG	64	-318.1361	822.0499	0.0001	0.9999
GTRIG	65	-317.8135	825.6270	0.0000	0.9999
JC69I	56	-335.9648	825.9297	0.0000	1.0000
JC69	55	-338.1356	826.5326	0.0000	1.0000
K2PI	57	-334.8841	827.5563	0.0000	1.0000
JC69G	56	-336.8343	827.6687	0.0000	1.0000
K2P	56	-336.9307	827.8614	0.0000	1.0000
K2PG	57	-335.6436	829.0753	0.0000	1.0000
JC69IG	57	-335.9647	829.7175	0.0000	1.0000
K2PIG	58	-334.8841	831.3948	0.0000	1.0000
SYM	60	-334.8385	839.1365	0.0000	1.0000
SYMI	61	-333.2607	839.9773	0.0000	1.0000
SYMG	61	-333.6687	840.7931	0.0000	1.0000
SYMIG	62	-333.2610	844.0289	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	60	-319.2556	758.5112	0.2529	0.2529
F81I	59	-320.6330	759.2660	0.1734	0.4263
HKYG	60	-319.9162	759.8324	0.1306	0.5570
HKY	59	-321.1819	760.3638	0.1002	0.6572
HKYIG	61	-319.2554	760.5107	0.0931	0.7502
F81G	59	-321.4604	760.9208	0.0758	0.8260
F81IG	60	-320.6330	761.2661	0.0638	0.8898
F81	58	-322.7495	761.4990	0.0568	0.9466
GTRI	64	-317.8133	763.6266	0.0196	0.9662
GTRG	64	-318.1361	764.2721	0.0142	0.9804
GTR	63	-319.2715	764.5431	0.0124	0.9928
GTRIG	65	-317.8135	765.6270	0.0072	1.0000
K2PI	57	-334.8841	783.7682	0.0000	1.0000
JC69I	56	-335.9648	783.9297	0.0000	1.0000
K2PG	57	-335.6436	785.2872	0.0000	1.0000
JC69G	56	-336.8343	785.6687	0.0000	1.0000
K2PIG	58	-334.8841	785.7681	0.0000	1.0000
K2P	56	-336.9307	785.8614	0.0000	1.0000
JC69IG	57	-335.9647	785.9294	0.0000	1.0000
JC69	55	-338.1356	786.2712	0.0000	1.0000
SYMI	61	-333.2607	788.5215	0.0000	1.0000
SYMG	61	-333.6687	789.3374	0.0000	1.0000
SYM	60	-334.8385	789.6770	0.0000	1.0000
SYMIG	62	-333.2610	790.5220	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
F81	58	-322.7495	955.3543	0.4085	0.4085
F81I	59	-320.6330	956.4637	0.2346	0.6431
HKY	59	-321.1819	957.5615	0.1355	0.7786
F81G	59	-321.4604	958.1185	0.1026	0.8811
HKYI	60	-319.2556	959.0512	0.0643	0.9455
HKYG	60	-319.9162	960.3724	0.0332	0.9787
F81IG	60	-320.6330	961.8061	0.0162	0.9949
HKYIG	61	-319.2554	964.3931	0.0045	0.9994
JC69	55	-338.1356	970.0996	0.0003	0.9996
JC69I	56	-335.9648	971.1004	0.0002	0.9998
JC69G	56	-336.8343	972.8394	0.0001	0.9998
K2P	56	-336.9307	973.0321	0.0001	0.9999
K2PI	57	-334.8841	974.2813	0.0000	0.9999
GTR	63	-319.2715	975.1101	0.0000	1.0000
K2PG	57	-335.6436	975.8003	0.0000	1.0000
JC69IG	57	-335.9647	976.4425	0.0000	1.0000
GTRI	64	-317.8133	977.5360	0.0000	1.0000
GTRG	64	-318.1361	978.1815	0.0000	1.0000
K2PIG	58	-334.8841	979.6235	0.0000	1.0000
GTRIG	65	-317.8135	982.8787	0.0000	1.0000
SYM	60	-334.8385	990.2171	0.0000	1.0000
SYMI	61	-333.2607	992.4039	0.0000	1.0000
SYMG	61	-333.6687	993.2197	0.0000	1.0000
SYMIG	62	-333.2610	997.7468	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (HKYI)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/2888_hg18.phylip = 1 - 209;
 Partition Dummy = 1:261-aggressive-phylip/2888_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best AICc model (F81I)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/2888_hg18.phylip = 1 - 209;
 Partition Dummy = 1:261-aggressive-phylip/2888_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=1 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/2888_hg18.phylip = 1 - 209;
 Partition Dummy = 1:261-aggressive-phylip/2888_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=1 rates=equal;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0);
END;


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

