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

Input data from file "261-aggressive-phylip/1882_hg18.phylip" (29 taxa, 170 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	-497.5822	1182.3204	0.5079	0.5079
HKYI	60	-497.9399	1183.0358	0.3552	0.8631
HKYIG	61	-497.4698	1186.9766	0.0495	0.9126
K2PG	57	-507.2334	1187.5026	0.0381	0.9506
K2PI	57	-507.4252	1187.8861	0.0314	0.9821
F81G	59	-504.7208	1191.8053	0.0044	0.9865
K2PIG	58	-507.1426	1191.9429	0.0041	0.9906
HKY	59	-504.8368	1192.0372	0.0039	0.9946
F81I	59	-505.0128	1192.3893	0.0033	0.9979
JC69G	56	-513.6105	1195.7165	0.0006	0.9985
JC69I	56	-513.8095	1196.1147	0.0005	0.9990
F81IG	60	-504.5998	1196.3555	0.0005	0.9995
GTRG	64	-495.3780	1197.9941	0.0002	0.9997
GTRI	64	-495.8438	1198.9257	0.0001	0.9998
K2P	56	-515.6784	1199.8524	0.0001	0.9999
JC69IG	57	-513.5059	1200.0475	0.0001	0.9999
F81	58	-512.4482	1202.5541	0.0000	1.0000
GTRIG	65	-495.2586	1203.0172	0.0000	1.0000
SYMG	61	-505.9526	1203.9422	0.0000	1.0000
SYMI	61	-506.0240	1204.0850	0.0000	1.0000
GTR	63	-502.3719	1206.8192	0.0000	1.0000
JC69	55	-522.0322	1208.0995	0.0000	1.0000
SYMIG	62	-505.9304	1208.8702	0.0000	1.0000
SYM	60	-514.4328	1216.0215	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-497.5822	1115.1644	0.4065	0.4065
HKYI	60	-497.9399	1115.8798	0.2843	0.6908
HKYIG	61	-497.4698	1116.9395	0.1673	0.8581
GTRG	64	-495.3780	1118.7560	0.0675	0.9256
GTRI	64	-495.8438	1119.6876	0.0424	0.9679
GTRIG	65	-495.2586	1120.5172	0.0280	0.9959
F81G	59	-504.7208	1127.4417	0.0009	0.9968
HKY	59	-504.8368	1127.6735	0.0008	0.9976
F81I	59	-505.0128	1128.0256	0.0007	0.9982
K2PG	57	-507.2334	1128.4669	0.0005	0.9987
K2PI	57	-507.4252	1128.8504	0.0004	0.9992
F81IG	60	-504.5998	1129.1996	0.0004	0.9995
K2PIG	58	-507.1426	1130.2853	0.0002	0.9997
GTR	63	-502.3719	1130.7438	0.0002	0.9999
SYMG	61	-505.9526	1133.9052	0.0000	0.9999
SYMI	61	-506.0240	1134.0480	0.0000	1.0000
SYMIG	62	-505.9304	1135.8608	0.0000	1.0000
JC69G	56	-513.6105	1139.2210	0.0000	1.0000
JC69I	56	-513.8095	1139.6191	0.0000	1.0000
F81	58	-512.4482	1140.8964	0.0000	1.0000
JC69IG	57	-513.5059	1141.0118	0.0000	1.0000
K2P	56	-515.6784	1143.3568	0.0000	1.0000
SYM	60	-514.4328	1148.8655	0.0000	1.0000
JC69	55	-522.0322	1154.0644	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-497.5822	1303.3123	0.4779	0.4779
HKYI	60	-497.9399	1304.0277	0.3341	0.8120
K2PG	57	-507.2334	1307.2074	0.0682	0.8802
K2PI	57	-507.4252	1307.5909	0.0563	0.9364
HKYIG	61	-497.4698	1308.2232	0.0410	0.9774
K2PIG	58	-507.1426	1312.1616	0.0057	0.9832
F81G	59	-504.7208	1312.4538	0.0049	0.9881
HKY	59	-504.8368	1312.6856	0.0044	0.9925
F81I	59	-505.0128	1313.0377	0.0037	0.9962
JC69G	56	-513.6105	1314.8257	0.0015	0.9977
JC69I	56	-513.8095	1315.2238	0.0012	0.9989
F81IG	60	-504.5998	1317.3475	0.0004	0.9994
K2P	56	-515.6784	1318.9616	0.0002	0.9996
GTRG	64	-495.3780	1319.4471	0.0001	0.9997
JC69IG	57	-513.5059	1319.7523	0.0001	0.9998
GTRI	64	-495.8438	1320.3787	0.0001	0.9999
F81	58	-512.4482	1322.7727	0.0000	1.0000
GTRIG	65	-495.2586	1324.3441	0.0000	1.0000
SYMG	61	-505.9526	1325.1889	0.0000	1.0000
SYMI	61	-506.0240	1325.3317	0.0000	1.0000
JC69	55	-522.0322	1326.5333	0.0000	1.0000
GTR	63	-502.3719	1328.2991	0.0000	1.0000
SYMIG	62	-505.9304	1330.2803	0.0000	1.0000
SYM	60	-514.4328	1337.0134	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/1882_hg18.phylip = 1 - 170;
 Partition Dummy = 1:261-aggressive-phylip/1882_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/1882_hg18.phylip = 1 - 170;
 Partition Dummy = 1:261-aggressive-phylip/1882_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/1882_hg18.phylip = 1 - 170;
 Partition Dummy = 1:261-aggressive-phylip/1882_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

