ISO/TC69 Meeting in Beijing and Has the rate of
Download
Report
Transcript ISO/TC69 Meeting in Beijing and Has the rate of
ASQ 1401 Section
El Paso, TEXAS
2009 January 14
Rudy Kittlitz
ISO/TC 69
Application of Statistical Methods
This Technical Committee [TC] has several
subcommittees [SC] and working groups [WG]
Terminology and symbols
Statistical interpretation of data
Applications of statistical methods in process
management
Acceptance sampling
Measurement methods and results
Six Sigma applications
th
30
Plenary Meeting
Beijing 2008 October 11-17
Approximately 75 delegates
From India, Germany, France, Denmark, United
Kingdom, etc.
US Technical Advisory Group [TAG]
One-day meetings in the spring and the fall
Rudy became involved in ISO/TC 69 with the 1995 March
meeting of the US TAG
His first international meeting was in 1996 June in
Stockholm and has attended each year
In 2001 Rudy was elected as Chair of the US TAG and has
continued
ANSI [American National Standards Institute] is the
official contact with ISO
All USA Delegates represent the USA and ANSI, not their
company, university, etc.
Meetings in Beijing
Arrived Friday night, 10/10
First meeting on Saturday
A special committee meeting on Sunday
Monday through Friday meetings
Took a tour on Tuesday and on Thursday
Left on Saturday, 10/18
Using Statistics To Answer The Question:
“Has The Rate of Category 3+ Atlantic Hurricanes Changed
Since 1940?”
This is not a talk on Global Warming
But a statistical assessment on whether or not the rate
(i.e., number of hurricanes per year) of category 3+
Atlantic hurricanes has changed since 1940
Cat 3+ are also known as “Major Hurricanes”
Over the past few years all sorts of statements about
“cycles of Atlantic hurricanes”, “increased intensity”, and
others
A simple statistical assessment should be able to answer
Figure 1. Atlantic Category 3+ Hurricanes By Year
Ref. http://weather.unisys.com/hurricane/atlantic
9
8
Category 3+
7
6
5
4
3
2
1
0
1940
1950
1960
1970
1980
Year
1990
2000
2010
Initial Comments About The Data
Some possible “up and down” for these 68 data points,
but has the rate changed?
Proper application of statistical analysis and Statistical
Process Control (SPC) should be able to answer this
question
Mean of the data is 2.588
Standard deviation is 1.863
Figure 2. Frequency Of Category 3+ Atlantic Hurricanes
Ref. Figure 1
20
Frequency
15
10
5
0
0
2
4
Category 3+
6
8
The Poisson Distribution
The Poisson is a candidate to describe this data
It is a count distribution
Only need to know the mean or the average
For small averages, the positive skew is evident
This is seen in Figure 2
The theoretical standard deviation is
Avg
Observed std dev = 1.893 vs theoretical std dev = 1.609
F-test indicates no significant difference [p = 0.234]
A SPC chart of Poisson data is the c-chart
Usual Calculations For The c-Chart
Upper Control Limit [UCL]
UCL = Avg + 3
UCL = 2.588 + 3
Avg
2.588
= 7.414
Lower Control Limit [LCL]
LCL = Avg – 3
LCL = 2.588 – 3
Avg
2.588
= - 2.238
No LCL since calculated LCL is below zero
These simple calculations ignore the skewness of the
Poisson
Improved Limits For Poisson Data
Article in Quality Engineering
Kittlitz, R. G. Jr. (2006). Calculating the (Almost) Exact Control
Limits for a C-Chart. Quality Engineering, 18:359-366.
The improved limits are based on a simple transformation of
the original data
2
It is c 0.25 3
Kittlitz, R. G. Jr. (2003). Transforming The Poisson Distribution
To Symmetry For SPC Applications and Other Statistical
Analysis. MS Capstone Project for University of Alabama at
Huntsville
For the hurricane data this transformation produces a
skewness of 0.47 vs original of 0.98
Improved Limits For Poisson Data Cont’d
UCL =
1
Avg
12
23
1
Avg
12
23
LCL =
2
16
3 Avg
3
2
16
3 Avg
3
32
32
3
4
1
4
Don’t let these equations scare you!
Programmable calculator performs the calculations
Improved UCL = 8.07
Improved LCL = No Lower Limit
Calculations produces a negative number inside bracket
Figure 3. Atlantic Category 3+ Hurricanes By Year
Ref. Figure 1
9
UB=8.07
8
Individual Value
7
6
5
4
_
X=2.588
3
2
1
0
No LB
1940
1950
1960
1970
1980
Year
1990
2000
2010
Initial Conclusions From Analysis
The typical “run-rules” for an SPC chart did not trigger
any signals
The 1950 data point of 8 “Cat 3+” hurricanes is close to
the limit, but the cumulative probability of 8 for an
average of 2.588 is 0.99856 which is less than the limit
of 0.99865
Unless the rate changes, we can expect 0 to 8 “Cat 3+”
hurricanes per year
Some Additional Analysis
Advances have been made to detect a shift in the mean of data
more efficiently/quicker and is an improvement over the
typical “counting rules”
Exponentially Weighted Moving Average [EWMA]
Cumulative Summation or CUSUM
An EWMA chart will be used to analyze the transformed data
The 2003 reference details how the transformed mean and
the transformed standard deviation can be calculated from
the original mean
Transformed mean = 1.92443
Transformed standard deviation = 0.777823
The EWMA chart is shown in Figure 4
Figure 4. EWMA Chart of Transformed Atlantic Category 3+ Hurricanes by Year
Ref. Figure 1
2.8
UCL=2.702
2.6
2.4
EWMA
2.2
_
_
X=1.924
2.0
1.8
1.6
1.4
1.2
LCL=1.147
1.0
1940
1950
1960
1970
1980
Year
1990
2000
Comments About Figure 4
Remember for an EWMA chart that “counting rules”
cannot be used since the points are not independent
Likewise the “wandering” of the points are typical of an
EWMA chart and no conclusions can be drawn from any
apparent “cycles”
The only valid signal is if a point exceeds the limits
No points exceed either the upper or lower limits
Thus, no significant change in the mean
Questions?