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

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

Minimum AIC  model: HKY
Minimum AICc model: HKY
Minimum BIC  model: K2P


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
HKY	59	-463.9568	1085.2470	0.3188	0.3188
K2P	56	-469.4910	1085.8673	0.2338	0.5527
HKYI	60	-463.1899	1087.2737	0.1157	0.6684
HKYG	60	-463.2147	1087.3232	0.1129	0.7813
K2PI	57	-468.7328	1087.7953	0.0892	0.8705
K2PG	57	-468.7547	1087.8390	0.0872	0.9577
HKYIG	61	-463.2214	1090.9372	0.0185	0.9763
K2PIG	58	-468.7607	1091.3335	0.0152	0.9915
SYM	60	-467.0927	1095.0792	0.0023	0.9938
GTR	63	-461.7349	1095.2881	0.0021	0.9959
SYMI	61	-466.6058	1097.7060	0.0006	0.9965
SYMG	61	-466.6391	1097.7725	0.0006	0.9972
GTRI	64	-461.1857	1097.9143	0.0006	0.9977
GTRG	64	-461.2238	1097.9905	0.0005	0.9983
F81	58	-472.3285	1098.4691	0.0004	0.9987
JC69	55	-477.5030	1098.4842	0.0004	0.9991
JC69I	56	-476.7486	1100.3823	0.0002	0.9993
JC69G	56	-476.7723	1100.4298	0.0002	0.9994
F81I	59	-471.5507	1100.4346	0.0002	0.9996
F81G	59	-471.5741	1100.4816	0.0002	0.9998
SYMIG	62	-466.6402	1101.4161	0.0001	0.9999
GTRIG	65	-461.2270	1101.7644	0.0001	0.9999
JC69IG	57	-476.7786	1103.8869	0.0000	1.0000
F81IG	60	-471.5736	1104.0412	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKY	59	-463.9568	1045.9137	0.2832	0.2832
HKYI	60	-463.1899	1046.3798	0.2244	0.5076
HKYG	60	-463.2147	1046.4293	0.2189	0.7265
HKYIG	61	-463.2214	1048.4428	0.0800	0.8065
GTR	63	-461.7349	1049.4699	0.0479	0.8543
GTRI	64	-461.1857	1050.3714	0.0305	0.8848
GTRG	64	-461.2238	1050.4476	0.0294	0.9142
K2P	56	-469.4910	1050.9821	0.0225	0.9366
K2PI	57	-468.7328	1051.4656	0.0176	0.9543
K2PG	57	-468.7547	1051.5093	0.0173	0.9715
GTRIG	65	-461.2270	1052.4540	0.0108	0.9823
K2PIG	58	-468.7607	1053.5214	0.0063	0.9886
SYM	60	-467.0927	1054.1853	0.0045	0.9931
SYMI	61	-466.6058	1055.2116	0.0027	0.9958
SYMG	61	-466.6391	1055.2781	0.0026	0.9985
SYMIG	62	-466.6402	1057.2805	0.0010	0.9994
F81	58	-472.3285	1060.6569	0.0002	0.9996
F81I	59	-471.5507	1061.1013	0.0001	0.9998
F81G	59	-471.5741	1061.1482	0.0001	0.9999
F81IG	60	-471.5736	1063.1473	0.0001	0.9999
JC69	55	-477.5030	1065.0059	0.0000	1.0000
JC69I	56	-476.7486	1065.4971	0.0000	1.0000
JC69G	56	-476.7723	1065.5446	0.0000	1.0000
JC69IG	57	-476.7786	1067.5572	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2P	56	-469.4910	1245.8978	0.7269	0.7269
K2PI	57	-468.7328	1249.8620	0.1001	0.8270
K2PG	57	-468.7547	1249.9057	0.0980	0.9250
HKY	59	-463.9568	1251.2714	0.0495	0.9745
HKYI	60	-463.1899	1255.2182	0.0069	0.9814
HKYG	60	-463.2147	1255.2676	0.0067	0.9881
K2PIG	58	-468.7607	1255.3984	0.0063	0.9944
JC69	55	-477.5030	1256.4411	0.0037	0.9981
JC69I	56	-476.7486	1260.4129	0.0005	0.9986
JC69G	56	-476.7723	1260.4603	0.0005	0.9991
HKYIG	61	-463.2214	1260.7618	0.0004	0.9996
F81	58	-472.3285	1262.5340	0.0002	0.9997
SYM	60	-467.0927	1263.0237	0.0001	0.9999
JC69IG	57	-476.7786	1265.9536	0.0000	0.9999
F81I	59	-471.5507	1266.4590	0.0000	0.9999
F81G	59	-471.5741	1266.5059	0.0000	1.0000
SYMI	61	-466.6058	1267.5306	0.0000	1.0000
SYMG	61	-466.6391	1267.5971	0.0000	1.0000
GTR	63	-461.7349	1268.7501	0.0000	1.0000
F81IG	60	-471.5736	1271.9856	0.0000	1.0000
SYMIG	62	-466.6402	1273.0801	0.0000	1.0000
GTRI	64	-461.1857	1273.1323	0.0000	1.0000
GTRG	64	-461.2238	1273.2085	0.0000	1.0000
GTRIG	65	-461.2270	1278.6956	0.0000	1.0000



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

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

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


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

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


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/11799_hg18.phylip = 1 - 240;
 Partition Dummy = 1:261-aggressive-phylip/11799_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0);
END;


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

