A Bootstrap Confidence Limit for Process Capability Indices
|
|
- Jesse Pearson
- 5 years ago
- Views:
Transcription
1 A ootstrap Confdence Lmt for Process Capablty Indces YANG Janfeng School of usness, Zhengzhou Unversty, P.R.Chna, Abstract The process capablty ndces are wdely used by qualty professonals as an estmate of process capablty. The lower confdence lmts of PCIs are dffcult to be estmated by parametrc methods for some non-normal dstrbuted processes. The non-parametrc but computer ntensve ootstrap technques are utlzed for these cases. The Percentle-t ootstrap (PT) method s used to estmate the confdence nterval for PCIs of these non-normal dstrbuted processes n the paper. A seres of smulatons usng dfferent dstrbutons, parameters and sample sze for the ndex C pk was undertaken to compare the performance and relablty of PT to S, P and CP. The results show that the PT method s better to estmate lower confdence lmts of PCIs for most non-normal dstrbuted processes. Keywords Process capablty ndex, ootstrap, Confdence lmt, non-normal dstrbuton 1 Introducton A seres of process capablty ndces (PCI) have receved substantal attenton snce Juran ntroduced the frst ndex C p. PCIs was wdely used n manufacture ndustry all over the world. PCIs are used to determne the capablty of a process. The most basc modfcaton s C p, C pk (Kane, 1986). The most basc PCIs can be defned as the followng: USL LSL C p = (1) 6σ USL µ µ LSL C p = mn, (2) 6σ 6σ where USL s the upper specfcaton lmt, LSL s the lower specfcaton lmt, µ s the process mean, σ s the process standard devaton, and T s target value. The mportant part of PCI s the confdence nterval because PCI s an estmated value. Chou et al. (1990) provded tables for constructng the 95% lower confdence lmts for C p, C pk. oyles (1991) provded an approxmate method for fndng lower confdence nterval. ut the calculaton of all these confdence nterval assumes a normal dstrbuted process. There s not a good parametrc method to confrm the confdence nterval of process capablty ndex for the non-normal dstrbuted processes. ootstrappng s a non-parametrc, but computer-ntensve estmaton method, ntroduced by Efron. ootstrap confdence ntervals may be obtaned wthout any assumpton of underlyng dstrbuton. Frankln and Wasserman (1992) ntroduced the ootstrap method to calculate confdence nterval of PCI and constructed the confdence nterval for some basc PCIs. Cho, et al. (1996), Han and Chou (2000), Tong (1998), Nam et al. (2001) and alamural (2002, 2003) dscussed ootstrap confdence nterval for non-normal dstrbutons and non-normal PCIs. ut the researchs show that some standard ootstrap methods are not precse for non-normal dstrbuted processes. The artcle reported the performance and relablty of Percentle-t ootstrap (PT) method (Lm, Shn and Km, 2004) for estmatng confdence nterval of PCIs n non-normal dstrbuted processes. 2 ootstrap Confdence Intervals of PCI ootstrap method s a computer-based method for estmatng the standard error of a summary statstc. The ootstrap estmate of standard error requres no theoretcal calculatons, and s avalable no matter how mathematcally complcated the summary statstc may be. Several types of ootstrap method for constructng confdence ntervals have been developed: the standard method (S), the percentle method (P), the bas-corrected percentle method (CP), and percentle-t method. 1201
2 Let x 1, x 2,, x n be a sample X taken from a process,.e., a sequence of n..d. random varables. A ootstrap sample X, denoted by x 1, x 2,, x n, s a sample of sze n drawn wth replacement from the orgnal sample X. There are a total of n n such resamples possble. In our case, these resamples would be used to calculate n n values of Ĉ. (C denotes PCIs.) Each of these would be an estmate of C and the entre collecton would consttute the complete ootstrap dstrbuton for Ĉ. ootstrap samplng s equvalent to samplng wth replacement form the emprcal probablty dstrbuton functon. Thus, the ootstrap dstrbutons of Ĉ s estmate of the dstrbuton of C. In practce, a rough mnmum of 1000 ootstrap resamples are usually suffcent to compute reasonably accurate confdence nterval estmates (Efron & Tbshran, 1986). Therefore, t s assumed that =1000 tmes bootstrap samples are taken and =1000 tmes ootstrap Ĉ are calculated and ordered from smallest to largest. The followeng s four possble constructons for the confdence nterval of PCIs usng ootstrap technques. 2.1 Standard ootstrap (S) Confdence Interval C ˆ From the =1000 ootstrap estmates, for =1,2,,, calculate the sample average 1 C = C ˆ and the sample standard devaton 1 2 s = ( Cˆ C ).The standard ootstrap c 1 = 1 confdence nterval for C s ˆ [, ˆ C Zα sc C+ Zαs c ] (3) 2.2 Percentle ootstrap (P) Confdence Interval From the ordered collecton of C, for =1, 2,, n, the percentle ootstrap confdence nterval for ˆ C s ˆ ˆ [ C ( α), C ((1 α ) )] (4) 2.3 ased-corrected Percentle ootstrap (CP) Confdence Interval It s possble that ootstrap dstrbutons obtaned usng only a sample of the complete ootstrap dstrbuton may be shfted hgher or lower than would be expected. Thus, a thrd method has been developed to correct for ths potental bas. Frst, usng the ordered dstrbuton of Ĉ, calculate the probablty ˆ ( ˆ 1 P0 = Pr C C). Then, Z 0 = Φ ( P 0 ). Thus, we can get the percentle of ordered dstrbuton F ( C ˆ ), P L = Φ( 2Z 0 Zα ) and P U = Φ( 2Z 0 + Zα ), where Φ ( ) s the standard normal cumulatve functon. Fnally, the CP confdence nterval s gven by ˆ ( ), ˆ [ C PL C ( PU )] (5) 2.4 Percentle-t ootstrap (PT) confdence nterval Cˆ ( Cˆ) ased on the ootstrap resample X, we compute Z =, = 1,2,...,, where Cˆ s the se estmate of C based on the ootstrap resample {X 1, X 2,,X }, and se s the estmated standard error of C. Thus, we can get a collecton Z, = 1,2,...,. Let Z ( ) denote the th smallest value among Z, = 1,2,...,. Then, the αth quantle of these Z values s estmated by the value t ˆα ( ) such that =1 I[ Z ( ) tˆ( α ) =α. Note that f α s not an nteger, the αth quantle s equal to the [(+1)α]th smallest value and the (1-α)th quantle s equal to the [(+1)-(+1)α]th smallest value among Z, = 1,2,...,. The 100(1-2α)% ootstrap confdence nterval for C s constructed as Cˆ + tˆ( α) s, Cˆ + tˆ(1 α) s ] (6) [ d d 1202
3 where s d denotes the standard error of Ĉ. However, there are no smple formulas to compute the standard error of Ĉ because of ts complcated structure. That s, t s not qute possble to compute ether estmate s d and c se for gven ootstrap resamples {X 1,X 2,,X }. Thus, we may use ootstrap s d. To estmate se, we frst compute a bootstrap estmate of standard error s nstead of for each ootstrap resample by utlzng two nested levels of ootstrap resamplng. To carry out such a resamplng, we generate 25 ootstrap samples {X 1,X 2,,X } for the purpose of estmatng the standard error. For each 25 ootstrap resample, we compute C n = 2 ˆ C j j= 1 Therefore, a 100(1-2α)% ootstrap confdence nterval for C s obtaned as ˆ ˆ [ C+ tˆ( α) s, C+ tˆ(1 α) ] c s c 3 The Smulaton for ootstrap Lower Confdence Lmts of PCI To compare the performance of ootstrap lower confdence lmts to those based on the assumpton of non-normal processes, a seres of smulatons was undertaken. In the smulaton, some dfferent dstrbutons and parameters were used to nvestgate four dfferent processes based on the assumpton that datum come from normal, skew and heavy-tal dstrbuted processes. For process capablty ndex, the lower confdence lmt s mportant not upper confdence lmt. So, we dscuss the ootstrap lower confdence lmt n the artcle. We choose normal, lognormal and t dstrbuton wth dfferent parameters to calculate the ootstrap lower confdence lmts of C pk n these smulatons. The value USL=3, LSL=-3 and T=0 was used for all smulatons. The sample sze n are 10, 30, 50 n each smulaton. We choose normal dstrbuton N(0,1) for normal processes, Lognormal(0,0.2) for slght skew processes, Lognormal(0,0.4) for severe skew processes and t 6 for heavy-taled processes. Then, there s C pk =1 for N(0,1). In order to make each dstrbuton have specfed value of µ=0 and σ=1, we transform three non-normal dstrbutons. Then, these transformed dstrbuton s gotten wth µ=0 and σ=1. Thus, these dfferent dstrbutons are comparable. Sngle smulaton was replcated N=1000 tmes. Thus, we were able to calculate the proporton of tmes the four ootstrap lower confdence lmts and the normally derved lower lmts were less than the correspondng ndex. The actual coverage proporton could be compared wth the expected value of 95%. These smulatons were undertaken by the MATLA software. The smulaton results are tabulated n Table 1-4. In these tables, Cov. denotes the coverage proporton. E ( C ) denotes the mean of 95%lower confdence lmts. E (Cˆ ) denotes the mean of estmates of C. ˆ C C ) denotes the mean of the wdth between Ĉ and C. A. Normal dstrbuted process The results show that, for a normal dstrbuted process wth n 30, all four ootstrap methods provde enough coverage proportons. ut S and PT provde wder confdence nterval than other two methods. In general, CP s better for n s larger. For small sample sze (), P and CP methods have lower coverage proportons. The two methods are not applcable to estmate ootstrap confdence nterval. For the value of ˆ C C ), PT method s better than S method for ths case. (Showed n Table1). Lognormal dstrbuted process For slght skew dstrbuted processes, S and PT methods provde stable and hgher coverage proportons, whether n s large or small. P and CP methods provde poor coverage proportons. Comparng the smulaton results of Lognormal(0,0.2) and Lognormal(0,0.4), the hger the skewness s, and se ˆ = n2 j= 1 ( Cˆ j n C 2 ) 2. (7) 1203
4 and the lower the coverage proporton s. For the dstrbuton Lognormal(0,0.4) wth hgh skewness, the coverage proporton of PT method s more than Therefore, PT lower confdence lmts are relable for those underlyng skew dstrbutons. (Showed n Table 2 and 3) Table 1 Converage Proporton for 95% Lower Table 2 Converage Proporton for 95% Lower Confdence Lmts of C pk Normal Process Confdence Lmts of C pk Lognormal(0,.2) Process n Method Cov. ) ˆ C C ) n Method Cov. ) ˆ C C ) S 0.981# P CP PT S P CP PT 0.975# S P CP PT 0.976# S (0.273) (0.279) (0.278) (0.359) (0.306) (0.336) P (0.242) (0.279) (0.069) (0.275) (0.306) (0.070) CP (0.227) (0.279) (0.082) (0.259) (0.306) (0.083) PT (0.293) (0.279) (0.264) (0.393) (0.306) (0.340) S (0.120) (0.138) (0.056) (0.153) (0.147) (0.073) P (0.123) (0.138) (0.034) (0.142) (0.147) (0.042) CP (0.122) (0.138) (0.041) (0.145) (0.147) (0.048) # PT (0.117) (0.138) (0.059) (0.146) (0.147) (0.069) S (0.095) (0.109) (0.031) (0.119) (0.119) (0.036) P (0.097) (0.109) (0.025) (0.117) (0.119) (0.027) CP (0.099) (0.109) (0.028) (0.120) (0.119) (0.031) # PT (0.094) (0.109) (0.034) (0.116) (0.119) (0.037) C. t dstrbuted process For the heavy-taled dstrbutons, P and CP methods provde poor coverage proportons, less than PT and S methods have approxmate coverage proportons for all sample sze. ut the PT method has smaller standard devaton of confdence lmts. Therefore, PT lower confdence lmts are relable for the heavy-taled processes. (Showed n Table 4) 4 Concluson In general, S and PT methods can provde relable confdence lmts (ntervals) for underlyng dstrbuted processes whatever the sample sze. P and CP methods have poor coverage proportons for most cases. The PT method provdes hgher coverage proporton wth small confdence ntervals, whether the sample sze s large or small. In the same tme, the method has wld applcablty for most non-normal dstrbuted processes. Therefore, we suggest, for the underlyng processes wth sample sze beng more than 30, the PT confdence lmt be used frstly. ut the PT method has ts shortcomng that ts calculatng program s very complex. When the twce resample sze s more than 50, t take a lot of tme to calculate. 1204
5 Acknowledgment Ths artcle were supported n part by the Natonal Scence Fund Chna Grant ( ) and Henan Soft Scence Grant ( ). Table 3 Converage Proporton for 95% Lower Table 4 Converage Proporton for 95% Lower Confdence Lmts of C pk Lognormal(0,.4) Process Confdence Lmts of C pk t 6 Process n Method Cov. ) ˆ C C ) n Method Cov. ) ˆ C C ) S P CP PT S P CP PT S P CP PT S (0.485) (0.375) (0.412) (0.437) (0.363) (0.393) P (0.347) (0.375) (0.076) (0.326) (0.363) (0.078) CP (0.326) (0.375) (0.088) (0.307) (0.363) (0.090) PT (0.500) (0.375) (0.388) (0.468) (0.363) (0.392) S (0.232) (0.203) (0.088) (0.192) (0.181) (0.083) P (0.209) (0.203) (0.044) (0.173) (0.181) (0.042) CP (0.209) (0.203) (0.049) (0.175) (0.181) (0.049) PT (0.225) (0.203) (0.079) (0.182) (0.181) (0.077) S (0.181) (0.157) (0.063) (0.173) (0.152) (0.064) P (0.167) (0.157) (0.038) (0.154) (0.152) (0.035) CP (0.169) (0.157) (0.043) (0.157) (0.152) (0.042) PT (0.178) (0.157) (0.059) (0.165) (0.152) (0.058) References [1] V. E. Kane. Process capablty ndces. Journal of Qualty Technology, 1986, 18(1): [2] Y. Chou, D.. Owen. On the dstrbuton of the estmated process capablty ndces. Communcaton n Statstcs Theory and Methods, 1989, 18: [3] R. A. oyles. The Taguch capablty ndex. Journal of Qualty Technology, 1991, 23(1): [4] L. A. Frankln, G. Wasserman. ootstrap lower confdence lmts for capablty ndces. Journal of Qualty Technology, 1992, 24(4): [5] K.C. Cho, K. H. Nam, D. H. Park. Estmaton of capablty ndex based on ootstrap method. Mcroelectron Relablty, 1996, 36(9): [6] J. H. Han, J. J. Cho, C. S. Leem. ootstrap confdence lmts for Wrght s C s. Communcaton n Statstcs - Theory and Method, 2000, 29: [7] L. I. Tong, J. P. Chen. Lower confdence lmts of process capablty ndces for non-normal process dstrbutons. Internal Journal Qualty & Relablty Management, 1998, 15(8): [8] K. H. Nam, D. K. Km, D. H. Park. Large-sample nterval estmators for process capablty ndces. Qualty Engneerng, 2001, 14(2): [9] S. alamural, M. Kalyanasundaram. ootstrap lower confdence lmts for the process capablty ndces C p, C pk and C pm. Internatonal Journal of Qualty & Relablty Management, 2002, 19(8):
6 [10] S. alamural. ootstrap confdence lmts for short-run capablty ndces. Qualty Engneerng, 2003, 15(4): [11] J. H. Lm, S. W. Shn, D. K. Km, et al. ootstrap confdence ntervals for steady-state avalablty. Asa-Pacfc Journal of Operatonal Research, 2004, 21(3): The Author can be contacted from Emal: 1206
Capability Analysis. Chapter 255. Introduction. Capability Analysis
Chapter 55 Introducton Ths procedure summarzes the performance of a process based on user-specfed specfcaton lmts. The observed performance as well as the performance relatve to the Normal dstrbuton are
More informationEconomic Design of Short-Run CSP-1 Plan Under Linear Inspection Cost
Tamkang Journal of Scence and Engneerng, Vol. 9, No 1, pp. 19 23 (2006) 19 Economc Desgn of Short-Run CSP-1 Plan Under Lnear Inspecton Cost Chung-Ho Chen 1 * and Chao-Yu Chou 2 1 Department of Industral
More informationTests for Two Correlations
PASS Sample Sze Software Chapter 805 Tests for Two Correlatons Introducton The correlaton coeffcent (or correlaton), ρ, s a popular parameter for descrbng the strength of the assocaton between two varables.
More informationInterval Estimation for a Linear Function of. Variances of Nonnormal Distributions. that Utilize the Kurtosis
Appled Mathematcal Scences, Vol. 7, 013, no. 99, 4909-4918 HIKARI Ltd, www.m-hkar.com http://dx.do.org/10.1988/ams.013.37366 Interval Estmaton for a Lnear Functon of Varances of Nonnormal Dstrbutons that
More information4. Greek Letters, Value-at-Risk
4 Greek Letters, Value-at-Rsk 4 Value-at-Rsk (Hull s, Chapter 8) Math443 W08, HM Zhu Outlne (Hull, Chap 8) What s Value at Rsk (VaR)? Hstorcal smulatons Monte Carlo smulatons Model based approach Varance-covarance
More informationReal Exchange Rate Fluctuations, Wage Stickiness and Markup Adjustments
Real Exchange Rate Fluctuatons, Wage Stckness and Markup Adjustments Yothn Jnjarak and Kanda Nakno Nanyang Technologcal Unversty and Purdue Unversty January 2009 Abstract Motvated by emprcal evdence on
More informationRandom Variables. b 2.
Random Varables Generally the object of an nvestgators nterest s not necessarly the acton n the sample space but rather some functon of t. Techncally a real valued functon or mappng whose doman s the sample
More informationLinear Combinations of Random Variables and Sampling (100 points)
Economcs 30330: Statstcs for Economcs Problem Set 6 Unversty of Notre Dame Instructor: Julo Garín Sprng 2012 Lnear Combnatons of Random Varables and Samplng 100 ponts 1. Four-part problem. Go get some
More information3: Central Limit Theorem, Systematic Errors
3: Central Lmt Theorem, Systematc Errors 1 Errors 1.1 Central Lmt Theorem Ths theorem s of prme mportance when measurng physcal quanttes because usually the mperfectons n the measurements are due to several
More informationAn Approximate E-Bayesian Estimation of Step-stress Accelerated Life Testing with Exponential Distribution
Send Orders for Reprnts to reprnts@benthamscenceae The Open Cybernetcs & Systemcs Journal, 25, 9, 729-733 729 Open Access An Approxmate E-Bayesan Estmaton of Step-stress Accelerated Lfe Testng wth Exponental
More informationoccurrence of a larger storm than our culvert or bridge is barely capable of handling? (what is The main question is: What is the possibility of
Module 8: Probablty and Statstcal Methods n Water Resources Engneerng Bob Ptt Unversty of Alabama Tuscaloosa, AL Flow data are avalable from numerous USGS operated flow recordng statons. Data s usually
More informationMode is the value which occurs most frequency. The mode may not exist, and even if it does, it may not be unique.
1.7.4 Mode Mode s the value whch occurs most frequency. The mode may not exst, and even f t does, t may not be unque. For ungrouped data, we smply count the largest frequency of the gven value. If all
More informationMgtOp 215 Chapter 13 Dr. Ahn
MgtOp 5 Chapter 3 Dr Ahn Consder two random varables X and Y wth,,, In order to study the relatonshp between the two random varables, we need a numercal measure that descrbes the relatonshp The covarance
More informationTests for Two Ordered Categorical Variables
Chapter 253 Tests for Two Ordered Categorcal Varables Introducton Ths module computes power and sample sze for tests of ordered categorcal data such as Lkert scale data. Assumng proportonal odds, such
More informationInformation Flow and Recovering the. Estimating the Moments of. Normality of Asset Returns
Estmatng the Moments of Informaton Flow and Recoverng the Normalty of Asset Returns Ané and Geman (Journal of Fnance, 2000) Revsted Anthony Murphy, Nuffeld College, Oxford Marwan Izzeldn, Unversty of Lecester
More informationChapter 3 Student Lecture Notes 3-1
Chapter 3 Student Lecture otes 3-1 Busness Statstcs: A Decson-Makng Approach 6 th Edton Chapter 3 Descrbng Data Usng umercal Measures 005 Prentce-Hall, Inc. Chap 3-1 Chapter Goals After completng ths chapter,
More informationMidterm Exam. Use the end of month price data for the S&P 500 index in the table below to answer the following questions.
Unversty of Washngton Summer 2001 Department of Economcs Erc Zvot Economcs 483 Mdterm Exam Ths s a closed book and closed note exam. However, you are allowed one page of handwrtten notes. Answer all questons
More informationClearing Notice SIX x-clear Ltd
Clearng Notce SIX x-clear Ltd 1.0 Overvew Changes to margn and default fund model arrangements SIX x-clear ( x-clear ) s closely montorng the CCP envronment n Europe as well as the needs of ts Members.
More informationInternational ejournals
Avalable onlne at www.nternatonalejournals.com ISSN 0976 1411 Internatonal ejournals Internatonal ejournal of Mathematcs and Engneerng 7 (010) 86-95 MODELING AND PREDICTING URBAN MALE POPULATION OF BANGLADESH:
More informationLikelihood Fits. Craig Blocker Brandeis August 23, 2004
Lkelhood Fts Crag Blocker Brandes August 23, 2004 Outlne I. What s the queston? II. Lkelhood Bascs III. Mathematcal Propertes IV. Uncertantes on Parameters V. Mscellaneous VI. Goodness of Ft VII. Comparson
More information/ Computational Genomics. Normalization
0-80 /02-70 Computatonal Genomcs Normalzaton Gene Expresson Analyss Model Computatonal nformaton fuson Bologcal regulatory networks Pattern Recognton Data Analyss clusterng, classfcaton normalzaton, mss.
More informationChapter 5 Student Lecture Notes 5-1
Chapter 5 Student Lecture Notes 5-1 Basc Busness Statstcs (9 th Edton) Chapter 5 Some Important Dscrete Probablty Dstrbutons 004 Prentce-Hall, Inc. Chap 5-1 Chapter Topcs The Probablty Dstrbuton of a Dscrete
More informationUNIVERSITY OF VICTORIA Midterm June 6, 2018 Solutions
UIVERSITY OF VICTORIA Mdterm June 6, 08 Solutons Econ 45 Summer A0 08 age AME: STUDET UMBER: V00 Course ame & o. Descrptve Statstcs and robablty Economcs 45 Secton(s) A0 CR: 3067 Instructor: Betty Johnson
More informationPhysics 4A. Error Analysis or Experimental Uncertainty. Error
Physcs 4A Error Analyss or Expermental Uncertanty Slde Slde 2 Slde 3 Slde 4 Slde 5 Slde 6 Slde 7 Slde 8 Slde 9 Slde 0 Slde Slde 2 Slde 3 Slde 4 Slde 5 Slde 6 Slde 7 Slde 8 Slde 9 Slde 20 Slde 2 Error n
More informationII. Random Variables. Variable Types. Variables Map Outcomes to Numbers
II. Random Varables Random varables operate n much the same way as the outcomes or events n some arbtrary sample space the dstncton s that random varables are smply outcomes that are represented numercally.
More informationIntroduction. Chapter 7 - An Introduction to Portfolio Management
Introducton In the next three chapters, we wll examne dfferent aspects of captal market theory, ncludng: Brngng rsk and return nto the pcture of nvestment management Markowtz optmzaton Modelng rsk and
More informationRaising Food Prices and Welfare Change: A Simple Calibration. Xiaohua Yu
Rasng Food Prces and Welfare Change: A Smple Calbraton Xaohua Yu Professor of Agrcultural Economcs Courant Research Centre Poverty, Equty and Growth Unversty of Göttngen CRC-PEG, Wlhelm-weber-Str. 2 3773
More informationThe Integration of the Israel Labour Force Survey with the National Insurance File
The Integraton of the Israel Labour Force Survey wth the Natonal Insurance Fle Natale SHLOMO Central Bureau of Statstcs Kanfey Nesharm St. 66, corner of Bach Street, Jerusalem Natales@cbs.gov.l Abstact:
More informationGraphical Methods for Survival Distribution Fitting
Graphcal Methods for Survval Dstrbuton Fttng In ths Chapter we dscuss the followng two graphcal methods for survval dstrbuton fttng: 1. Probablty Plot, 2. Cox-Snell Resdual Method. Probablty Plot: The
More informationOCR Statistics 1 Working with data. Section 2: Measures of location
OCR Statstcs 1 Workng wth data Secton 2: Measures of locaton Notes and Examples These notes have sub-sectons on: The medan Estmatng the medan from grouped data The mean Estmatng the mean from grouped data
More informationGlobal sensitivity analysis of credit risk portfolios
Global senstvty analyss of credt rsk portfolos D. Baur, J. Carbon & F. Campolongo European Commsson, Jont Research Centre, Italy Abstract Ths paper proposes the use of global senstvty analyss to evaluate
More informationMeasures of Spread IQR and Deviation. For exam X, calculate the mean, median and mode. For exam Y, calculate the mean, median and mode.
Part 4 Measures of Spread IQR and Devaton In Part we learned how the three measures of center offer dfferent ways of provdng us wth a sngle representatve value for a data set. However, consder the followng
More informationA Utilitarian Approach of the Rawls s Difference Principle
1 A Utltaran Approach of the Rawls s Dfference Prncple Hyeok Yong Kwon a,1, Hang Keun Ryu b,2 a Department of Poltcal Scence, Korea Unversty, Seoul, Korea, 136-701 b Department of Economcs, Chung Ang Unversty,
More informationCreating a zero coupon curve by bootstrapping with cubic splines.
MMA 708 Analytcal Fnance II Creatng a zero coupon curve by bootstrappng wth cubc splnes. erg Gryshkevych Professor: Jan R. M. Röman 0.2.200 Dvson of Appled Mathematcs chool of Educaton, Culture and Communcaton
More informationASSESSING GOODNESS OF FIT OF GENERALIZED LINEAR MODELS TO SPARSE DATA USING HIGHER ORDER MOMENT CORRECTIONS
ASSESSING GOODNESS OF FIT OF GENERALIZED LINEAR MODELS TO SPARSE DATA USING HIGHER ORDER MOMENT CORRECTIONS S. R. PAUL Department of Mathematcs & Statstcs, Unversty of Wndsor, Wndsor, ON N9B 3P4, Canada
More informationTCOM501 Networking: Theory & Fundamentals Final Examination Professor Yannis A. Korilis April 26, 2002
TO5 Networng: Theory & undamentals nal xamnaton Professor Yanns. orls prl, Problem [ ponts]: onsder a rng networ wth nodes,,,. In ths networ, a customer that completes servce at node exts the networ wth
More informationAnalysis of Variance and Design of Experiments-II
Analyss of Varance and Desgn of Experments-II MODULE VI LECTURE - 4 SPLIT-PLOT AND STRIP-PLOT DESIGNS Dr. Shalabh Department of Mathematcs & Statstcs Indan Insttute of Technology Kanpur An example to motvate
More informationOPERATIONS RESEARCH. Game Theory
OPERATIONS RESEARCH Chapter 2 Game Theory Prof. Bbhas C. Gr Department of Mathematcs Jadavpur Unversty Kolkata, Inda Emal: bcgr.umath@gmal.com 1.0 Introducton Game theory was developed for decson makng
More informationDesigning of Skip-Lot Sampling Plan of Type (Sksp-3) For Life Tests Based On Percentiles of Exponentiated Rayleigh Distribution
Internatonal Journal of Research n Advent Technology, Vol.6, No.8, August 28 Desgnng of Skp-Lot Samplng Plan of Type (Sksp-3) For Lfe Tests Based On Percentles of Exponentated Raylegh Dstrbuton Pradeepa
More informationInt. Statistical Inst.: Proc. 58th World Statistical Congress, 2011, Dublin (Session STS041) p The Max-CUSUM Chart
Int. Statstcal Inst.: Proc. 58th World Statstcal Congress, 1, Dubln (Sesson STS41) p.2996 The Max-CUSUM Chart Smley W. Cheng Department of Statstcs Unversty of Mantoba Wnnpeg, Mantoba Canada, R3T 2N2 smley_cheng@umantoba.ca
More informationarxiv: v1 [q-fin.pm] 13 Feb 2018
WHAT IS THE SHARPE RATIO, AND HOW CAN EVERYONE GET IT WRONG? arxv:1802.04413v1 [q-fn.pm] 13 Feb 2018 IGOR RIVIN Abstract. The Sharpe rato s the most wdely used rsk metrc n the quanttatve fnance communty
More informationCopyright 2017 by Taylor Enterprises, Inc., All Rights Reserved. Dr. Wayne A. Taylor
Taylor Enterprses, Inc. ormalzed Indvduals (I ) Chart Copyrght 07 by Taylor Enterprses, Inc., All Rghts Reserved. ormalzed Indvduals (I) Control Chart Dr. Wayne A. Taylor Abstract: The only commonly used
More informationNotes on experimental uncertainties and their propagation
Ed Eyler 003 otes on epermental uncertantes and ther propagaton These notes are not ntended as a complete set of lecture notes, but nstead as an enumeraton of some of the key statstcal deas needed to obtan
More informationUsing Cumulative Count of Conforming CCC-Chart to Study the Expansion of the Cement
IOSR Journal of Engneerng (IOSRJEN) e-issn: 225-32, p-issn: 2278-879, www.osrjen.org Volume 2, Issue (October 22), PP 5-6 Usng Cumulatve Count of Conformng CCC-Chart to Study the Expanson of the Cement
More informationSelf-controlled case series analyses: small sample performance
Self-controlled case seres analyses: small sample performance Patrck Musonda 1, Mouna N. Hocne 1,2, Heather J. Whtaker 1 and C. Paddy Farrngton 1 * 1 The Open Unversty, Mlton Keynes, MK7 6AA, UK 2 INSERM
More informationNew Distance Measures on Dual Hesitant Fuzzy Sets and Their Application in Pattern Recognition
Journal of Artfcal Intellgence Practce (206) : 8-3 Clausus Scentfc Press, Canada New Dstance Measures on Dual Hestant Fuzzy Sets and Ther Applcaton n Pattern Recognton L Xn a, Zhang Xaohong* b College
More information15-451/651: Design & Analysis of Algorithms January 22, 2019 Lecture #3: Amortized Analysis last changed: January 18, 2019
5-45/65: Desgn & Analyss of Algorthms January, 09 Lecture #3: Amortzed Analyss last changed: January 8, 09 Introducton In ths lecture we dscuss a useful form of analyss, called amortzed analyss, for problems
More informationA FRAMEWORK FOR PRIORITY CONTACT OF NON RESPONDENTS
A FRAMEWORK FOR PRIORITY CONTACT OF NON RESPONDENTS Rchard McKenze, Australan Bureau of Statstcs. 12p36 Exchange Plaza, GPO Box K881, Perth, WA 6001. rchard.mckenze@abs.gov.au ABSTRACT Busnesses whch have
More information>1 indicates country i has a comparative advantage in production of j; the greater the index, the stronger the advantage. RCA 1 ij
69 APPENDIX 1 RCA Indces In the followng we present some maor RCA ndces reported n the lterature. For addtonal varants and other RCA ndces, Memedovc (1994) and Vollrath (1991) provde more thorough revews.
More informationMaturity Effect on Risk Measure in a Ratings-Based Default-Mode Model
TU Braunschweg - Insttut für Wrtschaftswssenschaften Lehrstuhl Fnanzwrtschaft Maturty Effect on Rsk Measure n a Ratngs-Based Default-Mode Model Marc Gürtler and Drk Hethecker Fnancal Modellng Workshop
More information3/3/2014. CDS M Phil Econometrics. Vijayamohanan Pillai N. Truncated standard normal distribution for a = 0.5, 0, and 0.5. CDS Mphil Econometrics
Lmted Dependent Varable Models: Tobt an Plla N 1 CDS Mphl Econometrcs Introducton Lmted Dependent Varable Models: Truncaton and Censorng Maddala, G. 1983. Lmted Dependent and Qualtatve Varables n Econometrcs.
More informationHarmonised Labour Cost Index. Methodology
Harmonsed Labour Cost Index Methodology March 2013 Index 1 Introducton 3 2 Scope, coverage and reference perod 4 3 Defntons 5 4 Sources of nformaton 7 5 Formulae employed 9 6 Results obtaned 10 7 Seres
More informationCHAPTER 3: BAYESIAN DECISION THEORY
CHATER 3: BAYESIAN DECISION THEORY Decson makng under uncertanty 3 rogrammng computers to make nference from data requres nterdscplnary knowledge from statstcs and computer scence Knowledge of statstcs
More informationExamining the Validity of Credit Ratings Assigned to Credit Derivatives
Examnng the Valdty of redt atngs Assgned to redt Dervatves hh-we Lee Department of Fnance, Natonal Tape ollege of Busness No. 321, Sec. 1, h-nan d., Tape 100, Tawan heng-kun Kuo Department of Internatonal
More informationData Mining Linear and Logistic Regression
07/02/207 Data Mnng Lnear and Logstc Regresson Mchael L of 26 Regresson In statstcal modellng, regresson analyss s a statstcal process for estmatng the relatonshps among varables. Regresson models are
More informationComparison of Singular Spectrum Analysis and ARIMA
Int. Statstcal Inst.: Proc. 58th World Statstcal Congress, 0, Dubln (Sesson CPS009) p.99 Comparson of Sngular Spectrum Analss and ARIMA Models Zokae, Mohammad Shahd Behesht Unverst, Department of Statstcs
More informationUsing Conditional Heteroskedastic
ITRON S FORECASTING BROWN BAG SEMINAR Usng Condtonal Heteroskedastc Varance Models n Load Research Sample Desgn Dr. J. Stuart McMenamn March 6, 2012 Please Remember» Phones are Muted: In order to help
More informationCHAPTER 9 FUNCTIONAL FORMS OF REGRESSION MODELS
CHAPTER 9 FUNCTIONAL FORMS OF REGRESSION MODELS QUESTIONS 9.1. (a) In a log-log model the dependent and all explanatory varables are n the logarthmc form. (b) In the log-ln model the dependent varable
More informationMultifactor Term Structure Models
1 Multfactor Term Structure Models A. Lmtatons of One-Factor Models 1. Returns on bonds of all maturtes are perfectly correlated. 2. Term structure (and prces of every other dervatves) are unquely determned
More informationMeasurement of Dynamic Portfolio VaR Based on Mixed Vine Copula Model
Journal of Fnance and Accountng 207; 5(2): 80-86 http://www.scencepublshnggroup.com/j/jfa do: 0.648/j.jfa.2070502.2 ISSN: 2330-733 (Prnt); ISSN: 2330-7323 (Onlne) Measurement of Dynamc Portfolo VaR Based
More informationApplications of Myerson s Lemma
Applcatons of Myerson s Lemma Professor Greenwald 28-2-7 We apply Myerson s lemma to solve the sngle-good aucton, and the generalzaton n whch there are k dentcal copes of the good. Our objectve s welfare
More informationA Comparison of Statistical Methods in Interrupted Time Series Analysis to Estimate an Intervention Effect
Transport and Road Safety (TARS) Research Joanna Wang A Comparson of Statstcal Methods n Interrupted Tme Seres Analyss to Estmate an Interventon Effect Research Fellow at Transport & Road Safety (TARS)
More informationA Comparative Study of Mean-Variance and Mean Gini Portfolio Selection Using VaR and CVaR
Journal of Fnancal Rsk Management, 5, 4, 7-8 Publshed Onlne 5 n ScRes. http://www.scrp.org/journal/jfrm http://dx.do.org/.436/jfrm.5.47 A Comparatve Study of Mean-Varance and Mean Gn Portfolo Selecton
More informationAppendix for Solving Asset Pricing Models when the Price-Dividend Function is Analytic
Appendx for Solvng Asset Prcng Models when the Prce-Dvdend Functon s Analytc Ovdu L. Caln Yu Chen Thomas F. Cosmano and Alex A. Hmonas January 3, 5 Ths appendx provdes proofs of some results stated n our
More informationProbability distribution of multi-hop-distance in one-dimensional sensor networks q
Computer etworks (7) 77 79 www.elsever.com/locate/comnet Probablty dstrbuton of mult-hop-dstance n one-dmensonal sensor networks q Serdar Vural *, Eylem Ekc Department of Electrcal and Computer Engneerng,
More informationISyE 512 Chapter 9. CUSUM and EWMA Control Charts. Instructor: Prof. Kaibo Liu. Department of Industrial and Systems Engineering UW-Madison
ISyE 512 hapter 9 USUM and EWMA ontrol harts Instructor: Prof. Kabo Lu Department of Industral and Systems Engneerng UW-Madson Emal: klu8@wsc.edu Offce: Room 317 (Mechancal Engneerng Buldng) ISyE 512 Instructor:
More informationEXAMINATIONS OF THE HONG KONG STATISTICAL SOCIETY
EXAMINATIONS OF THE HONG KONG STATISTICAL SOCIETY HIGHER CERTIFICATE IN STATISTICS, 2013 MODULE 7 : Tme seres and ndex numbers Tme allowed: One and a half hours Canddates should answer THREE questons.
More informationMarket Opening and Stock Market Behavior: Taiwan s Experience
Internatonal Journal of Busness and Economcs, 00, Vol., No., 9-5 Maret Openng and Stoc Maret Behavor: Tawan s Experence Q L * Department of Economcs, Texas A&M Unversty, U.S.A. and Department of Economcs,
More informationAlternatives to Shewhart Charts
Alternatves to Shewhart Charts CUSUM & EWMA S Wongsa Overvew Revstng Shewhart Control Charts Cumulatve Sum (CUSUM) Control Chart Eponentally Weghted Movng Average (EWMA) Control Chart 2 Revstng Shewhart
More informationElton, Gruber, Brown and Goetzmann. Modern Portfolio Theory and Investment Analysis, 7th Edition. Solutions to Text Problems: Chapter 4
Elton, Gruber, Brown and Goetzmann Modern ortfolo Theory and Investment Analyss, 7th Edton Solutons to Text roblems: Chapter 4 Chapter 4: roblem 1 A. Expected return s the sum of each outcome tmes ts assocated
More informationUnderstanding price volatility in electricity markets
Proceedngs of the 33rd Hawa Internatonal Conference on System Scences - 2 Understandng prce volatlty n electrcty markets Fernando L. Alvarado, The Unversty of Wsconsn Rajesh Rajaraman, Chrstensen Assocates
More informationAn Application of Alternative Weighting Matrix Collapsing Approaches for Improving Sample Estimates
Secton on Survey Research Methods An Applcaton of Alternatve Weghtng Matrx Collapsng Approaches for Improvng Sample Estmates Lnda Tompkns 1, Jay J. Km 2 1 Centers for Dsease Control and Preventon, atonal
More informationWhich of the following provides the most reasonable approximation to the least squares regression line? (a) y=50+10x (b) Y=50+x (d) Y=1+50x
Whch of the followng provdes the most reasonable approxmaton to the least squares regresson lne? (a) y=50+10x (b) Y=50+x (c) Y=10+50x (d) Y=1+50x (e) Y=10+x In smple lnear regresson the model that s begn
More informationDomestic Savings and International Capital Flows
Domestc Savngs and Internatonal Captal Flows Martn Feldsten and Charles Horoka The Economc Journal, June 1980 Presented by Mchael Mbate and Chrstoph Schnke Introducton The 2 Vews of Internatonal Captal
More informationImpact of CDO Tranches on Economic Capital of Credit Portfolios
Impact of CDO Tranches on Economc Captal of Credt Portfolos Ym T. Lee Market & Investment Bankng UnCredt Group Moor House, 120 London Wall London, EC2Y 5ET KEYWORDS: Credt rsk, Collateralzaton Debt Oblgaton,
More informationarxiv:cond-mat/ v1 [cond-mat.other] 28 Nov 2004
arxv:cond-mat/0411699v1 [cond-mat.other] 28 Nov 2004 Estmatng Probabltes of Default for Low Default Portfolos Katja Pluto and Drk Tasche November 23, 2004 Abstract For credt rsk management purposes n general,
More informationThe Effects of Industrial Structure Change on Economic Growth in China Based on LMDI Decomposition Approach
216 Internatonal Conference on Mathematcal, Computatonal and Statstcal Scences and Engneerng (MCSSE 216) ISBN: 978-1-6595-96- he Effects of Industral Structure Change on Economc Growth n Chna Based on
More informationDOUBLE IMPACT. Credit Risk Assessment for Secured Loans. Jean-Paul Laurent ISFA Actuarial School University of Lyon & BNP Paribas
DOUBLE IMPACT Credt Rsk Assessment for Secured Loans Al Chabaane BNP Parbas Jean-Paul Laurent ISFA Actuaral School Unversty of Lyon & BNP Parbas Julen Salomon BNP Parbas julen.salomon@bnpparbas.com Abstract
More informationDiscounted Cash Flow (DCF) Analysis: What s Wrong With It And How To Fix It
Dscounted Cash Flow (DCF Analyss: What s Wrong Wth It And How To Fx It Arturo Cfuentes (* CREM Facultad de Economa y Negocos Unversdad de Chle June 2014 (* Jont effort wth Francsco Hawas; Depto. de Ingenera
More informationStochastic Generation of Daily Rainfall Data
Stochastc Generaton of Daly Ranfall Data Srkanthan, R. CRC for Catchment Hydrology, Bureau of Meteorology, Melbourne, Australa, E-Mal: r.srkanthan@bom.gov.au Keywords: Stochastc generaton; daly ranfall;
More informationIntroduction. Why One-Pass Statistics?
BERKELE RESEARCH GROUP Ths manuscrpt s program documentaton for three ways to calculate the mean, varance, skewness, kurtoss, covarance, correlaton, regresson parameters and other regresson statstcs. Although
More informationIncreasing the Accuracy of Option Pricing by Using Implied Parameters Related to Higher Moments. Dasheng Ji. and. B. Wade Brorsen*
Increasng the Accuracy of Opton Prcng by Usng Impled Parameters Related to Hgher Moments Dasheng J and B. Wade Brorsen* Paper presented at the CR-34 Conference on Appled Commodty Prce Analyss, orecastng,
More informationOption Pricing Variance Reduction Techniques Under the Levy Process
Appled and Computatonal Mathematcs 205; 4(3): 74-80 Publshed onlne June 8, 205 (http://www.scencepublshnggroup.com//acm) do: 0.648/.acm.2050403.20 ISS: 2328-5605 (Prnt); ISS: 2328-563 (Onlne) Opton Prcng
More informationiii) pay F P 0,T = S 0 e δt when stock has dividend yield δ.
Fnal s Wed May 7, 12:50-2:50 You are allowed 15 sheets of notes and a calculator The fnal s cumulatve, so you should know everythng on the frst 4 revews Ths materal not on those revews 184) Suppose S t
More informationModule Contact: Dr P Moffatt, ECO Copyright of the University of East Anglia Version 2
UNIVERSITY OF EAST ANGLIA School of Economcs Man Seres PG Examnaton 2012-13 FINANCIAL ECONOMETRICS ECO-M017 Tme allowed: 2 hours Answer ALL FOUR questons. Queston 1 carres a weght of 25%; Queston 2 carres
More informationAppendix - Normally Distributed Admissible Choices are Optimal
Appendx - Normally Dstrbuted Admssble Choces are Optmal James N. Bodurtha, Jr. McDonough School of Busness Georgetown Unversty and Q Shen Stafford Partners Aprl 994 latest revson September 00 Abstract
More informationSpurious Seasonal Patterns and Excess Smoothness in the BLS Local Area Unemployment Statistics
Spurous Seasonal Patterns and Excess Smoothness n the BLS Local Area Unemployment Statstcs Keth R. Phllps and Janguo Wang Federal Reserve Bank of Dallas Research Department Workng Paper 1305 September
More informationFinancial mathematics
Fnancal mathematcs Jean-Luc Bouchot jean-luc.bouchot@drexel.edu February 19, 2013 Warnng Ths s a work n progress. I can not ensure t to be mstake free at the moment. It s also lackng some nformaton. But
More informationMonetary Tightening Cycles and the Predictability of Economic Activity. by Tobias Adrian and Arturo Estrella * October 2006.
Monetary Tghtenng Cycles and the Predctablty of Economc Actvty by Tobas Adran and Arturo Estrella * October 2006 Abstract Ten out of thrteen monetary tghtenng cycles snce 1955 were followed by ncreases
More informationCalibration Methods: Regression & Correlation. Calibration Methods: Regression & Correlation
Calbraton Methods: Regresson & Correlaton Calbraton A seres of standards run (n replcate fashon) over a gven concentraton range. Standards Comprsed of analte(s) of nterest n a gven matr composton. Matr
More informationBootstrap and Permutation tests in ANOVA for directional data
strap and utaton tests n ANOVA for drectonal data Adelade Fgueredo Faculty of Economcs of Unversty of Porto and LIAAD-INESC TEC Porto - PORTUGAL Abstract. The problem of testng the null hypothess of a
More informationSolution of periodic review inventory model with general constrains
Soluton of perodc revew nventory model wth general constrans Soluton of perodc revew nventory model wth general constrans Prof Dr J Benkő SZIU Gödöllő Summary Reasons for presence of nventory (stock of
More informationTechnological inefficiency and the skewness of the error component in stochastic frontier analysis
Economcs Letters 77 (00) 101 107 www.elsever.com/ locate/ econbase Technologcal neffcency and the skewness of the error component n stochastc fronter analyss Martn A. Carree a,b, * a Erasmus Unversty Rotterdam,
More informationElton, Gruber, Brown, and Goetzmann. Modern Portfolio Theory and Investment Analysis, 7th Edition. Solutions to Text Problems: Chapter 9
Elton, Gruber, Brown, and Goetzmann Modern Portfolo Theory and Investment Analyss, 7th Edton Solutons to Text Problems: Chapter 9 Chapter 9: Problem In the table below, gven that the rskless rate equals
More informationA MODEL OF COMPETITION AMONG TELECOMMUNICATION SERVICE PROVIDERS BASED ON REPEATED GAME
A MODEL OF COMPETITION AMONG TELECOMMUNICATION SERVICE PROVIDERS BASED ON REPEATED GAME Vesna Radonć Đogatovć, Valentna Radočć Unversty of Belgrade Faculty of Transport and Traffc Engneerng Belgrade, Serba
More informationPASS Sample Size Software. :log
PASS Sample Sze Software Chapter 70 Probt Analyss Introducton Probt and lot analyss may be used for comparatve LD 50 studes for testn the effcacy of drus desned to prevent lethalty. Ths proram module presents
More informationBasket options and implied correlations: a closed form approach
Basket optons and mpled correlatons: a closed form approach Svetlana Borovkova Free Unversty of Amsterdam CFC conference, London, January 7-8, 007 Basket opton: opton whose underlyng s a basket (.e. a
More informationECONOMETRICS - FINAL EXAM, 3rd YEAR (GECO & GADE)
ECONOMETRICS - FINAL EXAM, 3rd YEAR (GECO & GADE) May 17, 2016 15:30 Frst famly name: Name: DNI/ID: Moble: Second famly Name: GECO/GADE: Instructor: E-mal: Queston 1 A B C Blank Queston 2 A B C Blank Queston
More informationSpatial Variations in Covariates on Marriage and Marital Fertility: Geographically Weighted Regression Analyses in Japan
Spatal Varatons n Covarates on Marrage and Martal Fertlty: Geographcally Weghted Regresson Analyses n Japan Kenj Kamata (Natonal Insttute of Populaton and Socal Securty Research) Abstract (134) To understand
More informationТеоретические основы и методология имитационного и комплексного моделирования
MONTE-CARLO STATISTICAL MODELLING METHOD USING FOR INVESTIGA- TION OF ECONOMIC AND SOCIAL SYSTEMS Vladmrs Jansons, Vtaljs Jurenoks, Konstantns Ddenko (Latva). THE COMMO SCHEME OF USI G OF TRADITIO AL METHOD
More information