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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	60	-485.9159	1149.9271	0.4683	0.4683
HKYI	60	-485.9293	1149.9538	0.4621	0.9304
HKYIG	61	-485.8647	1154.2415	0.0542	0.9845
K2PI	57	-496.7302	1158.7162	0.0058	0.9903
K2PG	57	-496.7376	1158.7311	0.0057	0.9961
GTRG	64	-482.7998	1161.7963	0.0012	0.9973
GTRI	64	-482.8532	1161.9032	0.0012	0.9985
K2PIG	58	-496.6591	1162.7870	0.0008	0.9992
HKY	59	-495.5315	1164.8111	0.0003	0.9995
F81G	59	-496.1754	1166.0989	0.0001	0.9997
F81I	59	-496.1929	1166.1339	0.0001	0.9998
GTRIG	65	-482.7337	1166.3766	0.0001	0.9999
SYMI	61	-493.4579	1169.4279	0.0000	0.9999
SYMG	61	-493.4599	1169.4318	0.0000	1.0000
F81IG	60	-496.0554	1170.2061	0.0000	1.0000
SYMIG	62	-493.3908	1173.7816	0.0000	1.0000
K2P	56	-506.8528	1174.8132	0.0000	1.0000
GTR	63	-491.9536	1175.4682	0.0000	1.0000
JC69I	56	-507.4427	1175.9932	0.0000	1.0000
JC69G	56	-507.4448	1175.9973	0.0000	1.0000
JC69IG	57	-507.3167	1179.8893	0.0000	1.0000
F81	58	-505.6111	1180.6909	0.0000	1.0000
SYM	60	-503.4267	1184.9486	0.0000	1.0000
JC69	55	-517.3675	1191.7579	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-485.9159	1091.8318	0.2991	0.2991
HKYI	60	-485.9293	1091.8586	0.2951	0.5943
GTRG	64	-482.7998	1093.5996	0.1236	0.7179
GTRI	64	-482.8532	1093.7064	0.1172	0.8350
HKYIG	61	-485.8647	1093.7295	0.1158	0.9508
GTRIG	65	-482.7337	1095.4675	0.0486	0.9994
K2PI	57	-496.7302	1107.4604	0.0001	0.9995
K2PG	57	-496.7376	1107.4752	0.0001	0.9997
SYMI	61	-493.4579	1108.9159	0.0001	0.9997
SYMG	61	-493.4599	1108.9198	0.0001	0.9998
HKY	59	-495.5315	1109.0630	0.0001	0.9998
K2PIG	58	-496.6591	1109.3182	0.0000	0.9999
GTR	63	-491.9536	1109.9072	0.0000	0.9999
F81G	59	-496.1754	1110.3508	0.0000	0.9999
F81I	59	-496.1929	1110.3858	0.0000	1.0000
SYMIG	62	-493.3908	1110.7816	0.0000	1.0000
F81IG	60	-496.0554	1112.1109	0.0000	1.0000
K2P	56	-506.8528	1125.7055	0.0000	1.0000
SYM	60	-503.4267	1126.8534	0.0000	1.0000
JC69I	56	-507.4427	1126.8855	0.0000	1.0000
JC69G	56	-507.4448	1126.8896	0.0000	1.0000
F81	58	-505.6111	1127.2222	0.0000	1.0000
JC69IG	57	-507.3167	1128.6335	0.0000	1.0000
JC69	55	-517.3675	1144.7350	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-485.9159	1285.6983	0.4601	0.4601
HKYI	60	-485.9293	1285.7251	0.4540	0.9140
HKYIG	61	-485.8647	1290.8271	0.0354	0.9495
K2PI	57	-496.7302	1291.6336	0.0237	0.9731
K2PG	57	-496.7376	1291.6484	0.0235	0.9966
K2PIG	58	-496.6591	1296.7225	0.0019	0.9985
HKY	59	-495.5315	1299.6984	0.0004	0.9989
GTRG	64	-482.7998	1300.3906	0.0003	0.9992
GTRI	64	-482.8532	1300.4974	0.0003	0.9995
F81G	59	-496.1754	1300.9862	0.0002	0.9997
F81I	59	-496.1929	1301.0212	0.0002	0.9999
GTRIG	65	-482.7337	1305.4895	0.0000	0.9999
F81IG	60	-496.0554	1305.9774	0.0000	0.9999
SYMI	61	-493.4579	1306.0135	0.0000	1.0000
SYMG	61	-493.4599	1306.0175	0.0000	1.0000
K2P	56	-506.8528	1306.6476	0.0000	1.0000
JC69I	56	-507.4427	1307.8276	0.0000	1.0000
JC69G	56	-507.4448	1307.8317	0.0000	1.0000
SYMIG	62	-493.3908	1311.1104	0.0000	1.0000
JC69IG	57	-507.3167	1312.8067	0.0000	1.0000
GTR	63	-491.9536	1313.4670	0.0000	1.0000
F81	58	-505.6111	1314.6265	0.0000	1.0000
SYM	60	-503.4267	1320.7199	0.0000	1.0000
JC69	55	-517.3675	1322.4460	0.0000	1.0000



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

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

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


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

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


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

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

