GINI MEAN DIFFERENCE AND EWMA CHARTS. Muhammad Riaz, Department of Statistics, Quaid-e-Azam University Islamabad,
|
|
- Egbert Elliott
- 5 years ago
- Views:
Transcription
1 GINI MEAN DIFFEENCE AND EWMA CHATS Muhammad iaz, Deparmen of Saisics, Quaid-e-Azam Universiy Islamabad, Pakisan. Saddam Akbar Abbasi, Deparmen of Saisics, Quaid-e-Azam Universiy Islamabad, Pakisan. 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
2 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 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
4 ˆ = + 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
5 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
6 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
7 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
8 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 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 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 G S D Lambda 8
9 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 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 G S D Lambda 9
10 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 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 G S D Lambda 10
11 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 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= Proporion of ouliers G S D Lambda 11
12 Fig 3.8: Proporion of ouliers deeced by EWMASM, GMA, SdM and EWMAMD chars for monioring process locaion parameer for C10N3 and n= Proporion of ouliers G S D 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 ) and EWMASM, GMA, SdM and EWMAMD for locaion parameer ( in Figures ). 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
13 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
14 eferences David, H.A. (1968). Gini s mean difference rediscovered. Biomerika, 55, 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 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 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 iaz, M. and Saghir, A. (2007). Monioring process variabiliy using Gini s mean difference, Quaniaive Managemen 4 (4), ocke, D. M. (1992). 41, pp 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:
15 Appendix Table 1: Conrol Char Consans of G -Chars n b
16 Table 2: Facors for G- Chars Probabiliy Limi n M M 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 M
17 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 Disribuion N(0,1) N(0,1) C05N3 C10N3 C05N5 C10N5 C05N3 C10N3 C05N5 C10N5 λ GVT SdV EWMAVD λ λ λ
18 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 Disribuion N(0,1) N(0,1) C05N3 C10N3 C05N5 C10N5 C05N3 C10N3 C05N5 C10N5 λ GMT SdM EWMAMD λ λ λ
A Robust Exponentially Weighted Moving Average Control Chart for the Process Mean
Journal of Modern Applied Saisical Mehods Volume 5 Issue Aricle --005 A Robus Exponenially Weighed Moving Average Conrol Char for he Process Mean Michael B. C. Khoo Universii Sains, Malaysia, mkbc@usm.my
More informationExponential Weighted Moving Average (EWMA) Chart Under The Assumption of Moderateness And Its 3 Control Limits
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,
More informationExponentially Weighted Moving Average (EWMA) Chart Based on Six Delta Initiatives
hps://doi.org/0.545/mjis.08.600 Exponenially Weighed Moving Average (EWMA) Char Based on Six Dela Iniiaives KALPESH S. TAILOR Deparmen of Saisics, M. K. Bhavnagar Universiy, Bhavnagar-36400 E-mail: kalpesh_lr@yahoo.co.in
More informationVehicle Arrival Models : Headway
Chaper 12 Vehicle Arrival Models : Headway 12.1 Inroducion Modelling arrival of vehicle a secion of road is an imporan sep in raffic flow modelling. I has imporan applicaion in raffic flow simulaion where
More informationRobust estimation based on the first- and third-moment restrictions of the power transformation model
h Inernaional Congress on Modelling and Simulaion, Adelaide, Ausralia, 6 December 3 www.mssanz.org.au/modsim3 Robus esimaion based on he firs- and hird-momen resricions of he power ransformaion Nawaa,
More informationSTRUCTURAL CHANGE IN TIME SERIES OF THE EXCHANGE RATES BETWEEN YEN-DOLLAR AND YEN-EURO IN
Inernaional Journal of Applied Economerics and Quaniaive Sudies. Vol.1-3(004) STRUCTURAL CHANGE IN TIME SERIES OF THE EXCHANGE RATES BETWEEN YEN-DOLLAR AND YEN-EURO IN 001-004 OBARA, Takashi * Absrac The
More informationCOMPUTATION OF THE PERFORMANCE OF SHEWHART CONTROL CHARTS. Pieter Mulder, Julian Morris and Elaine B. Martin
COMUTATION OF THE ERFORMANCE OF SHEWHART CONTROL CHARTS ieer Mulder, Julian Morris and Elaine B. Marin Cenre for rocess Analyics and Conrol Technology, School of Chemical Engineering and Advanced Maerials,
More informationAn Improved Adaptive CUSUM Control Chart for Monitoring Process Mean
An Improved Adapive CUSUM Conrol Char for Monioring Process Mean Jun Du School of Managemen Tianjin Universiy Tianjin 37, China Zhang Wu, Roger J. Jiao School of Mechanical and Aerospace Engineering Nanyang
More informationTHE CUSUM VERSUS MCUSUM MODIFIED CONTROL CHARTS WHEN APPLIED ON DIESEL ENGINES PARAMETERS CONTROL
Proceedings of he 6h Inernaional Conference on Mechanics and Maerials in Design, Ediors: J.F. Silva Gomes & S.A. Meguid, P.Delgada/Azores, 26-30 July 2015 PAPER REF: 5649 THE CUSUM VERSUS MCUSUM MODIFIED
More informationTesting the Random Walk Model. i.i.d. ( ) r
he random walk heory saes: esing he Random Walk Model µ ε () np = + np + Momen Condiions where where ε ~ i.i.d he idea here is o es direcly he resricions imposed by momen condiions. lnp lnp µ ( lnp lnp
More informationSmoothing. Backward smoother: At any give T, replace the observation yt by a combination of observations at & before T
Smoohing Consan process Separae signal & noise Smooh he daa: Backward smooher: A an give, replace he observaion b a combinaion of observaions a & before Simple smooher : replace he curren observaion wih
More informationJ. Martin van Zyl Department of Mathematical Statistics and Actuarial Science, University of the Free State, PO Box 339, Bloemfontein, South Africa
A weighed leas squares procedure o approximae leas absolue deviaion esimaion in ime series wih specific reference o infinie variance uni roo problems J. Marin van Zyl Deparmen of Mahemaical Saisics and
More informationResearch Report Statistical Research Unit Department of Economics University of Gothenburg
Research Repor Saisical Research Uni Deparmen of Economics Universiy of Gohenburg Sweden Hoelling s T Mehod in Mulivariae On-Line Surveillance. On he Delay of an Alarm E. Andersson Research Repor 008:3
More informationAn Adaptive Generalized Likelihood Ratio Control Chart for Detecting an Unknown Mean Pattern
An adapive GLR conrol char... 1/32 An Adapive Generalized Likelihood Raio Conrol Char for Deecing an Unknown Mean Paern GIOVANNA CAPIZZI and GUIDO MASAROTTO Deparmen of Saisical Sciences Universiy of Padua
More informationTesting for a Single Factor Model in the Multivariate State Space Framework
esing for a Single Facor Model in he Mulivariae Sae Space Framework Chen C.-Y. M. Chiba and M. Kobayashi Inernaional Graduae School of Social Sciences Yokohama Naional Universiy Japan Faculy of Economics
More informationKriging Models Predicting Atrazine Concentrations in Surface Water Draining Agricultural Watersheds
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 Kriging Models Predicing Arazine Concenraions in Surface Waer Draining Agriculural Waersheds Paul L. Mosquin, Jeremy Aldworh, Wenlin Chen Supplemenal Maerial Number
More informationExcel-Based Solution Method For The Optimal Policy Of The Hadley And Whittin s Exact Model With Arma Demand
Excel-Based Soluion Mehod For The Opimal Policy Of The Hadley And Whiin s Exac Model Wih Arma Demand Kal Nami School of Business and Economics Winson Salem Sae Universiy Winson Salem, NC 27110 Phone: (336)750-2338
More informationDiebold, Chapter 7. Francis X. Diebold, Elements of Forecasting, 4th Edition (Mason, Ohio: Cengage Learning, 2006). Chapter 7. Characterizing Cycles
Diebold, Chaper 7 Francis X. Diebold, Elemens of Forecasing, 4h Ediion (Mason, Ohio: Cengage Learning, 006). Chaper 7. Characerizing Cycles Afer compleing his reading you should be able o: Define covariance
More informationArticle from. Predictive Analytics and Futurism. July 2016 Issue 13
Aricle from Predicive Analyics and Fuurism July 6 Issue An Inroducion o Incremenal Learning By Qiang Wu and Dave Snell Machine learning provides useful ools for predicive analyics The ypical machine learning
More informationReliability of Technical Systems
eliabiliy of Technical Sysems Main Topics Inroducion, Key erms, framing he problem eliabiliy parameers: Failure ae, Failure Probabiliy, Availabiliy, ec. Some imporan reliabiliy disribuions Componen reliabiliy
More informationObject tracking: Using HMMs to estimate the geographical location of fish
Objec racking: Using HMMs o esimae he geographical locaion of fish 02433 - Hidden Markov Models Marin Wæver Pedersen, Henrik Madsen Course week 13 MWP, compiled June 8, 2011 Objecive: Locae fish from agging
More information20. Applications of the Genetic-Drift Model
0. Applicaions of he Geneic-Drif Model 1) Deermining he probabiliy of forming any paricular combinaion of genoypes in he nex generaion: Example: If he parenal allele frequencies are p 0 = 0.35 and q 0
More informationShiva Akhtarian MSc Student, Department of Computer Engineering and Information Technology, Payame Noor University, Iran
Curren Trends in Technology and Science ISSN : 79-055 8hSASTech 04 Symposium on Advances in Science & Technology-Commission-IV Mashhad, Iran A New for Sofware Reliabiliy Evaluaion Based on NHPP wih Imperfec
More informationA Nonparametric Multivariate Control Chart Based on. Data Depth. Abstract
A Nonparameric Mulivariae Conrol Char Based on Daa Deph Amor Messaoud, Claus Weihs and Franz Hering Deparmen of Saisics, Universiy of Dormund, 44221 Dormund, Germany email: messaoud@saisik.uni-dormund.de
More informationDEPARTMENT OF STATISTICS
A Tes for Mulivariae ARCH Effecs R. Sco Hacker and Abdulnasser Haemi-J 004: DEPARTMENT OF STATISTICS S-0 07 LUND SWEDEN A Tes for Mulivariae ARCH Effecs R. Sco Hacker Jönköping Inernaional Business School
More informationACE 562 Fall Lecture 5: The Simple Linear Regression Model: Sampling Properties of the Least Squares Estimators. by Professor Scott H.
ACE 56 Fall 005 Lecure 5: he Simple Linear Regression Model: Sampling Properies of he Leas Squares Esimaors by Professor Sco H. Irwin Required Reading: Griffihs, Hill and Judge. "Inference in he Simple
More informationEstimation Uncertainty
Esimaion Uncerainy The sample mean is an esimae of β = E(y +h ) The esimaion error is = + = T h y T b ( ) = = + = + = = = T T h T h e T y T y T b β β β Esimaion Variance Under classical condiions, where
More informationPade and Laguerre Approximations Applied. to the Active Queue Management Model. of Internet Protocol
Applied Mahemaical Sciences, Vol. 7, 013, no. 16, 663-673 HIKARI Ld, www.m-hikari.com hp://dx.doi.org/10.1988/ams.013.39499 Pade and Laguerre Approximaions Applied o he Acive Queue Managemen Model of Inerne
More informationEFFECT OF AUTOCORRELATION ON SPC CHART PERFORMANCE
3 rd Research/Exper Conference wih Inernaional Paricipaions QUALIT 3, Zenica, B&H, 13 and 14 ovember, 3. EFFECT OF AUTOCORRELATIO O SPC CHART PERFORMACE Vesna Bucevsa, Ph.D. Universiy S. Cyril and Mehodius,
More informationAir Traffic Forecast Empirical Research Based on the MCMC Method
Compuer and Informaion Science; Vol. 5, No. 5; 0 ISSN 93-8989 E-ISSN 93-8997 Published by Canadian Cener of Science and Educaion Air Traffic Forecas Empirical Research Based on he MCMC Mehod Jian-bo Wang,
More informationOBJECTIVES OF TIME SERIES ANALYSIS
OBJECTIVES OF TIME SERIES ANALYSIS Undersanding he dynamic or imedependen srucure of he observaions of a single series (univariae analysis) Forecasing of fuure observaions Asceraining he leading, lagging
More informationStructural Break Detection in Time Series Models
Srucural Break Deecion in Time Series Models Richard A. Davis Thomas Lee Gabriel Rodriguez-Yam Colorado Sae Universiy (hp://www.sa.colosae.edu/~rdavis/lecures) This research suppored in par by: Naional
More informationA Generalized Poisson-Akash Distribution: Properties and Applications
Inernaional Journal of Saisics and Applicaions 08, 8(5): 49-58 DOI: 059/jsaisics0808050 A Generalized Poisson-Akash Disribuion: Properies and Applicaions Rama Shanker,*, Kamlesh Kumar Shukla, Tekie Asehun
More informationComparing Means: t-tests for One Sample & Two Related Samples
Comparing Means: -Tess for One Sample & Two Relaed Samples Using he z-tes: Assumpions -Tess for One Sample & Two Relaed Samples The z-es (of a sample mean agains a populaion mean) is based on he assumpion
More informationLecture 3: Exponential Smoothing
NATCOR: Forecasing & Predicive Analyics Lecure 3: Exponenial Smoohing John Boylan Lancaser Cenre for Forecasing Deparmen of Managemen Science Mehods and Models Forecasing Mehod A (numerical) procedure
More informationEcon Autocorrelation. Sanjaya DeSilva
Econ 39 - Auocorrelaion Sanjaya DeSilva Ocober 3, 008 1 Definiion Auocorrelaion (or serial correlaion) occurs when he error erm of one observaion is correlaed wih he error erm of any oher observaion. This
More informationComputer Simulates the Effect of Internal Restriction on Residuals in Linear Regression Model with First-order Autoregressive Procedures
MPRA Munich Personal RePEc Archive Compuer Simulaes he Effec of Inernal Resricion on Residuals in Linear Regression Model wih Firs-order Auoregressive Procedures Mei-Yu Lee Deparmen of Applied Finance,
More informationStructural Break Detection in Time Series Models
Srucural Break Deecion in Time Series Models Richard A. Davis Thomas Lee Gabriel Rodriguez-Yam Colorado Sae Universiy (hp://www.sa.colosae.edu/~rdavis/lecures) This research suppored in par by an IBM faculy
More informationAppendix to Creating Work Breaks From Available Idleness
Appendix o Creaing Work Breaks From Available Idleness Xu Sun and Ward Whi Deparmen of Indusrial Engineering and Operaions Research, Columbia Universiy, New York, NY, 127; {xs2235,ww24}@columbia.edu Sepember
More informationEcon107 Applied Econometrics Topic 7: Multicollinearity (Studenmund, Chapter 8)
I. Definiions and Problems A. Perfec Mulicollineariy Econ7 Applied Economerics Topic 7: Mulicollineariy (Sudenmund, Chaper 8) Definiion: Perfec mulicollineariy exiss in a following K-variable regression
More informationStructural Break Detection in Time Series Models
Srucural Break Deecion in Time Series Models Richard A. Davis Thomas Lee Gabriel Rodriguez-Yam Colorado Sae Universiy (hp://www.sa.colosae.edu/~rdavis/lecures) This research suppored in par by an IBM faculy
More informationDynamic Probability Control Limits for Risk-Adjusted Bernoulli Cumulative Sum Charts
Dynamic Probabiliy Conrol Limis for Risk-Adjused Bernoulli Cumulaive Sum Chars Xiang Zhang Disseraion submied o he faculy of he Virginia Polyechnic Insiue and Sae Universiy in parial fulfillmen of he requiremens
More informationStructural Break Detection in Time Series Models
Srucural Break Deecion in Time Series Models Richard A. Davis Thomas Lee Gabriel Rodriguez-Yam Colorado Sae Universiy (hp://www.sa.colosae.edu/~rdavis/lecures) This research suppored in par by an IBM faculy
More informationOn Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature
On Measuring Pro-Poor Growh 1. On Various Ways of Measuring Pro-Poor Growh: A Shor eview of he Lieraure During he pas en years or so here have been various suggesions concerning he way one should check
More informationA new flexible Weibull distribution
Communicaions for Saisical Applicaions and Mehods 2016, Vol. 23, No. 5, 399 409 hp://dx.doi.org/10.5351/csam.2016.23.5.399 Prin ISSN 2287-7843 / Online ISSN 2383-4757 A new flexible Weibull disribuion
More informationChapter 11. Heteroskedasticity The Nature of Heteroskedasticity. In Chapter 3 we introduced the linear model (11.1.1)
Chaper 11 Heeroskedasiciy 11.1 The Naure of Heeroskedasiciy In Chaper 3 we inroduced he linear model y = β+β x (11.1.1) 1 o explain household expendiure on food (y) as a funcion of household income (x).
More informationA note on spurious regressions between stationary series
A noe on spurious regressions beween saionary series Auhor Su, Jen-Je Published 008 Journal Tile Applied Economics Leers DOI hps://doi.org/10.1080/13504850601018106 Copyrigh Saemen 008 Rouledge. This is
More information1. Introduction. Rawid Banchuin
011 Inernaional Conerence on Inormaion and Elecronics Engineering IPCSIT vol.6 (011 (011 IACSIT Press, Singapore Process Induced Random Variaion Models o Nanoscale MOS Perormance: Eicien ool or he nanoscale
More informationChristos Papadimitriou & Luca Trevisan November 22, 2016
U.C. Bereley CS170: Algorihms Handou LN-11-22 Chrisos Papadimiriou & Luca Trevisan November 22, 2016 Sreaming algorihms In his lecure and he nex one we sudy memory-efficien algorihms ha process a sream
More informationRegression with Time Series Data
Regression wih Time Series Daa y = β 0 + β 1 x 1 +...+ β k x k + u Serial Correlaion and Heeroskedasiciy Time Series - Serial Correlaion and Heeroskedasiciy 1 Serially Correlaed Errors: Consequences Wih
More informationInventory Control of Perishable Items in a Two-Echelon Supply Chain
Journal of Indusrial Engineering, Universiy of ehran, Special Issue,, PP. 69-77 69 Invenory Conrol of Perishable Iems in a wo-echelon Supply Chain Fariborz Jolai *, Elmira Gheisariha and Farnaz Nojavan
More informationGeorey E. Hinton. University oftoronto. Technical Report CRG-TR February 22, Abstract
Parameer Esimaion for Linear Dynamical Sysems Zoubin Ghahramani Georey E. Hinon Deparmen of Compuer Science Universiy oftorono 6 King's College Road Torono, Canada M5S A4 Email: zoubin@cs.orono.edu Technical
More informationProcess Monitoring versus Process Adjustment
CHAPTER 1 Process Monioring versus Process Adjusmen I can be said ha he birh of modern saisical process conrol Ž SPC. ook place when Waler A. Shewhar, a physicis and saisician working for Bell Laboraories,
More informationFailure of the work-hamiltonian connection for free energy calculations. Abstract
Failure of he work-hamilonian connecion for free energy calculaions Jose M. G. Vilar 1 and J. Miguel Rubi 1 Compuaional Biology Program, Memorial Sloan-Keering Cancer Cener, 175 York Avenue, New York,
More informationA unit root test based on smooth transitions and nonlinear adjustment
MPRA Munich Personal RePEc Archive A uni roo es based on smooh ransiions and nonlinear adjusmen Aycan Hepsag Isanbul Universiy 5 Ocober 2017 Online a hps://mpra.ub.uni-muenchen.de/81788/ MPRA Paper No.
More informationA Specification Test for Linear Dynamic Stochastic General Equilibrium Models
Journal of Saisical and Economeric Mehods, vol.1, no.2, 2012, 65-70 ISSN: 2241-0384 (prin), 2241-0376 (online) Scienpress Ld, 2012 A Specificaion Tes for Linear Dynamic Sochasic General Equilibrium Models
More informationNavneet Saini, Mayank Goyal, Vishal Bansal (2013); Term Project AML310; Indian Institute of Technology Delhi
Creep in Viscoelasic Subsances Numerical mehods o calculae he coefficiens of he Prony equaion using creep es daa and Herediary Inegrals Mehod Navnee Saini, Mayank Goyal, Vishal Bansal (23); Term Projec
More informationStructural Break Detection in Time Series
Srucural Break Deecion in Time Series Richard A. Davis Thomas ee Gabriel Rodriguez-Yam Colorado Sae Universiy (hp://www.sa.colosae.edu/~rdavis/lecures This research suppored in par by an IBM faculy award.
More informationTime series Decomposition method
Time series Decomposiion mehod A ime series is described using a mulifacor model such as = f (rend, cyclical, seasonal, error) = f (T, C, S, e) Long- Iner-mediaed Seasonal Irregular erm erm effec, effec,
More informationSPC Procedures for Monitoring Autocorrelated Processes
Qualiy Technology & Quaniaive Managemen Vol. 4, No. 4, pp. 501-540, 007 QTQM ICAQM 007 SPC Procedures for Monioring Auocorrelaed Processes S. Psarakis and G. E. A. Papaleonida Ahens Universiy of Economics
More informationBOOTSTRAP PREDICTION INTERVALS FOR TIME SERIES MODELS WITH HETROSCEDASTIC ERRORS. Department of Statistics, Islamia College, Peshawar, KP, Pakistan 2
Pak. J. Sais. 017 Vol. 33(1), 1-13 BOOTSTRAP PREDICTIO ITERVAS FOR TIME SERIES MODES WITH HETROSCEDASTIC ERRORS Amjad Ali 1, Sajjad Ahmad Khan, Alamgir 3 Umair Khalil and Dos Muhammad Khan 1 Deparmen of
More informationSTATE-SPACE MODELLING. A mass balance across the tank gives:
B. Lennox and N.F. Thornhill, 9, Sae Space Modelling, IChemE Process Managemen and Conrol Subjec Group Newsleer STE-SPACE MODELLING Inroducion: Over he pas decade or so here has been an ever increasing
More information5 The fitting methods used in the normalization of DSD
The fiing mehods used in he normalizaion of DSD.1 Inroducion Sempere-Torres e al. 1994 presened a general formulaion for he DSD ha was able o reproduce and inerpre all previous sudies of DSD. The mehodology
More informationACE 564 Spring Lecture 7. Extensions of The Multiple Regression Model: Dummy Independent Variables. by Professor Scott H.
ACE 564 Spring 2006 Lecure 7 Exensions of The Muliple Regression Model: Dumm Independen Variables b Professor Sco H. Irwin Readings: Griffihs, Hill and Judge. "Dumm Variables and Varing Coefficien Models
More information13.3 Term structure models
13.3 Term srucure models 13.3.1 Expecaions hypohesis model - Simples "model" a) shor rae b) expecaions o ge oher prices Resul: y () = 1 h +1 δ = φ( δ)+ε +1 f () = E (y +1) (1) =δ + φ( δ) f (3) = E (y +)
More informationEvaluation of Mean Time to System Failure of a Repairable 3-out-of-4 System with Online Preventive Maintenance
American Journal of Applied Mahemaics and Saisics, 0, Vol., No., 9- Available online a hp://pubs.sciepub.com/ajams/// Science and Educaion Publishing DOI:0.69/ajams--- Evaluaion of Mean Time o Sysem Failure
More informationDistribution of Least Squares
Disribuion of Leas Squares In classic regression, if he errors are iid normal, and independen of he regressors, hen he leas squares esimaes have an exac normal disribuion, no jus asympoic his is no rue
More informationMonitoring and data filtering II. Dynamic Linear Models
Ouline Monioring and daa filering II. Dynamic Linear Models (Wes and Harrison, chaper 2 Updaing equaions: Kalman Filer Discoun facor as an aid o choose W Incorporae exernal informaion: Inervenion General
More informationThe General Linear Test in the Ridge Regression
ommunicaions for Saisical Applicaions Mehods 2014, Vol. 21, No. 4, 297 307 DOI: hp://dx.doi.org/10.5351/sam.2014.21.4.297 Prin ISSN 2287-7843 / Online ISSN 2383-4757 The General Linear Tes in he Ridge
More informationState-Space Models. Initialization, Estimation and Smoothing of the Kalman Filter
Sae-Space Models Iniializaion, Esimaion and Smoohing of he Kalman Filer Iniializaion of he Kalman Filer The Kalman filer shows how o updae pas predicors and he corresponding predicion error variances when
More informationACE 562 Fall Lecture 8: The Simple Linear Regression Model: R 2, Reporting the Results and Prediction. by Professor Scott H.
ACE 56 Fall 5 Lecure 8: The Simple Linear Regression Model: R, Reporing he Resuls and Predicion by Professor Sco H. Irwin Required Readings: Griffihs, Hill and Judge. "Explaining Variaion in he Dependen
More informationPENALIZED LEAST SQUARES AND PENALIZED LIKELIHOOD
PENALIZED LEAST SQUARES AND PENALIZED LIKELIHOOD HAN XIAO 1. Penalized Leas Squares Lasso solves he following opimizaion problem, ˆβ lasso = arg max β R p+1 1 N y i β 0 N x ij β j β j (1.1) for some 0.
More informationTime Series Models for Growth of Urban Population in SAARC Countries
Advances in Managemen & Applied Economics, vol., no.1, 01, 109-119 ISSN: 179-7544 (prin version), 179-755 (online) Inernaional Scienific Press, 01 Time Series Models for Growh of Urban Populaion in SAARC
More informationRandom Walk with Anti-Correlated Steps
Random Walk wih Ani-Correlaed Seps John Noga Dirk Wagner 2 Absrac We conjecure he expeced value of random walks wih ani-correlaed seps o be exacly. We suppor his conjecure wih 2 plausibiliy argumens and
More informationDistribution of Estimates
Disribuion of Esimaes From Economerics (40) Linear Regression Model Assume (y,x ) is iid and E(x e )0 Esimaion Consisency y α + βx + he esimaes approach he rue values as he sample size increases Esimaion
More informationStochastic Model for Cancer Cell Growth through Single Forward Mutation
Journal of Modern Applied Saisical Mehods Volume 16 Issue 1 Aricle 31 5-1-2017 Sochasic Model for Cancer Cell Growh hrough Single Forward Muaion Jayabharahiraj Jayabalan Pondicherry Universiy, jayabharahi8@gmail.com
More informationI. Return Calculations (20 pts, 4 points each)
Universiy of Washingon Spring 015 Deparmen of Economics Eric Zivo Econ 44 Miderm Exam Soluions This is a closed book and closed noe exam. However, you are allowed one page of noes (8.5 by 11 or A4 double-sided)
More informationStatistical procedures searching for suspicious measurements and their sequences in measurement series. Ing. Marek Brabec, PhD CHMU, SZU, UI Praha
Saisical procedures searching for suspicious measuremens and heir sequences in measuremen series Ing. Marek Brabec, PhD CHMU, SZU, UI Praha 8h AQ EIONET meeing Oslo, November 2003 Goals: Propose procedures
More informationChapter 4. Location-Scale-Based Parametric Distributions. William Q. Meeker and Luis A. Escobar Iowa State University and Louisiana State University
Chaper 4 Locaion-Scale-Based Parameric Disribuions William Q. Meeker and Luis A. Escobar Iowa Sae Universiy and Louisiana Sae Universiy Copyrigh 1998-2008 W. Q. Meeker and L. A. Escobar. Based on he auhors
More informationArticle Sequential Change-Point Detection via Online Convex Optimization. and Huan Xu
enropy Aricle Sequenial Change-Poin Deecion via Online Convex Opimizaion Yang Cao, Liyan Xie, Yao Xie * ID and Huan Xu H. Milon Sewar School of Indusrial and Sysems Engineering, Georgia Insiue of Technology,
More informationGMM - Generalized Method of Moments
GMM - Generalized Mehod of Momens Conens GMM esimaion, shor inroducion 2 GMM inuiion: Maching momens 2 3 General overview of GMM esimaion. 3 3. Weighing marix...........................................
More informationDepartment of Economics East Carolina University Greenville, NC Phone: Fax:
March 3, 999 Time Series Evidence on Wheher Adjusmen o Long-Run Equilibrium is Asymmeric Philip Rohman Eas Carolina Universiy Absrac The Enders and Granger (998) uni-roo es agains saionary alernaives wih
More informationSolutions to Odd Number Exercises in Chapter 6
1 Soluions o Odd Number Exercises in 6.1 R y eˆ 1.7151 y 6.3 From eˆ ( T K) ˆ R 1 1 SST SST SST (1 R ) 55.36(1.7911) we have, ˆ 6.414 T K ( ) 6.5 y ye ye y e 1 1 Consider he erms e and xe b b x e y e b
More informationA New Unit Root Test against Asymmetric ESTAR Nonlinearity with Smooth Breaks
Iran. Econ. Rev. Vol., No., 08. pp. 5-6 A New Uni Roo es agains Asymmeric ESAR Nonlineariy wih Smooh Breaks Omid Ranjbar*, sangyao Chang, Zahra (Mila) Elmi 3, Chien-Chiang Lee 4 Received: December 7, 06
More informationSemi-Competing Risks on A Trivariate Weibull Survival Model
Semi-Compeing Risks on A Trivariae Weibull Survival Model Jenq-Daw Lee Graduae Insiue of Poliical Economy Naional Cheng Kung Universiy Tainan Taiwan 70101 ROC Cheng K. Lee Loss Forecasing Home Loans &
More informationA Bayesian Approach to Spectral Analysis
Chirped Signals A Bayesian Approach o Specral Analysis Chirped signals are oscillaing signals wih ime variable frequencies, usually wih a linear variaion of frequency wih ime. E.g. f() = A cos(ω + α 2
More information4.1 Other Interpretations of Ridge Regression
CHAPTER 4 FURTHER RIDGE THEORY 4. Oher Inerpreaions of Ridge Regression In his secion we will presen hree inerpreaions for he use of ridge regression. The firs one is analogous o Hoerl and Kennard reasoning
More informationModule 2 F c i k c s la l w a s o s f dif di fusi s o i n
Module Fick s laws of diffusion Fick s laws of diffusion and hin film soluion Adolf Fick (1855) proposed: d J α d d d J (mole/m s) flu (m /s) diffusion coefficien and (mole/m 3 ) concenraion of ions, aoms
More informationMath 10B: Mock Mid II. April 13, 2016
Name: Soluions Mah 10B: Mock Mid II April 13, 016 1. ( poins) Sae, wih jusificaion, wheher he following saemens are rue or false. (a) If a 3 3 marix A saisfies A 3 A = 0, hen i canno be inverible. True.
More informationPerformance of Stochastically Intermittent Sensors in Detecting a Target Traveling between Two Areas
American Journal of Operaions Research, 06, 6, 99- Published Online March 06 in SciRes. hp://www.scirp.org/journal/ajor hp://dx.doi.org/0.436/ajor.06.60 Performance of Sochasically Inermien Sensors in
More informationStationary Time Series
3-Jul-3 Time Series Analysis Assoc. Prof. Dr. Sevap Kesel July 03 Saionary Time Series Sricly saionary process: If he oin dis. of is he same as he oin dis. of ( X,... X n) ( X h,... X nh) Weakly Saionary
More informationA Point Optimal Test for the Null of Near Integration. A. Aznar and M. I. Ayuda 1. University of Zaragoza
A Poin Opimal es for he Null of Near Inegraion A. Aznar and M. I. Ayuda Universiy of Zaragoza he objecive of his paper is o derive a poin opimal es for he null hypohesis of near inegraion (PONI-es). We
More informationWATER LEVEL TRACKING WITH CONDENSATION ALGORITHM
WATER LEVEL TRACKING WITH CONDENSATION ALGORITHM Shinsuke KOBAYASHI, Shogo MURAMATSU, Hisakazu KIKUCHI, Masahiro IWAHASHI Dep. of Elecrical and Elecronic Eng., Niigaa Universiy, 8050 2-no-cho Igarashi,
More informationProbabilistic Robotics
Probabilisic Roboics Bayes Filer Implemenaions Gaussian filers Bayes Filer Reminder Predicion bel p u bel d Correcion bel η p z bel Gaussians : ~ π e p N p - Univariae / / : ~ μ μ μ e p Ν p d π Mulivariae
More informationImproved Approximate Solutions for Nonlinear Evolutions Equations in Mathematical Physics Using the Reduced Differential Transform Method
Journal of Applied Mahemaics & Bioinformaics, vol., no., 01, 1-14 ISSN: 179-660 (prin), 179-699 (online) Scienpress Ld, 01 Improved Approimae Soluions for Nonlinear Evoluions Equaions in Mahemaical Physics
More informationNature Neuroscience: doi: /nn Supplementary Figure 1. Spike-count autocorrelations in time.
Supplemenary Figure 1 Spike-coun auocorrelaions in ime. Normalized auocorrelaion marices are shown for each area in a daase. The marix shows he mean correlaion of he spike coun in each ime bin wih he spike
More informationMATHEMATICAL DESCRIPTION OF THEORETICAL METHODS OF RESERVE ECONOMY OF CONSIGNMENT STORES
MAHEMAICAL DESCIPION OF HEOEICAL MEHODS OF ESEVE ECONOMY OF CONSIGNMEN SOES Péer elek, József Cselényi, György Demeer Universiy of Miskolc, Deparmen of Maerials Handling and Logisics Absrac: Opimizaion
More informationChapter 5. Heterocedastic Models. Introduction to time series (2008) 1
Chaper 5 Heerocedasic Models Inroducion o ime series (2008) 1 Chaper 5. Conens. 5.1. The ARCH model. 5.2. The GARCH model. 5.3. The exponenial GARCH model. 5.4. The CHARMA model. 5.5. Random coefficien
More informationThe Arcsine Distribution
The Arcsine Disribuion Chris H. Rycrof Ocober 6, 006 A common heme of he class has been ha he saisics of single walker are ofen very differen from hose of an ensemble of walkers. On he firs homework, we
More informationThe Simple Linear Regression Model: Reporting the Results and Choosing the Functional Form
Chaper 6 The Simple Linear Regression Model: Reporing he Resuls and Choosing he Funcional Form To complee he analysis of he simple linear regression model, in his chaper we will consider how o measure
More information