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

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

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

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

Vehicle Arrival Models : Headway

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

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

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

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

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

Development of a new metrological model for measuring of the water surface evaporation Tovmach L. Tovmach Yr. Abstract Introduction

20. Applications of the Genetic-Drift Model

A Dynamic Model of Economic Fluctuations

Appendix to Creating Work Breaks From Available Idleness

CONFIDENCE LIMITS AND THEIR ROBUSTNESS

School and Workshop on Market Microstructure: Design, Efficiency and Statistical Regularities March 2011

18 Biological models with discrete time

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

Outline. lse-logo. Outline. Outline. 1 Wald Test. 2 The Likelihood Ratio Test. 3 Lagrange Multiplier Tests

Internet Traffic Modeling for Efficient Network Research Management Prof. Zhili Sun, UniS Zhiyong Liu, CATR

An Improved Adaptive CUSUM Control Chart for Monitoring Process Mean

2. Nonlinear Conservation Law Equations

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

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

OBJECTIVES OF TIME SERIES ANALYSIS

Introduction to Probability and Statistics Slides 4 Chapter 4

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

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

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

Bias in Conditional and Unconditional Fixed Effects Logit Estimation: a Correction * Tom Coupé

Air Traffic Forecast Empirical Research Based on the MCMC Method

The General Linear Test in the Ridge Regression

The Arcsine Distribution

Types of Exponential Smoothing Methods. Simple Exponential Smoothing. Simple Exponential Smoothing

IMPROVED HYPERBOLIC FUNCTION METHOD AND EXACT SOLUTIONS FOR VARIABLE COEFFICIENT BENJAMIN-BONA-MAHONY-BURGERS EQUATION

A unit root test based on smooth transitions and nonlinear adjustment

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

MODULE 3 FUNCTION OF A RANDOM VARIABLE AND ITS DISTRIBUTION LECTURES PROBABILITY DISTRIBUTION OF A FUNCTION OF A RANDOM VARIABLE

Effects of Coordinate Curvature on Integration

The equation to any straight line can be expressed in the form:

DEPARTMENT OF STATISTICS

R t. C t P t. + u t. C t = αp t + βr t + v t. + β + w t

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

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

Christos Papadimitriou & Luca Trevisan November 22, 2016

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

FITTING OF A PARTIALLY REPARAMETERIZED GOMPERTZ MODEL TO BROILER DATA

SPC Procedures for Monitoring Autocorrelated Processes

Continuous Time. Time-Domain System Analysis. Impulse Response. Impulse Response. Impulse Response. Impulse Response. ( t) + b 0.

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

A First Course on Kinetics and Reaction Engineering. Class 19 on Unit 18

Accurate RMS Calculations for Periodic Signals by. Trapezoidal Rule with the Least Data Amount

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010

ACE 562 Fall Lecture 4: Simple Linear Regression Model: Specification and Estimation. by Professor Scott H. Irwin

Regression with Time Series Data

A Specification Test for Linear Dynamic Stochastic General Equilibrium Models

Introduction D P. r = constant discount rate, g = Gordon Model (1962): constant dividend growth rate.

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

Reliability of Technical Systems

How to Deal with Structural Breaks in Practical Cointegration Analysis

On a Discrete-In-Time Order Level Inventory Model for Items with Random Deterioration

EE 435. Lecture 31. Absolute and Relative Accuracy DAC Design. The String DAC

Nonlinearity Test on Time Series Data

Some Basic Information about M-S-D Systems

Two Coupled Oscillators / Normal Modes

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

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

5.2. The Natural Logarithm. Solution

Estimation Uncertainty

EE 435. Lecture 35. Absolute and Relative Accuracy DAC Design. The String DAC

6. COMPUTATION OF CENTILES AND Z-SCORES FOR VELOCITIES BASED ON WEIGHT, LENGTH AND HEAD CIRCUMFERENCE

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

1. Diagnostic (Misspeci cation) Tests: Testing the Assumptions

A Note on the Equivalence of Fractional Relaxation Equations to Differential Equations with Varying Coefficients

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

Average Number of Lattice Points in a Disk

Kriging Models Predicting Atrazine Concentrations in Surface Water Draining Agricultural Watersheds

Maintenance Models. Prof. Robert C. Leachman IEOR 130, Methods of Manufacturing Improvement Spring, 2011

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

Measurement Error 1: Consequences Page 1. Definitions. For two variables, X and Y, the following hold: Expectation, or Mean, of X.

Lecture 33: November 29

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

On the Separation Theorem of Stochastic Systems in the Case Of Continuous Observation Channels with Memory

1 Differential Equation Investigations using Customizable

EE 315 Notes. Gürdal Arslan CLASS 1. (Sections ) What is a signal?

Summer Term Albert-Ludwigs-Universität Freiburg Empirische Forschung und Okonometrie. Time Series Analysis

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

Research Report Statistical Research Unit Department of Economics University of Gothenburg

Physics 127b: Statistical Mechanics. Fokker-Planck Equation. Time Evolution

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

Econ Autocorrelation. Sanjaya DeSilva

Final Spring 2007

not to be republished NCERT MATHEMATICAL MODELLING Appendix 2 A.2.1 Introduction A.2.2 Why Mathematical Modelling?

Chapter 15. Time Series: Descriptive Analyses, Models, and Forecasting

IMPLICIT AND INVERSE FUNCTION THEOREMS PAUL SCHRIMPF 1 OCTOBER 25, 2013

Chapter 2. Models, Censoring, and Likelihood for Failure-Time Data

Distribution of Estimates

1 Evaluating Chromatograms

Lecture Notes 2. The Hilbert Space Approach to Time Series

Process Monitoring versus Process Adjustment

Ordinary Differential Equations

Transcription:

DOI: 0.545/mjis.07.5009 Exponenial Weighed Moving Average (EWMA) Char Under The Assumpion of Moderaeness And Is 3 Conrol Limis KALPESH S TAILOR Assisan Professor, Deparmen of Saisics, M. K. Bhavnagar Universiy, Bhavnagar-36400 Email: kalpesh_lr@yahoo.co.in Received: Ocober, 06 Revised: November 9, 06 Acceped: December 5, 06 Published online: March 05, 07 The Auhor(s) 06. This aricle is published wih open access a www.chikara.edu.in/publicaions Absrac Moderae disribuion is a very good alernaive of normal disribuion proposed by Naik V.D and Desai J.M. [4], which has mean deviaion as scale parameer raher han he sandard deviaion. Mean deviaion (δ) is a very good alernaive of sandard deviaion (σ) as mean deviaion is considered o be he mos inuiively and raionally defined measure of dispersion. This fac can be very useful in he field of qualiy conrol o consruc he conrol limis of he conrol chars. On he basis of his fac Naik V.D. and Tailor K.S. [5] have proposed 3δ conrol limis. In 3δ conrol limis, he upper and lower conrol limis are se a 3δ disance from he cenral line where δ is he mean deviaion of sampling disribuion of he saisic being used for consrucing he conrol char. In his paper i has been assumed ha he underlying disribuion of he variable of ineres follows moderae disribuion proposed by Naik V.D and Desai J.M. [4] and 3δ conrol limis of exponenial weighed moving average char are derived. Also an empirical sudy is carried ou o illusrae he use hese chars. Keywords Mean deviaion, Moderae disribuion, Exponenial weighed moving average, 3δ conrol limis.. INTRODUCTION A fundamenal assumpion in he developmen of conrol chars for variables is ha he underlying disribuion of he concerned qualiy characerisic is normal. The normal disribuion is one of he mos imporan disribuions in he saisical inference in which mean ( µ ) and sandard deviaion (σ) are Mahemaical Journal of Inerdisciplinary Sciences Vol-5, No-, March 07 pp. 9

Tailor, KS he parameers. Naik V.D and Desai J.M. [4] have suggesed an alernaive of normal disribuion, which is called moderae disribuion. In moderae disribuion mean ( µ ) and mean deviaion (δ) are he pivoal parameers and, hey have properies similar o normal disribuion. Naik V.D. and Tailor K.S. [5] have proposed he concep of 3δ conrol limis on he basis of moderae disribuion. Under his rule, he upper and lower conrol limis are se a 3δ disance from he cenral line where δ is he mean deviaion of sampling disribuion of he saisic being used for consrucing he conrol char. Thus in he proposed conrol chars, under he moderaeness assumpion, he conrol limis for any saisic T should be deermined as follows. Cenral line (CL) Expeced value of T μ Lower Conrol Limi (LCL) Mean of T -3δ T µ -3δ T Upper Conrol Limi (UCL) Mean of T + 3δ T µ + 3δ T Where μ is mean of T and δ T is he mean deviaion of T. I is found ha since δ provides exac average disance from mean and σ provides only an approximae average disance, 3δ limis can be considered o be more raional as compared o 3σ limis. Naik and Tailor have derived 3δ conrol limis of X -char, R-char, s (sample sandard deviaion) char and d (sample mean deviaion) char. Tailor has also suggesed 3δ conrol limis of moving average and moving range chars. Thus, in his paper i is assumed ha he underlying disribuion of he concerned variable follows moderae disribuion and 3δ conrol limis for exponenial weighed moving average is derived. An empirical sudy is also carried ou o illusrae he use of he char.. EXPONENTIAL WEIGHTED MOVING AVERAGE (EWMA) CHART The concep of EWMA char was inroduced by Robers S.W [7]. The exponenially weighed moving average char is a ype of moving mean char in which an exponenially weighed mean is calculaed each ime a new resul becomes available. The EWMA conrol char is a very good alernaive o he Shewhar char, when we are ineresed in deecing small shifs.

New weighed mean ( new resul) + ((- ) previous resul ), where is he smoohing consan. I has a value beween 0 and, many saisician use 0., bu choice of has o be lef o he judgmen of he qualiy conrol specialis, he smaller he value of, he greaer he influence of he hisorical daa. The EWMA char is much more effecive han moving average char for deecing small shifs. If i is imporan o recognize small shifs early in he process, hen he value of should be small. If, he EWMA char reduces iself o he usual X -char. This has been used by some organizaions, paricularly in he process indusries, as he basis of new conrol (performance) char sysems. Grea care mus be aken when using hese sysems since hey do no show changes in variabiliy very well and he basis for weighing daa is ofen eiher quesionable or arbirary. Exponenial Weighed Moving Average (EWMA) Char Under The Assumpion of Moderaeness And Is 3 Conrol Limis 3. (3δ) CONTROL LIMITS FOR EXPONENTIAL WEIGHTED MOVING AVERAGE (EWMA) CHART Suppose a measurable qualiy characerisic of he produc is denoed by X. Suppose ha m samples, each of size n, are drawn a more or less regular inerval of ime from he producion processes. These samples are known as subgroups, and for each of hese subgroups he values of exponenial weighed moving mean X are obained. Le he disribuion of he variable X be moderae wih mean µ and mean deviaion δ, hen, as proved by Naik V.D and Desai J.M. [4], he disribuion of X is also moderae wih mean µ and δ mean deviaion. Furher, if he disribuion of X is no moderae, and he n number of unis in each subgroup is 4 or more, hen on he basis of cenral limi heorem for moderae disribuion, i can be said ha X follows moderae disribuion. The EWMA funcion is defined as, Z xi + ( - ) Zi -, where 0< If he individual X are independen random variables wih variance σ n, hen he variance of Z is defined as σ σ α - - ( - ) 3

Tailor, KS Therefore σ σ α - - ( - ) σ - - ( - ) α () Since we are assuming moderaeness, he mean error of Z is defined as δ π δ α - - - ( ) () Thus, he 3δ- conrol limis of EWMA char can be deermined as follows Cenral line (C.L) E( X ) Lower conrol limi (L.C.L) E( X )- X (3) 3δ X -3 π δ - n - ( - α) δ X -5. 384 α - - ( - ) (4) Upper conrol limi (U.C.L) E( X )+ 3δ π δ + - - - X 3 ( α) n (5) δ X + 5. 384 α - - - ( ) Where X and δ are ypically esimaed from preliminary daa as sample mean and sample mean deviaion. As α is small and if increases, he effec of saring value soon dissipaes and he mean error converges o is asympoic value. 4

i.e δ π δ - The conrol limis for EWMA char are usually based on he asympoic mean deviaion of he saisic. Hence asympoic 3δ-conrol limis for his char can be derived as following way, Cenral line (C.L) E ( X ) Exponenial Weighed Moving Average (EWMA) Char Under The Assumpion of Moderaeness And Is 3 Conrol Limis Lower conrol limi (L.C.L) E( X )- X (6) 3δ X -3 X -5. 384 π δ - δ - (7) Upper conrol limi (U.C.L) E( X )+ 3δ X + 3 X + 5. 384 π δ - δ - (8) 4. AN EMPIRICAL STUDY FOR EWMA CHART: To illusrae he use of EWMA conrol scheme, we use a se of simulaed observaions aken from Lucas J. M and Crosier R.B []. The daa, ogeher wih he corresponding EWMA values, are shown in Table. The arge value is aken o be 0, so he process is in conrol for he firs 0 observaions. The mean level was shifed upward by approximaely one sandard deviaion for he las nine observaions. The parameers of he EWMA are chosen o be α 0.5, δ, n Asympoic 3-conrol limis for EWMA char are obained by 5

Tailor, KS Table : Observed value EWMA + (- ) 0-0.0.0 0.5-0.5 0.063 3 0.0 0.047 4-0.8-0.65 5-0.8-0.34 6 -. -0.543 7.5-0.03 8-0.6-0.74 9.0 0.9 0-0.9-0.35. 0.98 0.5 0.74 3.6 0.855 4 0.7 0.87 5. 0.887 6.0.66 7.4.4 8.9.393 9 0.8.45 LCL -.00, CL 0 and UCL.00 Similarly, asympoic 3σ-conrol limis for EWMA char are obained by LCL -.34, CL 0 and UCL.34 Now, EWMA char for moderaeness and normaliy assumpions can be consruced as follows. From figure, i can be seen ha under moderaeness assumpion wih 3δ -conrol limis, EWMA char is under he saisical conrol, while under normaliy assumpion wih 3σ-conrol limis, i is ou of conrol as one sample poin falls ouside he UCL. The poin which shows ou of conrol siuaion in EWMA char under normaliy assumpion shows under conrol siuaion in EWMA char under moderaeness assumpion. 6

Exponenial Weighed Moving Average (EWMA) Char Under The Assumpion of Moderaeness And Is 3 Conrol Limis Figure : I can also be seen from figure, ha using asympoic conrol limis raher han he ime varying limis, makes he EWMA char much less sensiive o process shifs in he firs few observaions. This could be a significan problem if a large shif occurs early, or if afer an ou of conrol condiion he process is no properly rese. 5. SUMMARY This paper derives he EWMA char under he assumpion of moderaeness. The 3δ - conrol limis are derived for he EWMA char wih ime varying conrol limis. Also he asympoic 3δ - conrol limis are derived for he same char. An empirical sudy is carried ou o illusrae his char. A comparaive sudy is carried ou for he EWMA char under moderaeness assumpion and under normaliy assumpion and i is found ha EWMA char under moderaeness assumpion perform beer han EWMA char under normaliy assumpion. Hence, i is suggesed ha EWMA char under moderaeness assumpion wih 3δ - conrol limis should be used for effecive performance of he char. 6. APPENDIX (a) Moderae Disribuion Suppose he p.d.f. of a disribuion of a random variable X is defined as, 7

Tailor, KS X -µ - π δ f( x) e,- < X <, δ > 0 πδ Then, he random variable X may be said o be following moderae disribuion wih parameers μ and δ and may be denoed as X M( µδ, ). I can be proved ha, (i) f ( x ) - (ii) Mean E(x) μ (iii) Mean deviaion E X -π δ (vi) Sandard deviaion π δ () x (v) M.G.F M e π π+ δ 4 (vi) f ( µ - X) f µ + x ( ) (b) 3σ-conrol limis of EWMA char C.L. X U.C.L. L.C.L. X + 3σ - - ( - α) X - - - - 3σ ( α) REFERENCES [] Huner J. S. (986) The Exponenially Weighed Moving Average, Journal of Qualiy Technology, 8, 03 0. [] Lucas J.M. and Crosier R.B. (98) Fas Iniial Response for CUSUM Qualiy Conrol Schemes, Technomerics, 4, 99 05 [3] Lucas J.M. and Saccucci M.S. (990) Exponenially Weighed Moving average Schemes, Properies and Enhancemens, Technomerics, 3, 9 [4] Naik V.D and Desai J.M. (05) Moderae Disribuion: A modified normal disribuion which has Mean as locaion parameer and Mean Deviaion as scale parameer, VNSGU Journal of Science and Technology, Vol.4, No. 56 70 8

[5] Naik V.D and Tailor K.S. (05) On performance of and R-chars under he assumpion of moderaeness raher han normaliy and wih 3 conrol limis raher han 3 conrol limis, VNSGU Journal of Science and Technology, Vol.4, No., 43 55 [6] Kalpesh S. Tailor (06) Moving average and moving range chars under he assumpion of moderaeness and is 3 conrol limis [7] Robers S.W. (959) Conrol char Tess Based on Geomeric Moving Average Chars. Technomerics, Vol.-, No.-3, pp.39 50 [8] Tailor K.S. and Naik V.D. (06) Mean deviaion () based conrol limis of SQC chars for sample sandard deviaion(s) and sample mean deviaion (d) and heir performance analysis under 3 conrol limis agains 3 conrol limis, VNSGU Journal of Science and Technology, (Acceped for Publicaion). Exponenial Weighed Moving Average (EWMA) Char Under The Assumpion of Moderaeness And Is 3 Conrol Limis 9