Outliers detection in the statistical accuracy test of a pK_a prediction

摘要

The regression diagnostics algorithm REGDIA in S-Plus is introduced to examine the accuracy of pK _a predicted with four programs: PALLAS, MARVIN, PERRIN and SYBYL. On basis of a statistical analysis of residuals, outlier diagnostics are proposed. Residual analysis of the ADSTAT program is based on examining goodness-of-fit via graphical diagnostics of 15 exploratory data analysis plots, such as bar plots, box-and-whisker plots, dot plots, midsum plots, symmetry plots, kurtosis plots, differential quantile plots, quantile-box plots, frequency polygons, histograms, quantile plots, quantile-quantile plots, rankit plots, scatter plots, and autocorrelation plots. Outliers in pK _a relate to molecules which are poorly characterized by the considered pK _a program. Of the seven most efficient diagnostic plots (the Williams graph, Graph of predicted residuals, Pregibon graph, Gray L–R graph, Index graph of Atkinson measure, Index graph of diagonal elements of the hat matrix and Rankit Q–Q graph of jackknife residuals) the Williams graph was selected to give the most reliable detection of outliers. The six statistical characteristics, , and s in pK _a units, successfully examine the specimen of 25 acids and bases of a Perrin’s data set classifying four pK _a prediction algorithms. The highest values and the lowest value of MEP and s and the most negative AIC have been found for PERRIN algorithm of pK _a prediction so this algorithm achieves the best predictive power and the most accurate results. The proposed accuracy test of the REGDIA program can also be extended to test other predicted values, as log P, log D, aqueous solubility or some physicochemical properties. Keywords pK _a prediction - Dissociation constants - Outliers - Residuals - Goodness-of-fit - Williams graph