GINI MEAN DIFFERENCE AND EWMA CHARTS. Muhammad Riaz, Department of Statistics, Quaid-e-Azam University Islamabad,

Similar documents
A Robust Exponentially Weighted Moving Average Control Chart for the Process Mean

Exponential Weighted Moving Average (EWMA) Chart Under The Assumption of Moderateness And Its 3 Control Limits

Exponentially Weighted Moving Average (EWMA) Chart Based on Six Delta Initiatives

Vehicle Arrival Models : Headway

Robust estimation based on the first- and third-moment restrictions of the power transformation model

STRUCTURAL CHANGE IN TIME SERIES OF THE EXCHANGE RATES BETWEEN YEN-DOLLAR AND YEN-EURO IN

COMPUTATION OF THE PERFORMANCE OF SHEWHART CONTROL CHARTS. Pieter Mulder, Julian Morris and Elaine B. Martin

An Improved Adaptive CUSUM Control Chart for Monitoring Process Mean

THE CUSUM VERSUS MCUSUM MODIFIED CONTROL CHARTS WHEN APPLIED ON DIESEL ENGINES PARAMETERS CONTROL

Testing the Random Walk Model. i.i.d. ( ) r

Smoothing. Backward smoother: At any give T, replace the observation yt by a combination of observations at & before T

J. Martin van Zyl Department of Mathematical Statistics and Actuarial Science, University of the Free State, PO Box 339, Bloemfontein, South Africa

Research Report Statistical Research Unit Department of Economics University of Gothenburg

An Adaptive Generalized Likelihood Ratio Control Chart for Detecting an Unknown Mean Pattern

Testing for a Single Factor Model in the Multivariate State Space Framework

Kriging Models Predicting Atrazine Concentrations in Surface Water Draining Agricultural Watersheds

Excel-Based Solution Method For The Optimal Policy Of The Hadley And Whittin s Exact Model With Arma Demand

Diebold, Chapter 7. Francis X. Diebold, Elements of Forecasting, 4th Edition (Mason, Ohio: Cengage Learning, 2006). Chapter 7. Characterizing Cycles

Article from. Predictive Analytics and Futurism. July 2016 Issue 13

Reliability of Technical Systems

Object tracking: Using HMMs to estimate the geographical location of fish

20. Applications of the Genetic-Drift Model

Shiva Akhtarian MSc Student, Department of Computer Engineering and Information Technology, Payame Noor University, Iran

A Nonparametric Multivariate Control Chart Based on. Data Depth. Abstract

DEPARTMENT OF STATISTICS

ACE 562 Fall Lecture 5: The Simple Linear Regression Model: Sampling Properties of the Least Squares Estimators. by Professor Scott H.

Estimation Uncertainty

Pade and Laguerre Approximations Applied. to the Active Queue Management Model. of Internet Protocol

EFFECT OF AUTOCORRELATION ON SPC CHART PERFORMANCE

Air Traffic Forecast Empirical Research Based on the MCMC Method

OBJECTIVES OF TIME SERIES ANALYSIS

Structural Break Detection in Time Series Models

A Generalized Poisson-Akash Distribution: Properties and Applications

Comparing Means: t-tests for One Sample & Two Related Samples

Lecture 3: Exponential Smoothing

Econ Autocorrelation. Sanjaya DeSilva

Computer Simulates the Effect of Internal Restriction on Residuals in Linear Regression Model with First-order Autoregressive Procedures

Structural Break Detection in Time Series Models

Appendix to Creating Work Breaks From Available Idleness

Econ107 Applied Econometrics Topic 7: Multicollinearity (Studenmund, Chapter 8)

Structural Break Detection in Time Series Models

Dynamic Probability Control Limits for Risk-Adjusted Bernoulli Cumulative Sum Charts

Structural Break Detection in Time Series Models

On Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature

A new flexible Weibull distribution

Chapter 11. Heteroskedasticity The Nature of Heteroskedasticity. In Chapter 3 we introduced the linear model (11.1.1)

A note on spurious regressions between stationary series

1. Introduction. Rawid Banchuin

Christos Papadimitriou & Luca Trevisan November 22, 2016

Regression with Time Series Data

Inventory Control of Perishable Items in a Two-Echelon Supply Chain

Georey E. Hinton. University oftoronto. Technical Report CRG-TR February 22, Abstract

Process Monitoring versus Process Adjustment

Failure of the work-hamiltonian connection for free energy calculations. Abstract

A unit root test based on smooth transitions and nonlinear adjustment

A Specification Test for Linear Dynamic Stochastic General Equilibrium Models

Navneet Saini, Mayank Goyal, Vishal Bansal (2013); Term Project AML310; Indian Institute of Technology Delhi

Structural Break Detection in Time Series

Time series Decomposition method

SPC Procedures for Monitoring Autocorrelated Processes

BOOTSTRAP PREDICTION INTERVALS FOR TIME SERIES MODELS WITH HETROSCEDASTIC ERRORS. Department of Statistics, Islamia College, Peshawar, KP, Pakistan 2

STATE-SPACE MODELLING. A mass balance across the tank gives:

5 The fitting methods used in the normalization of DSD

ACE 564 Spring Lecture 7. Extensions of The Multiple Regression Model: Dummy Independent Variables. by Professor Scott H.

13.3 Term structure models

Evaluation of Mean Time to System Failure of a Repairable 3-out-of-4 System with Online Preventive Maintenance

Distribution of Least Squares

Monitoring and data filtering II. Dynamic Linear Models

The General Linear Test in the Ridge Regression

State-Space Models. Initialization, Estimation and Smoothing of the Kalman Filter

ACE 562 Fall Lecture 8: The Simple Linear Regression Model: R 2, Reporting the Results and Prediction. by Professor Scott H.

PENALIZED LEAST SQUARES AND PENALIZED LIKELIHOOD

Time Series Models for Growth of Urban Population in SAARC Countries

Random Walk with Anti-Correlated Steps

Distribution of Estimates

Stochastic Model for Cancer Cell Growth through Single Forward Mutation

I. Return Calculations (20 pts, 4 points each)

Statistical procedures searching for suspicious measurements and their sequences in measurement series. Ing. Marek Brabec, PhD CHMU, SZU, UI Praha

Chapter 4. Location-Scale-Based Parametric Distributions. William Q. Meeker and Luis A. Escobar Iowa State University and Louisiana State University

Article Sequential Change-Point Detection via Online Convex Optimization. and Huan Xu

GMM - Generalized Method of Moments

Department of Economics East Carolina University Greenville, NC Phone: Fax:

Solutions to Odd Number Exercises in Chapter 6

A New Unit Root Test against Asymmetric ESTAR Nonlinearity with Smooth Breaks

Semi-Competing Risks on A Trivariate Weibull Survival Model

A Bayesian Approach to Spectral Analysis

4.1 Other Interpretations of Ridge Regression

Module 2 F c i k c s la l w a s o s f dif di fusi s o i n

Math 10B: Mock Mid II. April 13, 2016

Performance of Stochastically Intermittent Sensors in Detecting a Target Traveling between Two Areas

Stationary Time Series

A Point Optimal Test for the Null of Near Integration. A. Aznar and M. I. Ayuda 1. University of Zaragoza

WATER LEVEL TRACKING WITH CONDENSATION ALGORITHM

Probabilistic Robotics

Improved Approximate Solutions for Nonlinear Evolutions Equations in Mathematical Physics Using the Reduced Differential Transform Method

Nature Neuroscience: doi: /nn Supplementary Figure 1. Spike-count autocorrelations in time.

MATHEMATICAL DESCRIPTION OF THEORETICAL METHODS OF RESERVE ECONOMY OF CONSIGNMENT STORES

Chapter 5. Heterocedastic Models. Introduction to time series (2008) 1

The Arcsine Distribution

The Simple Linear Regression Model: Reporting the Results and Choosing the Functional Form

Transcription:

GINI MEAN DIFFEENCE AND EWMA CHATS Muhammad iaz, Deparmen of Saisics, Quaid-e-Azam Universiy Islamabad, Pakisan. E-Mail: riaz76qau@yahoo.com Saddam Akbar Abbasi, Deparmen of Saisics, Quaid-e-Azam Universiy Islamabad, Pakisan. E-Mail: saddamabbasi@yahoo.com Absrac: For an improved monioring of process parameers, an EWMA conrol charing scheme is proposed. The proposal inroduces Gini mean difference based EWMA conrol scheme for locaion and scale parameers using asympoic limis. The comparison of he proposed schemes are made wih he EWMA conrol schemes based on sample range, sample sandard deviaion S and Downon s esimaor used by Khoo (2004). Key Words: Conrol Chars; Downon s esimaor; EWMA; Gini mean difference; locaion and scale parameers; Normaliy; Ouliers. 1. Inroducion In saisical Process Conrol ool ki, conrol char is he mos commonly used and he mos powerful ool for process monioring. Process parameers (e.g. locaion and scale) need a imely deecion of unfavorable variaions which may cause deerioraion in qualiy of he producs. EWMA and CUSUM conrol srucures address smaller shifs whereas Shewhar conrol schemes address larger shifs in process parameers. In EWMA, CUSUM and Shewhar seup, robus conrol schemes are always desirable. 1

There are many robus conrol chars in qualiy conrol lieraure e.g. ocke (1992) proposed robus EWMA based on sample range and sample mean, Khoo (2004) proposed proposed downon s based EWMA schemes, Maravelakis e al. (2005) examines robusness o non normaliy of he EWMA Conrol Chars for he spread parameer where as Soumbos and Sullivan (2002) examines robusness o non normaliy of he mulivariae EWMA conrol char ec. Ng and Case (1989) proposed wo chars EWMAS and EWMASM for monioring spread and locaion parameers respecively based on sample range.. EWMAS char is proposed by ploing he saisic following asympoic limis. r$ ( 1 ) $ = λ + λ r 1 agains he d3 λ UCL = 1+ 3 d2 2 λ CL = d3 λ LCL = 1 3 d 2 λ 2 (1.1) where λ (lying beween 0 and 1) is he weigh assigned o each observaion, d and d 2 3 are well known conrol char coefficiens available in qualiy conrol lieraure (e.g. see Mongomery (2001)) as a funcion of sample size (n) and is he average of sample s esimaed from iniially seleced samples assumed o come from an in-conrol process. EWMASM char was proposed by ploing he saisic Z λx ( 1 λ) 1 following conrol limis: ˆ = + Z ˆ agains he 2

3 λ UCL = X + d 2 2 n λ CL Zˆ = o = X 3 λ LCL = X d n 2 λ 2 (1.2) where X is he grand mean of average sample X s esimaed from iniially seleced samples assumed o come from an in-conrol process. Khoo (2004) proposed wo chars namely EWMAVD and EWMAMD chars for monioring spread and locaion parameers respecively based on Downon s esimaor. EWMAVD char was proposed by ploing he saisic ˆ λ ( λ) he following conrol limis: Y 1 ˆ = S + Y agains -1 ˆ 2 λ UCL = Yo + 3ˆ σˆ = c4 3 ( 1 c Y + 4) σ 2 λ CL Yˆ c = o = 4σ ˆ 2 λ 3ˆ LCL = Yo + σ ˆ = max 0, c4 3 ( 1 c Y 4) σ 2 λ (1.3) where c 4 is a well known conrol char coefficien available in qualiy conrol lieraure * (e.g. see Mongomery (2001)) as a funcion of sample size (n), σ is he Downon s Esimaor defined as: n * 2 π 1 σ = n n 1 + i= 1 2 ( ) i ( n 1) X() (1.4) i and σ * is he average of sample σ s esimaed from iniially seleced samples assumed o come from an in-conrol process. 3

ˆ = + 1 Z agains he EWMAMD char was proposed by ploing he saisic Z λx ( λ) $ 1 following conrol limis ˆ 3 λ UCL = Zo + 3 ˆ σzˆ = X + σ n 2 λ CL Zˆ = o = X ˆ 3 λ LCL = Zo 3 ˆ σzˆ = X σ n 2 λ (1.5) The similar EWMA conrol srucures based on sample saisic S for monioring process spread and locaion parameers are defined respecively in erms of SdV and SdM chars as: SdV char is defined by ploing he saisic ˆ λ ( λ) following conrol limis: Y 1 ˆ = S + Y agains he -1 UCL = c + 3 λ 2 ( 1 c4 ) S 4 2 λ c4 CL = S 2 λ ( 1 c4 ) S LCL = c4 3 2 λ c4 (1.6) where S is he average of sample S s esimaed from iniially seleced samples assumed o come from an in-conrol process. SdM char is defined by ploing he saisic Z λx ( λ) Z $ 1 conrol limis ˆ = + 1 agains he following 4

3 λ S UCL = X + n 2 λ c4 CL Zˆ = o = X 3 λ S LCL = X n 2 λ c 4 (1.7) The following secion inroduces EWMA conrol srucures for monioring locaion and scale parameers based on Gini mean difference using asympoic limis. 2. The Proposed Chars The proposals are based on Gini mean difference (G), which was used by Gini (1912), as a measure of variabiliy. Le x1, x2, x3,..., x n be an i.i.d random sample, hen G is defined as: G = n n j= 1 i= 1 j i x n( n 1) / 2 i x j (2.1) Gini (1912) had shown superioriy of sample saisic G over he sample saisic as an esimae ofσ, in erms of efficiency and robusness agains non normaliy. David (1968) showed ha π 2 G serves as an unbiased esimae ofσ. iaz and Saghir (2007) proposed a Shewhar ype conrol char namely G char for robus monioring of process scale parameerσ. They defined he wo random variables K and M as: π K = G, (2.2) 2 M K =, (2.3) σ 5

and developed he design srucure of G-char in erms of coefficiens b 3 and he quanile poins of M. These quaniies (i.e. Tables 1 and 2 respecively. b 3 and quanile poins of M) are given in Appendix We propose here he use of K in EWMA ype paern for improved monioring of variabiliy and locaion parameers. The GVA char for monioring process variabiliy is proposed as: ˆ 2 λ UCL = Yo + 3ˆ σ ˆ = c4 3 ( 1 c Y + 4) K 2 λ CL Yˆ c K, = o = 4 ˆ 2 λ 3ˆ LCL = Yo σ ˆ = max 0, c4 3 ( 1 c Y 4) K 2 λ (2.4) where K he average of sample K s esimaed from iniially seleced samples assumed o come from an in-conrol process. Afer consrucing he conrol srucure of GVA char we plo he EWMA saisic ( λ) Yˆ = λs + 1 Yˆ following Khoo (2004). Any -1 Ŷ falling ouside he conrol limis is an indicaion of ou of conrol siuaion wih respec o scale parameer. Similarly an asympoic limi EWMA scheme namely GMA char is proposed for monioring he locaion parameer as: ˆ 3 λ UCL = Zo + 3 ˆ σ ˆ = X + Z K n 2 λ CL Zˆ = o = X, ˆ 3 λ LCL = Zo 3 ˆ σ Zˆ = X K n 2 λ (2.5) 6

Afer consrucing he conrol srucure for GMA char, we plo he EWMA saisic ( λ) Z $ 1 Zˆ = λx + 1 following Khoo (2004). Any Z ˆ falling ouside he conrol limis is an indicaion of ou of conrol siuaion wih respec o locaion parameer. 3. Comparisons In his secion, comparison of he proposed EWMA charing schemes is made wih ha of Khoo (2004) proposed scheme and he well known and S based EWMA schemes for locaion and spread parameers. For comparison purposes, he abiliy o deec ouliers is used as a performance measure. Following Khoo (2004), we have chosen m=10 subgroups of sizes n = 5 and 10 by considering he following wo siuaions: a) Un-conaminaed disribuion where all he samples are generaed from he sandard normal disribuion i.e. N (0, 1). b) Conaminaed disribuions where ouliers are inroduced in he daa in four differen combinaions namely C05N3 (i.e. a siuaion where 95% observaion come from N (0, 1) and 5% from N (0, 9) ), C10N3 (i.e. 90% from N (0, 1) & 10% N (0, 9)), C05N5 (i.e. 95% from N (0, 1) & 5% from N (0, 25)) and C10N5 (i.e. 90% from N (0, 1) & 10% from N (0, 25)). The subgroups (m=10) of sizes n=5 and 10 are generaed for boh he above menioned siuaions (a) and (b), asympoic limis are compued for all he EWMA conrol srucures based on, S, G and * σ for spread and locaion parameers. Afer he compuaion of he conrol limis, he number of subgroups for which he EWMA saisics fall ouside heir respecive conrol limis, are couned. The whole procedure is simulaed 10,000 imes and he proporion of subgroups ha fall ouside he limis is 7

calculaed for all he EWMA conrol srucures. The resuling proporions are repored in Appendix Tables 3 and 4 for spread and locaion parameers respecively. These simulaion resuls are in close agreemen wih hose of Khoo (2004) provided for EWMAVD, EWMAS, EWMAMD and EWMASM, which shows he validiy of our rouines used o obain he simulaion resuls. For ease in comparison among differen conrol schemes under discussion, graphs of proporion of ouliers agains values of λ are ploed. The graphs of columns iled N (0, 1) for n=5 and 10, CO5N3 for n=5 and C10N3 for n=10 (from Appendix Table 3) are provided here in he following Figures 3.1-3.4 for scale parameer. The graphs of columns iled N (0, 1) for n=5, C1ON3 for n=5, C0N5 for n=10 and C10N3 for n=10 (from Appendix Table 4) are provided here in he following Figures 3.5-3.8 for locaion parameer. Fig 3.1: Proporion of ouliers deeced by EWMAS, GVA, SdV and EWMAVD for monioring process scale parameer for N(0,1) and n=5 Proporion of ouliers 0.0035 0.0030 0.0025 0.0020 0.0015 0.0010 0.0005 0.0000 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 G S D Lambda 8

Fig 3.2: Proporion of ouliers deeced by EWMAS, GVA, SdV and EWMAVD for monioring process scale parameer for N (0, 1) and n=10 Proporion of ouliers 0.0035 0.0030 0.0025 0.0020 0.0015 0.0010 0.0005 0.0000 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 G S D Lambda Fig 3.3: Proporion of ouliers deeced by EWMAS, GVA, SdV and EWMAVD for monioring process scale parameer for C05N3 and n=5 Proporion of ouliers 0.0350 0.0300 0.0250 0.0200 0.0150 0.0100 0.0050 0.0000 G S D 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 Lambda 9

Fig 3.4: Proporion of ouliers deeced by EWMAS, GVA, SdV and EWMAVD for monioring process scale parameer for C10N3 and n=10 Proporion of ouliers 0.0700 0.0600 0.0500 0.0400 0.0300 0.0200 0.0100 0.0000 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 G S D Lambda Fig 3.5: Proporion of ouliers deeced by EWMASM, GMA, SdM and EWMAMD chars for monioring process locaion parameer for N(0,1) and n=5 Proporion of ouliers 0.0035 0.0030 0.0025 0.0020 0.0015 0.0010 0.0005 0.0000 G S D 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 Lambda 10

Fig 3.6: Proporion of ouliers deeced by EWMASM, GMA, SdM and EWMAMD chars for monioring process locaion parameer for C10N3 and n=5 Proporion of ouliers 0.0080 0.0070 0.0060 0.0050 0.0040 0.0030 0.0020 0.0010 0.0000 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 Lambda G S D Fig 3.7: Proporion of ouliers deeced by EWMASM, GMA, SdM and EWMAMD chars for monioring process locaion parameer for C05N5 and n=10 0.0140 Proporion of ouliers 0.0120 0.0100 0.0080 0.0060 0.0040 0.0020 G S D 0.0000 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 Lambda 11

Fig 3.8: Proporion of ouliers deeced by EWMASM, GMA, SdM and EWMAMD chars for monioring process locaion parameer for C10N3 and n=10 0.0070 Proporion of ouliers 0.0060 0.0050 0.0040 0.0030 0.0020 0.0010 G S D 0.0000 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 Lambda The bars in he above chars represen he proporion of ouliers deeced using differen EWMA conrol schemes. The bars referred o as, G, S and D respecively represen he proporion of ouliers deeced by EWMAS, GVA, SdV and EWMAVD for scale parameer (in Figures 3.1-3.4) and EWMASM, GMA, SdM and EWMAMD for locaion parameer ( in Figures 3.5-3.8). From he above graphs i is clear ha: i) for in conrol siuaions (i.e. unconaminaed disribuions), he proporion of ouliers deeced (i.e. false alarm rae) by he proposed GVA do no exceed he respecive false alarm raes of EWMAS, EWMAVD and SdV chars. The same is he case wih GMA char as compared o EWMASM, EWMAMD and SdM chars; ii) for ou of conrol siuaions (i.e. conaminaed disribuions), he proporion of ouliers deeced by he proposed GVA char is higher han hose of EWMAS and SdV chars and almos equal o ha of 12

EWMAVD char. Similarly he proporion of ouliers deeced by he proposed GMA char is higher han hose of EWMASM and SdM chars, and almos equal o ha of EWMAMD char. All his is due o he fac ha he limis of he proposed chars based on G are less influenced by he presence of ouliers, which makes he proposed srucures of GVA and GMA chars more capable of deecing ouliers. Hence he proposed GVA char is superior o EWMAS and SdV chars and almos equally efficien o EWMAVD char for scale parameer, in erms of ouliers deecion. Also he proposed GMA char is superior o EWMASM and SdM chars and almos equally efficien o EWMAMD char for locaion parameer, in erms of ouliers deecion. 4 Conclusion The proposed GVA and GMA chars are more efficien, in erms of oulier deecion, han and S based EWMA design srucures and almos equally efficien o he Khoo (2004) Downon s esimaor based EWMA design srucures. 13

eferences David, H.A. (1968). Gini s mean difference rediscovered. Biomerika, 55, 573-575. Gini, C. (1912). Variabilia e Muabilia, Bologna. Khoo, M.B.C. (2004). Some Conrol Chars for he Process Mean and Variance based on Downon s Esimaor, Inernaional Engineering Managemen Conference 2004. Maravelakis, P.E., Panareos, J. and Psarakis, S. (2005). An Examinaion of he obusness o Non Normaliy of he EWMA Conrol Chars for he Dispersion, Communicaions in Saisics - Simulaion and Compuaion, Volume 34, Issue 4 Ocober 2005, pages 1069 1079. Mongomery, D. C. (2001). Inroducion o Saisical Qualiy Conrol, 4h ediion. Ng, C. H. and Case, K. E. (1989). Developmen and evaluaion of conrol char using exponenially weighed moving averages, Journal of Qualiy Technology, vol. 21, pp. 242-250. iaz, M. and Saghir, A. (2007). Monioring process variabiliy using Gini s mean difference, Quaniaive Managemen 4 (4), 439-455. ocke, D. M. (1992). 41, pp. 97-104. X Q and Q chars: robus conrol chars, The Saisician, vol. Soumbos,Z.G., Sullivan, J. H. (2002). obusness o non normaliy of he mulivariae EWMA conrol char, Journal of Qualiy Technology. 34:260-276. 14

Appendix Table 1: Conrol Char Consans of G -Chars n b 3 2 0.7345 3 0.5280 4 0.4222 5 0.3641 6 0.3244 7 0.2962 8 0.2754 9 0.2545 10 0.2436 11 0.2301 12 0.2185 13 0.2095 14 0.2014 15 0.1939 16 0.1873 17 0.1818 18 0.1779 19 0.1649 20 0.1601 15

Table 2: Facors for G- Chars Probabiliy Limi n M 0.001 M 0.005 M 0.01 M 0.05 M 0.10 M 0.20 M 0.25 M 0.50 M 0.75 M 0.80 M 0.90 M 0.95 M 0.99 M 0.995 M 0.999 2 0.002 0.009 0.018 0.076 0.154 0.308 0.391 0.843 1.444 1.615 2.086 2.478 3.333 3.644 4.510 3 0.035 0.084 0.117 0.262 0.368 0.537 0.608 0.940 1.343 1.442 1.721 1.954 2.449 2.649 3.081 4 0.092 0.165 0.202 0.358 0.466 0.614 0.676 0.948 1.255 1.339 1.565 1.750 2.151 2.286 2.700 5 0.164 0.257 0.300 0.454 0.552 0.682 0.736 0.972 1.233 1.300 1.493 1.656 1.959 2.064 2.311 6 0.229 0.302 0.352 0.509 0.598 0.725 0.774 0.985 1.210 1.267 1.436 1.574 1.852 1.923 2.160 7 0.299 0.365 0.405 0.543 0.626 0.742 0.787 0.986 1.195 1.249 1.394 1.517 1.753 1.834 2.010 8 0.303 0.3902 0.434 0.572 0.647 0.759 0.805 0.988 1.183 1.232 1.362 1.464 1.697 1.781 1.996 9 0.349 0.421 0.464 0.603 0.678 0.783 0.824 0.992 1.166 1.211 1.338 1.437 1.641 1.732 1.899 10 0.377 0.460 0.501 0.624 0.694 0.794 0.831 0.994 1.161 1.206 1.324 1.419 1.627 1.709 1.868 11 0.395 0.476 0.508 0.637 0.707 0.798 0.834 0.996 1.148 1.187 1.295 1.393 1.576 1.646 1.801 12 0.427 0.501 0.539 0.652 0.725 0.813 0.847 0.991 1.142 1.178 1.286 1.375 1.552 1.623 1.763 13 0.445 0.525 0.567 0.674 0.735 0.819 0.854 0.995 1.141 1.176 1.277 1.359 1.520 1.610 1.750 14 0.473 0.538 0.577 0.685 0.750 0.832 0.865 0.998 1.135 1.171 1.262 1.340 1.513 1.592 1.769 15 0.483 0.542 0.582 0.693 0.756 0.835 0.867 0.997 1.131 1.162 1.254 1.336 1.471 1.534 1.696 16 0.507 0.561 0.605 0.711 0.770 0.843 0.872 0.995 1.127 1.160 1.252 1.324 1.467 1.517 1.630 17 0.513 0.576 0.618 0.719 0.781 0.848 0.876 0.995 1.120 1.153 1.245 1.313 1.453 1.511 1.623 18 0.521 0.589 0.621 0.724 0.784 0.853 0.879 0.996 1.119 1.147 1.227 1.304 1.442 1.508 1.609 19 0.527 0.602 0.639 0.728 0.789 0.856 0.880 0.997 1.118 1.140 1.210 1.293 1.428 1.483 1.591 20 0.538 0.615 0.643 0.739 0.791 0.861 0.885 0.998 1.116 1.139 1.209 1.280 1.410 1.459 1.580 16

Table 3: Proporion of ou-of-conrol for EWMAS, GVT, SdV and EWMAVD Chars for monioring process scale parameer Sae of Conrol EWMAS n In Conrol Siuaion ou of Conrol Siuaion 5 10 5 10 Disribuion N(0,1) N(0,1) C05N3 C10N3 C05N5 C10N5 C05N3 C10N3 C05N5 C10N5 λ 0.1 0.0005 0.0004 0.0045 0.0080 0.0187 0.0259 0.0112 0.0133 0.0489 0.0501 0.2 0.0007 0.0006 0.0057 0.0093 0.0225 0.0290 0.0119 0.0188 0.0585 0.0668 0.3 0.0009 0.0008 0.0112 0.0127 0.0320 0.0412 0.0192 0.0243 0.0692 0.0807 0.4 0.0009 0.0008 0.0156 0.0190 0.0402 0.0501 0.0218 0.0319 0.0869 0.0957 0.5 0.0012 0.0009 0.0163 0.0217 0.0477 0.0582 0.0295 0.0380 0.0899 0.1018 0.6 0.0012 0.0020 0.0211 0.0252 0.0541 0.0645 0.0386 0.0385 0.0902 0.1082 0.7 0.0019 0.0022 0.0218 0.0288 0.0547 0.0738 0.0419 0.0481 0.0918 0.1096 0.8 0.0022 0.0025 0.0230 0.0317 0.0563 0.0740 0.0428 0.0515 0.0935 0.1098 0.9 0.0023 0.0025 0.0230 0.0372 0.0575 0.0748 0.0430 0.0517 0.0958 0.1101 1 0.0032 0.0030 0.0246 0.0385 0.0621 0.0798 0.0440 0.0536 0.1015 0.1130 GVT SdV EWMAVD λ λ λ 0.1 0.0004 0.0003 0.0069 0.0104 0.0268 0.0384 0.0137 0.0179 0.0841 0.0991 0.2 0.0005 0.0004 0.0080 0.0106 0.0342 0.0465 0.0161 0.0285 0.0977 0.1323 0.3 0.0005 0.0004 0.0123 0.0171 0.0453 0.0619 0.0251 0.0350 0.1015 0.1383 0.4 0.0006 0.0005 0.0170 0.0229 0.0549 0.0693 0.0295 0.0422 0.1047 0.1402 0.5 0.0011 0.0009 0.0191 0.0279 0.0601 0.0762 0.0322 0.0458 0.1104 0.1409 0.6 0.0011 0.0011 0.0219 0.0308 0.0634 0.0789 0.0394 0.0485 0.1124 0.1427 0.7 0.0016 0.0014 0.0223 0.0331 0.0635 0.0914 0.0424 0.0533 0.1143 0.1474 0.8 0.0018 0.0018 0.0246 0.0369 0.0640 0.0917 0.0434 0.0547 0.1159 0.1484 0.9 0.0022 0.0018 0.0267 0.0411 0.0645 0.0923 0.0439 0.0554 0.1161 0.1488 1 0.0028 0.0020 0.0308 0.0425 0.0677 0.0931 0.0451 0.0623 0.1177 0.1507 0.1 0.0004 0.0003 0.0056 0.0081 0.0201 0.0279 0.0087 0.0108 0.0464 0.0489 0.2 0.0004 0.0004 0.0062 0.0092 0.0243 0.0315 0.0094 0.0159 0.0561 0.0657 0.3 0.0005 0.0004 0.0118 0.0130 0.0342 0.0440 0.0148 0.0202 0.0656 0.0795 0.4 0.0006 0.0005 0.0158 0.0191 0.0401 0.0528 0.0166 0.0265 0.0811 0.0959 0.5 0.0011 0.0007 0.0168 0.0217 0.0455 0.0613 0.0217 0.0328 0.0841 0.1011 0.6 0.0011 0.0011 0.0209 0.0253 0.0528 0.0674 0.0306 0.0346 0.0897 0.1013 0.7 0.0017 0.0014 0.0210 0.0288 0.0541 0.0766 0.0359 0.0397 0.0909 0.1084 0.8 0.0020 0.0019 0.0226 0.0325 0.0574 0.0768 0.0361 0.0429 0.0931 0.1105 0.9 0.0022 0.0020 0.0230 0.0375 0.0579 0.0771 0.0365 0.0441 0.0952 0.1121 1 0.0028 0.0020 0.0283 0.0387 0.0627 0.0821 0.0381 0.0461 0.0972 0.1140 0.1 0.0004 0.0003 0.0069 0.0103 0.0267 0.0384 0.0137 0.0178 0.0839 0.0989 0.2 0.0006 0.0004 0.0079 0.0106 0.0340 0.0463 0.0160 0.0285 0.0975 0.1321 0.3 0.0008 0.0006 0.0121 0.0171 0.0453 0.0618 0.0250 0.0350 0.1014 0.1381 0.4 0.0008 0.0008 0.0169 0.0229 0.0547 0.0692 0.0293 0.0422 0.1046 0.1401 0.5 0.0011 0.0013 0.0189 0.0278 0.0597 0.0761 0.0321 0.0457 0.1102 0.1408 0.6 0.0011 0.0015 0.0219 0.0307 0.0633 0.0787 0.0394 0.0483 0.1123 0.1425 0.7 0.0016 0.0016 0.0223 0.0331 0.0635 0.0912 0.0423 0.0532 0.1143 0.1473 0.8 0.0021 0.0018 0.0244 0.0367 0.0639 0.0912 0.0433 0.0547 0.1157 0.1483 0.9 0.0022 0.0020 0.0267 0.0410 0.0643 0.0922 0.0439 0.0553 0.1160 0.1488 1 0.0028 0.0020 0.0307 0.0425 0.0676 0.0931 0.0451 0.0622 0.1176 0.1505 17

Table 4: Proporion of ou-of-conrol for EWMASM, GMT, SdM and EWMAMD Chars for monioring process locaion parameer Sae of Conrol EWMASM n In Conrol Siuaion ou of Conrol Siuaion 5 10 5 10 Disribuion N(0,1) N(0,1) C05N3 C10N3 C05N5 C10N5 C05N3 C10N3 C05N5 C10N5 λ 0.1 0.0005 0.0004 0.0008 0.0011 0.0020 0.0021 0.0003 0.0005 0.0001 0.0004 0.2 0.0006 0.0005 0.0010 0.0012 0.0021 0.0023 0.0004 0.0007 0.0010 0.0010 0.3 0.0008 0.0005 0.0010 0.0013 0.0025 0.0033 0.0006 0.0007 0.0013 0.0013 0.4 0.0008 0.0007 0.0016 0.0021 0.0039 0.0046 0.0008 0.0009 0.0016 0.0019 0.5 0.0011 0.0008 0.0023 0.0023 0.0062 0.0068 0.0010 0.0011 0.0023 0.0026 0.6 0.0018 0.0011 0.0029 0.0039 0.0076 0.0087 0.0014 0.0021 0.0035 0.0036 0.7 0.0021 0.0014 0.0037 0.0048 0.0085 0.0112 0.0021 0.0023 0.0047 0.0049 0.8 0.0024 0.0016 0.0043 0.0051 0.0115 0.0120 0.0024 0.0027 0.0053 0.0054 0.9 0.0027 0.0021 0.0050 0.0063 0.0114 0.0142 0.0028 0.0030 0.0065 0.0065 1 0.0029 0.0023 0.0057 0.0067 0.0141 0.0157 0.0030 0.0036 0.0066 0.0066 GMT SdM EWMAMD λ λ λ 0.1 0.0005 0.0003 0.0009 0.0014 0.0025 0.0027 0.0005 0.0009 0.0022 0.0024 0.2 0.0006 0.0005 0.0011 0.0014 0.0025 0.0028 0.0007 0.0011 0.0023 0.0026 0.3 0.0007 0.0005 0.0012 0.0015 0.0033 0.0044 0.0009 0.0011 0.0045 0.0034 0.4 0.0008 0.0006 0.0017 0.0024 0.0050 0.0057 0.0014 0.0016 0.0037 0.0048 0.5 0.0012 0.0008 0.0026 0.0027 0.0077 0.0087 0.0016 0.0020 0.0050 0.0068 0.6 0.0016 0.0010 0.0031 0.0044 0.0093 0.0114 0.0023 0.0035 0.0073 0.0091 0.7 0.0020 0.0013 0.0039 0.0054 0.0108 0.0141 0.0033 0.0039 0.0089 0.0115 0.8 0.0023 0.0016 0.0048 0.0058 0.0138 0.0149 0.0034 0.0045 0.0107 0.0131 0.9 0.0028 0.0020 0.0056 0.0071 0.0137 0.0177 0.0042 0.0047 0.0120 0.0135 1 0.0027 0.0022 0.0060 0.0074 0.0169 0.0194 0.0043 0.0058 0.0128 0.0143 0.1 0.0005 0.0004 0.0008 0.0011 0.0020 0.0020 0.0004 0.0006 0.0010 0.0011 0.2 0.0006 0.0004 0.0010 0.0012 0.0020 0.0021 0.0005 0.0009 0.0016 0.0016 0.3 0.0007 0.0006 0.0011 0.0013 0.0024 0.0032 0.0008 0.0012 0.0021 0.0023 0.4 0.0008 0.0007 0.0015 0.0021 0.0038 0.0043 0.0011 0.0017 0.0025 0.0024 0.5 0.0011 0.0010 0.0023 0.0023 0.0060 0.0067 0.0015 0.0020 0.0029 0.0030 0.6 0.0017 0.0015 0.0029 0.0390 0.0073 0.0850 0.0016 0.0022 0.0031 0.0031 0.7 0.0019 0.0016 0.0036 0.0048 0.0085 0.0109 0.0020 0.0023 0.0048 0.0050 0.8 0.0024 0.0016 0.0043 0.0052 0.0112 0.0116 0.0023 0.0026 0.0051 0.0052 0.9 0.0028 0.0020 0.0050 0.0061 0.0113 0.0139 0.0025 0.0028 0.0063 0.0063 1 0.0027 0.0023 0.0056 0.0066 0.0139 0.0155 0.0029 0.0034 0.0064 0.0065 0.1 0.0005 0.0003 0.0009 0.0014 0.0025 0.0027 0.0005 0.0009 0.0022 0.0024 0.2 0.0006 0.0005 0.0011 0.0014 0.0025 0.0028 0.0007 0.0011 0.0022 0.0026 0.3 0.0007 0.0005 0.0012 0.0015 0.0033 0.0044 0.0009 0.0011 0.0027 0.0034 0.4 0.0008 0.0006 0.0017 0.0024 0.0050 0.0057 0.0014 0.0016 0.0037 0.0048 0.5 0.0012 0.0008 0.0026 0.0027 0.0077 0.0086 0.0016 0.0020 0.0050 0.0068 0.6 0.0016 0.0010 0.0031 0.0045 0.0093 0.0114 0.0023 0.0035 0.0073 0.0091 0.7 0.0020 0.0013 0.0039 0.0054 0.0108 0.0141 0.0033 0.0039 0.0089 0.0115 0.8 0.0024 0.0016 0.0047 0.0058 0.0138 0.0149 0.0034 0.0045 0.0107 0.0131 0.9 0.0027 0.0020 0.0055 0.0071 0.0137 0.0176 0.0042 0.0047 0.0120 0.0135 1 0.0028 0.0022 0.0060 0.0074 0.0168 0.0194 0.0043 0.0058 0.0128 0.0143 18