International Journal of Data Mining & Knowledge Management Process (IJDKP) Vol.3, No.4, July 2013

IDENTIFICATION OF OUTLIERS IN OXAZOLINES
AND OXAZOLES HIGH DIMENSION MOLECULAR
DESCRIPTOR DATASET USING PRINCIPAL
COMPONENT OUTLIER DETECTION ALGORITHM
AND COMPARATIVE NUMERICAL STUDY OF OTHER
ROBUST ESTIMATORS
Doreswamy1 and Chanabasayya .M. Vastrad2
1

Department of Computer Science, Mangalore University,
Mangalagangotri-574 199, Karnataka, India
1

Doreswamyh@yahoo.com
channu.vastrad@gmail.com

2

ABSTRACT
From the past decade outlier detection has been in use. Detection of outliers is an emerging topic and is
having robust applications in medical sciences and pharmaceutical sciences. Outlier detection is used to
detect anomalous behaviour of data. Typical problems in Bioinformatics can be addressed by outlier
detection. A computationally fast method for detecting outliers is shown, that is particularly effective in
high dimensions. PrCmpOut algorithm make use of simple properties of principal components to detect
outliers in the transformed space, leading to significant computational advantages for high dimensional
data. This procedure requires considerably less computational time than existing methods for outlier
detection. The properties of this estimator (Outlier error rate (FN), Non-Outlier error rate(FP) and
computational costs) are analyzed and compared with those of other robust estimators described in the
literature through simulation studies. Numerical evidence based Oxazolines and Oxazoles molecular
descriptor dataset shows that the proposed method performs well in a variety of situations of practical
interest. It is thus a valuable companion to the existing outlier detection methods.

KEYWORDS
Robust distance, Mahalanobis distance, median absolute deviation, Principal Components Analysis,
Projection-Persuit ,

1. INTRODUCTION
Accurate detection of outliers plays an important role in statistical analysis. If classical statistical
models are randomly applied to data containing outliers, the results can be deceptive at best. In
addition , An outlier is an observation that lies an abnormal distance from other values in a data
from a dataset and their identification is the main purpose of the investigation. Classical methods
based on the mean and covariance matrix are not often able to detect all the multivariate outliers
DOI : 10.5121/ijdkp.2013.3405

75

International Journal of Data Mining & Knowledge Management Process (IJDKP) Vol.3, No.4, July 2013

in a descriptor dataset due to the masking effect [1] , with the result that methods based on
classical measures are inappropriate for general use unless it is certain that outliers are not
present. Erroneous data are usually found in several situations, and so robust methods that detect
or down weight outliers are important methods for statisticians. The objective of this examination
is to provide an detection of outliers, prior to whatever modelling process is visualized.
Sometimes detection of outliers is the primary purpose of the analysis, other times the outliers
need to be removed or down weighted prior to fitting non robust models. We do not distinguish
between the various reasons for outlier detection, we simply aim to inform the analyst of
observations that are considerably different from the majority. Our procedures are therefore
exploratory, and applicable to a wide variety of settings.
Most methods with a high resistance to outliers are computationally exhaustive; not accordingly,
the availability of cheap computing resources has enabled this field to develop significantly in
recent years.
Among other proposals, there currently exist a wide variety of statistical models ranging from
regression to principal components [2] that can include outliers without being excessively
influenced, as well as several algorithms that explicitly focus on outlier detection.
There are several applications where multi-dimensional outlier identification and/or robust
estimation are important. The field of Computational drug discovery, for instance, has recently
received a lot of concentration from statisticians (e.g. the project Bioconductor,
http://www.bioconductor.org). Improvements in computing power have allowed pharmacists to
record and store extremely large databases of information. Such information likely to contain a
fair amount of large errors, however, so robust methods is needed to prevent these errors from
influencing the statistical model. Undoubtedly, algorithms that take a long time to compute are
not perfect or even practical for such large data sets.
Also, there is a further difficulty encountered in molecular descriptor data. The number of
dimensions is typically several orders of importance larger than the number of observations,
leading to a singular covariance matrix, so the majority of statistical methods cannot be applied in
the usual way. As will be discussed later, this situation can be solved through singular value
decomposition but it does require special attention. It can thus be seen that there are a number of
important applications in which current robust statistical models are impractical.

2. MATERIALS AND ALGORITHAMS
2.1 The Data Set
The molecular descriptors of 100 Oxazolines and Oxazoles derivatives [30-31] based H37Rv
inhibitors analyzed. These molecular descriptors are generated using Padel-Descriptor tool [32].
The dataset covers a diverse set of molecular descriptors with a wide range of inhibitory activities
against H37Rv.

2.2 A Brief Overview of Outlier Detection
There are two basic ways to outlier detection – distance-based methods, and projection pursuit.
Distance-based methods aim to detect outliers by computing a measure of how far a particular
point is from the centre of the data. The common measure of “outlyingness” for a data point
x ∈ ℝ , i = 1, … . . , n, is a robust version of the Mahalanobis distance,
76

International Journal of Data Mining & Knowledge Management Process (IJDKP) Vol.3, No.4, July 2013

 =  −     − 

1

Where  is a robust measure of location of the descriptor data set  and  is a robust estimate of
the covariance matrix. Problems faced by distance-based methods include (i) obtaining a reliable
estimate of
in addition to (ii) how large RD should be before a point is classified as outlying.
This focuses the devoted connection between outlier detection and robust estimation – the latter is
required as part of the prior. Retrieving good robust estimators of  and  are essential for
distance-based outlier detection methods. It is then crucial to find a metric (based on  and C) set
apart outliers from regular points. The final separation boundary commonly depends on userspecified penalties for misclassification of outliers as well as regular points(inliers).
2.2.1 Robust Estimation as Main Goal
A simple robust estimate of location is the coordinatewise median. This estimator is not
orthogonally equivariant (does not transform correctly under orthogonal transformations), but if
this property is important, the L1 median should be used alternatively, expressed as
-

! = "#$%&' +‖ − ‖ ,
(∈) *

2

.

where ‖·‖ stands for the Euclidean norm. The L1 median has maximal breakdown point, and a
fast algorithm for its computation is given in [3] .
A simple robust estimate of scale is the MAD (median absolute deviation), expressed for a
dataset {  ,….,- }⊂ ℝ as
34 , … . , -  = 1.4826 · %89 ;: − %89  ;. 3
:



More complex estimators of location and scale are presented by the class of S-estimators [4] ,
expressed as the vector  and positive definite symmetric matrix  that satisfy

min||

?. @.

-

1
+ A9 ⁄B  = DE ,
'
.

9 =  −     −  ,

4

where ρ· is a non-decreasing function on G0, ∞ , and
and bE are tuning constants that can
be as one chosen to support particular breakdown properties. It is commonly easier to work with
ψ = ∂ρ⁄∂d after all
has a root where ρ has minimum.
Distance based algorithms that follow robust estimation as a main goal – without explicit outlier
detection – contain the OGK estimate [5] the minimum volume ellipsoid (MVE) and minimum
covariance determinant (MCD)[6-7]

MCD tries to find the covariance matrix of minimum determinant containing at least h data
points, where h determines the robustness of the estimator; it should be at least n + p − 1⁄2. The
77

International Journal of Data Mining & Knowledge Management Process (IJDKP) Vol.3, No.4, July 2013

MCD and MVE are instances of S-estimators with non-differentiable ρd  because ρd  is either
0 or 1. MCD shows good performance on data sets with low dimension but on larger data sets the
computational stress can be restrictive – the accurate solution requires a combinational search. In
the latter case good starting points need to be acquired, accommodating an approximately correct
method. Equivariant procedures of acquiring these starting points, however, are based on
subsampling methods and the number of subsamples needed to acquire an tolerable level of
accuracy increases swiftly with dimension. Rousseeuw and van Driessen (1999) developed a
faster version of MCD which was a considerable advancement, but is still quite computationally
exhaustive. The OGK estimator [5] is located on pair wise robust estimates of the covariance.
Gnanadesikan and Kettenring (1972) computed a robust covariance estimate for two variables 
and
based on the uniqueness
1
QR, S = T + SU − T − SU ,
4

5

where σ is a robust estimate of the variance. The matrix build from these pair wise estimates will
not necessarily be positive semi definite, so Maronna and Zamar (2002) carried out by
performing an eigen decomposition of this matrix. After all the variables in eigenvector space are
orthogonal, the covariances are zero and it is enough to get robust variance estimates of the data
projected onto each eigenvector direction. The eigenvalues are then restored with these robust
variances, and the eigenvector conversion is applied in reverse to give in a positive semidefinite
robust covariance matrix. If the original data matrix is robustly scaled (every component divided
by its robust variance), the OGK will be scale invariant. This method can be iterated, despite
Maronna and Zamar (2002) find this is not consistently better. Maronna and Zamar (2002) in
addition to find that using weighted estimates is to some extent better, in which case the
observations are weighted according to their robust distances 9 as scaled by the robust covariance
matrix. They concern a weighting function of the form W9 < 9E  where W· is the indicator
function and 9E is taken to be
9E =

YZU [%899 , … . . , 9U 
,
YZU 0.5

6

where χU β is the -quantile of the χU distribution. Observations thus admit full weight except
that their robust distance d > dE , in which case they admit zero weight. Maronna and Zamar
(2002) note that the robust distances
can be swiftly calculated in the eigenvector space without
the requirement for matrix inversion because the _ components are orthogonal in this space. That
is,


U

zb − μdZb f
d = + `
h , i = 1, … … , n,
σdZb f
b.

where zb are the data in the space of eigenvectors , Zb are the
robust location estimate and T is a robust variance estimate.

7
components in this space,

is

78

International Journal of Data Mining & Knowledge Management Process (IJDKP) Vol.3, No.4, July 2013

2.2.2 Explicit Outlier Detection
Continuing robust estimation to outlier detection needs some knowledge of the distribution of
robust distances If
pursues a multivariate normal distribution, the squared classic Mahalanobis
distance (based upon the sample mean and covariance matrix) pursues a χU distribution [8]. In

addition to, if robust estimators
and
are applied to a large data set in which the non-outliers
are normally distributed. Hardin and Rocke (2005) discovered that the squared distances could be
expressed by a scaled - distribution. After all, for non-normal data, it is not clear how the
outlier boundary should be driven to give optimal classification rates. These regards form the
basis for the use of dE by Maronna and Zamar (2002). The complete change of equation (6) helps
the distribution of d ′s be like that of χU for non-normal original data, heading to better results for
the cutoff value than simply χU β.

Hopeful algorithms that focus on detection of outliers contain rocke-Estimator [25], sfastEstimator[26],M-Estimator[27],MVE-Estimator[28],NNC-Estimator
[29]
,BACON[9],
PCDist[10-11] ,sign1[12-13] and sign2[14], Outlier identification using robust (mahalanobis)
distances based on robust multivariate location and covariance matrix [15]. BACON and robust
multivariate location and covariance matrix are distance-based and in an appropriate, direct the
larger part of computational effort toward obtaining robust estimators
and . BACON begins
with a small subset of observations believed to be outlier-free, to which it iteratively adds points
of the current subset. One reason that
that have a small Mahalanobis distance based on T and
builds MCD unreliable for high
is that its contamination bias evolves very swiftly with p [16].
Robust multivariate location and covariance matrix aims to reduce the computational burden by
subdividing the data into cells and running MCD on each cell, i.e. reducing the number of
observations that MCD performs on, with the same number of dimensions. It then associates the
results from each cell to yield a beginning point for an S-estimator [17] that deals a complex
minimization problem to yields a robust estimate of the covariance matrix . S-estimators can
occasionally converge to an inaccurate local solution, so a good beginning point is needed.
Despite, awaiting on MCD in the first stage confines Robust multivariate location and covariance
matrix from investigating large data sets, particularly those of high dimension. It would appear
that methods based on combinatorial search and alternative there of obtains an inherent inability
to investigate large data sets.
2.2.3 Projection Pursuit
In adverse to distance-based procedures are projection pursuit methods [18] which can
equivalently be applied to robust estimation as a main goal or carries towards explicit outlier
detection. The fundamental purpose of projection pursuit methods is to find suitable projections
of the data in which the outliers are readily credible and can thus be downweighted to turn out a
robust estimator, which in turn can be used to detect the outliers. Because they do not consider
the data to begin from a particular distribution but only search for useful projections, projection
pursuit methods are not altered by non-normality and can be mainly applied in diverse data
situations. The penalty for such independence comes in the form of increased computational
stress, because it is not clear which projections should be tested; an exact method would require
that all attainable directions be tested. The primal equivariant robust estimator having a high
breakdown point in arbitrary dimension was the Stahel-Donoho estimator [19-20]. Computer
approximation based on instructions from random subsamples was developed by Stahel (1981),
but without any doubt a large amount of time is necessary to obtain acceptable results. Even
79

International Journal of Data Mining & Knowledge Management Process (IJDKP) Vol.3, No.4, July 2013

though projection pursuit algorithms have the advantage of being appropriate in different data
situations, their computational difficulties seem horrible.

2.3 The High Dimensional Situation
High dimensional data present several problems to classical statistical analysis. As earlier
discussed, computation time increases more rapidly with _ than with '. For combinatorial and
projection pursuit algorithms, this increase is of enough importance to put in question the
practicability of such methods for high dimensional molecular descriptor data. In between the
speedy distance-based methods, computation times of algorithms increase linearly with ' and
cubically with _. This indicates that for very high dimensional molecular descriptor data, the
computational stress of inverting the scatter matrix is nontrivial. This is particularly evident in
iterative methods which require many iterations to converge, since the covariance matrix is
inverted on each iteration. Thus, while the Mahalanobis distance is a very useful metric for
finding correlated multivariate outliers, it is expensive to calculate. Every other methods of
detection of outliers fare even worse, however, usually give up either computational time or
detection accuracy. The MCD is a good case of this in that the precise solution is very accurate
but infeasible to compute for all but small molecular descriptor data sets, whereas a faster
solution can be obtained if
random subsampling is used to produce an approximate solution. It will be investigated in results
and discussion whether the subsampling version of MCD is competitive regarding both accuracy
and computation time. Projection pursuit methods including the Stahel-Donoho estimator have
computation times that increase very rapidly in higher dimensions, and are often at least an order
of magnitude slower than distance based methods since their search for appropriate projections is
an naturally time-consuming task. Thus even if the Mahalanobis distance may be computationally
difficult due to the matrix inversion step, the robust version – RD, as defined in equation (1) – is
an accurate metric for outlier detection and could well be more computationally fair than other
ways.
This is appropriate to several biological applications where the data commonly have orders of
weight more dimensions than observations. This is also the typical situation in chemometrics,
which led to the development of Partial Least Squares (PLS) [21] with other methods.Because the
covariance matrix is singular the robust Mahalanobis distance cannot be calculated. This is not as
big a problem as initially occurs, after all, because the data can be transformed via singular value
decomposition to an equivalent space of dimension ' − 1 [2] and the analysis carried in the same
way as _ < '. However, this situation needs special attention and most outlier algorithms have to
be modified to way this type of high-dimensional molecular descriptor data.
High dimensional molecular descriptor data have various interesting geometrical properties,
discussed in [22]. One such property that is particularly appropriate to outlier detection, is that
high dimensional data points lie near the surface of an enlarging sphere. For example , if ‖‖ is
the norm of  =  , … . . , -  drawn from a normal distribution with zero mean and identity
covariance matrix, then, for ample _ we have
‖‖
_

=

l∑Z:.d: f
_

U

⟶ 1,
80

International Journal of Data Mining & Knowledge Management Process (IJDKP) Vol.3, No.4, July 2013

because the summation includes a YZU distribution In this manner, if the outliers have even a
slightly dissimilar covariance structure from the inliers(non-outliers), they will lie on a different
sphere. This does not help low dimensional outlier detection, but if an algorithm is efficient of
processing high dimensional molecular descriptor data, it should not be too hard to discover the
different spheres of the outliers and inliers.
Principal components are a well known method of dimension reduction, that also advice an way
to detecting high dimensional outliers. Recall that principal components are those directions that
maximize the variance along each component, directed to the circumstance of orthogonality.
Because outliers increase the variance along their respective directions, it appears instinctive that
outliers will come into more visible in principal component space than the original data space; i.e.
, Minimum some of the directions of maximum variance are hopefully to be those that enable the
outliers to “tie-up” more. Exploring for outliers in principal component space should at least, till,
not be any worse than searching for them in the original data space. If the data originally exist in
a high-dimensional space, many of these dimensions likely do not provide significant additional
information and are irrelevant. Principal components thus pick out a handful of highly
informative components (relative to the total number of components), thereby performing a high
degree of dimension reduction and making the data set much more computationally manageable
without losing a lot of information.
For high dimension molecular descriptor data, a ample portion of the smaller principal
components are actually noise [23]. Particularly if _ ≪ ' , the larger part of principal components
will actually be noise and will not add to the total variance. By considering only those principal
components that constitute some agreed level of the total variance, the number of components can
be extensively reduced so that only those components that are truly meaningful are kept. Almost
we found good results using a level of 99%. It can be argued this yields similar results to
transforming the data via SVD to a dimension less than the minimum of ' and _. Thus, in place
of imposing a level of improvement to the variance such as 99%, it would also be likely to choose
the ' − 1 (or fewer) components with the largest variance.
As marked in equation (7) in the OGK method [5] , after dividing by the MAD, the Euclidean
distance in principal components space is for that reason similar to a robust Mahalanobis
distance, because the off-diagonal elements of the scatter matrix are zero. Mahalanobis distance,
because the off-diagonal elements of the scatter matrix are zero.

Hence, it is not essential to invert a p ⨯ p matrix when computing a measure of outlyingness for
every point (i.e. the robust Mahalanobis distance), on the other hand slightly to divide (or
“standardize”) each principal component by its specific variance element. Because eigenvector
decomposition has computational complexity _q to matrix inversion, doing the robust distance
computations in principal component space is not more time-consuming than in common data
space. If this transformation benefits the outliers become more visible and reduces the number of
iterations required to detect them, the result will be a net savings in computational time.

It can be seen that the above approaches are based on simple basic properties of principal
components; this is additional proof of how principal components continue to present appealing
properties to both theoretical and applied statisticians.

81

International Journal of Data Mining & Knowledge Management Process (IJDKP) Vol.3, No.4, July 2013

2.4 Details of the Proposed PrCmpOut Method
The method we propose consists of two basic parts: First one aims to detect location outliers, and
a next that aims to detect scatter outliers. Scatter outliers obtains a different scatter matrix than
the rest of the data, while location outliers are described by a different location parameter. To
begin, it is useful to robustly rescale or sphere each component using the coordinatewise median
and the MAD, presenting to
∗
x,b
=

x,b − meddx,b , … . . , xtb f
MADdx,b , … … , xt,b f

,

j = 1, … … , p .

8

Dimensions with a MAD of zero should be one or the other excluded , or alternatively another
∗
, we calculate a weighted
scale measure has to be used. Beginning with the rescaled data x,b
covariance matrix, from which we compute the eigen values and eigen vectors and hence a semirobust principal components breakdown. We absorb only those eigen vectors or eigen values that
provide to at least 99% of the total variance; call this new dimension _∗ . The remaining
components are generally useless noise and d only deal to complicate any fundamental structure.
For the case _ ≫ ', this also clarifies the singularity problem because _∗ < '. For the _∗ ⨯ _∗
matrix R of eigenvectors we thus get the matrix of principal components as
y =  ∗ R,

9

∗
Rescale these principal components by the median
where  ∗ is the matrix with the elements x,b
and the MAD very much alike to equation (8),

∗
z,b
=

z,b − meddz,b , … . . , ztb f
MADdz,b , … … , zt,b f

,

j = 1, … … , p∗ .

10

Collect y ∗ for the second stage of the algorithm. After the above pre-processing steps, the
location outlier stage is initiated by computing the absolute value of a robust kurtosis measure for
each component as stated in :


∗
∗
∗
}~:
− %89d~:
, … . . , ~f 
1
− 3| ,
{: = | +
∗
∗ 
'
, … . . , ~34d~:
f
.
-

 = 1, … … . , _∗

11

We make use of the absolute value because very much like to Multivariate outlier detection and
robust covariance matrix estimation [24] , both small and large values of the kurtosis coefficient
can be characteristics of outliers. This allows us to assign weights to each component according
to how likely we think it is to disclose the outliers. We use respective weights wb ⁄∑ w to present
a well known scale 0 ≤ wb ≤ 1. If no outliers are there in a given component, we assume the
principal components to be nearly normally distributed similar to the original data, producing a
kurtosis near to zero. Because the being of outliers is likely to originate the kurtosis to become
different than zero, we weight each of the p∗ dimensions proportional to the absolute value of its
kurtosis coefficient. Assigning equal weights to all components (meanwhile the computation of
82

International Journal of Data Mining & Knowledge Management Process (IJDKP) Vol.3, No.4, July 2013

robust Mahalanobis distances) reduce the strength of the discriminatory power because if outliers
undoubtedly stick out in one component, the information in this component will be weakened
unless it is given higher weight. Specifically in principal component space, outliers are more
likely to be clearly visible in one specific component than slightly visible in several components,
so it is important to assign this component higher weight. Because the components are
uncorrelated, we compute a robust Mahalanobis distance employing the distance from the median
(as scaled by the MAD), weighting each component according to the relative weights wb ⁄∑ w
with the kurtosis measure wb expressed in equation (11). The kurtosis measure expressed in
equation (11) helps to guarantee that important information included in a particular component is
not diluted by components which do not separate the outliers.
To complete the first stage of the algorithm, we want to determine how large the robust
Mahalanobis distance should be to get an accurate classification between outliers and nonoutliers. The kurtosis weighting strategy destroys any similarity to a χU∗ distribution that might
have been introduce , so it is not possible use a χU∗ quantile as a separation barrier. Although ,
resembling to Maronna and Zamar (2002) and to equation (6), we got that transforming the robust
distances {  as stated in
9 =  ·

lχU∗ ,E.

%89 , … … . , - 

# & = 1, … … , '

12

helped the experimental distances {d  to have the same median as the assumed distances and thus
bring the earlier somewhat near to χU∗ , where χU∗ ,E. is the χU∗ 50th quantile. We make use of
the adapted biweight function [25] to assign weights to every observation and use these weights
as a measure of outlyingness. The adapted biweight fits into the general scheme of S-estimators
defined by equation (4) and is similar to Tukey’s biweight function other than that ψ begins from
0 at some point 3 from the origin. That is, observations closer than the scaled distance 3 to
location estimate accepts full weight of 1. The ψ function of the adapted biweight is thus given by
9,

U
9−3 U

ψd, c, M = `1 −  B  h ,

0,



which corresponds to the weighting function
{9; B, 3 =







`1 − 

1,

U U

9−3
 h ,
B
0,

0≤9<3

3 ≤9 ≤3+B ,
9 >3+B

0≤9<3

3 ≤ 9 ≤ 3+B ,
9 >3+B

13

14

83

International Journal of Data Mining & Knowledge Management Process (IJDKP) Vol.3, No.4, July 2013

Precisely assigning known non-outliers full weights of one while assigning known outliers
weights of zero results in increased capabilities of the estimators (supported the classifications are
correct), and is also computationally faster. Between these peaks is a subset of points that receive
weights similar to the usual biweight function. To support a high level of robustness, it is best to
be moderate in assigning weights of one, since if any outliers enter the course of action with a
weight of one (or close to it) that will make the other outliers harder to detect due to the masking
effect. Remember that since the principal components have been scaled by the median and MAD
these robust distances measure a weighted distance from the median (using transformed units of
the MAD). We found good experimental values assigning a weight of one to the 1⁄3 of points
acquiring the smallest robust distances. At the other end of the weighting pattern we assign zero
weight to points with 9 > B , where
B = %899 , … . . , 9-  + 2.5 · 349 , … . . , 9-  ,

15

corresponding nearly to classical outlier barriers. Similar to equation (14), the weights for each
observation are calculated by the interpreted biweight function as stated in

{ =

0,




U U

9 ≥ B

9 − 3
 h , 3 < 9 < B ,
`1 − 
B−3


1,
9 ≤ 3


16

where & = 1, … … , ' and 3 is the 33 q  of the quantile of the distances {9 , … . . , 9U . Alternate
weighing strategies were tested ; the benefit of the interpreted biweight is that it allows a subset
of points (that we are quite definite are non-outliers) to be given full weight, while another subset
of points that is likely to contain outliers can be given weights of zero, thereby excluding
unacceptable influence by potential outliers, and a smooth weighting curve for in between points
The weights. {w  from equation (16) are saved; we will use them again at the end of the
algorithm. The second stage of our algorithm is similar to the first except that we don’t use the
kurtosis weighting layout. Principal components addresses on those directions that have large
variance, so it is possibly not unexpected that we find good results searching for scatter outliers in
the semi robust principal component space expressed at the beginning of this section. That is, we
search for outliers in the space determined by y ∗ from equation (10). Earlier , computing the
Euclidian norm for data in principal component space is similar to the Mahalanobis distance in
the original data space, other than that it is faster to compute.


Because the distribution of these distances has not been modified like it was through the kurtosis
weighting layout and assuming we start with normally distributed non-outliers, transforming the
robust distances as before via equation (12) results in a distribution that is fairly close to χU∗ . In
establishing the interpreted biweight as in equation (16), then, acceptable results can be acquired
by assigning 3U equal to the χU∗ 25th quantile and B U equal to the χU∗ 99th quantile. This
distribution is certainly not exactly equal to χU∗ , so there are chances when visual examination of
these distances could head to a better boundary than this automated algorithm. Denote the
weights calculated in this way, {U , & = 1, … . . , '.
84

International Journal of Data Mining & Knowledge Management Process (IJDKP) Vol.3, No.4, July 2013

Lastly we combine the weights from these two steps to compute eventual weights w , i = 1, … . . , n,
as stated in.
w =

w + swU + s
,
1 + sU

17

Where usually the scaling constant ? = 0.25. The justification for presenting ? is that frequently
too many non-outliers accept a weight of 0 in only one of the two steps; setting ? ≠ 0 helps to
certify that the last weight { = 0 only if both steps assign a low weight. Outliers are then
classified as points that have weight w < 0.25. These values means that if one of the two weights
w or wU equals one, the other must be less than 0.0625 for the point x to be classified an
outlier. Or, if w = wU , then this common value must be less than 0.375 for x to be classified as
outlying.
We will hereafter call this algorithm as PrCmpOut. It is beneficial to recap the algorithm in brief.
Stage 1: Detection of location outliers
a) Robustly sphere the data as stated in equation (8). Compute the sample covariance
matrix of the transformed data ∗ .

b) Compute a principal component break down of the semi robust covariance matrix from
Step a and retain only those _∗ eigenvectors whose eigenvalues provide to at least 99% of the
total variance. Robustly sphere the transformed data as in equation (10).

c) Compute the robust kurtosis weights for each component as in equation (11), and thus
weighted norms for the sphered data from Step b. Because the data have been scaled by the
MAD, these Euclidean norms in principal component space are similar to robust Mahalanobis
distances. Convert these distances as stated in equation (12).

d) Select weights { for every robust distance according the adapted biweight in equation

(16), with 3 equal to the 33 q  quantile of the distances {9 , … . . , 9-  and B = %899 , … . . , 9-  +
2.5 · 349 , … . . , 9-  .
Stage 2: Detection of scatter outliers
e) Use the same semi-robust principal component break down computed in Step b and
compute the (unweighted) Euclidean norms of the data in principal component space. Transform
as stated in equation (12) to produce a set of distances for use in Step f.
f) Select weights wU for every robust distance as stated in the adapted biweight in equation
(16) with B U equal to the χU∗ 99th quantile and 3U equal to the χU∗ 25th quantile.

Combining Stage 1 and Stage 2: Use the weights from Steps d and f to select final weights for
all observations as stated in equation (17).

85

International Journal of Data Mining & Knowledge Management Process (IJDKP) Vol.3, No.4, July 2013

2.5 Preliminary Investigation of developed method
As a preliminary step to developing a new outlier detection method, we briefly examined and
compared existing methods to determine possible areas of improvement We selected nine
algorithms that appeared to have good potential for finding outliers: rocke-Estimator [25], sfastEstimator[26],M-Estimator[27],MVE-Estimator[28],NNC-Estimator [29] BACON[9],PCDist[1011], sign1[12-13] and sign2[14]. Somewhat similar to our method, the Sign procedure is also
based on a type of robust principal component analysis. It obtains robust estimates of location
and spread based upon projecting the data onto a sphere. In this way, the effects of outlying
observations are limited since they are placed on the boundary of the ellipsoid and the resulting
mean and covariance matrix are robust. Standard principal components can thus be carried out on
the sphered data without undue influence by any single point (or small subset of points). We
considered attributes p = 10, 20 ,30,40 and critical values α= 0.05,0.1,0.15,0.2 , so for each
parameter combination we carried out 14 simulations with n = 100 observations.
Frequently, results of this type are presented in a 2 × 2 table showing success and failure in
identifying outliers and non-outliers, which we henceforth call inliers. In this paper, we present
the percentage false negatives(FN) followed by the percentage false positives(FP) in the same
cell; we henceforth refer to these respective values as the outlier error rate and inlier error
rate(non outlier error rate). We propose these names because the outlier error rate specifies the
percentage errors recorded within the group of true outliers, and similarly for the inlier error rate.
True outliers among 100 observations are 10,16,18,22,23,25,27,29, 30,47,66,70,72,80,84 90,99
and 100. This is more compact than a series of 2 × 2 table and allows us to easily focus on
bringing the error rates down to zero.
Table 1. Sample 2 × 2 table illustrating the notation used in this paper
Predicted Outliers

Predicted Inliers

True Outliers

a

b

True Inliers

c

d

Q@&8# ### #"@8  =

D

"+D

, W'&8# ### #"@8  =

B

B+9

3. RESULTS AND DISCUSSION
Even if this algorithm was designed basically for computational performance at high dimension,
we compare its performance against other outlier algorithms in low dimension, because a high
dimensional comparison is not achievable. In the following we examine a variety of outlier
configurations in molecular descriptor data. Examination of Table 2 unfolds that PrcCmpOut
performs well at identifying outliers (low false positives), although it has a higher percentage of
false negatives than most of the methods. PrCmpOut has the lowest percentage of false negatives,
often by a tolerable margin, and is a competitive outlier detection method.

86

International Journal of Data Mining & Knowledge Management Process (IJDKP) Vol.3, No.4, July 2013

Table 2 Outlyingness measures: average percentage of outliers that were not identified and
average percentage of regular observations that were declared outliers

We present simulation results in which the dimension was increased from _ = 10 to _ = 40,
based on the mean of 16 simulations at each level. In contrast to the previous simulation
experiment in dimension _ = 10 , in this case we were not able to examine the performance of
the other algorithms since they were not computationally feasible for these dimensions. The
number of observations was held constant at ' = 100, as was the number of outliers at 18. With
increasing dimension PrCmpOut can identify almost all outliers. None of the known methods
experience much success in identifying outliers for small dimensions because geometrically, the
outliers are not very different from the non-outliers. However, as dimension increases, it can be
seen how the outliers separate from the non-outliers and become easier to detect. At p =
20 dimensions, barely more than half of the outliers can be detected, at p = 30 dimensions almost
90% are detected, and at _ = 40 dimensions more than 95% of the outliers are detected.

Figure 1. Average outlier error rate(left) and Average non-outlier error rate(right) for varying
fraction of Outliers
87

International Journal of Data Mining & Knowledge Management Process (IJDKP) Vol.3, No.4, July 2013

The results for the PrCmpOut and rest of the estimators are presented in Figure 1. For an outlier
fraction of 10% all estimators except PCDist perform excellent in terms of outlier error rate (FN)
and detect outliers independently of the percentage of outliers as seen from the left panel of
Figure 1. The average percentage of non-outliers that were declared outliers (FP) differ and
PrCmpOut performs best, followed closely by Sign2 ,PCDist (below 10). With more than 20%
comes next for the rest of the estimators. MVE declares somewhat more than 30% of regular
observations as outliers.Sign1 performs worst with the average non-outlier error rate increasing
with increase of the fraction of non outliers as an outliers.In terms of the outlier error rate
PrCmpOut performs best (11%) followed by NNC, SIGN1, BACON and MVE Estimators error
rate between 13% and 34% and PCDist (89%) is again last.
Exploratory data analysis is often used to get some understanding of the data at hand, with one
important aspect being the possible occurrence of outliers. Sometimes the detection of these
outliers is the goal of the data analysis, more often however they must be identified and dealt with
in order to assure the validity of inferential methods. Most outlier detection methods are not very
robust when the multinormality assumption is violated, and in particular when outliers are
present. Robust distance plots are commonly used to identify outliers. To demonstrate this idea,
we first compute the robust distances based on the sample mean vector and the sample covariance
matrix. Points which have distances larger than χU, are usually viewed as potential outliers,
and so we will label such points accordingly. Figure 2 shows the distance plots of all the methods.

Figure 2. Robust distance plots for the Oxazolines and Oxazoles molecular descriptor data set.
The red points are according to the robust distances outliers.

88

International Journal of Data Mining & Knowledge Management Process (IJDKP) Vol.3, No.4, July 2013

It is important to consider the computational performance of the different outlier detection
algorithms in case of large Oxazolines and Oxazoles molecular descriptor dataset with n = 100.
To evaluate and compare the computational times a simulation experiment was carried out. The
experiment was performed on a Intel core i3 with 6Gb RAM running Windows 7 Professional.
All computations were performed in R i386 2.15.3. The rocke, M,sfast, MVE, NNC, BACON,
PCDist, Sign1, Sign2 and Proposed algorithm PrCmpOut algorithms from were used. Fastest is
PrCmpOut, followed closely by Sign1,Sign2,PCDist ,NNC and BACON. Slowest is sfast-Est
followed by rocke-Est, M-Est and MVE. These computation times are presented graphically in
Figure 3 using a seconds scale

Figure 3. Comparing the computation time of outlier detection algorithms in the presence
of outliers
A fast multivariate outlier detection method is particularly useful in the field of bioinformatics
where hundreds or even thousands of molecular descriptors need to be analyzed. Here we will
focus only on outlier detection among the molecular descriptors, clearly a high dimensional data
set. First, columns with MAD equal to zero were removed, with the remaining columns
investigated for outliers. Figure 4 shows the results from PrCmpOut; The intermediate graphs
provide more insight into the workflow of the method: the kurtosis weights of Step 3 of the
PrCmpOut algorithm are shown in the upper left panel, together with the weight boundaries
described in Step 4, leading to the weights in the upper right panel. Similarly, the distances from
Step 5 and the weights from Step 6 are respectively shown in the left and right panel of the
second row. The lower left panel shows the combined weights (Phase 1 and 2 combined) together
with the outlier boundary 0.25, which results in 0/1 weights (lower right panel). The observations
are clearly visible as multivariate outliers in the intermediate steps of the algorithm

89

International Journal of Data Mining & Knowledge Management Process (IJDKP) Vol.3, No.4, July 2013

Figure 4. The panels show the intermediate steps of the PrCmpOut algorithm (distances and
weights) for analyzing Oxazolines and Oxazoles molecular descriptor dataset.
We compare the performance of PrCmpOut on this data set with the Sign2 method; the other
algorithms are not feasible due to the high dimensionality. Figure 5 shows the distances (left) and
weights (right) as calculated by the Sign2 method.

90

International Journal of Data Mining & Knowledge Management Process (IJDKP) Vol.3, No.4, July 2013

Figure 9. Distances and weights for analyzing the Oxazolines and Oxazoles molecular
descriptor dataset with the Sign method.
A possible explanation for the difficulty experienced by the Sign2 method is a masking effect for
PCA. It is evident that PrCmpOut has better performance than the Sign2 method, which is also
evident from the results in Table 2. We infer that PrCmpOut is a competitive outlier detection
algorithm regarding detection accuracy as well as computation time.

4. CONCLUSIONS
PrCmpOut is a method for detecting outliers in multivariate data that utilizes inherent properties
of principal components decomposition. It demonstrates very good performance for high
dimensional data and through the use of a robust kurtosis measure. In this paper we tested several
approaches for identifying outliers in Oxazolines and Oxazoles molecular descriptor dataset. Two
aspects seem to be of major importance: the computation time and the accuracy of the outlier
detection method. For the latter we used the fraction of false negatives(FN) - outliers that were
not identified - and the fraction of false positives(FP) - non-outliers that were declared as outliers.
It is very fast to compute and can easily handle high dimensions. Thus, it can be extended to
fields such as bioinformatics and data mining where computational feasibility of statistical
routines has usually been a limiting factor. At lower dimensions, it still produces competitive
results when compared to well-known outlier detection methods.

ACKNOLDGEMENTS
We gratefully thank to the Department of Computer Science Mangalore University, Mangalore
India for technical support of this research.

91

International Journal of Data Mining & Knowledge Management Process (IJDKP) Vol.3, No.4, July 2013

REFERENCES
[1]

[2]
[3]
[4]
[5]
[6]

[7]
[8]
[9]
[10]
[11]
[12]
[13]

[14]
[15]
[16]
[17]
[18]
[19]
[20]
[21]
[22]
[23]
[24]
[25]
[26]

Ben-Gal I., “Outlier detection” , In: Maimon O. and Rockach L. (Eds.) Data Mining and Knowledge
Discovery Handbook: A Complete Guide for Practitioners and Researchers," Kluwer Academic
Publishers, 2005, ISBN 0-387-24435-2.
Hubert, M., Rousseeuw, P., Branden, K. V., “Robpca: A new approach to robust principal
components analysis”,Technometrics 2005.47, 64–79
H¨ossjer, O., Croux, C., “Generalizing univariate signed rank statistics for testing and estimating a
multivariate location parameter. “ ,Journal of Nonparametric Statistics 4 1995. 293–308.
Maronna, R., Martin, R., Yohai, V., “Robust Statistics: Theory and Methods.” ,2006. John Wiley.
Maronna, R., Zamar, R., “Robust estimates of location and dispersion for high-dimensional data
sets.”,Technometrics 2002. 44 (4), 307–317.
Rousseeuw, P., In: Grossmann, W., Pflug, G., Vincze, I., Wertz, W. (Eds.), “Multivariate estimation
with high breakdown point. “Mathematical Statistics and Applications. Vol. B. Reidel Publishing
Company, pp. 1985. 283–297.
Rousseeuw, P., van Driessen, K., “A fast algorithm for the minimum covariance determinant
estimator.”, Technometrics 1999. 41, 212–223.
Johnson, R., Wichern, D. W “Applied Multivariate Statistical Analysis”, 4th Edition. , 1998. PrenticeHall.
Billor, N., Hadi, A., Velleman, P., “BACON: Blocked adaptive computationally-efficient outlier
nominators.”, Computational Statistics and Data Analysis 2000. 34, 279–2980 .
Peter Filzmoser,, Valentin Todorov, “Robust Tools for the Imperfect World”, Elsevier September 25,
2012 http://www.statistik.tuwien.ac.at/public/filz/papers/2012IS.pdf
A.D. Shieh and Y.S. Hung,“Detecting Outlier Samples in Microarray Data”, Statistical Applications
in Genetics and Molecular Biology (2009)Vol. 8
Jerome H. Friedman ,“Exploratory Projection Pursuit”, Journal Of The American Statistical
Association Vol. 82, No. 397(Mar., 1987) , 249-266.
Robert Serfling and Satyaki Mazumder ,” Computationally Easy Outlier Detection via Projection
Pursuit with Finitely Many Directions” ,
http://www.utdallas.edu/~serfling/papers/Serfling_Mazumder_final_version.pdf
Locantore, N., Marron, J., Simpson, D., Tripoli, N., Zhang, J., Cohen, K.,. “Robust principal
components for functional data.” Test 1999 8, 1–73.
P. J. Rousseeuw and B. C. Van Zomeren. “Unmasking multivariate outliers and leverage points.”,
Journal of the American Statistical Association. Vol. 85(411), (1990) pp. 633-651.
Adrover, J., Yohai, V. ,”Projection estimates of multivariate location.” ,The Annals of Statistics 30,
2002 1760–1781.
Davies, P.,“Asymptotic behaviour of s-estimates of multivariate location parameters and dispersion
matrices.“,The Annals of Statistics 15 (3), 1987. 1269–1292.
Huber, P.,” Projection pursuit.“, Annals of Statistics 1985 13 (2), 435–475,
Donoho, D., ”Breakdown properties of multivariate location estimators.”, 1982. Ph.D. thesis,
Harvard University.
Stahel, W., “Breakdown of covariance estimators.” , Research Report 1981. 31, Fachgruppe
f¨urStatistik, E.T.H. Z¨urich.
Tenenhaus, M., Vinzi, V. E., Chatelin, Y.-M., Lauro, C., “Pls path modeling”, 48 (1), 2005. 159–205.
Hall, P., Marron, J., Neeman, A., “Geometric representation of high dimension, low sample size
data”, Journal of the Royal Statistical Society, Series B: Statistical Methodology 67, 2005. 427–444
Jackson, J., “A User’s Guide to Principal Components. “,1991.,Wiley & Sons, New York.
Pe˜na, D., Prieto, F., “Multivariate outlier detection and robust covariance matrix estimation.”,
Technometrics 43 (3) 2001., 286–310.
Rocke, D., ,”Robustness properties of s-estimators of multivariate location and shape in high
dimension.”, Annals of Statistics 24 (3), 1996. 1327–1345.
C. Croux, C. Dehon, A.Yadine ,“ON THE OPTIMALITY OF MULTIVARIATE S-ESTIMATORS”,
Scandinavian Journal of Statistics, 38(2), 332-341
92

International Journal of Data Mining & Knowledge Management Process (IJDKP) Vol.3, No.4, July 2013
[27] Leonard A. Stefanski and Dennis D. Boos ,“Calculus of M-Estimation” ,The American Statistician,
February 2002. Vol 56. No. 1.
[28] Christophe Croux, Gentiane Haesbroeck, “Location Adjustment for the Minimum Volume Ellipsoid
Estimator” , Journal of Statistical Computation and Simulation , 72, 585-596
[29] Croux, C., and Van Aelst, C., "Nearest Neighbor Variance Estimation (NNVE): Robust Covariance
Estimation via Nearest Neighbor Cleaning" by Wang, N., and Raftery, A.E., The Journal of the
American Statistical Association, 97, (2002) 1006-1009.
[30] Andrew J. Phillips, Yoshikazu Uto, Peter Wipf, Michael J. Reno, and David R. Williams, “Synthesis
of Functionalized Oxazolines and Oxazoles with DAST and Deoxo-Fluor” Organic Letters 2000 Vol
2 ,No.8 1165-1168
[31] Moraski GC, Chang M, Villegas-Estrada A, Franzblau SG, Möllmann U, Miller MJ.,”Structureactivity relationship of new anti-tuberculosis agents derived from oxazoline and oxazole benzyl
esters”, Eur J Med Chem. 2010 May;45(5):1703-16. doi: 10.1016/j.ejmech.2009.12.074. Epub 2010
Jan 14.
[32] “Padel-Descriptor” http://padel.nus.edu.sg/software/padeldescriptor/

Authors
Doreswamy received B.Sc degree in Computer Science and M.Sc Degree in Computer
Science from University of Mysore in 1993 and 1995 respectively. Ph.D degree in
Computer Science from Mangalore University in the year 2007. After completion of
his Post-Graduation Degree, he subsequently joined and served as Lecturer in
Computer Science at St. Joseph’s College, Bangalore from 1996-1999.Then he has
elevated to the position Reader in Computer Science at Mangalore University in year
2003. He was the Chairman of the Department of Post-Graduate Studies and research
in computer science from 2003-2005 and from 2009-2008 and served at varies
capacities in Mangalore University at present he is the Chairman of Board of Studies and Associate
Professor in Computer Science of Mangalore University. His areas of Research interests include Data
Mining and Knowledge Discovery, Artificial Intelligence and Expert Systems, Bioinformatics, Molecular
modelling and simulation, Computational Intelligence, Nanotechnology, Image Processing and Pattern
recognition. He has been granted a Major Research project entitled “Scientific Knowledge Discovery
Systems (SKDS) for Advanced Engineering Materials Design Applications” from the funding agency
University Grant Commission, New Delhi, India. He has been published about 30 contributed peer
reviewed Papers at national/International Journal and Conferences. He received SHIKSHA RATTAN
PURASKAR for his outstanding achievements in
the year 2009 and
RASTRIYA VIDYA
SARASWATHI AWARD for outstanding achievement in chosen field of activity in the year 2010.
Chanabasayya M. Vastrad received B.E. degree and M.Tech. degree in the year 2001 and
2006 respectively. Currently working towards his Ph.D Degree in Computer Science and
Technology under the guidance of Dr. Doreswamy in the Department of Post-Graduate
Studies and Research in Computer Science , Mangalore University.

93

