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

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

Minimum AIC  model: GTR
Minimum AICc model: GTR
Minimum BIC  model: HKY


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
GTR	63	-439.0412	1042.6661	0.5832	0.5832
HKY	59	-447.4079	1046.0552	0.1071	0.6903
GTRI	64	-439.0407	1046.0814	0.1057	0.7960
GTRG	64	-439.0416	1046.0832	0.1056	0.9017
F81	58	-450.5433	1049.0680	0.0237	0.9254
HKYI	60	-447.4074	1049.3432	0.0207	0.9461
HKYG	60	-447.4086	1049.3455	0.0207	0.9668
GTRIG	65	-439.0435	1049.5362	0.0188	0.9856
F81I	59	-450.5429	1052.3253	0.0047	0.9902
F81G	59	-450.5440	1052.3275	0.0047	0.9949
HKYIG	61	-447.4092	1052.6668	0.0039	0.9988
F81IG	60	-450.5447	1055.6178	0.0009	0.9997
SYM	60	-452.1133	1058.7550	0.0002	0.9999
SYMI	61	-452.1127	1062.0738	0.0000	0.9999
SYMG	61	-452.1140	1062.0763	0.0000	1.0000
K2P	56	-461.8286	1065.2128	0.0000	1.0000
SYMIG	62	-452.1151	1065.4302	0.0000	1.0000
JC69	55	-464.6503	1067.6877	0.0000	1.0000
K2PI	57	-461.8281	1068.4097	0.0000	1.0000
K2PG	57	-461.8291	1068.4117	0.0000	1.0000
JC69I	56	-464.6498	1070.8552	0.0000	1.0000
JC69G	56	-464.6509	1070.8574	0.0000	1.0000
K2PIG	58	-461.8296	1071.6406	0.0000	1.0000
JC69IG	57	-464.6516	1074.0566	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTR	63	-439.0412	1004.0824	0.5270	0.5270
GTRI	64	-439.0407	1006.0814	0.1940	0.7210
GTRG	64	-439.0416	1006.0832	0.1938	0.9148
GTRIG	65	-439.0435	1008.0869	0.0712	0.9860
HKY	59	-447.4079	1012.8158	0.0067	0.9927
HKYI	60	-447.4074	1014.8149	0.0025	0.9951
HKYG	60	-447.4086	1014.8172	0.0025	0.9976
HKYIG	61	-447.4092	1016.8185	0.0009	0.9985
F81	58	-450.5433	1017.0867	0.0008	0.9993
F81I	59	-450.5429	1019.0859	0.0003	0.9996
F81G	59	-450.5440	1019.0880	0.0003	0.9999
F81IG	60	-450.5447	1021.0895	0.0001	1.0000
SYM	60	-452.1133	1024.2267	0.0000	1.0000
SYMI	61	-452.1127	1026.2254	0.0000	1.0000
SYMG	61	-452.1140	1026.2279	0.0000	1.0000
SYMIG	62	-452.1151	1028.2302	0.0000	1.0000
K2P	56	-461.8286	1035.6572	0.0000	1.0000
K2PI	57	-461.8281	1037.6562	0.0000	1.0000
K2PG	57	-461.8291	1037.6582	0.0000	1.0000
JC69	55	-464.6503	1039.3006	0.0000	1.0000
K2PIG	58	-461.8296	1039.6593	0.0000	1.0000
JC69I	56	-464.6498	1041.2996	0.0000	1.0000
JC69G	56	-464.6509	1041.3018	0.0000	1.0000
JC69IG	57	-464.6516	1043.3032	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKY	59	-447.4079	1225.7746	0.4991	0.4991
F81	58	-450.5433	1226.4360	0.3585	0.8576
HKYI	60	-447.4074	1231.3832	0.0302	0.8878
HKYG	60	-447.4086	1231.3855	0.0302	0.9180
GTR	63	-439.0412	1231.4791	0.0288	0.9468
F81I	59	-450.5429	1232.0447	0.0217	0.9685
F81G	59	-450.5440	1232.0469	0.0217	0.9902
HKYIG	61	-447.4092	1236.9963	0.0018	0.9920
GTRI	64	-439.0407	1237.0876	0.0017	0.9938
GTRG	64	-439.0416	1237.0894	0.0017	0.9955
F81IG	60	-450.5447	1237.6578	0.0013	0.9968
K2P	56	-461.8286	1237.7876	0.0012	0.9981
JC69	55	-464.6503	1237.8216	0.0012	0.9993
SYM	60	-452.1133	1240.7950	0.0003	0.9996
GTRIG	65	-439.0435	1242.7026	0.0001	0.9997
K2PI	57	-461.8281	1243.3961	0.0001	0.9997
K2PG	57	-461.8291	1243.3981	0.0001	0.9998
JC69I	56	-464.6498	1243.4300	0.0001	0.9999
JC69G	56	-464.6509	1243.4323	0.0001	1.0000
SYMI	61	-452.1127	1246.4032	0.0000	1.0000
SYMG	61	-452.1140	1246.4057	0.0000	1.0000
K2PIG	58	-461.8296	1249.0086	0.0000	1.0000
JC69IG	57	-464.6516	1249.0431	0.0000	1.0000
SYMIG	62	-452.1151	1252.0175	0.0000	1.0000



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

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

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


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

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


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

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


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

