Download V. Control Charts

V. Control Charts
A. Overview
Consider an injection molding process for a pen barrel.
The goal of this process: To produce barrels whose true mean outside diameter
is μ0.
Consider a plot of the mean, μ(t), over time.
What are our options?
• Caveat emptor. “Let the buyer beware.”
Basic attitude: The customer will buy whatever we make.
Therefore, let the process mean drift!
• Let μ(t) wander and perform 100% inspection expansion.
Basic idea: Inspect out the bad barrels.
Problems:
• Very expensive in terms of inspector's time.
• Lots of wasted barrels.
• Even 100% inspection is not 100% effective!
• Use a statistical sampling plan.
Basic idea: Use statistics to suggest ways to cut down the inspection costs.
Advantage: Generally cheaper than 100% inspection.
Problems:
• Lots of wasted barrels.
• Definitely not 100% effective.
• Statistical Process Control
Basic idea: Constantly monitor the process; adjust the process once a drift is
detected.
Advantage: Stops problems before they become serious.
Drawback: Does not directly deal with the sources of problems.
• Create processes which do not produce unacceptable product.
Basic idea: Identify sources of problems; set up safeguards to prevent their
occurrence.
Primary statistical tool: Experimental Designs (Chapter 7)
Consider monitoring a process.
Recall, our goal is to keep μ(t) at some target value, μ0.
Note: μ0 may or may not be the “ideal” value for our characteristic of interest.
The first goal in any quality improvement program is to get the process stable.
Only after the process is stable, can we begin to move μ0 to its “ideal” value.
When μ(t) = μ0, we say that the process is in control.
The process is assumed to remain in control until acted upon by an assignable
cause.
The assignable cause shifts μ(t) to some other value μ1 and the process is said to be
out of control.
Our goal, then, is to develop a method for monitoring μ(t).
Problem: Do we ever see μ(t)?
No, but we do see estimates of μ(t)!
Therefore, what should we do?
Assume that we are concerned with the mean of the characteristic of
interest, and that the characteristics variance is known.
What seems reasonable?
Step 1: Take samples at regular intervals.
Step 2: Perform hypothesis tests of the form
H0: μ = μ0
H1: μ ≠ μ0
process is in control
process is out of control
Step 3: The appropriate test statistic
y  0
Z
/ n
Step 4: The critical region
We reject if H0
| Z | z / 2
What does this critical region mean?
We conclude the process is out of control if
Z   z / 2
or
Z  z / 2
Equivalently, we act as if the process is in control (we have insufficient evidence
to conclude the process is out of control) if
 z / 2  Z  z / 2
y  0
 z / 2 
 z / 2
/ n
 z / 2   / n  y   0  z / 2   / n




 0  z / 2   / n   y   0  z / 2   / n 
As a result, we conclude that the process is out of control if either




y  0  z / 2   / n
or
y  0  z / 2   / n
We call
•
•
0  z / 2   / n 
0  z / 2   / n 
the upper control limit (UCL) and
the lower control limit (LCL).
What about  ?
Traditionally, z  3 , which corresponds to an  =.0027.
/2
We summarize this sequence of tests by a graphical procedure.
Any observed sample mean above or below the control limits is considered
evidence that the process is out of control.
If a “out of control” condition signals, then we must search for the assignable
cause, correct it, and restart the chart.
Note: It is possible for a “false alarm” to occur.
False alarm: The chart indicates that the process is out of control when it is not.
What is P(false alarm)?

Let N be the number of samples until the process signals an out of control
condition.
If the process is in control and we use z  3 ,
/2
1
1
E(N )  
 370.4
 0.0027
It is important to note: A control chart will signal on out of control situation at
some point in time even if the process remains in control the entire time.
z / 2  3
was chosen to insure that the process will signal very infrequently when
it is in control.
On the other hand, if μ1 is significantly different from μ0, then the chart will signal
quickly.
How quickly depends upon
1   0
/ n
The general form of our control charts is
UCL   0  3 ˆ
LCL   0  3 ˆ
CL   0
 θ represents the characteristic of interest we wish to monitor.
We shall consider the cases of:
• sample mean,
• sample range,
• sample proportion
• sample “counts”.
 θ0 represents the target value for θ.
   represents the standard error for our estimator of θ.
ˆ
B.
X - and R-Charts
Consider monitoring the ash contents for a pencil lead process.
Suppose we take a sample of 5 lots each shift.
1. The R-Chart
First, consider monitoring the within sample variability.
Thus, at each time period, we need to estimate σ.
Let yij be the jth observation (j = 1, …, n) taken at the ith interval (i = 1,2, …, ).
Let Ri be the sample range of these n observations.
Ri  yi ( n )  yi (1)
 largest - smallest
We note that Ri is a measure of the “spread” or variability at the ith times interval.
Also, it is easy to calculate.
If our data follow a normal distribution, then Ri is an estimate of a multiple of σ.
Typically, we do not know the true in control value for σ.
Problem: How should we estimate it?
Solution: A “base period”.
We use an initial base period of m samples.
m may be as small as 20 intervals.
Personally, I recommend 40-60 intervals.
We then base our control limits for the sample ranges on the average sample
range from this base period, R , given by
1 m
R   Ri
m i 1
The control limits for the range are:
UCL  D4 R
LCL  D3 R
where D3 and D4 are appropriate constants. (see the text).
Note: for n ≤ 6, D3 = 0.
For our example, the data follow.
Let the first 20 samples form the base period.
Example: the data are:
Let the first 20 samples
form the base period.
OBS
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
ASH1
41.3
42.4
41.0
42.9
43.4
42.1
43.2
42.3
42.6
42.6
42.0
42.3
43.2
42.0
43.7
42.7
40.3
41.0
41.0
40.7
41.9
41.3
41.7
40.4
42.0
42.5
42.9
42.6
42.3
42.1
ASH2
42.1
41.7
41.6
41.7
43.0
42.6
42.7
42.5
42.5
40.2
41.3
42.4
40.7
42.2
42.6
42.5
42.6
42.3
40.3
40.8
41.8
41.7
42.7
42.4
42.8
42.6
41.8
42.2
41.6
42.2
ASH3
41.5
40.6
41.8
41.5
42.7
42.4
40.9
43.0
42.5
42.0
40.3
41.7
41.0
42.0
43.7
42.1
42.6
42.2
40.5
40.5
41.7
42.8
41.0
42.9
42.3
42.0
42.4
42.3
41.6
41.8
ASH4
41.5
40.2
41.2
42.5
42.2
43.1
40.8
42.5
42.2
41.4
41.7
42.5
41.4
42.0
41.3
42.2
42.1
42.3
40.1
41.5
41.6
42.4
40.4
42.1
42.6
42.0
42.4
42.3
42.1
42.2
ASH5
41.9
43.3
41.6
40.3
42.3
42.8
42.3
42.2
42.3
42.1
41.2
42.5
42.0
42.5
42.7
42.3
42.4
40.9
40.5
42.3
41.4
41.5
42.3
42.7
42.5
42.6
42.5
42.3
43.0
43.3
MEAN
41.66
41.64
41.44
41.78
42.72
42.60
41.98
42.50
42.42
41.66
41.34
42.28
41.66
42.14
42.80
42.36
42.00
41.74
40.48
41.16
41.68
41.94
41.62
42.10
42.44
42.34
42.40
42.34
42.12
42.32
RANGE
0.8
3.1
0.8
2.6
1.2
1.0
2.4
0.8
0.4
2.4
1.7
0.8
2.5
0.5
2.4
0.6
2.3
1.4
0.9
1.8
0.5
1.5
2.3
2.5
0.8
0.6
1.1
0.4
1.4
1.5
In this example,
n=5
D4 = 2.115
D3 = 0
The control limits are:
UCL  (2.115) R  (2.115)(1.52)  3.21
LCL  0
The resulting chart is:
2. The
The
X
X -Chart
-Chart simply refers to a control chart for sample means or averages.
Let y i be the sample mean of the n observations taken at the ith time interval;
thus,
1 n
yi   yij
n j 1
Often, both μ0 and σ0 are unknown.
In such situations, we use a base period of m samples to estimate it by
1 m
̂ 0  y   yi
m i 1
Traditionally, people have used
R  (1 / m) Ri to estimate σ.
2
Actually a better estimate of σ is based on s ,
1 m 2
s   si
m i 1
2
2
where the si is the sample standard deviation for the ith sample.
See the text for details.
It can be shown that
UCL  y  A2 R
which is an estimate of
0 
3
n
and
LCL  y  A2 R
which is an estimate of
3
0 
n
where A2 is a constant which comes from the Table 5 in the Appendix.
Note: A2 depends upon the number of observations in each sample.
The center line for this chart is y .
The control chart thus will have the form:
y  A2 R
y
y  A2 R
________________________________
________________________________
________________________________
As an example, again consider the pencil lead ash data set.
The observations are ash contents for lots of pencil lead.
Five lots constitute a sub-group or sample.
Again, we shall use the first 20 samples as our base period.
1 20
1
y   yi  (838.36)  41.918
20 i 1
20
1 20
1
R   Ri  (30.4)  1.52
20 i 1
20
Since n=5, A2 = 0.577, and A2 R  0.877 .
Thus, our control limits are
UCL  y  A2 R  41.918  0.877  42.794
LCL  y  A2 R  41.918  0.877  41.041
The resulting chart is:
3. Comments or
X
and R charts:
1. Typically both charts are run simultaneously.
Recall, a basis assumption of the X chart is that σ2 is constant.
Thus, the
X chart assumes that the R chart is in control.
2. In general the R chart is more stable than the
X
chart.
Generally, the mean tends to change more frequently then the variance.
3. In some cases, the R chart will signal a shift in the mean before the
X
chart.
This can occur if the shift in the mean happens during a particular sample.
In such a case, the Range, which is sensitive to “outliers”, may pick up the shift
while X may not.
4. The R-chart assumes that the data come from a normal distribution (quite
restrictive!).
5. The
X -chart assumes the Central Limit Theorem.
We should use graphical techniques such as the stem-and-leaf, boxplot and the
normal probability plot on the observations in the base period to determine how
reasonable these assumptions are.
Once again, consider the pencil lead ash data.
The stem-and-leaf for the base period follows.
N = 100
Median = 42.1
Quartiles = 41.3, 42.5
Decimal point is at the colon
40 : 1223333
40 : 5556778899
41 : 00002233344
41 : 555566777789
42 : 000000111112222223333333334444
42 : 555555555666666777789
43 : 0012234
43 : 77
The boxplot follows.
The normal probability plot follows.
C. The np-chart
In class, simulate a process which produces good/bad data.
A prime example is Deming's Red Bead Experiment, available from several
sources (check the periodical Quality Progress published by the American Society
for Quality).
The basic idea: The experiment uses a container filled with 500 red and 2500
white beads (20% red beads).
Using a special paddle with 50 holes in it, students sample from the population.
Someone records the number of red beads on each paddle drawn.
An inexpensive variant is to use red and white poker chips.
In class, construct the appropriate control chart.
• yi represents the number of red beads on the ith paddle drawn.
• p is the true proportion of red beads which should be drawn.
Note: As a first approximation, p is 0.2, but as Deming pointed out every time
he conducted the experiment, the true value of p is actually something different,
close but different.
The red and white beads are actually different sizes, so there is no reason to
expect that the true proportion is exactly .2.
Deming's paddle actually had a true proportion of something less than 0.2,
based on the many times he conducted the experiment!
• n is the number of beads drawn each time. Usually, n = 50.
The number of red beads per paddle follows a binomial distribution; thus,
• E(yi) = np, and
• var(yi) = np(1-p).
If p is unknown, then we can use a base period of m samples to estimate p by
1
p y
n
where
1 m
y   yi
m i 1
The appropriate control limits are
UCL  np  3 np (1  p )
LCL  np  3 np (1  p )
CL  np
D. The c-Chart
Consider items which we evaluated in terms of the number of defects present.
Let ci represent the “count” for the number of defects found in
the ith sample of n items inspected.
Note: n may be one.
Assume that n is fixed and constant for all samples.
What is an appropriate distribution for ci?
Poisson
Let λ be the expected number of defects per samples.
What is a reasonable procedure for monitoring ci?
A sequence of hypothesis tests of the form
H0: λ = λ0
Ha: λ ≠ λ0.
The Test Statistic:
Zi 
ci  0
0
Critical Region
Reject H0 if |Zi| > 3
This leads to the following control limits
UCL  0  3 0
LCL  0  3 0
Problem: Do we know
0 ?
No, therefore we must estimate it.
Again use a base period of m samples.
1 m
̂0  c   ci
m i 1
Therefore,
UCL  c  3 c
LCL  c  3 c
CL  c
The following data represent the number of contaminating particles on silicon
wafers after a rinsing operation.
Base Period
7 4 9 9 2 10 3 6 6 5
5 7 5 7 3 4 8 4 5 5
m
 c  114
i 1
Next 40 samples
9 8 8 8 13 10 6 10 11 3
11 13 9 11 13 15 6 10 11 12
12 2 4 7 2 4 7 6 7 4
6 4 6 6 8 5 6 9 3 6
i
For our data, m = 20 and
 c  114; thus,
m
i 1
i
1 m
114
c   ci 
 5.7
m i 1
20
The resulting control limits are
UCL  c  3 c  5.7  3 5.7  12.86
LCL  c  3 c  5.7  3 5.7  1.46  0
The control chart:
E. CUSUM and EWMA Charts
Standard control charts use only the information in the current subgroup.
This leads to them having problems detecting small shifts. Another
approach is to use information from successive groups.
Let’s start with the CUSUM (cumulative sum) chart. Consider the
hypotheses
H 0 :   0
H a :   0
Let Si be the CUSUM statistic for the ith observation where
Si  max 0, Si 1  ( zi  d )
Si-1 is the value of the CUSUM statistic from previous observation,
Zi=(yi-μ0)/σ, d is an appropriate constant.
If Si>0, the current information suggests the process might be out-of-control. How
large is large? We need a threshold, so if Si>h, we will claim the process has shifted.
We typically choose h=5 along with S0=0 and d=0.5.
Consider testing: H 0 :   10
then we say the process is out-of control if Si>5.
H a :   10
We estimate the standard deviation from our data as
MR
ˆ 
 1.461 , then
d2

 yi  10 
S i  max 0, S i 1  

 1.461 

The data and the CUSUM values are shown below:
y 9.45 7.99 9.29 11.66 12.16 10.18 8.04 11.46 9.20 10.34 9.03 11.47 10.51 9.40 10.08
Si
0
0
0
.636 1.614 1.238 0
.499
0
0
0
.506 .355
0
0
We see from the CUSUM chart that the process is in control.
CUSUM Chart
5
h=5
4
Si
3
2
1
0
1
2
3
4
5
6
7
8
Day
9
10
11
12
13
14
15
The CUSUM chart is typically used for one-sided tests, but can be
used for two-sided alternatives by plotting the values we discussed
and comparing them to h=5 and plotting
Si  min0, Si 1  ( zi  d )
using h=-5 as the threshold.
An alternative time-weighted chart is the EWMA (exponentially
weighted moving average) chart. This chart is typically used for twosided alternatives. The EWMA statistic is
Z i  X i  (1   ) Z i 1
where Xi is ith observation, Zi-1 is value of the EWMA statistic for the
previous observation and 0 < Φ ≤ 1 is weighting constant. We will let
Z0=μ0.
The control limits for the EWMA are
UCL    k
0
LCL    k
0

(2   )

(2   )
[1  (1   ) ]
2i
[1  (1   ) ]
2i
with the center line at μ0. Using k=2.7 and Φ=0.1 is approximately equal to the
CUSUM with h=5 and d=0.5.
Consider the previous example:
H 0 :   10
H a :   10
We see from the EWMA chart that the process is in control.
EWMA Chart
11.0
10.886
Zi
10.5
10.0
9.5
9.114
9.0
1
2
3
4
5
6
7
8
Day
9
10
11
12
13
14
15