CUSUM

From Wikipedia, the free encyclopedia
Jump to: navigation, search
For the Roman town of Cusum, see Petrovaradin.
CUSUM chart
Originally proposed by E. S. Page
Process observations
Rational subgroup size n = 1
Measurement type Cumulative sum of a quality characteristic
Quality characteristic type Variables data
Underlying distribution Normal distribution
Performance
Size of shift to detect ≤ 1.5σ
Process variation chart
Not applicable
Process mean chart
Center line The target value, T, of the quality characteristic
Upper control limit C_i^+ = max \lbrack 0, x_i - \left ( T + K \right ) + C_{i - 1}^+\rbrack
Lower control limit C_i^- = max \lbrack 0, \left ( T - K \right ) - x_i + C_{i - 1}^-\rbrack
Plotted statistic C_i = \sum_{j=1}^i \bar x_j - T

In statistical quality control, the CUSUM (or cumulative sum control chart) is a sequential analysis technique developed by E. S. Page of the University of Cambridge. It is typically used for monitoring change detection.[1] CUSUM was announced in Biometrika, in 1954, a few years after the publication of Wald's SPRT algorithm.[2]

Page referred to a "quality number" \theta, by which he meant a parameter of the probability distribution; for example, the mean. He devised CUSUM as a method to determine changes in it, and proposed a criterion for deciding when to take corrective action. When the CUSUM method is applied to changes in mean, it can be used for step detection of a time series.

A few years later, George Alfred Barnard developed a visualization method, the V-mask chart, to detect both increases and decreases in \theta.[3]

Method[edit]

As its name implies, CUSUM involves the calculation of a cumulative sum (which is what makes it "sequential"). Samples from a process x_n are assigned weights \omega_n, and summed as follows:

S_0=0
S_{n+1}=\max(0, S_n+x_n-\omega_n)

When the value of S exceeds a certain threshold value, a change in value has been found. The above formula only detects changes in the positive direction. When negative changes need to be found as well, the min operation should be used instead of the max operation, and this time a change has been found when the value of S is below the (negative) value of the threshold value.

Page did not explicitly say that \omega represents the likelihood function, but this is common usage.

Note that this differs from SPRT by always using zero function as the lower "holding barrier" rather than a lower "holding barrier".[1] Also, CUSUM does not require the use of the likelihood function.

As a means of assessing CUSUM's performance, Page defined the average run length (A.R.L.) metric; "the expected number of articles sampled before action is taken." He further wrote:[2]

When the quality of the output is satisfactory the A.R.L. is a measure of the expense incurred by the scheme when it gives false alarms, i.e., Type I errors (Neyman & Pearson, 1936[4]). On the other hand, for constant poor quality the A.R.L. measures the delay and thus the amount of scrap produced before the rectifying action is taken, i.e., Type II errors.

Example[edit]

The following example shows 15 observations of a process with a mean value of X 0 and a standard deviation of 0.5. It can be seen as the value of Z is never greater than 3, so other control charts should not be detected as a failure while using the Cusum 17 that shows the value of SH is greater than 4.

Comments

Performance

Variants[edit]

Cumulative observed-minus-expected plots[1] are a related method.

References[edit]

  1. ^ a b c Grigg et al.; Farewell, VT; Spiegelhalter, DJ (2003). "The Use of Risk-Adjusted CUSUM and RSPRT Charts for Monitoring in Medical Contexts". Statistical Methods in Medical Research 12 (2): 147–170. doi:10.1177/096228020301200205. PMID 12665208. 
  2. ^ a b Page, E. S. (June 1954). "Continuous Inspection Scheme". Biometrika 41 (1/2): 100–115. doi:10.1093/biomet/41.1-2.100. JSTOR 2333009. 
  3. ^ Barnard, G.A. (1959). "Control charts and stochastic processes". Journal of the Royal Statistical Society. B (Methodological) (21, number 2): 239–71. JSTOR 2983801. 
  4. ^ "Sufficient statistics and uniformly most powerful tests of statistical hypotheses". Statistical Research Memoirs I: 113–137. 

Further reading[edit]

External links[edit]