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

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

Minimum AIC  model: HKYIG
Minimum AICc model: K2PG
Minimum BIC  model: K2PG


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
K2PG	57	-530.5594	1228.4415	0.4168	0.4168
K2PI	57	-530.8455	1229.0135	0.3131	0.7300
K2PIG	58	-529.3983	1230.4389	0.1535	0.8835
HKYG	60	-526.0229	1232.5416	0.0537	0.9372
HKYI	60	-526.3230	1233.1419	0.0397	0.9769
HKYIG	61	-524.8931	1234.8196	0.0172	0.9941
SYMG	61	-527.2376	1239.5086	0.0016	0.9957
SYMI	61	-527.5765	1240.1862	0.0012	0.9969
JC69G	56	-538.9192	1240.9104	0.0008	0.9977
JC69I	56	-539.2143	1241.5006	0.0006	0.9983
SYMIG	62	-525.9977	1241.6424	0.0006	0.9989
JC69IG	57	-537.5889	1242.5003	0.0004	0.9993
GTRG	64	-522.0397	1243.1905	0.0003	0.9995
GTRI	64	-522.2536	1243.6183	0.0002	0.9997
F81G	59	-534.6317	1245.2961	0.0001	0.9998
F81I	59	-534.9001	1245.8330	0.0001	0.9999
GTRIG	65	-520.9942	1245.9540	0.0001	1.0000
F81IG	60	-533.3329	1247.1616	0.0000	1.0000
K2P	56	-544.8859	1252.8437	0.0000	1.0000
HKY	59	-540.1793	1256.3914	0.0000	1.0000
SYM	60	-542.1828	1264.8615	0.0000	1.0000
JC69	55	-553.0564	1265.0017	0.0000	1.0000
GTR	63	-537.2556	1268.8503	0.0000	1.0000
F81	58	-548.7473	1269.1369	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYIG	61	-524.8931	1171.7862	0.1714	0.1714
GTRIG	65	-520.9942	1171.9885	0.1549	0.3263
HKYG	60	-526.0229	1172.0458	0.1505	0.4768
GTRG	64	-522.0397	1172.0794	0.1480	0.6249
GTRI	64	-522.2536	1172.5071	0.1195	0.7444
HKYI	60	-526.3230	1172.6460	0.1115	0.8559
K2PIG	58	-529.3983	1174.7966	0.0380	0.8939
K2PG	57	-530.5594	1175.1189	0.0324	0.9263
K2PI	57	-530.8455	1175.6909	0.0243	0.9507
SYMIG	62	-525.9977	1175.9953	0.0209	0.9716
SYMG	61	-527.2376	1176.4753	0.0164	0.9880
SYMI	61	-527.5765	1177.1529	0.0117	0.9997
F81IG	60	-533.3329	1186.6657	0.0001	0.9998
F81G	59	-534.6317	1187.2633	0.0001	0.9999
F81I	59	-534.9001	1187.8002	0.0001	0.9999
JC69IG	57	-537.5889	1189.1778	0.0000	1.0000
JC69G	56	-538.9192	1189.8384	0.0000	1.0000
JC69I	56	-539.2143	1190.4286	0.0000	1.0000
HKY	59	-540.1793	1198.3586	0.0000	1.0000
GTR	63	-537.2556	1200.5113	0.0000	1.0000
K2P	56	-544.8859	1201.7717	0.0000	1.0000
SYM	60	-542.1828	1204.3656	0.0000	1.0000
F81	58	-548.7473	1213.4946	0.0000	1.0000
JC69	55	-553.0564	1216.1128	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-530.5594	1357.7473	0.4828	0.4828
K2PI	57	-530.8455	1358.3193	0.3627	0.8455
K2PIG	58	-529.3983	1360.6290	0.1143	0.9598
HKYG	60	-526.0229	1364.2862	0.0184	0.9782
HKYI	60	-526.3230	1364.8864	0.0136	0.9917
HKYIG	61	-524.8931	1367.2306	0.0042	0.9960
JC69G	56	-538.9192	1369.2628	0.0015	0.9975
JC69I	56	-539.2143	1369.8530	0.0011	0.9986
JC69IG	57	-537.5889	1371.8061	0.0004	0.9990
SYMG	61	-527.2376	1371.9197	0.0004	0.9995
SYMI	61	-527.5765	1372.5973	0.0003	0.9997
SYMIG	62	-525.9977	1374.6438	0.0001	0.9998
F81G	59	-534.6317	1376.2997	0.0000	0.9999
F81I	59	-534.9001	1376.8366	0.0000	0.9999
GTRG	64	-522.0397	1377.1358	0.0000	1.0000
GTRI	64	-522.2536	1377.5636	0.0000	1.0000
F81IG	60	-533.3329	1378.9061	0.0000	1.0000
GTRIG	65	-520.9942	1380.2489	0.0000	1.0000
K2P	56	-544.8859	1381.1961	0.0000	1.0000
HKY	59	-540.1793	1387.3950	0.0000	1.0000
JC69	55	-553.0564	1392.3332	0.0000	1.0000
SYM	60	-542.1828	1396.6060	0.0000	1.0000
F81	58	-548.7473	1399.3270	0.0000	1.0000
GTR	63	-537.2556	1402.3637	0.0000	1.0000



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

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

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


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/1620_hg18.phylip = 1 - 182;
 Partition Dummy = 1:261-aggressive-phylip/1620_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/1620_hg18.phylip = 1 - 182;
 Partition Dummy = 1:261-aggressive-phylip/1620_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

