Attribute Statistical Process Control Chart
Attribute Statistical Process Control Chart
-1-
Soo King Lim
-2-
Soo King Lim
p0 (1 p0 )
LCL or UCL = p0 Z / 2 (1)
n
p0 (1 p0 )
For α = 0.0026, the control limits are p0 3 . The central line is p0 .
n
If p0 value is not specified, it must be estimated from the data. Let’s assume
that there are m sample batches, each with sample size n are collected and that the
total number of defective part in these m sample batches is d then the estimate of the
proportion of defectives per sample batch of sample size n is p = d/(mn), p̂0 and p .
Then the control limits are
p (1 p )
LCL p 3 (2)
n
CL p (3)
p (1 p )
UCL p 3 (4)
n
-3-
Soo King Lim
Example 1
The readings in Fig. 1 are the number of defective parts in 18 sample batches, each
containing a total of 50 items i.e. n = 50. Establish the control limits.
Solution
The total number of defective parts in the 18 sample batches collected is 229. Each
sample batch has 50 items. Thus, the average fraction of defective in these 18 sample
batches is p = 229/(18x50) = 0.254. The limits using average fraction of defective
p (1 p ) 0.254(1 0.254)
as the estimate of p0 are p 3 ; 0.254 3 , which are 0.0693
n 50
and 0.439 and the center line is 0.254.
In this example, all the p values are within the calculated control limits. Hence,
these limits can be used for monitoring p. If one or more p value falls outside the
limits, then these values have to be removed and a new p is calculated. This
procedure has to be repeated until all the p values used in the estimation of p are
within the control limits.
-4-
Soo King Lim
The test statistic used is the number of defective part x in a sample batch of size
n, which is denoted as np. The control limits of this chart are obtained by multiplying
the lower and upper control limits of the p chart by n, which are
p (1 p0 )
LCL or UCL n p0 3 0 for α = 0.0026 (5)
n
In equation (5), p0 can be replaced by its estimate, which is p . The center line
CL is CL = np0. This is called np chart.
Example 2
Using data shown in Fig. 1, establish the control limits for number of defective part
in the lot and plot the control chart.
Solution
The control limits are calculated using p = 0.254 and n = 50, and equation
p (1 p0 ) 0.254(1 0.254)
LCL or UCL n p0 3 0 , which is 500.254 3 after
n 50
substituting p0 value and n value.
It yields lower control limit LCL equal to 3.45 and upper control limit UCL equal
to 21.95.
The statistical control chart is shown in Fig. 2. It may sound ridiculous to have
fractional control limits since the data is non parametric type. Nevertheless, it is
acceptable for academic purpose.
-5-
Soo King Lim
The test hypothesis is; null hypothesis H0: mean number of defect per piece or
some units or per batch = c0 versus alternative hypothesis H1: mean number of defect
per piece or some units or product c0.
The test statistic is the number of defect per the appropriate unit, which is c.
The expected value and the standard deviation with assumption that c obeys a
Poisson distribution are E(c) = c0; standard deviation of c = c0. Poisson distribution
are E(c) = c0 and standard deviation of c = c0. Hence, the control limits c0 ± Zα/2 c0
becomes equation (6) taking α is α = 0.0026.
-6-
Soo King Lim
c0 ± 3 c0 (6)
The center line is c0. If c0 cannot be specified, it can be estimated by the mean
number of defect from one or more sample batches, which is denoted by c then, the
control limits become c 3 c .The center line is located at c .
Example 3
Data shown in Fig. 3 are the number of defect in 10 sample batches, where each
batch contains 15 items. Establish the control limits of number of defect in batch
and plot the control chart.
Number of
Batch number i
Defect in Batch c
1 18
2 12
3 14
4 9
5 16
6 10
7 14
8 15
9 19
10 17
Figure 3: Data of example 3
Solution
Assuming that the process is in control when these observations are collected, the
calculated limits of the c chart are c = (18 + 12 + 14 + 9 + 16 + 10 + 14 + 15 + 19
+ 17)/10 = 12.4.
The control limits are LCL = 12.4 - 3 12.4 = 1.84 and UCL = 12.4+ 3 12.4 = 22.96.
The center line CL is 12.4. All the c values used in the estimation are within the
control limits as is shown in Fig. 4.
-7-
Soo King Lim
The test statistic used is the average number of defect per unit U. The expected
value E(U) and standard deviation of U are
E(U) = U0 (7)
U = U0 / n (8)
where U = c/n; E(c) = c0; Var(c) = c0; E(U) = c0/n = U0; Var (U) = Var(c)/n2 = c0/n2
= nU0/n2 = U0/n. n is the sample size. Hence, the lower and upper control limits of
the U chart are
-8-
Soo King Lim
Example 4
The data set shown in Fig. 5 contains the number of defect in eight sample batches,
where each batch contains 15 items. Assuming that the process is in control when
these observations are collected, calculate the control limits of the U statistical
process control chart and plot the control chart.
c i
U i 1
= 0.86
8x15
Figure 5: Data of example 4
Solution
The average number of defect U is U = [18 + 12 + 7 + 9 + 16 + 12 + 14 + 15]/(8x15)
= 0.86.
The control limits are U 3 U / n = 0.86 3 0.86 / 15 , which are 0.14 and 1.58. The
center line CL is 0.86. The statistical control chart is shown in Fig. 6.
-9-
Soo King Lim
Example 5
The data of the number of defect found in eight sample batches with unequal sample
batch size are given in Fig. 7. Calculate the control limits of the U chart and plot the
control chart.
Solution
In this example, the sample sizes are not equal. Hence each sample batch will have
its own control limits. This is because the control limits U 3 U / n are functions of
the sample size. Grand mean U is calculated using all observations. Thus, U is U
= [5 + 6 + 2 + 7 + 9 + 15 + 7 +14]/[2 + 3 + 2 + 5 + 8 + + 3 + 7] = 65/40 = 1.63.
The control limits are equal to U 3 U / ni , where the subscript i denotes batch i.
The center line CL is 1.63.
For i = 1, the control limits are 1.63 3 1.63 / 2 , which are -1.08 and 4.34 respectively.
Since negative value is not permitted, the control limits for i = 1 are lower control
- 10 -
Soo King Lim
limit LCL = 0, center line CL = 1.63, and upper control limit UCL = 4.34. The rest
of control limits for i value batches are shown in Fig. 8.
Average
Batch Sample Size
Number of LCLi CL UCLi
Number i ni
Defect Ui
1 2 2.5 0.00 1.63 4.34
2 3 2.0 0.00 1.63 3.84
3 2 1.0 0.00 1.63 4.34
4 5 1.4 0.00 1.63 3.34
8 8 1.1 0.28 1.63 2.98
6 10 1.5 0.42 1.63 2.84
7 3 2.3 0.00 1.63 3.84
8 7 2.0 0.18 1.63 3.08
Figure 8: Control limits of U chart for i = 1, 2, 3, 4, 5, 6, 7, and 8
The control limits of eight sample batches are given in Fig. 8. It can be seen that the
average numbers of defect Ui’s for all sample batches are within their respective
control limits. The U control chart is shown in Fig. 9.
- 11 -
Soo King Lim
p(1 p )
p3 . The probability q for m out of n samples from the lot is defective is
n
given by
q p
m
(11)
ARL = 1 / p
m
(12)
Example 6
Let’s take the earlier illustration of Shewart control chart that sample size is n = 6,
the probability of the plotted point outside the control limit is p = 1 -
Pr 5.44 Z 0.55 = 1 - (0.7088 - 0) = 0.2912 and the ARL is 3.43. Calculate the
ARL of p-control chart using this data for 2 out of 6 and 6 out of 6 samples are
defective. Comment your results.
Solution
The average run length of the p-control chart for 2 out of 6 samples are defective is
ARL = 1 / p = 1 / 0.29122 = 11.79.
m
The average run length of the p-control chart for 6 out of 6 samples are defective is
ARL = 1 / p = 1 / 0.29126 = 1,640.
m
For 2 out of 6 samples are defective case, it would most probable occur once
between 11 and 12 sample batches.
For 6 out of 6 samples are defective case, it would most probable occur once every
1,640 sample batches.
- 12 -
Soo King Lim
A P
Average run length ................................................. 10 p chart ....................................................................... 3
C U
c control chart ........................................................... 5 U control chart .......................................................... 7
N
np chart ..................................................................... 4
- 13 -