# Patent

## Title
Method and system for semi-supervised anomaly detection with feed-forward neural network for high-dimensional sensor data

## Abstract

One embodiment can provide a system for detecting anomaly for high-dimensional sensor data associated with one or more machines. During operation, the system can obtain sensor data from a set of sensor associated with one or machines, generate a first set of outputs by using a set of clustering models learned in parallel from the unlabeled sensor data and user-provided partial label information, generate a second set of outputs by using a set of feed-forward neural network (FNN) models learned in parallel from the first set of outputs and the unlabeled sensor data, and determine whether an anomaly is present in the operation of the one or more machines based on the second set of outputs and a user-specified threshold.

## Background
BACKGROUND
Field
This disclosure is generally related to a system and method for anomaly detection in operation of industrial machines or systems. More specifically, this disclosure is related to detection of anomaly in high-dimensional sensor data.
Related Art
With the increasing use of Internet of Things (IoT) enabled equipments in numerous industrial applications, e.g., digital manufacturing, a large amount of sensor data is readily available. In particular, in many Cyber Physical System (CPS) applications, sensor data is continuously generated from IoT equipped devices. Detecting anomaly events using the available sensor data is a valuable task for planning, diagnostics, and prognostics. However, in CPS applications, labeled datasets are rarely available or expensive to be obtained. In other words, there is no apriori knowledge about which points are normal or abnormal. As a result, few or no labels are available for sensor data in real-world applications.
Unlike other conventional data, e.g., text, image, and voice, acquiring labeled sensor data is difficult in many CPS applications. Another reason for the difficulty of obtaining labeled sensor data is that the complexity of an application-specific system allows only few experts in the application to assign labels to the sensor data. Furthermore, an abnormal event rarely occurs compared with a normal event which causes a large label-imbalance problem for training data sets.
The above-mentioned problems make it difficult to use supervised learning approaches for anomaly detection in CPS applications. Therefore, to make anomaly detection methods scale well to real-world applications, the anomaly detection method needs to learn models from the unlabeled sensor data and user-provided partial label information.
SUMMARY
According to one embodiment of the present invention, a system for detecting anomalies for high-dimensional senor data associated with one or more machines is provided. During operation, the system can obtain sensor data from a set of sensors associated with one or machines, generate a first set of outputs by using a set of clustering models learned in parallel from the unlabeled sensor data and user-provided partial label information, generate a second set of outputs by using a set of feed-forward neural network (FNN) models learned in parallel from the first set of outputs and the unlabeled sensor data, and determine whether an anomaly is present in the operation of the one or more machines based on the second set of outputs and a user-specified threshold
In a variation on this embodiment, pre-processing the sensor data prior to generating the first set of outputs.
In a variation on this embodiment, the first set of outputs include a set of labels.
In a variation on this embodiment, each label in the set of labels is: a normal label; or an abnormal label.
In a variation on this embodiment, at least one clustering model in the set of clustering models includes a Gaussian Mixture Model (GMM).
In a variation on this embodiment, the user-provided partial label information includes a set of normal labels.
In a further variation on this embodiment, generating the first set of outputs further comprises: training the set of clustering models with the unlabeled sensor data to generate an initial set of outputs, wherein each output in the initial set of outputs is associated with each clustering model in the set of clustering models; computing a set of weights associated with the set of clustering models by using the user-provided partial label information and the initial set of outputs; using the set of weights and the initial set of outputs to compute a set of abnormal label probabilities; and generating the first set of outputs by applying a Bernoulli random sample generator to the set of abnormal label probabilities.
In a further variation on this embodiment, determining whether an anomaly occurs further comprises: computing an average for each sample across the second set of outputs associated the set of FNN models; and detecting an anomaly when the average for each sample is greater than the user-specified threshold.



BRIEF DESCRIPTION OF THE FIGURES
FIG. 1 shows an exemplary semi-supervised anomaly detection data flow architecture, according to one embodiment of the present invention.
FIG. 2 shows an exemplary semi-supervised anomaly detection system architecture, according to one embodiment of the present invention.
FIG. 3A shows the pseudocode for different computation steps in a Semi-Supervised Expectation Maximization algorithm for a Gaussian Mixture Model, according to one embodiment of the present invention.
FIG. 3B shows the pseudocode for Semi-Supervised Expectation Maximization algorithm, according to one embodiment of the present invention.
FIG. 4A presents a flowchart illustrating an exemplary process for performing semi-supervised anomaly detection with a feed-forward neural network, according to one embodiment of the present invention.
FIG. 4B presents a continuation of flowchart in FIG. 4A, illustrating an exemplary process for performing semi-supervised anomaly detection with feed-forward neural network, according to one embodiment of the present invention.
FIG. 5 shows experimental results, according to one embodiment of the present invention.
FIG. 6 illustrates an exemplary computer system that facilitates a semi-supervised anomaly-detection system, according to one embodiment of the present invention.
FIG. 7 illustrates an exemplary apparatus that facilitates a semi-supervised anomaly-detection system, according to one embodiment of the present invention.



In the figures, like reference numerals refer to the same figure elements.
DETAILED DESCRIPTION
The following description is presented to enable any person skilled in the art to make and use the embodiments, and is provided in the context of a particular application and its requirements. Various modifications to the disclosed embodiments will be readily apparent to those skilled in the art, and the general principles defined herein may be applied to other embodiments and applications without departing from the spirit and scope of the present disclosure. Thus, the present invention is not limited to the embodiments shown, but is to be accorded the widest scope consistent with the principles and features disclosed herein.
Overview
The embodiments described herein solve the technical problem of detecting anomalies in operation of industrial machines. More specifically, an anomaly-monitoring-and-detection system can obtain unlabeled sensor data from sensors embedded in the machines and preprocess the sensor data. The system then incrementally learns anomaly models from unlabeled high-dimensional sensor data with a few normal labels provided by users.
In some embodiments of the present invention, the system uses a mixture of generative and discriminative model for generating labels and building a classifier, respectively. The generative model can re-sample labels to train the discriminative model for performing a binary classification of normal or abnormal states. In other words, the system can generate initial labels by using an ensemble Gaussian Mixture models (GMMs) learned from the unlabeled sensor data. The system can also generate probabilistic (soft) labels using the GMM ensemble, weighted by partial label information provided by users, to train an ensemble of Feed Forward Neural Networks (FNN). The system can then detect anomalies in the operation of the machines by applying a user-defined threshold to average of the FNNs outputs.
Semi-Supervised Anomaly Detection System Architecture
FIG. 1 shows an exemplary semi-supervised anomaly detection data flow architecture, according to one embodiment of the present invention. The semi-supervised anomaly detection data flow architecture, can include a sensor data collector 102, a sensor data pre-processor 104, a sensor database 106, an user input 108 including partial normal labels, a label generator 110, a model estimator 112, and an abnormal label detector 114.
Sensor data collector 102 collects sensor data from a set of sensors embedded in one or more machines. Sensor data pre-processor 104 can pre-process the collected sensor data by normalizing the collected sensor data to have zero-mean and unit-variance. Time samples of the pre-processed sensor data xt can be stored in database 106 for further processing. In some embodiments of the present invention, the system can send the normalized sensor data to another computer system implementing label generator 110 to generate estimated labels V. The normalized sensor data includes unlabeled sensor data X.
Label generator 110 can build a set of clustering models using the unlabeled sensor data X to generate an initial set of label outputs. In one embodiment of the present invention, a Gaussian mixture model (GMM) is used as a clustering model. The initial set of label outputs can include normal labels and/or abnormal labels. Label generator 110, can in addition to building a set of clustering models, estimate a Bernoulli probability for each sample using the initial set of label outputs generated from the clustering models and partial normal labels Yh provided by user input 108. The Bernoulli probability can be interpreted as a probabilistic soft label of a sample. Label generator 110 can then generate a set of estimated labels V.
Model estimator 112 can build a set of FNN models by using the unlabeled sensor data X and the set of estimated labels V. Model estimator 112 can generate a set of predicted label outputs, with each predicted label sample for all outputs denoted as zt. The system can then perform anomaly detection by using anomaly detector 114. Anomaly detector 114 can detect an anomaly by computing an average of the set of predicted label outputs for each sample. The system then applies user-specified threshold to this computed average to detect an anomaly.
FIG. 2 shows an exemplary semi-supervised anomaly detection system architecture, according to one embodiment of the present invention. In this example, a semi-supervised anomaly detection system 200 provides architectural details of label generator 110, model estimator 112 and abnormal detector 114 shown in FIG. 1. Semi-supervised anomaly detection system 200 can include a database 202, a GMM model estimator 204, a W estimator 206, an user input module 208, a Q estimator, 210, a Bernoulli random sample generation module 212, a FNN model estimator 214, an anomaly label prediction module 216, and an abnormal event detection module 218.
The following simplified notations for a matrix are used for describing the operation of different elements in the semi-supervised anomaly detection system 200. For an n×m matrix, A=[aij]nm, an ith row vector and a jth column vector are denoted by Ai and Aj respectively, where 1≤i≤n and 1≤j≤m.
Database 202 stores the pre-processed sensor data collected from a set of sensors embedded in one or more machines. The system can pre-process the collected sensor data by normalizing the collected sensor data to have zero-mean and unit-variance. Let the set of sensors corresponding to p feature sensors with n samples each represent certain unknown states of a system of interest. In other words, the sensor data include unlabeled sensor data. Let X denote a normalized feature matrix defined by X=[xij]np, where xij represents the ith sample of the jth feature sensor for 1≤i≤n. The column vector X.j=[x1j, x2j, . . . , xnj] represents data from jth feature sensor. Since all the sensor data are normalized, their variance is Var[X.j]=1 and mean is E[X.j]=0.
GMM model estimator 204 trains M randomly initialized clustering models, e.g., GMMs, in parallel from a set of sub-matrix X to generate an initial label hypothesis independently for each model. The system builds the GMMs by using the following formulations. Let x denote a random variable and a multivariate Gaussian distribution of the random variable be given by: N(x|μk, Σk) with a 1×p mean vector μk and a p×p covariance matrix Σk for kth cluster. Accordingly, a probability density function of the GMM with K centroids can be determined by:



 
  
   
    
     
      g
      ⁡
      (
      x
      )
     
     =
     
      
       ∑
       
        1
        ≤
        k
        ≤
        K
       
      
      
       
        π
        k
       
       ⁢
       
        N
        (
        
         
          x
          |
          
           μ
           k
          
         
         ,
         
          Σ
          k
         
        
        )
       
      
     
    
   
   
    
     (
     1
     )
    
   
  
 



where πk represents a weight probability with 0≤πk≤1 and Σkπk=1. Let Θ=[(πk, μk, Σk)]k=1, . . . , K denote a GMM parameter. The system can train the GMM models to compute the GMM parameter, Θ, by using training data X and user-provided normal labels, Yh. The GMMs can be trained by using a Semi-Supervised Expectation Maximization (SSEM) algorithm. The SSEM algorithm includes a user input update step and represents a modified version of a standard Expectation Maximization algorithm. FIG. 3A shows the pseudocode for different computation steps in a Semi-Supervised Expectation Maximization algorithm for a Gaussian Mixture Model, according to one embodiment of the present invention. FIG. 3B shows the pseudocode for Semi-Supervised Expectation Maximization algorithm, according to one embodiment of the present invention.

In one embodiment of the present invention, the system can estimate the state (i.e., the label) based on the observed sensor measurements x* by computing k*:



 
  
   
    
     
      k
      *
     
     =
     
      arg
       
      
       max
       k
      
       
      
       π
       k
      
      ⁢
      
       N
       (
       
        
         
          x
          *
         
         |
         
          μ
          k
         
        
        ,
        
         Σ
         k
        
       
       )
      
     
    
   
   
    
     (
     2
     )
    
   
  
 



For example, when the system performs anomaly detection, the number of clusters can be set to 2, i.e., K=2. A cluster with a smaller weight matrix, πk, can be assigned to an anomaly class, e.g., cluster 1 can be an anomaly class if π1<π2.

Given X and Yh, the GMM parameter, Θ, may randomly converge to one of many local optimal parameters, thereby resulting in a large number of possible GMM hypothesis. The system finds M hypothesis or estimated labels {um|m=1, 2, . . . , M}, by training M GMMs in parallel, where um denotes labels from mth GMM.
W estimator 206 computes weights W for each GMM's label um. The W estimator can be used for evaluating a reliability of each GMM by using user-provided normal labels Yh, which can be provided by a user via user input module 208. Let y denote ground truth labels, where each sample yn∈{0,1}. Let Yh denote the index of normal labels in y provided by users, where yn=0,∀n∈Yh. The system can compute the weights W as follows:



 
  
   
    
     
      w
      m
      ′
     
     ⁢
     
      
       ∑
       
        n
        ∈
        
         Y
         h
        
       
      
      
       
        I
        ⁡
        (
        
         
          y
          n
         
         ,
         
          u
          
           m
           ⁢
           n
          
         
        
        )
       
       /
       
        
        
         Y
         h
        
        
       
      
     
    
   
   
    
     (
     3
     )
    
   
  
 




 
  
   
    
     W
     =
     
      
       [
       
        
         w
         1
         ′
        
        ,
        
         w
         2
         ′
        
        ,
        ...
            
        ,
        
         w
         M
         ′
        
       
       ]
      
      /
      
       
        ∑
        m
       
       
        w
        m
        ′
       
      
     
    
   
   
    
     (
     4
     )
    
   
  
 



where l(x,y) is an identity function with output 1 if x=y, and 0 otherwise; and ∥x∥ is the cardinality of x.

Q estimator 210 computes abnormal probability q by computing a weighted average of labels {u1, u2, . . . , uM} by applying the weights W. In one embodiment, the system can compute the abnormal probability q=(q1, q2, . . . , qN) as:

q=(UWT)°1h  (5)

where U is a matrix of GMM labels {u1, u2, . . . , uM}, ° is a Hadamard product, and 1h is a N×1 binary matrix, where 1h[n]=0 for n∈Yh and 1 otherwise.

Bernoulli random sample generation module 212 can generate a set of Bernoulli random labels {v1, v2, . . . , vL} from the abnormal probability q:

vnl ̃Bern(qn)  (6)

In one embodiment, the Bernoulli random labels can be used as training labels for nth sample in lth FNN model in the FNN model estimator 214.

FNN model estimator 214 can use a back propagation algorithm to create an ensemble of FNN models trained from labeled sensor data (X,vl) where vl represents the Bernoulli random labels for the lth FNN model. The system can train L FNN models in parallel with their respective training data set (X,vl) and observed sensor data xt to generate a set of predicted labels {zt1, zt2, . . . , ztL}.
Anomaly label prediction module 216 can output a final abnormal score {circumflex over (z)}t by computing an average of the FNN model estimator 214 outputs {zt1, zt2, . . . , ztL}, where 0≤{circumflex over (z)}t≤1. Abnormal event detection module 218 can detect abnormal events from each sample, which is defined by:



 
  
   
    
     
      s
      t
     
     =
     
      {
      
       
        
         
          0
          ,
         
        
        
         normal
        
        
         
          
           if
           ⁢
              
           
            
             z
             ^
            
            t
           
          
          ≤
          
           z
           th
          
         
        
       
       
        
         
          1
          ,
         
        
        
         abnormal
        
        
         
          
           if
           ⁢
              
           
            
             z
             ^
            
            t
           
          
          >
          
           z
           th
          
         
        
       
      
     
    
   
   
    
     (
     15
     )
    
   
  
 



where st represents a state of each sample and 0≤zth≤1 is a user-defined threshold.

FIG. 4A presents a flowchart illustrating an exemplary process for performing semi-supervised anomaly detection with a feed-forward neural network, according to one embodiment of the present invention. During operation, the system may obtain unlabeled high-dimensional sensor data stored in a sensor database (operation 402). The sensor data is associated with a plurality of sensors embedded in one or more machines in a factory layout.
The system can pre-process the unlabeled sensor data to provide normalized unlabeled sensor data with zero-mean and unit-variance (operation 404). Subsequent to performing sensor data pre-processing, the system can build a set of clustering models, e.g., GMMs. In other words, the system can use the unlabeled sensor data to train a set of M clustering models to generate an initial set of normal and/or abnormal labels independently (operation 406).
In one embodiment of the present invention, the system can compute weights for each clustering model by using the initial set of labels and user-provided partial normal labels (operation 408). These weights represent reliability of each GMM model in the semi-supervised anomaly detection process. The system then uses these weights to compute an abnormal label probability (operation 410).
FIG. 4B presents a continuation of flowchart in FIG. 4A, illustrating an exemplary process for performing semi-supervised anomaly detection with a feed-forward neural network, according to one embodiment of the present invention. The system subsequently generates the Bernoulli random labels by using the abnormal label probability (operation 412). Each of these Bernoulli random labels can be interpreted as a probabilistic soft label of a sample.
The system then uses the Bernoulli random labels and the unlabeled sensor data to build an ensemble of FNN models. Net the system outputs a set of predicted labels for the ensemble of the FNN models (operation 414). Finally, a collective decision is made for an anomaly by computing an average of the FNN outputs. The system then applies a user-defined threshold to the computed average value for each sample to detect an anomaly (operation 416).
FIG. 5 shows experimental results, according to one embodiment of the present invention. For evaluating the performance of the anomaly detection process presented in FIG. 4A and FIG. 4B, two different experiments have been performed on 5 different data sets. In the first experiment, 5 GMM ensemble and one FNN were generated without normal label inputs. In the second experiment, the 5 GMM ensemble and one FNN were generated with 5% normal label input. In both the experiments an Area Under Curve (AUC) metric has been used to evaluate the performance of the anomaly detection process. In both the experiments an average AUC has been computed for the final anomaly score {circumflex over (z)}t and for an estimated label q generated by the 5 GMM ensemble. FIG. 5 shows that by using a normal label input of 5% a significant AUC gain over the experiment without normal label input can be obtained. Further, a comparison of the AUC values for the final anomaly score {circumflex over (z)}t and the estimated label q, depicts that there is a significant AUC gain by using the FNN ensemble. FIG. 5 also shows that the anomaly detection system and method of the present invention performs better than a conventional anomaly detection method iForest (Isolation Forest is a model-based method for isolating anomalies).
Exemplary Computer System and Apparatus
FIG. 6 illustrates an exemplary computer system that facilitates a semi-supervised anomaly-detection system, according to one embodiment of the present invention. In this example, a computer system 600 can include a processor 602, a memory 604, and a storage device 606. Computer system 600 can be coupled to a display device 610, a keyboard 612, and a pointing device 614, and can also be coupled via one or more network interfaces to network 408. Storage device 606 can store an operating system 616, and an anomaly detection system 618.
Anomaly-detection system 618 can include instructions, which when executed by computer system 600 can cause computer system 400 to perform methods and/or processes described in this disclosure. Anomaly detection system 618 can also include instructions for receiving sensor data associated with the one or more sensors (sensor data receiving module 620), instructions for performing sensor data pre-processing (sensor data pre-processing module 622). Furthermore, anomaly detection system 618 can include instructions for generating labels using an ensemble of GMMs (label generation module 624), and instructions for building an ensemble of FNN models (model estimation module 626). Anomaly detection system 618 can further include instructions for detecting anomaly events (anomaly detection module 628).
FIG. 7 illustrates an exemplary apparatus that facilitates a semi-supervised anomaly detection system, according to one embodiment of the present invention. In this example, an apparatus 700 can include a plurality of units or apparatuses that may communicate with one another via a wired or wireless communication channel. Apparatus 700 may be realized using one or more integrated circuits, and may include fewer or more units or apparatuses than those shown in FIG. 7. Further, apparatus 700 may be integrated in a computer system, or realized as a separate device that is capable of communicating with other computer systems and/or devices. Specifically, apparatus 700 can include units 702-712, which perform functions or operations similar to modules 620-628 of computer system 600 of FIG. 6. Apparatus 700 can include: a sensor data receiving unit 702, a sensor data pre-processing unit 704, a label generation unit 706, a model estimation unit 708, and an anomaly detection unit 710. Apparatus 700 can further include a communication unit 712.
In general, embodiments of the present invention provide a method and system to detect anomaly for high-dimensional and unlabeled sensor data using a semi-supervised anomaly detecting with FNN. Detecting anomalies in sensor date associated with factory machines has been used as an example. In practice, this solution is not limited to detecting anomalies associated with factory machines. It can also be used to detect anomalies in other types of equipment or machinery.
The methods and processes described in the detailed description section can be embodied as code and/or data, which can be stored in a computer-readable storage medium as described above. When a computer system reads and executes the code and/or data stored on the computer-readable storage medium, the computer system performs the methods and processes embodied as data structures and code and stored within the computer-readable storage medium.
Furthermore, the methods and processes described above can be included in hardware modules or apparatus. The hardware modules or apparatus can include, but are not limited to, application-specific integrated circuit (ASIC) chips, field-programmable gate arrays (FPGAs), dedicated or shared processors that execute a particular software module or a piece of code at a particular time, and other programmable-logic devices now known or later developed. When the hardware modules or apparatus are activated, they perform the methods and processes included within them.
The foregoing descriptions of embodiments of the present invention have been presented for purposes of illustration and description only. They are not intended to be exhaustive or to limit the present invention to the forms disclosed. Accordingly, many modifications and variations will be apparent to practitioners skilled in the art. Additionally, the above disclosure is not intended to limit the present invention. The scope of the present invention is defined by the appended claims.

## Claims

1. A method for detecting an anomaly in operation of one or more machines, the method comprising:
obtaining a plurality of unlabeled sensor data samples from one or more sensors associated with the one or more machines;
training a plurality of randomly initialized clustering models in parallel using the unlabeled sensor data and user-provided partial label information including a set of normal labels to generate a set of estimated labels, wherein
training each randomly initialized clustering model comprises evaluating a reliability of the randomly initialized clustering model, wherein
evaluating the reliability of each randomly initialized clustering model comprises computing a weight associated with the randomly initialized clustering model based on a matching rate between estimated labels outputted by the randomly initialized clustering model and user-provided ground-truth labels, and wherein
the estimated labels are determined by calculating




 
  
   
    k
    *
   
   =
   
    arg
    ⁢
      
    
     
      max
        
     
     k
    
    ⁢
    
     π
     k
    
    ⁢
    
     N
     ⁡
     (
     
      
       x
       *
      
      ⁢
      
       
        ❘
        "\[LeftBracketingBar]"
       
       
        
         μ
         k
        
        ,
        
         ∑
         k
        
       
      
     
     )
    
   
  
  ,
 



with x* being the unlabeled sensor data;
computing, for each unlabeled sensor data sample, an abnormal probability based on the estimated labels and the weight associated with each randomly initialized clustering model;
applying a random sample generator to generate multiple sets of labeled training samples based on the abnormal probability;
training a set of feed-forward neural network (FNN) models in parallel, wherein a respective FNN model of the set of trained FNN models is trained using a corresponding set of labeled training samples;
obtaining, for an observed sensor data sample, a set of predicted labels outputted by the set of trained FNN models, wherein each trained FNN model outputs a predicted label;
computing an average of the set of predicted labels outputted by the set of trained FNN models; and
determining whether an anomaly is present in the operation of the one or more machines based on whether the average of the set of predicted labels is greater than a user-specified threshold.

2. The method of claim 1, further comprising pre-processing the sensor data samples prior to generating the set of estimated labels.

3. The method of claim 1, wherein each label in the set of estimated labels is:
a normal label; or
an abnormal label.

4. The method of claim 1, wherein at least one randomly initialized clustering model in the plurality of randomly initialized clustering models includes a Gaussian Mixture Model (GMM).

5. The method of claim 1, wherein computing the abnormal probability further comprises computing a weighted average of labels outputted by the plurality of randomly initialized cluttering models, and
wherein applying the random sample generator to generate the multiple sets of labeled training samples further comprises applying a Bernoulli random sample generator to generate the multiple sets of labeled training samples based on the abnormal probability.

6. An apparatus for detecting an anomaly in operation of one or more machines, comprising: one or more processors;
a set of sensors embedded in the one or more machines; and memory storing instructions that, when executed by the one or more processors, cause the apparatus to:
obtain a plurality of unlabeled sensor data samples from one or more sensors associated with the one or more machines;
train a plurality of randomly initialized clustering models in parallel using the unlabeled sensor data samples and user-provided partial label information including a set of normal labels to generate a set of estimated labels, wherein
training each randomly initialized clustering model comprises evaluating a reliability of the randomly initialized clustering model, wherein
evaluating the reliability of each randomly initialized clustering model comprises computing a weight associated with the randomly initialized clustering model based on a matching rate between estimated labels outputted by the randomly initialized clustering model and user-provided ground-truth labels, and wherein
the estimated labels are determined by calculating




 
  
   
    k
    *
   
   =
   
    arg
    ⁢
      
    
     
      max
        
     
     k
    
    ⁢
    
     π
     k
    
    ⁢
    
     N
     ⁡
     (
     
      
       x
       *
      
      ⁢
      
       
        ❘
        "\[LeftBracketingBar]"
       
       
        
         μ
         k
        
        ,
        
         ∑
         k
        
       
      
     
     )
    
   
  
  ,
 



with x* being the unlabeled sensor data;
computing, for each unlabeled sensor data sample, an abnormal probability based on the estimated labels and the weight associated with each randomly initialized clustering model;
apply a random sample generator to generate multiple sets of labeled training samples based on the abnormal probability;
train a set of feed-forward neural network (FNN) models in parallel, wherein a respective FNN model of the set of trained FNN models is trained using a corresponding set of labeled training samples;
obtain, for an observed sensor data sample, a set of predicted labels outputted by the set of trained FNN models, wherein each trained FNN model outputs a predicted label;
compute an average of the set of predicted labels outputted by the set of trained FNN models; and
determine whether an anomaly is present in the operation of the one or more machines based on whether the average of the set of predicted labels is greater than a user-specified threshold.

7. The apparatus of claim 6, wherein obtaining the sensor data samples further comprises: pre-processing the sensor data samples.

8. The apparatus of claim 6, wherein each label in the set of estimated labels is:
a normal label; or
an abnormal label.

9. The apparatus of claim 6, wherein at least one randomly initialized clustering model in the plurality of randomly initialized clustering models includes a Gaussian Mixture Model (GMM).

10. The apparatus of claim 6, wherein computing the abnormal probability further comprises computing a weighted average of labels outputted by the plurality of randomly initialized cluttering models, and wherein applying the random sample generator to generate the multiple sets of labeled training samples further comprises applying a Bernoulli random sample generator to generate the multiple sets of labeled training samples based on the abnormal probability.

11. A non-transitory computer-readable storage medium storing instructions that when executed by a computer cause the computer to perform a method for detecting an anomaly in operation of one or more machines, the method comprising:
obtaining a plurality of unlabeled sensor data sample from one or more sensors associated with the one or more machines;
training a plurality of randomly initialized clustering models in parallel using the unlabeled sensor data and user-provided partial label information including a set of normal labels to generate a set of estimated labels, wherein
training each randomly initialized clustering model comprises evaluating a reliability of the randomly initialized clustering model, wherein
evaluating the reliability of each randomly initialized clustering model comprises computing a weight associated with the randomly initialized clustering model based on a matching rate between estimated labels outputted by the randomly initialized clustering model and user-provided ground-truth labels, and wherein
the estimated labels are determined by calculating




 
  
   
    k
    *
   
   =
   
    arg
    ⁢
      
    
     
      max
        
     
     k
    
    ⁢
    
     π
     k
    
    ⁢
    
     N
     ⁡
     (
     
      
       x
       *
      
      ⁢
      
       
        ❘
        "\[LeftBracketingBar]"
       
       
        
         μ
         k
        
        ,
        
         ∑
         k
        
       
      
     
     )
    
   
  
  ,
 



with x* being the unlabeled sensor data;
computing, for each unlabeled sensor data sample, an abnormal probability based on the estimated labels and the weight associated with each randomly initialized clustering model;
applying a random sample generator to generate multiple sets of labeled training samples based on the plurality of unlabeled sensor data samples and the set of estimated labels;
training a set of feed-forward neural network (FNN) models in parallel, wherein a respective FNN model of the set of trained FNN models is trained using a corresponding set of labeled training samples;
obtaining, for an observed sensor data sample, a set of predicted labels outputted by the set of trained FNN models, wherein each trained FNN model outputs a predicted label;
computing an average of the set of predicted labels outputted by the set of trained FNN models; and
determining whether an anomaly is present in the operation of the one or more machines based on whether the average of the set of predicted labels is greater than a user-specified threshold.

12. The non-transitory computer-readable storage medium claim 11, wherein at least one randomly initialized clustering model in the plurality of randomly initialized clustering models includes a Gaussian Mixture Model (GMM).

13. The non-transitory computer-readable storage medium claim 11, wherein each label in the set of estimated labels is one of:
a normal label; and
an abnormal label.

14. The non-transitory computer-readable storage medium claim 11, wherein computing the abnormal probability further comprises computing a weighted average of labels outputted by the plurality of randomly initialized cluttering models, and wherein applying the random sample generator to generate the multiple sets of labeled training samples further comprises applying a Bernoulli random sample generator to generate the multiple sets of labeled training samples based on the abnormal probability.