The embodiments herein generally relate to binary classification, and more particularly to systems and methods for constructing generalized one-class support vector machines with jointly optimized hyperparameters thereof.
It is seen that in many real-life applications, only single labeled training class is available to learn and classify, whereas testing phase may involve unknown number of classes. This is a typical problem faced when performing smart analytics, where positive examples (non-anomalous samples) are provided during learning and at the testing phase, negative or anomalous samples are required to be separated out.
With wide-spread adoption of Internet of Things (IoT), data explosion has rendered active domain expert involvement a massively costly affair. Labeling and annotation of training datasets by domain experts in many applications are not feasible in many scenarios. For example, in medical applications, anomaly detection is seemingly the most vital analytics decision. Yet, computational methodology cannot be completely dependent on labeling efforts made by the experts. Many bio-medical applications require sensor signal analysis. For example, Electrocardiogram (ECG) and phonocardiogram (PCG) are required to identify cardio-vascular abnormality. However, in most cases clinically normal training of PCG, ECG signals are much higher than clinically abnormal PCG, ECG signals. Such imbalance problem renders the classification task difficult. Thus, clinical decision making in data-driven computational methods is a challenging task due to scarcity of negative examples.
Embodiments of the present disclosure present technological improvements as solutions to one or more of the above-mentioned technical problems recognized by the inventors in conventional systems.
In an aspect, there is provided a processor implemented method comprising: jointly optimizing hyperparameters (i) kernel co-efficient γ and (ii) rejection rate hyperparameter ν, corresponding to a maximum performance max of a one-class support vector machine (OC-SVM), wherein
max is identified from a matrix
of combinational values of the hyperparameters; and obtaining an optimal non-linear decision boundary based on the jointly optimized hyperparameters (γopt and νopt) for binary classification.
In another aspect, there is provided a system comprising: one or more data storage devices operatively coupled to the one or more processors and configured to store instructions configured for execution by the one or more processors to: jointly optimize hyperparameters (i) kernel co-efficient γ and (ii) rejection rate hyperparameter ν, corresponding to a maximum performance max of a one-class support vector machine (OC-SVM), wherein
max is identified from a matrix
of combinational values of the hyperparameters; and obtain an optimal non-linear decision boundary based on the jointly optimized hyperparameters (γopt and νopt) for binary classification.
In yet another aspect, there is provided a computer program product comprising a non-transitory computer readable medium having a computer readable program embodied therein, wherein the computer readable program, when executed on a computing device, causes the computing device to: jointly optimize hyperparameters (i) kernel co-efficient γ and (ii) rejection rate hyperparameter ν, corresponding to a maximum performance max of a one-class support vector machine (OC-SVM), wherein
max is identified from a matrix
of combinational values of the hyperparameters; and obtain an optimal non-linear decision boundary based on the jointly optimized hyperparameters (γopt and νopt) for binary classification.
In an embodiment of the present disclosure, the one or more hardware processors are further configured to jointly optimize the hyperparameters by: (i) eliminating outliers in the matrix to obtain a steadiness matrix
steady; (ii) computing a steadiness parameter
steady based on maximum performance and standard deviation associated with the steadiness matrix
steady; (iii) diversifying the steadiness parameter
steadyby forming a plurality of matrices
new representing a plurality of regions comprising the steadiness matrix
steady for analyzing new steadiness parameter
newsteady corresponding to each of the four matrices
new; (iv) iteratively computing new steadiness parameter
newsteady based on step (ii) until a stopping criterion is satisfied, wherein the stopping criterion is a ratio of the steadiness parameter
steady and the new steadiness parameter
newsteady less than or equal to ϵ, wherein ϵ represents a deviation coefficient tending to 1; (v) selecting a new steadiness matrix
newsteady corresponding to the new steadiness parameter
newsteady that meets the stopping criterion; and (vi) determining a pair
opt of optimal kernel co-efficient γopt and optimal rejection rate hyperparameter νopt corresponding to a maximum performance element of selected new steadiness matrix
newsteady.
It is to be understood that both the foregoing general description and the following detailed description are exemplary and explanatory only and are not restrictive of the embodiments of the present disclosure, as claimed.
The embodiments herein will be better understood from the following detailed description with reference to the drawings, in which:
It should be appreciated by those skilled in the art that any block diagram herein represent conceptual views of illustrative systems embodying the principles of the present subject matter. Similarly, it will be appreciated that any flow charts, flow diagrams, state transition diagrams, pseudo code, and the like represent various processes which may be substantially represented in computer readable medium and so executed by a computing device or processor, whether or not such computing device or processor is explicitly shown.
Exemplary embodiments are described with reference to the accompanying drawings. In the figures, the left-most digit(s) of a reference number identifies the figure in which the reference number first appears. Wherever convenient, the same reference numbers are used throughout the drawings to refer to the same or like parts. While examples and features of disclosed principles are described herein, modifications, adaptations, and other implementations are possible without departing from the spirit and scope of the disclosed embodiments. It is intended that the following detailed description be considered as exemplary only, with the true scope and spirit being indicated by the following claims.
Before setting forth the detailed explanation, it is noted that all of the discussion below, regardless of the particular implementation being described, is exemplary in nature, rather than limiting.
Traditional classification problem discriminates binary classes. In reality, there exist problems where the discrimination is between one class from any other class(es), particularly, when only a positive example set is available for reliable training. De-corruption of physiological signals like phonocardiogram (PCG) is one of such application cases. For ensuring clinical inference from such signals without human intervention, appropriate de-corruption is of utmost importance. It may be noted that corrupted PCG signals have very less information for appropriate clinical analytics. In fact, contaminated signals may render high mi-classification which is undesirable for clinical inference purposes. Noise or corruption as part of signal space as well as class noise leads to poor classification outcome independent of the power of machine learning techniques employed.
In order to identify noisy PCG signals, learner models have to be adequately trained by both clean and noisy examples. However, training sets mainly contain clean signals only and it is impractical to train with corrupted signals as corrupted signal universe can rarely be captured through a smaller set of negative examples. In such a scenario, optimal decision boundary construction plays a major role in reducing misclassification errors. Accordingly, systems and methods of the present disclosure kernel-optimize one-class support vector machine (OC-SVM) by finding optimal kernel co-efficient γopt and rejection rate hyperparameter νopt.
Referring now to the drawings, and more particularly to
The I/O interface device(s) 106 can include a variety of software and hardware interfaces, for example, a web interface, a graphical user interface, and the like and can facilitate multiple communications within a wide variety of networks N/W and protocol types, including wired networks, for example, LAN, cable, etc., and wireless networks, such as WLAN, cellular, or satellite. In an embodiment, the I/O interface device(s) can include one or more ports for connecting a number of devices to one another or to another server.
The memory 102 may include any computer-readable medium known in the art including, for example, volatile memory, such as static random access memory (SRAM) and dynamic random access memory (DRAM), and/or non-volatile memory, such as read only memory (ROM), erasable programmable ROM, flash memories, hard disks, optical disks, and magnetic tapes. In an embodiment, one or more modules (not shown) of the system 100 can be stored in the memory 102.
In an embodiment, the system 100 comprises one or more data storage devices or memory 102 operatively coupled to the one or more processors 104 and is configured to store instructions configured for execution of steps of the method 200 by the one or more processors 104.
Variance of Radial Basis Function (RBF) is determined by the kernel co-efficient γ. For larger γ, narrower would be the kernel and the corresponding hypersurface is spiky, which means it is zero almost everywhere except at the support vectors, while low value of γ corresponds to larger RBF bandwidth and the hypersurface is very flat. The rejection rate hyperparameter ν represents a lower bound on the support vector and higher bound on number of outliers. Accordingly, OC-SVM has high sensitivity with both γ and ν. Isolated attempts of optimizing performances of the two hyperparameters in the past have resulted in sub-optimal performance of the OC-SVM. The methods and systems of the present disclosure facilitate joint optimization of the two hyperparameters for robust detection of anomalous events to solve the class imbalance problem. It ensures that the decision boundary of the OC-SVM is not biased for positive examples or known classes, thereby enabling construction of a generic application independent classifier.
The present disclosure aims to find optimal hyperparameters of OC-SVM such that a smooth non-linear decision boundary can be formed. The OC-SVM is very sensitive to the hyperparameters and thus hyperparameters impact the decision making process of the classifier. Let, co be the vector perpendicular to the decision boundary of the OC-SVM and p be a bias that parameterizes a hypersurface in the feature space and k be the associated kernel function and
be a reproducing kernel Hilbert space. It is assumed that κ(χ, ·) be bounded for χϵχ, where χ={χi}, i=1,2, . . . n, χiϵ
d, χ be the input training vector space. The constraint objective function of OC-SVM as known is:
subject to:
ωTϕ(χi)≥ρ−ξi, ξi≥0 (2)
where a non-linear function ϕ:χ→ maps vector χ from input vector space χ to feature space
and ξi is the ith positive slack variable that penalizes the objective function, but allows few of the points to lie on the other side of the decision boundary ƒ(χ). ϕ(χi)=[ƒ1(χi), ƒ2(χi), . . . ,]T containing the features ƒi(χi) from
and rejection rate hyperparameter νϵ(0,1). ν is a trade-off parameter that has a significant role in anomaly detection. It is to be noted that without an optimal choice of ν, number of anomaly detection may suffer; ν is an upper bound on the fraction of training vectors external to the constructed decision boundary. The decision of anomaly detection by the decision function ƒ may be sub-optimal optimal and misclassification rate may be high, when
it means that when ν is not close to what is supposed to be the optimal ν(νoptimal), then either high false positive error or high false negative error are reported. If,
some sure anomalies are rejected and when
some normal/non-anomalous samples are detected as anomalies.
Let, α be the Lagrangian multipliers. The solution of equation 1, 2 is equivalent to Lagrange duality form and its Wolfe dual representation form:
From equation (3), the decision function may be represented as:
ƒ(χ)=sgn(Σi=1nαiϕ(χi),ϕ(χ)
−ρ) (4)
Using “kernel trick” of kernel κ that satisfies Mercer's conditions:
κ(χi,χj)=ϕ(χi), ϕ(χj)
(5)
wherein κ is a dot product to fit into the transformed feature space for providing maximum-margin hyperplane.
In the present disclosure, Radial Basis Function (RBF) kernel has been chosen, where the samples are linearly dependent in the feature space
γ denotes the kernel bandwidth. For larger γ, the narrower would be the kernel and the corresponding hypersurface is spiky, which means it is zero in almost everywhere except at the support vectors, while low value of γ corresponds to larger RBF bandwidth and the hypersurface is very flat. It is conclude that in order to achieve optimal OC-SVM decision boundary, dual optimization of γ and ν is required. Subsequently, it is demonstrated that optimal or near-optimal hyperparameters (γopt,νopt) would sufficiently induce less misclassification error.
Let χ={χtrain, χlearn} be divided into disjoint sets χtrain, χlearn for training and model learning respectively, where
The present disclosure provides joint optimization of the hyperparameters γ,ν to find opt that is steady with non-spurious maximum consistent performance over the learning set χlearn.
The range of ν=(0,1) and typically the range of γ:(γlow=2−4, γhigh=25). Let, Δν and Δγ be the quantization level of ν,γ respectively,
be the number of total instances of ν,γ respectively and ρ be the performance parameter which may be F1-score, accuracy, sensitivity, specificity, geometric mean of sensitivity and specificity, the more common F1-score being considered in the present disclosure. ={ρij}i,j=1i=θ
of dimension θν×θγ is formed, where each element corresponds to the performance of OC-SVM on χlearn in terms of ρ. The objective of the present disclosure is to find the ν,γ corresponding to
max=maxiϵν,jϵγ(ρij) such that
max is not inconsistent with its vicinity. In accordance with the present disclosure, the divide-conquer algorithm to find the
max and corresponding
opt={γoptimal,νoptimal} is explained hereinafter.
In accordance with an embodiment of the present disclosure, the one or more processors 104 are configured to jointly optimize, at step 202, hyperparameters (i) kernel co-efficient γ and (ii) rejection rate hyperparameterν, corresponding to a maximum performance max a one-class support vector machine (OC-SVM), wherein
max is identified from a matrix
of combinational values of the hyperparameters.
Input: ={ρij}i,j=1i=θ
consists of all the possible combination of γ, ν within the range with the dimension θν×θγ.
In an embodiment, the step 202 of jointly optimizing the hyperparameters firstly comprises eliminating outliers in the matrix to obtain a steadiness matrix
steady (step 202a). As represented in the high-level flow chart of
, n are to be tuned, where
, n are the distance and density parameters respectively, where
is defined the furthest distance for which a point is density-reachable and n is the minimum number of points required to form a density cluster.
In the present disclosure, is chosen as
=min ((max(θν, θγ))1/3, ζ:s.t.ζ! is nearest to max(θν, θγ)). The consideration of
stems from the fact that
should vary in inverse polynomial for smaller number of elements in the cluster and in negative exponential for larger number of elements in the density cluster.
In accordance with the present disclosure, =┌3σ┐, σ=standard deviation (
), which is determined intuitively considering the quasi-homogeneity in the
pattern.
Then a steadiness parameter steady is computed based on maximum performance and standard deviation associated with the steadiness matrix
steady (step 202b).
ensures that the maximum performance parameter reported in the
steady is simultaneously not too deviated from the mean performance of
steady.
Once the steadiness parameter steady is determined, it is diversified by forming a plurality of matrices
new representing a plurality of regions comprising the steadiness matrix
steady for analyzing new steadiness parameter
newsteady corresponding to each of the plurality of matrices
new (step 202c). For instance, four
new matrices may be formed from the steadiness matrix
steady as represented in the high-level flow chart of
In accordance with the present disclosure, for each of the four new matrices, new steadiness parameter
newsteady is iteratively computed as explained above until a stopping criterion is satisfied, wherein the stopping criterion is a ratio of the steadiness parameter
steady and the new steadiness parameter
newsteady less than or equal to ϵ, wherein ϵ represents a deviation coefficient tending to 1 (step 202d).
In accordance with an embodiment of the present disclosure, a new steadiness matrix newsteady corresponding to the new steadiness parameter
newsteady that meets the stopping criterion is selected (step 202e). A pair
opt of optimal hyperparameters kernel co-efficient γopt and optimal rejection rate hyperparameter νopt corresponding to a maximum performance element of the selected new steadiness matrix
newsteady is then determined (step 202f). Incremental optimization may be performed to derive
opt={γoptimal, νoptimal}. F1—score may be considered as the performance of merit for the optimization. The objective of discovering
opt for a given unbalanced training set is to find that {γ, ν} which is in that center of white region of a heat-map visualization, where in both X-Y direction, boundary of the white region is maximum, which is probable through joint optimization of γ, ν The method of the present disclosure for finding
opt={γoptimal, νoptimal} nearly converges to the center of the white region (equivalent to consistently maximum performance).
opt={γoptimal, νoptimal} derived in accordance with the present disclosure outclasses anomaly detection performance over the other methods known in the art.
Thus in accordance with the present disclosure, the joint optimization of hyperparameters of one-class classifiers can augment clinical utility of automated cardiac condition screening. In the absence of negative examples, one class classification is an effective method to tackle class imbalance and improve clinical decision making outcomes. However, the optimality of parameters of the one-class learner kernel plays a deterministic role on the performance of the learner model of the OC-SVM. Hyperparameters, viz., kernel co-efficient γ and rejection rate hyperparameterν are responsible for the OC-SVM to form a non-linear boundary with training vectors (positive examples). Optimizing the hyperparameters such that steadiness is achieved ensures that one-off incident of excellent performance that can result in over-fitting can be ignored, thereby making the OC-SVM decision boundary smooth and unperturbed by outlier training examples. Since the methods of the present disclosure are independent of the class of the class of training set, the present disclosure facilitates constructing generalized one-class support vector machines.
The written description describes the subject matter herein to enable any person skilled in the art to make and use the embodiments of the present disclosure. The scope of the subject matter embodiments defined here may include other modifications that occur to those skilled in the art. Such other modifications are intended to be within the scope if they have similar elements that do not differ from the literal language of the claims or if they include equivalent elements with insubstantial differences from the literal language.
The scope of the subject matter embodiments defined here may include other modifications that occur to those skilled in the art. Such other modifications are intended to be within the scope if they have similar elements that do not differ from the literal language of the claims or if they include equivalent elements with insubstantial differences from the literal language.
It is, however to be understood that the scope of the protection is extended to such a program and in addition to a computer-readable means having a message therein; such computer-readable storage means contain program-code means for implementation of one or more steps of the method, when the program runs on a server or mobile device or any suitable programmable device. The hardware device can be any kind of device which can be programmed including e.g. any kind of computer like a server or a personal computer, or the like, or any combination thereof. The device may also include means which could be e.g. hardware means like e.g. an application-specific integrated circuit (ASIC), a field-programmable gate array (FPGA), or a combination of hardware and software means, e.g. an ASIC and an FPGA, or at least one microprocessor and at least one memory with software modules located therein. Thus, the means can include both hardware means and software means. The method embodiments described herein could be implemented in hardware and software. The device may also include software means. Alternatively, the embodiments of the present disclosure may be implemented on different hardware devices, e.g. using a plurality of CPUs.
The embodiments herein can comprise hardware and software elements. The embodiments that are implemented in software include but are not limited to, firmware, resident software, microcode, etc. The functions performed by various modules comprising the system of the present disclosure and described herein may be implemented in other modules or combinations of other modules. For the purposes of this description, a computer-usable or computer readable medium can be any apparatus that can comprise, store, communicate, propagate, or transport the program for use by or in connection with the instruction execution system, apparatus, or device. The various modules described herein may be implemented as software and/or hardware modules and may be stored in any type of non-transitory computer readable medium or other storage device. Some non-limiting examples of non-transitory computer-readable media include CDs, DVDs, BLU-RAY, flash memory, and hard disk drives.
Further, although process steps, method steps, techniques or the like may be described in a sequential order, such processes, methods and techniques may be configured to work in alternate orders. In other words, any sequence or order of steps that may be described does not necessarily indicate a requirement that the steps be performed in that order. The steps of processes described herein may be performed in any order practical. Further, some steps may be performed simultaneously.
The illustrated steps are set out to explain the exemplary embodiments shown, and it should be anticipated that ongoing technological development will change the manner in which particular functions are performed. These examples are presented herein for purposes of illustration, and not limitation. Further, the boundaries of the functional building blocks have been arbitrarily defined herein for the convenience of the description. Alternative boundaries can be defined so long as the specified functions and relationships thereof are appropriately performed. Alternatives (including equivalents, extensions, variations, deviations, etc., of those described herein) will be apparent to persons skilled in the relevant art(s) based on the teachings contained herein. Such alternatives fall within the scope and spirit of the disclosed embodiments. Also, the words “comprising,” “having,” “containing,” and “including,” and other similar forms are intended to be equivalent in meaning and be open ended in that an item or items following any one of these words is not meant to be an exhaustive listing of such item or items, or meant to be limited to only the listed item or items. It must also be noted that as used herein and in the appended claims, the singular forms “a,” “an,” and “the” include plural references unless the context clearly dictates otherwise.
It is intended that the disclosure and examples be considered as exemplary only, with a true scope and spirit of disclosed embodiments being indicated by the following claims.
Number | Date | Country | Kind |
---|---|---|---|
201721028487 | Aug 2017 | IN | national |
This U.S. patent application Ser. No. claims priority under 35 U.S.C. § 119 to: Indian Patent Application No. 201721028487, filed on 10 Aug. 2017. The entire contents of the aforementioned application are incorporated herein by reference.