Abstract:Many industrial and engineering processes monitored as times series have smooth trends that indicate normal behavior and occasionally anomalous patterns that can indicate a problem. This kind of behavior can be modeled by a smooth trend such as a spline or Gaussian process and a disruption based on a sparser representation. Our approach is to expand the process signal into two sets of basis functions: one set uses $L_2$ penalties on the coefficients and the other set uses $L_1$ penalties to control sparsity. From a frequentist perspective, this results in a hybrid smoother that combines cubic smoothing splines and the LASSO, and as a Bayesian hierarchical model (BHM), this is equivalent to priors giving a Gaussian process and a Laplace distribution for anomaly coefficients. For the hybrid smoother we propose two new ways of determining the penalty parameters that use effective degrees of freedom and contrast this with the BHM that uses loosely informative inverse gamma priors. Several reformulations are used to make sampling the BHM posterior more efficient including some novel features in orthogonalizing and regularizing the model basis functions. This methodology is motivated by a substantive application, monitoring the water treatment process for the Denver Metropolitan area. We also test the methods with a Monte Carlo study designed around the kind of anomalies expected in this application. Both the hybrid smoother and the full BHM give comparable results with small false positive and false negative rates. Besides being successful in the water treatment application, this work can be easily extended to other Gaussian process models and other features that represent process disruptions.

Consistent second-order discrete kernel smoothing using dispersed Conway-Maxwell-Poisson kernels

Accurate Smoothing Methods for State Estimation of Continuous-Discrete Nonlinear Dynamic Systems

Approximate Bayesian Computation for Smoothing

A Useful Distribution for Fitting Discrete Data: Revival of the Conway–Maxwell–Poisson Distribution

Fast Bivariate Penalized Splines: the Sandwich Smoother

Smoothed Hierarchical Dirichlet Process: A Non-Parametric Approach to Constraint Measures

Estimating Mixture of Gaussian Processes by Kernel Smoothing

Distributed Estimation and Inference for Semi-parametric Binary Response Models

Hybrid Smoothing for Anomaly Detection in Time Series

Kernel Smoothing for Nested Estimation with Application to Portfolio Risk Measurement.

Efficient and fast spline-backfitted kernel smoothing of additive models

Fast Bayesian Inference for Spatial Mean-Parameterized Conway–Maxwell–Poisson Models

Convergence of sequential quasi-Monte Carlo smoothing algorithms

Efficient estimation for a smoothing thin plate spline in a two-dimensional space

Minimax Nonparametric Two-sample Test under Smoothing

Nonlinear state space smoothing using the conditional particle filter

Efficient Computation of Smoothing Splines Via Adaptive Basis Sampling

A Tree-Based Semi-Varying Coefficient Model for the COM-Poisson Distribution

Smooth Simultaneous Confidence Bands for Cumulative Distribution Functions

A statistical and computational theory for robust and sparse Kalman smoothing

Beyond Smoothness: Incorporating Low-Rank Analysis into Nonparametric Density Estimation