statistical quality control simple applications of statistics in tqm
TRANSCRIPT
Statistical Quality Control
Simple applications of Statistics in TQM
2
Introduction
• The production processes are not perfect!
• Which means that the output of these processes will not be perfect – correct and deterministic.
• Successive runs of the same production process will produce non-identical parts.
• Alternately, seemingly similar runs of the production process will vary, by some degree, and impart the variation into the some product characteristics.
• Because of these variations in the products, we need probabilistic models and robust statistical techniques to analyze quality of such products.
3
Introduction
• No matter how carefully a production process is controlled, these quality measurements will vary from item to item, and there will be a probability distribution associated with the population of such measurements.
• If all important sources of variations are under control in a production process, then the slight variations among the quality measurements usually cause no serious problems.
• Such a process should produce the same distribution of quality measurements no matter when it is sampled, thus this is a “stable system.”
4
Introduction
• Objective of quality control is to develop a scheme for sampling a process, making a quality measurement of interest on sample items, and then making a decision as to whether or not the process is in the stable state, or “in control.”
• If the sample data suggests that the process is “out of control,” a cause is for the abnormality is sought.
• A common method for making these decisions involves the use of control charts.
• These are very important and widely used techniques in industry, and everyone in the industry, even if not directly related to quality control, should be aware of these.
5
SPC
Statistical process control
• Methodology for monitoring a process to identify special causes of variation and signaling the need to take corrective action.
• When special causes are present, the system said to be statistically out of control.
• If the variations are due to common causes alone, the process is said to be in statistical control.
• SPC relies heavily on control charts.
6
Quality control measurements
• Attributes – A performance characteristics that is either present or absent in the product or service under consideration.
• Examples: Order is either complete or incomplete; an invoice can have one, two, or more errors.
• Attributes data are discrete and tell whether the characteristics conforms to specifications.
• Attributes measurements typically represented as proportions or rates. e.g. rate of errors per opportunity.
• Typically measured by “Go-No Go” gauges.
7
Quality control measurements
• Variable – Continuous data that is concerned with degree of conformance to specifications.
• Generally expressed with statistical measures such as averages and standard deviations.
• Sophisticated instruments (caliper) used.
• In statistical sense, attributes inspection less efficient than variable inspection.
• Attribute data requires larger sample than variable inspection to obtain same amount of statistical information.
• Most quality characteristics in service industry are attributes.
8
The chart
• A typical quality control plan requires sampling one or more items from a production process periodically, and making the appropriate quality measurements.
• Usually more than one items are measured each time to increase accuracy and measure variability.
• The chart helps the quality control person decide whether the center (or average, or the location of central tendency) of the measurement has shifted.
x
x
9
The chart
• Suppose that n observations are to be made at each time in which the process is checked.
• Let Xi be the ith observation (i = 1,…,n) at the specified time point.
• Let be the average of n observations at time j. If
for a process in control, then should be approximately normally distributed with
jX
2 and ii XVarXE
./ and 2 nXVarXE jj
jX
x
10
The chart
• Now since has approximately a normal distribution, we can find the interval that will have a probability of 1- α of containing .
• This interval is
• If the mean and variation of the population (the µ and σ) is known then,
we can set the lower control limit at
and the upper control limit at
jX
n
z 2/
,2/
n
z
jX
,2/
n
z
x
11
The chart
• If a control chart is being constructed for a new process, the µ and σ (the population parameters) will not be known and hence must be estimated from the sample data.
• For establishing control limits, it is generally recommended that at least 30 time points be sampled before the control limits are calculated.
• For each of the k samples, we compute the sample mean, the sample variance, and the range .
• Note that the range of a sample is the difference between the largest and smallest value in the sample.
),,( jjj RSX
x
12
The chart
deviation. std. population ofestimator unbiasednot is thisHowever,
.`1
`
- quantities ` of average calculate wefashion,similar In the
.1
:calculatecan we
means sample from ofestimator unbiasedan form To
. than measurebetter a is 1
`
: thatfoundbeen hasIt
1
2
2
1
222
k
jj
j
k
jj
j
jjj
Sk
S
S
Xk
X
X
SSn
nS
x
13
The chart
tables.ref. from obtained isfor which value,3
where
`, aslimit control thecan write weNow,
.3 have wecesignifican of level 99.74% For,
` :becomes oflimit control theSo,
.for estimator unbiasedan is `
Thus,
tables.ref. from constant aby dividingby unbiased made becan It
2
1
1
2/
2
2/2/
2
2
ncA
SAX
z
nc
SzX
nz
c
S
c
x
14
The chart
• Computation of the mean standard deviation can be avoided by using the range data instead of calculating the adjusted standard deviation.
tables.ref. std. from obtained is Here, is! ratio theHowever,
estimator. unbiasedan not isit , with iprelationsh a has Though
.1
22
1
dd
R
σR
Rk
Rk
jj
x
15
The chart
ce.significan of level 99.74% indicating ,3
where
,or
:becomeslimit control new theThus,
2
2
2
2
2/
ndA
RAXnd
RzX
x
16
The r-Chart
• Deciding whether the center of the distribution of quality measurements has shifted up or down may not be enough.
• It is frequently of interest to decide if the variability of the process measurements has significantly increased or decreased.
• A process that suddenly starts turning out highly variable products could cause severe problems in the operations.
• Since we already have established the fact that there is a relationship between the ranges (Rj) and σ, it would be natural to base our control chart for variability on the Rj’s.
• A control chart for variability could be based on the adjusted std dev also (Sj), but use of ranges provides nearly as much accuracy for much less computation.
17
The r-Chart
• Using the same argument as control chart for central tendency (x-bar), one could argue that almost all the Rj’s should be within three standard deviations of the mean of Rj. Now
32/232/2
32
223
2
:is then,interval, confidence The
tables.ref. stdin found are and where
;
;
dzddzd
dd
dRVar
dRE
j
j
18
The r-Chart
• As in the previous chart, if σ is not specified, it must be estimated from the data. The best estimator of σ based on the range is , and the estimator of the control then becomes:
2/ dR
.3 ce,significan of level 99.74%at again, Once
1or ;
2/
2
32/32/2
2
z
d
dzRdzd
d
R
19
The r-Chart
tables.lstatisticain found are and of Values
.,
becomeslimit control the,31 and ;31 Letting,
43
43
42
33
2
3
DD
DRDR
Dd
dD
d
d
20
Control charts for variable data
• Two basic types: x-bar-chart and R-chart – one used to monitor centering of the process, the other for variation.
• Construction and establishing statistical control:
RAxLCLRDLCL
RAxUCLRDUCL
k
xx
k
RR
xR
xR
k
ii
k
ii
23
24
11
;
;
;
21
Interpreting patterns in control charts
• General rules to determine whether a process is in control:1. No points outside the control limits.2. The number of points above and below the center line are
about the same.3. Points seem to fall randomly above and below the center
line.4. Most points are near center line, and only a few are close to
control limits.
• Basic assumption: Central Limit Theorem.
22
Interpreting patterns in control charts
• One point outside control limits: measurement or calculation error, power surge, a broken tool,
incomplete operation.• Sudden shift in process average: new operator or inspector, new machine setting.• Cycles: operator rotation or fatigue at the end of shift, different gauges,
seasonal effects such as temperature and humidity.• Trends: x-bar-chart – learning effect, dirt or chip buildup, tool wear,
aging of equipment; R-chart (increasing trend) – gradual decline in material quality; R-chart (decreasing trend) – improved skills, better materials.
23
Interpreting patterns in control charts
• Hugging the center line: sample taken over various machines canceling out the variation
within the sample.• Hugging the control limits: sample taken over various machines not canceling out the
variation within the sample.• Instability: difficult to identify causes. Typically, over-adjustment of
machine.
• Always, R-chart analysis before the x-bar-chart analysis.
24
Interpreting patterns in control charts
Downward trend in R-chart…
25
Interpreting patterns in control charts
….causes smaller variation in x-bar-chart.
26
Process capability
• Particularly process characteristics study: Performed over a period of time under actual operating conditions to capture the variations in material and operators.
• Here we differentiate between the control limits and specified tolerance limits (also called specification limits).
• Control limits are obtained from the data itself. Probability of finding similar data within the calculated control limits can be found.
• Specification limits are specified by product/process designers.
• Process capability study checks the probability of getting product dimension within the specified limits.
27
Process capability
• Note that we are interested in finding individual product dimensions within limits (and not the sample means).
• We need to use the estimate for population variance and not variance of sample mean.
• From the previous analyses, we have:
n
d
R
c
S
:mean sample of dev. Std.
or `
:std.dev. Population22
28
Process capability
• Let the specified tolerance limits be (a, b).• To determine that the process of capable of producing parts
within these limits, we calculate the probability of finding parts outside these limits. Hence, calculate –
• We calculate the probability area to the right of z1 and left of z2 from the standard normal table.
• Addition of these two probabilities gives us the probability that a part will not meet the specifications.
ax
zxb
z
21 and
29
Process capability
%1.31ion specificat meet thenot doesproduct that prob. Total
0.1590.341-0.5 is ofleft the toProb.
0.1520.348-0.5 is ofright the toProb.
0.137.0
0.1637.16
027.137.0
37.1675.16
.37.08407.0
31.0.8407.0,5For
31.0` ;37.16
.75.160,.16let example, numerical theFrom
2
1
2
1
2
z
z
z
z
cn
Sx
ba
30
Process capability
%5.30ion specificat meet thenot doesproduct that prob. Total
0.1560.3438-0.5 is ofleft the toProb.
0.1490.3508-0.5 is ofright the toProb.
013.1365.0
0.1637.16
041.1365.0
37.1675.16
.365.0326.2
85.0.326.2,5For
85.0 ;37.16
.75.160,.16let example, numerical theFrom
2
1
2
1
2
z
z
z
z
dn
Rx
ba
31
Process capability
• PCI – Process Capability Index
.337.026.1548.17
0.1675.16
y,Alternatel
26.1537.0*337.16 ;48.1737.0*337.16
.337.037.0*6
0.1675.16
6
:example numerical For the
.6
p
p
p
C
LSUS
abC
LTLUTLC
32
The p-Chart
• The control charts we looked at so far were applicable to quality measurements that possessed continuous probability distribution.
• This sampling is commonly referred to as “sampling by variables.”
• In many cases we merely want to assess whether or not a certain item is defective.
• We can observe a number of defective items from a particular sample of series of samples. This is referred as “sampling by attributes.”
33
The p-Chart
• Suppose a series of k independent samples, each of size n, is selected for a particular process.
• Let p denote the proportion of defective items in the population (total production for a certain time period) for a process in control.
• Let Xi denote the number of defective items in the ith sample.
• Then Xi is a binomial variable, assuming random sampling from a large lot, with mean E(Xi) = np, and Var(Xi) = np(1-p).
• We typically consider the fraction of defective items in a sample (Xi/n) rather than the observed number of defectives.
34
The p-Chart
.size sample total
observed defective ofnumber total
: ofestimator unbiased The
samples. from data theusing
limits, control theestimatemust weunknown, is Since
./)1(/ and / since,
)1(3
interval in theor mean, their of
deviations standard ree within thlie should fractions sample All
1
nk
Xp
p
k
p
nppnXVarpnXEn
ppp
k
ii
ii
35
The p-Chart
• Thus the estimated control limits are given by:
.)1(
3n
ppp
36
The c-Chart
• In many quality control problems the particular items being subjected to inspection may have more than one defect.
• We may wish to count # of defects instead of merely classifying at item as to whether or not it is defective.
• If Ci denotes the # of defects observed in the ith inspected item, we can safely assume that Ci has a Poisson distribution.
• Let this Poisson distribution have a mean of λ for a process in control.
• For us to be 99.74% sure, almost all the Ci’s should fall within three standard deviations of the mean if the process is in control.
37
The c-Chart
.3
:become then limits control estimated The
.1
:is ofestimator unbiased then theinspected, are items If
.3
:are limits control the, and Since
1
CC
Ck
C
k
CVarCE
k
ii
ii
38
Acceptance sampling by attributes
• Inspections in which a unit of product is classified simply as defective or non-defective is called inspection by attributes.
• Lots of items, either raw materials or finished products, are sold by a producer to consumer with some guarantee for quality.
• We determine quality, here, by the proportion p of defective items in the lot.
• To check this characteristics of the lot, the consumer will sample some of the items, test them, and observe the number of defectives.
39
Acceptance sampling by attributes
• Generally, some defectives are allowed because defective-free lots may be too expensive for the consumer to purchase.
• But if the number of defectives is too large, the consumer will reject the lot and return it to the producer.
• Sampling is generally used because the cost of inspecting the entire lot may be too high or the inspection process may be destructive.
• Before sampling inspection takes place for a lot, the consumer must have in mind a proportion p0 of defectives that will be acceptable.
• Thus, if the true proportion of defectives p is no greater than p0, the consumer wants to accept the lot, and reject otherwise.
40
Acceptance sampling by attributes
• This maximum proportion of defectives satisfactory to the consumer is called the Acceptable Quality Level (AQL).
• So, we have a hypothesis-testing problem. • We are testing null hypothesis of the true defectives
proportions being less than or equal to the acceptable proportion; against an alternative hypothesis that the true proportion is actually greater than the acceptable one.
.:against : 00 ppHppH ao
41
Acceptance sampling by attributes
• A random sample of n items is selected from the lot, and the number of defectives Y is observed.
• The decision concerning rejection or non-rejection of Ho (and consequently rejecting or accepting the lot) will be based on the observed value of Y.
• The probability of Type I error (α) in this problem is the probability that the lot is rejected by the consumer when, in fact, the proportion of defectives is satisfactory to the consumer.
• This is referred to as the producer’s risk.
42
Acceptance sampling by attributes
• The value of α calculated at p = p0 is the upper limit to the proportion of good lots rejected by the sampling plan being considered.
• The probability of Type II error β calculated for some proportion of defective p1, where p1 > p0, represents the probability that an unsatisfactory lot will be accepted by the consumer.
• This is referred to as the consumer’s risk.
• The value of β calculated at p1 = p is the upper limit to the proportion of bad lots accepted by the sampling plan being considered, for all .1pp
43
Acceptance sampling by attributes
• The null hypothesis will be rejected if the observed value of Y is larger than some constant a.
• Since the lot is accepted if Y is less than or equal to a, the constant a is called the acceptance number.
• In this context, the significance level α is given by:
)(lot when theAccepting
when rejectingNot 1
,
when Reject
0
00
0
00
APppP
ppHP
ppaYP
ppHP
44
Acceptance sampling by attributes
• Since p0 may not be known precisely, it is of interest to see how P(A) behaves as a function of true p, for a given n and a.
• A plot of these probabilities is called an operating characteristic curve (OC curve).
45
Acceptance sampling by attributes
Example of OC curve plotting:• Suppose that one sampling plan calls for (n = 10, a = 1) while
another calls for (n = 25, a = 3). Plot the OC curve for both of them.
• If the plant using these plans can operate with 30% defective raw materials, considering the price, but cannot operate efficiently if the proportion gets close to 40%, which plan do you recommend?
Solution: We need to calculate the P(A) for various values of p. • For each plan, P(A) is the binomial probability of finding at
the most a number of defects in n trials.
46
OC curve plotting
• For the first sampling plan, we have:
• For the second plan:
p 0 0.1 0.2 0.3 0.4 0.6
P(A) 1 0.736 0.376 0.15 0.046 0.002
p 0 0.1 0.2 0.3 0.4
P(A) 1 0.764 0.234 0.033 0.002
47
OC curve plotting
48
OC curve plotting
• Note that OC curve for first plan (n = 10, a = 1) drops slowly and does not get appreciably small until p is in the neighborhood of 0.4.
• For p around 0.3 this plan still has fairly high probability of accepting the lot.
• The other curve for the second plan (n = 25, a = 3) drops much more rapidly and falls to a small P(A) at p = 0.3.
• Hence, the second plan would be better for inspection if the plant can afford to sample n = 25 items out of each lot before making a decision.
49
Deming’s kp rule
• Objective is to minimize the average total cost of inspection of incoming materials and final products for processes that are stable.
• It calls for 0% or 100% inspection.• It can be shown that, if the process is stable, the distribution of
the nonconforming items in the sample is independent of the distribution of nonconforming items in the remainder of the lot.
• Items are submitted in lots, and a random sample is drawn from that lot to make a decision regarding the entire lot.
50
Deming’s kp rule
Assumptions
• Inspection process is completely reliable.
• All items are inspected prior to moving forward to the next customer in the process.
• The vendor provides the buyer with an extra supply of items in order to replace any nonconforming items that are found. The cost is assumed to included in the contract.
51
Deming’s kp rule
Notations• p: average proportion of nonconforming items in the lots• k1: cost of initial inspection of an item• k2: cost of repair or reassembly due to the usage of a
nonconforming item• k: average cost to find a conforming item from the additional
supply to replace a detected nonconforming item = k1/(1-p)• xi =1 if item i is nonconforming; 0 otherwise.
52
Deming’s kp rule
• The unit cost of initial inspection and replacement, if an item is nonconforming, is given by:
• Unit cost to repair and replace a nonconforming item:
• These costs components are mutually exclusive (if one is present, other is zero.)
inspectednot item if 0
inspected item if 11 i
ikxkC i
inspected item if 0
inspectednot item if )( 22 i
ixkkC i
53
Deming’s kp rule
• If the item is inspected, the total cost is (k1+kxi). If the item is not inspected, the total cost is (k2+k)xi.
• It can be extended to the whole lot, where the proportion of nonconforming items is p.
• So, the average total cost per item if items are inspected is (k1+kp). And the average total cost per item if items are not inspected is (k2+k)p.
• The break-even point for p is obtained by equating the average total costs.
• We get: p = k1/k2.
54
Deming’s kp rule
• If k1/k2 is greater than p, conduct no inspection. This situation occurs if the proportion of nonconforming items is low, cost of inspection is high, and cost of repairing a nonconformity is low.
• If k1/k2 is less than p, conduct 100% inspection. In this situation, it is quite expensive for a nonconforming item to be allowed into production.
• If k1/k2 is equal to p, either no or 100% inspection may be conducted. Usually, if the estimated value of p is not very reliable, 100% inspection is conducted.
55
Critique of the kp rule
• Estimating proportion p requires sampling. Contradiction: for sampling, initially some inspection would have to be carried out.
• Proportion p may not be stationary over a long period of time. Once again, sampling may be needed to verify.
• Estimating the costs is not easy. Linearity of costs may not hold.
• What 100% inspection policy, if it involves destructive testing?
56
Acceptance sampling by variables
• At times the characteristics under study involves a measurement on each sampled item.
• In these cases, one could base the decision to accept or reject on the measurements themselves, rather than simply on the number of defective items.
• Defective items would mean those items having measurements that do not meet the standards.
• When looking at actual measurements obtained from variables method, the experimenter may gain some insights into degree of nonconformance and may be able to quickly suggest methods of improvements.
57
Acceptance sampling by variables
• For variables, arriving at the correct sample size and correct acceptance factor is slightly tedious and depends on the non-standard statistical tables.
• US Military standards (MIL-STD-414) is one such standard for calculating the same.