2 - Statistical Data Analysis and Quality Management
2 - Statistical Data Analysis and Quality Management
Quality Management
IE519 _Summer semester
2023_GJU
Ground Rules
• Learn! Apply!!
• Participate - ask questions.
• Take notes!
• Help each other.
• No PC multi-tasking - email, web surfing, etc. during class.
• Mobile phone on silent/vibrate mode.
• Be on time.
• What else?
Homework and Attendance
Homework Attendance:
• 100% completion required for • Please notify instructor prior to class if
credit. you cannot attend.
• Submit hard copy to instructor. • Unexcused absences are not acceptable.
• This icon will appear • Please sign-in at each class.
whenever a homework is • If you miss a class it is your responsibility
assigned: to catch-up by reading the notes, doing
the in-class exercises/examples, and
homework.
• Show exercises and examples to
instructor.
Outline
• Module 1: Introduction
Assign HW#1
• Secondary Factors:
5. Sample Size – Constant or Varying
6. Matching – Matchable or Matchable with Delta
7. Entities – Single or Multiple.
8. Sidedness – One Sided or Two Sided limits
1. Data Type
• Continuous
• Values of continuous variables fall at any point along a continuum. Theoretically may take on
any value in a range.
• Continuous variables always have numeric values.
• Examples: Temperature, Viscosity, Thickness, Contact Resistance, Offset, Coplanarity
• Binomial
• Values of discrete variable only have two possible values.
• Pass/Fail, Good/Bad
• Examples: Yield, Loss Codes, Test Bins, Recovery, Retest
2. Monitoring Level
• Group Level
• Monitoring lot or group-level summary statistics.
• OOC decisions are not made until the entire sample is collected and available.
• Unit Level
• Monitoring individual “unit” (package, die, bond, etc.) level data points within
a product.
• OOC decisions are made after the data for each “unit” is collected.
• “Real Time” monitoring
3. Summary Statistics
• Process Mean
• Monitoring to detect changes in the location of the process mean.
• Process Variance
• Monitoring to detect changes in the process variation.
Note: Binomial data required different charts discussed late in the semester
4. Detection Objective
• Equipment/Entity/Tool Level
• Monitoring each link, tool, entity, or sub-entity independently to detect
instability.
• In contrast to the differential below.
• Differential
• Monitoring the relative performance of links, tools, entities, or sub-entities to
only detect differences in performance.
Entity/Tool Level Definition
• Need to carefully define the “entity/tool level” for all control charts.
• Most processes or links have multiple tools, entities, or sub-entities.
• For example:
Primary Types of Control Charts
• For each combination of these primary factors, there is a
recommended control chart.
• Varying
• The within lot sample size varies for each lot and/or entity.
6. Matching
• Matchable
• Entities produce “identical” output.
• Multiple Entities
• Monitoring multiple entities.
• Typically, after development, all processes have multiple entities.
• The BKM requirement is CE! control limits across a fleet of tools or entities, or sub-
entities within a tool.
• Note: In this course, the single entity/tool case is covered primarily for illustration and the focus is on setting control limits
across multiple entities.
8. Sidedness
• Two Sided Control Limits
• For process means, if both LCL and UCL limits exist, use two-sided control
limits.
• Examples: Placement offset, thickness, weight
one sided
Voids 0 no mean
upper
High Speed one sided
infinite no mean
Shear lower
Coplanarity two sided unknown no mean
Control Limits
• Statistical control limits are set at 3 standard deviations from the
centerline.
𝐿𝐶𝐿 = 𝐶𝐿 − 3𝜎
𝑈𝐶𝐿 = 𝐶𝐿 + 3𝜎
Too Tight: Response would result Just Right! Minimize Tampering Too Loose: No Response would
in Tampering and Ignoring result in Ignoring
False Alarm Rates
The normal distribution can be used to estimate false alarm rates under
stable “common cause” conditions.
• 2 sigma limits will falsely trigger 5% (50 out of 1000) of the time.
• 3 sigma limits will falsely trigger 0.3% (3 out of 1000) of the time.
Parameter degrades over Tool wear causes the performance to degrade. Engineering
time Parts are then replenished.
Tool to Tool differences that Metrology tools not able to be matched TTD (Tool to Tool
cannot be eliminated across sites due to incoming facilities Delta)
differences.
Setting Control Limits - How
• Setting control limits for a new parameter:
• Set on first tool Initial Limits:
• Start with legacy limits if applicable.
• Otherwise, leave empty.
• Collect data, evaluate excursions, and define responses.
• Refined limits:
• Minimum of 16 runs.
• 2-3 weeks to capture noise (recommended).
• Use methods in the remainder of this lesson.
Computing Control Limits
• Three methods are available in JMP for control limits:
• Moving Range (MR) Method: “Local” Sigma Method
• Levey Jennings (LJ) Method: “Total” Sigma Method
• When evaluating different methods, ask:
• “How well do the limits separate common cause from special cause variation?”
• The most appropriate limits are the ones that best contain common cause variation and best
identify special cause variation.
𝑀𝑅 6.66
𝜎= = = 5.9
1.128 1.128
Note: 1.128 is a constant.
𝑀𝑅
𝐿𝐶𝐿 = 𝑋ത − 3 = 99.68 − 3 ∗ 5.9 = 81.97
1.128
𝑀𝑅
𝑈𝐶𝐿 = 𝑋ത + 3 = 99.68 + 3 ∗ 5.9 = 117.39
1.128
• Note: The computation of MR requires that the data is in process time order.
Example: Setting Limits using Minitab
Thickness measurements made following a plating step.
• Response: Thickness
• Sampling Plan:
– Single measurement per lot (n=1)
– 30 lots
• But for an unstable process, the control limits can be very different.
MR vs. LJ Methods
Moving Range Method Levey Jennings Method
• Since outliers will inflate estimates of the mean and sigma, robust
estimates are used:
• Replace mean with median.
• Replace sigma with pseudo sigma based on percentile and median.
Outlier Screening (cont)
• Calculate Upper Screen and Lower Screen values:
• Exclude runs from control limit computations if above Upper Screen or below
Lower Screen.
• Excluded runs are still included in computing OOC rates.
Example: Effect of Outliers
Thickness data with one outlier.
MR limits using all data:
𝑀𝑅
𝜎=
1.128
𝑀𝑅50
𝜎=
0.954
Note: 0.954 is a constant.
Median is less influenced by outliers
Example: Mean vs Median MR
𝑀𝑅 8.94
𝜎𝑀𝑅 = = =7.92
1.128 1.128
𝑀𝑅50 6.68
𝜎𝑀𝑅50 = = =7.00
0.954 0.954
Exercise : Dispense Weight
• The dataset DispenseWeight contains dispense weight measurements from 45
lots.
• The process target is 15 mg.
• Note: Be sure to specify the target in SOS.
• Note: Use the default of k=4 for outlier screening.
• Example:
• Thickness measurements made to monitor a plating
process.
• Five thickness measurements per lot.
• n =5
𝑋ത Chart: Moving Range
• The same moving range formulas for an
individual chart are applied to the column of Lot
Means:
𝐶𝐿 = 𝑋ത = 99.07
𝑀𝑅𝐿𝑜𝑡 𝑀𝑒𝑎𝑛
𝐿𝐶𝐿 = 𝑋ത − 3 = 99.07 − 3 ∗ 6.78 = 78.73
1.128
𝑀𝑅𝐿𝑜𝑡 𝑀𝑒𝑎𝑛
ത
𝑈𝐶𝐿 = 𝑋 + 3 = 99.07 + 3 ∗ 6.78 = 119.41
1.128
𝑀𝑅𝐿𝑜𝑡 𝑀𝑒𝑎𝑛𝑠
• Where = 7.65/1.128 = 6.78 an estimate of
1.128
s.
s Chart: Moving Range Lot Std Dev(Thickness) Moving Range
𝐶𝐿 = 𝑋ത𝐿𝑜𝑡𝑆𝑡𝑑𝐷𝑒𝑣 = 1.89 10
11
2.75
1.21
0.8
1.55
12 1.67 0.47
13 1.64 0.04
14 1.94 0.3
𝑀𝑅𝐿𝑜𝑡𝑆𝑡𝑑𝐷𝑒𝑣 15 1.46 0.48
1.128 17
18
2.4
0.58
0.7
1.82
𝑀𝑅𝐿𝑜𝑡𝑆𝑡𝑑𝐷𝑒𝑣 19 2.93 2.35
1.128 21
22
3.37
0.85
2.21
2.52
23 1.81 0.96
24 2.49 0.67
25 2 0.49
𝑀𝑅
• Where 𝐿𝑜𝑡𝑆𝑡𝑑𝐷𝑒𝑣 = 1.02/1.128 = 0.9 an estimate of s
26 3.76 1.76
27 2.19 1.56
1.128 28 1.65 0.54
for the control chart. 29 1.44 0.21
30 2.11 0.67
Mean 1.89 1.02
Std Dev 0.77
s Chart Lower Control Limits
• In practice, s charts typically use only an UCL since the objective is to
detect an increase in variation.
• If the computed LCL is < 0, then clearly it is not logical to implement a
LCL.
• In some cases, if the computed LCL > 0, LCL may be implemented to
detect:
• Metrology or data integrity problems, or
• Actual process improvement.
Control Limit BKM Summary
• defaults to start.
• ±3 sigma using mean moving range with ±4 sigma outlier screening.
• Evaluate the separation of special and common cause variation.
• Fine tune accordingly.
• Rounding for convenience.
• Consider other sigma estimates.
• Levey Jennings, Median moving range, different k multiplier for outlier screening
• Manual adjustment.
Example: 𝑋ത and s chart
Thickness measurements made following a plating step.
• Sampling Plan:
• Response: Thickness – 5 measurements per lot (n=5)
– 30 lots
In JMP:
• Open the data set: Thickness003
• Use Minitab- Stat- Control Chart
• Fill in the dialog box.
• Click OK
• Does the UCL appear to distinguish between common and special cause?
• Note: Computed LCL < 0 is displayed but not used.
ഥ and s Charts
Summary: 𝑿
Three graphs are required for setting statistical control limits for lot level data:
• Variability chart of raw data
• 𝑋ത Chart
• s Chart
Exercise : Trace Width
• The dataset TraceWidth. contains trace width data from 36 lots.
• The sample size is 4 units per lot.
• The target for the process is 100 um.
Thickness
105
100 Avg=99.68
95
90
85
80 LCL=81.98
12
15
18
21
24
27
30
3
9
Lot
Online and Offline Stability
• Let’s review a simple case to illustrate the potential difference in the
online and offline assessment of stability.
Is This Process Stable? (Online)
• How many points are out of
control relative to the online
limits?
• 0
• What is the %OOC rate?
• 0%
Recall:
• Online Control Limits: Used on the factory floor.
• Offline Control Limits: Calculated by JMP or Catalyst.
• Appears random.
• Results from common cause variation.
• Centered around the center line.
• Variation is contained within control limits (0% OOC).
Trend Rules
• The different types of instabilities:
• Outliers
• Drifts
• Shifts
• Chaos (Any combination of the above)
Trend rules are used in offline analysis to help identify these patterns of
instability.
Types of Instability: Terminology
• There are 4 typical terms used to describe patterns in an unstable process.
1. Outliers
2. Shift (upward)
Types of Instability: Terminology
3. Drift (upward)
4. Chaos
WECO Trend Rules: Zones
WECO Trend Rules
• The following trend rules are used for offline analysis.
Rule Uses Example
#1. A single point is Detects very large, UCL
outside either sudden shifts, excursions,
control limit. or outliers. CL
LCL
centerline.
LCL
WECO Trend Rules
Rule Uses Example
#5. 2 out of 3 Detects large shifts.
UCL
consecutive points
+2s
are more than 2
CL
sigma from the CL,
-2s
on the same side. LCL
LCL
WECO Trend Rules (Tests)
Test 1 One point beyond Zone A.
Test 2 Nine points in a row in a single (upper or lower) side of Zone C or beyond.
Two out of three points in a row in Zone A or beyond and the point itself is in
Test 5
Zone A or beyond.
Four out of five points in a row in Zone B or beyond and the point itself is in
Test 6
Zone B or beyond.
Test 7 Fifteen points in a row in Zone C, above and below the center line.
Test 8 Eight points in a row on both sides of the center line with none in Zones C.
Trend Rule BKM
• All online control charts use rule #1.
• Offline control charts use rule #1 and WECO (Western Electric) rules #2, 5 and 6
to highlight instability.
• Offline %OOC computation is based only on rule #1.
• Note: WECO rules aren’t applied when engineering limits are used.
• On an exception basis, online charts may implement selected WECO rules to
more quickly detect shifts and drifts.
• Experience has found that WECO rules are rarely required online if the control limits are statistically
valid.
• EWMA charts (discussed in another lesson) are often a better approach for quickly detecting shifts and
drifts.
Checkpoint : Trends
• Assess the stability in each offline control chart. Trend rules #1, 2, 5, and 6 have been
applied. Describe the nature of the instability where applicable.
B.
Checkpoint 2.1: Trends (cont)
• Assess the stability in each offline control chart. Trend rules #1, 2, 5, and 6 have been
applied. Describe the nature of the instability where applicable.
Control Chart Interpretation
• Unstable.
C.
– Single o
• Unstable.
D.
– Chaos.
Example: Laser Scribe Placement
• Laser Scribe Process
• A laser is used to blast away (ablate) material in the streets of a wafer prior to saw.
• Key Parameter: Laser scribe placement.
• Measured in the x direction and y direction (horizontal and vertical streets); on both sides of
the scribe (x1 and x2).
• Note: Leave the SOS window open for assessing control limit validity in Lesson 3.
Example: Analysis
What is the online %OOC rate?
Conclusions
• 0% for both charts.