Learn R Programming


title: "Multi-State Adaptive-Dynamic Process Monitoring" author: Gabriel Odom, Ben Barnard, and Melissa Innerst

Overview

We create this package, mvMonitoring, from the foundation laid by Kazor et al (2016). This package is designed to make simulation of multi-state multivariate process monitoring statistics easy and straightforward, as well as streamlining the online process monitoring component.

Installation from CRAN

Install the stable version of this package via

install.packages("mvMonitoring")

Installation of Development Version

Make sure you have the latest version of the devtools package, and pull the package from GitHub.

devtools::install_github("gabrielodom/mvMonitoring")

Load the library after installation by

library(mvMonitoring)

Examples

These are the examples shown in the help files for the mspProcessData(), mspTrain(), mspMonitor(), and mspWarning() functions.

# Generate one week's worth of normal operating (NOC) data recorded at the one-
# minute level
nrml <- mspProcessData(faults = "NOC")
# The state values are recorded in the first column.
n <- nrow(nrml)

# Calculate the training summary, but save five observations for monitoring.
# This function will treat the first 3 days as in control (IC), and then update
# the training window each day.
trainResults_ls <- mspTrain(
  data = nrml[1:(n - 5), -1],
  labelVector = nrml[1:(n - 5), 1],
  trainObs = 4320
)

# While training, we included 1 lag (the default), so we will also lag the
# observations we will test.
testObs <- nrml[(n - 6):n, -1]
testObs <- xts:::lag.xts(testObs, 0:1)
testObs <- testObs[-1,]
testObs <- cbind(nrml[(n - 5):n, 1], testObs)

# Run the monitoring function.
dataAndFlags <- mspMonitor(
  observations = testObs[, -1],
  labelVector = testObs[, 1],
  trainingSummary = trainResults_ls$TrainingSpecs
)

# Alarm check the last row of the matrix returned by the mspMonitor function
mspWarning(dataAndFlags)

Paper Graphics

The R code to build and save the simulation graphics from the paper are in the inst/mspGraphsGrid.R file.

Acknowledgements

This work is supported by the King Abdullah University of Science and Technology (KAUST) Office of Sponsored Research (OSR) under Award No: OSR-2015-CRG4-2582 and by the National Science Foundation PFI:BIC Award No: 1632227.

Copy Link

Version

Install

install.packages('mvMonitoring')

Monthly Downloads

185

Version

0.2.4

License

GPL-2

Issues

Pull Requests

Stars

Forks

Maintainer

Gabriel Odom

Last Published

November 21st, 2023

Functions in mvMonitoring (0.2.4)

oneDay_clean

Real Process Data for Testing
rotate3D

Three-Dimensional Rotation Matrix
threshold

Non-parametric Threshold Estimation
tenDay_clean

Real Process Data for Training
mspTrain

Multi-State Adaptive-Dynamic Process Training
quantile.density

Extract Quantiles from 'density' Objects
processNOCdata

Simulate NOC Observations from a Single-State or Multi-State Process
rotateScale3D

Three-Dimensional Rotation and Scaling Matrix
mvMonitoring

A Package for Multivariate Statistical Process Monitoring
mspSPEPlot

Squared Prediction Error Contribution Plots
mspSubset

Multi-State Subsetting
processMonitor

Adaptive Process Training
dataStateSwitch

Alternate Observations in a Data Frame over States
fault3A_xts

Process Data under a System Signal Amplification
faultFilter

Process Fault Filtering
faultDetect

Process Fault Detection
mspMonitor

Real-Time Process Monitoring Function
faultSwitch

Induce the Specified Fault on NOC Observations
mspContributionPlot

Contribution Plots
mspProcessData

Simulate Normal or Fault Observations from a Single-State or Multi-State Process
mspWarning

Process Alarms
pca

PCA for Data Scatter Matrix
mspT2Plot

T-Squared Contribution Plots
fault1A_xts

Process Data under a System Shift Fault
normal_switch_xts

Process Data under Normal Conditions
fault2A_xts

Process Data under a System Drift Fault