Statistical Quality Control
Statistical Quality Control
Basic Concept
W. Edward Deming advocated the implementation of a statistical quality management
approach.
His philosophy behind this approach is ‘reduce variation’- fundamental to
the principle of continuous improvement and
the achievement of consistency, reliability, and uniformity.
It helps in trustworthiness, competitive position, and success.
Statistical Quality Control
Statistics: data sufficient enough to obtain a reliable result.
Quality: relative term and can be defined as totality of features and characteristics of a
product or service that bear on its ability to satisfy stated or implied need
(ISO).
Control: The operational techniques and activities (a system for measuring and checking)
used to fulfill the requirements for quality. It
incorporates a feedback mechanism system to explore the causes of poor quality
or unsatisfactory performance and
takes corrective actions.
also suggests when to inspect, how often to inspect, and how much to inspect.
Basic Concept ….Cont’d
Statistical Quality Control
A quality control system using statistical techniques to control quality by
performing
inspection,
testing and
analysis
to conclude whether the product is as stated or designed quality standard.
Relying on the probability theory, SQC
evaluates batch quality and
controls the quality of processes or products
It makes the inspection more reliable and less costly.
The basis of the measurement is the performance indicator, either individual,
group or departmental calculated over time (hourly, daily, or weekly).
These performance measures are plotted on a chart.
Pattern obtained from plotting these measures are basis of taking
appropriate actions so that
The process variation in minimized and
Major problems are prevented in future.
The timing and type of, and responsibility for, these actions depends on
whether the causes of variation is controlled or uncontrolled
Basic Concept ….Cont’d
Statistical Quality Control
In repetitive manufacture of a product, even with refined machinery,
skilled operator, and selected material, variations are inevitable in the
quality of units produced due to interactions of various causes.
Variation may be due to
Common or random causes of variation (as a result of normal
variation in material, method, and so on that causes natural variation
in product or process quality) resulting in stable pattern of variation.
Special causes (changes in men, machine, materials or tools, jigs and
fixture and so on) resulting in a shift from the stable pattern of
variation.
SQC assists in timely identification and elimination of the problem with
an object of reducing variations in process or product.
x
the accuracy or centering of a process.
i
The Mean- simply the average of a set of data
i 1
Sum of all the measurements/data divided by the
number of observations. x
The Median- simply the value of middle item if the
data are arranged in ascending or descending order.
n
Applies directly if the number in the series is odd. K
It lies between two middle numbers if the number
of the series is even.
The Mode- value that repeat itself maximum number
X
j 1
j
of times in the series.
Shape of Distribution of Observed Data K
A measure of distribution of data
Normal or bell shaped
Skewed
Distribution of Data
Also a measure of quality
characteristics.
Symmetric distribution - same
number of data are observed above
and below the mean.
This is what we see only when
normal variation is present in the
data
Skewed distribution – a
disproportionate number of data are
observed either above or below the
mean.
Mean and median fall at different
points in the distribution
Centre of gravity is shifted to
oneside or other.
Traditional Statistical Tools …cont’d
Measure of Precision or Spread
Reveals the extent to which numerical data
tend to spread about the mean value.
The Range- the simplest possible measure
of dispersion.
Difference between largest and smallest
observations in a set of data.
o Depends on sample size and it tends
to increase as sample size increases.
o Remains the same despite changes
in values lying between two extreme
values.
x
n 2
of the values from the mean.
i X
i 1
Small values >> data are closely σ
clustered around the mean n 1
Large values >> data are spread out
around the mean.
Statistical Process Control
Process Control
Refers to procedures or techniques adopted to evaluate, maintain and
improve the quality standard in various stages of manufacture.
A process is considered satisfactory as long as it produces items within
designed specification.
Process should be continuously monitored to ensure that the
process behaves as it is expected.
Salient features of process control
Controling the process at the right level and variability.
Detecting the deviation as quickly as possible so as to take
immediate corrective actions.
Ultimate aim is not only to detect trouble, but also to find out the
cause.
Developing an efficient information system in order to establish an
efficient system of process control.
Statistical Process Control
Statistical Process Control (SPC)
Statistical evaluation of the output of a process during production.
Goal is to make the process stable over time and then keep it stable unless the
planned changes are made.
Statistical description of stability requires that ‘pattern of variation’ remains
stable over time, not that there be no variation in the variable measured.
In statistical process control language:
A process that is in control has only common or random cause variation -
an inherent variability of the system.
When the normal functioning of the prosess is disturbed by some
unpredictable events, special cause variation is added to common cause
variation.
characteristics
Measured
One or more data points above an UCL or
below a LCL mark statistically significant
changes in the process
A process is in control if
Time period
No sample points outside limits
Most points near process average
About equal number of points above and below centerline
Points appear randomly distributed
A process is assumed to be out of control if
Rule 1: A single point plots outside the control limits;
Rule 2: Two out of three consecutive points fall outside the two sigma warning limits on
the same side of the center line;
Rule 3: Four out of five consecutive points fall beyond the 1 sigma limit on the same
side of the center line;
Rule 4: Nine or more consecutive points fall to one side of the center line;
Rule 5: There is a run of six or more consecutive points steadily increasing or
decreasing
Control Chart …Cont’d
Setting Control Limits
Type I error
Concluding a process is not in control when it actually is.
Type II error
Concluding a process is in control when it is not.
/2 /2
Mean
Type I error
In control No Error
(producers risk)
Type II Error
Out of control No error
(consumers risk)
Control Chart …Cont’d
UCL = μ + kσ
CL = μ
LCL = μ – kσ
where
μ is the mean of the variable
σ is the standard deviation of the variable
UCL=upper control limit; LCL = lower control limit;
CL = center line.
k is the distance of the control limits from the center line,
expressed in terms of standard deviation units.
When k is set to 3, we speak of 3-sigma control charts.
Historically, k = 3 has become an accepted standard in
industry.
Control Chart …Cont’d
Sample Size
Attribute charts require larger sample sizes
50 to 100 parts in a sample
Variable charts require smaller sample sizes
2 to 10 parts in a sample
Control Chart …Cont’d
CL =μ i X i
LWL = μ – 2/3 kσx
X
m
LCL or LAL = μ – kσx
X
n
Control Chart …Cont’d
Mean control charts Interpret the R-chart first:
Used to detect the If R-chart is in control -> interpret the
variations in mean of X-bar chart -> (i) if in control:
a process.
X-bar chart the process is in control; (ii) if out
Range control charts
of control: the process average
Used to detect the changes in is out of control
dispersion or variability If R-chart is out of control: the
of a process process variation is out of
R chart
control -> investigate the
System can show acceptable central cause; no need to interpret the X-bar
chart
tendencies but unacceptable variability
or
System can show acceptable variability
but unacceptable central tendencies
Control Chart
Plot the sample mean in the
sequence from which it was
generated and interpret the
pattern in the control chart.
Control Chart …Cont’d
Second Method for X-bar Chart using Range and A2 factor
Use this method when standard deviation for the process
distribution is unknown.
Control limits solution:
k
R i
R i 1
;
k
R R
or ;& x
d2 dn n
Center line and 3-sigma
control Fomulas:
CLx X
R
UCLx X 3 X A2 R
d2 n
R
LCLx X 3 X A2 R
d2 n
Control Chart …Cont’d
Calculate limits for X-bar chart using sample range
(x i x j) 2
S
j 1
j
Sj i 1
&S
n 1 k
Control Chart …Cont’d
UCL =
3 sigma = x + A2R
Zone A
= 2
2 sigma = x + 3 (A2R)
Zone B
= 1
1 sigma = x + 3 (A2R)
Zone C
Process =
x
average
Zone C
=
1 sigma = x - 1 (A2R)
3
Zone B
=
2 sigma = x - 2 (A2R)
3
Zone A
=
LCL 3 sigma = x - A2R
| | | | | | | | | | | | |
1 2 3 4 5 6 7 8 9 10 11 12 13
Sample number
Interpreting the Patterns …cont’d
Performing a Pattern Test
5.10 –
5.08 –
UCL = 5.08
5.06 –
5.04 –
5.02 –
x= = 5.01
Mean 5.00 –
4.98 –
x- bar
Chart 4.96 –
Example 4.94 –
LCL = 4.94
(cont.) 4.92 –
| | | | | | | | | |
1 2 3 4 5 6 7 8 9 10
Sample number
Interpreting the Patterns …cont’d
Performing a Pattern Test
1 4.98 B — B
2 5.00 B U C
3 4.95 B D A
4 4.96 B D A
5 4.99 B U C
6 5.01 — U C
7 5.02 A U C
8 5.05 A U B
9 5.08 A U A
10 5.03 A D B
Control Chart …Cont’d
p(1 p )
UCL p 3
ni
p (1 p )
LCL p 3
ni
p (1 p )
n i
UCL p 3 n Where
i 1
n m
p (1 p )
LCL p 3
n
Control Chart : c chart
No standard given
Average number of nonconformities per sample unit is found from the sample
observation and is denoted by c-bar.
The center line and control limits are:
CLc c
UCLc c 3 c
LCLc c 3 c
If lower control limit is found to be less than zero, it is converted to zero.
Standard given
if the specified target for the number of nonconformities per sample unit be co..
The center line and control limits are then calculated from:
CLc co
UCLc co 3 co
LCLo co 3 co
Control Chart: c chart …Cont’d
SAMPLE DEFECTS 21
190 18
c= = 12.67
Number of defects
c = 12.67
1 12 15
15
2 8
12
UCL = c + zc 9
3 16
= 12.67 + 3 12.67
6
: :
3 LCL = 1.99
= 23.35 2 4 6 8 10 12 14 16
Sample number
: : LCL = c + zc
15 15 = 12.67 - 3 12.67
190 = 1.99
Process Capability Analysis
Process Variability Lower Upper
Specification Specification
Natural variation in the process.
Tolerance or specification
Range of acceptable values established
design engineers or product design specialist. 6σ
For the product to be considered acceptable, its quality characteristics must fall
within this preset range.
Process Capability
Process variability relative to specification.
Relationship between the process variability and the tolerance can be
formulized by the consideration of standard deviation, σ of the process.
In order to manufacture product within the specification,
the distance between the upper specification limit (USL) and the lower
specification limit (LSL) i.e. (USL-LSL) or 2T must be equal to or greater than
width of the process variability defined by the control limits, i.e. 6σ
Process Capability Analysis …Cont’d
Relationship between the process width and specification
The relationship between (USL - LSL) and 6σ results in three levels of
situations:
Design
Specifications
Process
Process variation is large relative to the specifications.
A large percentage of the product will fall outside the specification.
Process is not capable of meeting specifications all the time.
The process cannot be considered capable regardless of the
process centering.
Process variation must be reduced drastically
Process Capability Analysis …Cont’d
Design
Specifications
Process
Design
Specifications
Process
Process variation is small relative to the specifications.
The process mean can shift about without causing the process to
degrade its capability
Process is capable of meeting specifications all the time.
This will reduce the defects per million (DPM), reduce the cost of
quality (COQ), and hence increase profitability.
Using only
the Cp
measure
would lead
to an
incorrect
conclusion
Process Capability Analysis …Cont’d
CpK Index
Another measure of process capability
This measure accounts for the location of the process mean and is used when
the process mean is not at the target value, which is assumed to be the
halfway between the specification limits.
The process capability of each half of the normal distribution is computed and
minimum of the two is used. i.e.
where
USL LSL µ = the mean of the process
Cpk min( , ) σ = the standard deviation of the
3 3 process
Process Capability Analysis …Cont’d
Net weight specification = 9.0 oz 0.5 oz Net weight specification = 9.0 oz 0.5 oz
Process mean = 8.80 oz Process mean = 8.80 oz
Process standard deviation = 0.12 oz Process standard deviation = 0.12 oz
=
x - lower specification limit
upper specification limit - ,
3
lower specification limit Cpk = minimum =
Cp = upper specification limit - x
6 3
The Cp value leads us to conclude that the process is capable. Whereas, the C pk
value is less than 1, revealing that process is not capable.
Reasons for the difference in measures is that process is not centered on
specification range.
Process Capability Analysis …Cont’d
Moonlight Jeans
Moonlight Jeans store receives a shipment of 300 pairs of jeans from its
warehouse. It is common practice for the store to sample 5% of the total
received. The acceptance number under any and all circumstances for
Moonlight Jeans is 10. Of the 15 pairs of jeans observed, 2 were defective.
What conclusion should the store manager come to based on this
information?
The store manager has just found out that the clerk who inspected the
samples made a huge mistake… The actual number of defective
pairs of jeans sampled was 12.
What type of risk is involved with the error made at
Moonlight Jeans?
How might this error affect the store and their customers?
Acceptance Sampling …Cont’d
Random Sampling
Important
Units selected for inspection from lot must be chosen at random
Should be representative of all units in a lot.
Watch for Salting
Vendor may put “good” units on top layer of lot knowing a lax
inspector might only sample from the top layer.
Suggested Technique
Assign a number to each unit, or use location of unit in lot.
Generate / pick a random number for each unit / location in lot.
Sort on the random number – reordering the lot / location pairs.
Select first (or last) n items to make sample.
Acceptance Sampling …Cont’d
Acceptance Sampling Plan
Specifies the lot size, sample size, number of samples and
acceptance/rejection criteria for lot sentencing.
Unless it is mentioned following convention is practiced
Sampling is performed without replacement
Sampling is a simple random sample
Each item in the lot has equal probability of being in the sample.
Designing Sampling Plan for Attributes
Simpliest sample plan
Based on binomial distribution (if sample is less than 20 units otherwise use
poisson’s distribution)
Requires large sample size
Sampling plan involves
Single sampling plan
Double sampling plan,
Multiple sampling plan
Acceptance Sampling …Cont’d
Single Sampling Plan
Quality characteristic is an attribute, i.e., conforming or nonconforming
Define:
N: lot size
n: sample size, and
c: Acceptance number
Procedure:
Take a sample of size n and inspect each of the items drawn
If d ≤ c, accept lot; else reject,
d: number of defective items in sample
Acceptance Sampling …Cont’d
Double Sampling Plan
Define:
n1: sample size of the first sample
c1: acceptance number for the first sample
n2: sample size of the second sample
c2: acceptance number for the second sample
Procedure:
Take an initial sample of size n1
If number of defective items, d1 ≤ c1, accept the lot
If number of defective items, d1 > c2, reject the lot
If c1 < d1 ≤ c2, take second sample of size, n2
If the combined number of defective items (d1+d2) ≤ c2, accept the
lot; otherwise reject the lot
Acceptance Sampling …Cont’d
Characterizing attribute sampling plan
Typically four graphs are used to characterize a sampling plan.
Operating Characteristic (OC) curve
The probability of acceptance for a given quality level.
1 20 e 1.2
prob( x 0 | p 0.02) 0.3012
0!
1 21e 1.2
prob( x 1 | p 0.02) 0.3614
1!
1 22 e 1.2
prob( x 2 | p 0.02) 0.2169
2!
1 23 e 1.2
prob( x 3 | p 0.02) 0.3012
3!
prob( x 3 | p 0.02) 0.3012 0.3614 0.2169 0.0867
0.9662 = probability of accepting the lot.
Acceptance Sampling …Cont’d
How to Construct OC Curve
OC curve is plot of Pa vs p
Pa = P (Accepting Lot | true proportion of defective is p)
p = lot fraction defective
Suppose, n = 89 and c = 2.
For each value of p, compute the probability of acceptance using a Binomial or
Poisson approximation as given in the table:
p = fraction defective in Pa = P [Accepting Probability of Acceptance, Pa
lot Lot]
1.0
0.005 0.9897
0.8 n=89
0.010 0.9397 Pa
0.6 c=2
0.4
0.015 0.8502
0.2
0.020 0.7366 0.0
0.025 0.6153 0.00 0.02 0.04 0.06 0.08 0.10
Lot fraction defective, p
0.030 0.4985
0.035 0.3936
Acceptance Sampling …Cont’d
Ideal OC Curve
This characterizes the ideal sampling plan
that discriminates perfectly between good
Keep whole
shipment
and bad shipments
If the supplier’s process average 100 –
Whole Shipment)
nonconforming is below the AQL, the
consumer will accept all the shipped 75 – Return whole
lots. shipment
If the supplier’s process average 50 –
nonconforming is above the AQL, the Cut-Off
consumer will reject all the shipped lots. 25 –
P(Accept
Both α and β are zero
0 |– | | | | | | | | | |
it is obtainable by 100% inspection IF 0 10 20 30 40 50 60 70 80 90 100
inspection are error free. % Defective in Lot
ideal OC curve is unobtainable in
practice
Acceptance Sampling …Cont’d
Effect of n on OC Curve
Precision with which a sampling plan differentiates between good and bad lots
increases as the sample size increases
Increasing n (with c proportional) approaches the ideal OC curve.
100.0%
Probability of acceptance, Pa
n= 50, c=1
80.0%
n=100, c=2
60.0%
n=200, c=4
40.0%
20.0%
0.0%
0.0% 2.0% 4.0% 6.0% 8.0% 10.0% 12.0% 14.0%
Percent nonconform ing, p
Acceptance Sampling …Cont’d
Effect of c on OC Curve
Changing acceptance number, c, does not dramatically change slope of OC
curve.
Increasing c (with n constant) approaches the ideal OC curve.
Operating Characteristic Curve
100.0%
Probability of acceptance, Pa
20.0%
0.0%
0.0% 2.0% 4.0% 6.0% 8.0% 10.0% 12.0% 14.0%
Percent nonconform ing, p
Acceptance Sampling …Cont’d
Balancing the producer and consumer risks
The true value of p is unknown
Even if p known, sampling involves randomness, and we can still
reject a lot even if p < c/n
accept a lot even if p > c/n
The OC curve gives an indication of the values of α and β.
Producer’s risk- α
Producer wants as many lots accepted by consumer as possible so
Producer “makes sure” the process produces a level of fraction defective
equal to or less than:
p1 = AQL = Acceptable Quality Level
α is the probability that a good lot will be rejected by the consumer even
though the lot really has a fraction defective ≤ p1
Lot is rejected given that process
P
has an acceptable quality level
P
Lot is rejected p AQL
Acceptance Sampling …Cont’d
Consumer’s risk- β
Consumer wants to make sure that no bad lots are accepted
Consumer says, “I will not accept a lot if percent defective is greater than or
equal to p2”
p2 = LTPD = Lot Tolerance Percent Defective
β is the probability a bad lot is accepted by the consumer when the lot really
has a fraction defective ≥ p2
Lot accepted given that lot
P
has unacceptable quality level
P
Lot accepted p LTPD
Example: Suppose
p1 = 0.01, α = 0.05, p2 = 0.06, and β = 0.10.
Find the acceptance sampling plan.
Acceptance Sampling …Cont’d
Acceptance Sampling …Cont’d
Designing a Sampling Plan using Tables
For a given producer’s and consumer’s risks, various tables have been
developed for constructing single and double sampling plans.
Excerpt From a Sampling Plan Table with Producers Risk = 0.05 and Consumers Risk = 0.10
C LTPD/AQL n(AQL) n(LTPD)
0 44.89 0.052 2.334
1 10.946 .355 3.886
2 6.509 .818 5.324
3 4.89 1.366 6.68
4 4.057 1.97 7.992
5 3.549 2.613 9.274
6 3.206 3.286 10.535
7 2.957 3.981 11.772
8 2.768 4.695 12.996
9 2.618 5.426 14.205
For double sampling plan, use Grubbs’ Tables (Table 10-6 and
Table 10-7: Amitava Mitra, Fundamentals of Quality Control and
Improvement, 2nd ed., Pages: 445-446)
Acceptance Sampling …Cont’d
Rectifying Sampling
All known defective units replaced with good ones, that is,
If lot is rejected, replace all bad units in lot
If lot is accepted, just replace the bad units in sample
Such sampling program is known as rectifying inspection program
Since such inspection activity affects the final quality of the outgoing product,
two questions come to mind :
How many items are inspected on average after rectifying inspection?
What is the average outgoing quality after rectifying inspection?
Acceptance Sampling …Cont’d
Average Outgoing Quality: AOQ
Quality that results from application of rectifying inspection.
The long-run ratio of expected number of defectives to expected number of
items successfully passing through the inspection plan
Average value obtained over long sequence of lots from process with fraction
defective p:
Pa p N n
AOQ
N
N - Lot size, n = # units in sample that, after inspection, contains no
defectives.
If the lot is accepted: number of defective units in the lot
# units
fraction
p N n remaining
defective in lot
Acceptance Sampling …Cont’d
Average Outgoing Quality: AOQ
Expected number of defective
units:
Lot # defective
Pa p N n Prob
accepted units in lot
Averag e fractio n n o n co n fo rm in g , o u tg o in g lo ts
7.0%
AOQL
possible value of AOQ over all 6.0%
(AOQL). 3.0%
revised. 0.0%
0.0% 20.0% 40.0% 60.0% 80.0% 100.0%
Pe rcent nonconform ing, p
Acceptance Sampling …Cont’d
Average Total Inspected: ATI
Total number of inspection per lot
required for sampling plan.
If the lot is fully conforming,
p=0.0 (Pa=1.0), then no lot will Average Total Inspection Curve
be rejected and hence, inspect
160.0
only the sample.
If the lot is totally 140.0
0.0
p = fraction defectives; 0.0% 20.0% 40.0% 60.0% 80.0% 100.0%
Percent nonconform ing, p
Pa = Probability of
accepting a lot.
Design of Experiment
Why Experiment?
To increase knowledge about a process
Design
Speed Octane Tire Pressure Matrix
55 85 30
60 85 30
55 90 30
60 90 30
55 85 35
60 85 35
55 90 35 A treatment
60 90 35
combination
Design of Experiment ….Cont’d
What makes up an experiment?
Response Variable(s)
Factors
Randomization
Repetition and Replication
Response Variable
The variable that is measured and the object of the characterization or
optimization (the Y).
Defining the response variable can be difficult.
Often selected due to ease of measurement.
Factors
A variable which is controlled or varied in a systematic way during the
experiment (the X)
Tested at 2 or more levels to observe its effect on the response variable(s)
(Ys)
Design of Experiment ….Cont’d
Randomization
Randomization can be done in several ways :
Run the treatment combinations in random order
Assign experimental units to treatment combinations randomly
An experimental unit is the entity to which a specific treatment
combination is applied
Advantage of randomization is to “average out” the effects of extraneous
factors (called noise)
that may be present
but are not controlled or measured during the experiment
Effect sizes determine which factors have the most significant impact on the
results.
ANOVA determine the significance of each factor based on these effect
calculations.
Design of Experiment ….Cont’d
Main Effects Plots:
Effects plots are a quick and efficient way to
visualize effect size.
The grand mean is plotted as a horizontal
line.
The average result is represented by dots
for each factor level.
The Y axis is always the same for each factor
in Main Effects Plots.
Factors with steeper slopes have larger
effects and thus
larger impacts on the results.
This graph shows that
A+ has a higher mean fill height than A-
B+ and C+ also have higher means than B-
and C- respectively, and
The effect size of factor A, Carbonation, is
larger than the other factor effects.
Design of Experiment ….Cont’d
Interaction Effects:
It is often critical to identify how multiple factors
interact in affecting the results.
An interaction occurs when one factor affects
the results differently depending on a second
factor.
To find the AB interaction effect, first calculate the
average result for each of the four level
combinations of A and B:
Calculate the average when factors A and B
are both at the low level (-4 + -1) / 2 = -2.5
Calculate the mean when factors A and B are
both at the high level (5 + 11) / 2 = 8
Also calculate the average result for each of
the levels of AC and BC.
Design of Experiment ….Cont’d
Interaction Effect plots:
Interaction plots are used to determine the
effect size of interactions.
For example, the AB plot below shows that
the effect of B is larger when A is 12%.
AB plot also shows that when the data is
restricted to A+, the B effect is more
steeper than when we restrict our data to
A-.
The bottom right plot shows the interaction (or
2-way effects) of all three factors:
When the lines are parallel, no interaction
effects between the factors.
The more different the slopes, the more
influence the interaction effect has on the
results.
This graph shows that the AB
interaction effect is the largest.