Using the IML Procedure to Examine the Efficacy of a New Control Charting Technique
|
|
- Jemimah Elliott
- 5 years ago
- Views:
Transcription
1 Paper Usig the IML Procedure to Examie the Efficacy of a New Cotrol Chartig Techique Austi Brow, M.S., Uiversity of Norther Colorado; Bryce Whitehead, M.S., Uiversity of Norther Colorado ABSTRACT For may years, cotrol charts have bee utilized to moitor processes, improve quality, ad icrease profitability. However, the body of literature tilts overwhelmigly toward charts moitorig ormally distributed processes. I practice, the uderlyig distributio of a process may ot follow a ormal distributio, ad may of those techiques may ot be most effective. Mukherjee, McCracke, ad Chakraborti (2015) suggested three cotrol charts for simultaeous moitorig of the locatio ad scale parameters for processes followig the shifted expoetial distributio. This study examies their proposed, Shifted Expoetial Maximum Likelihood Estimator-Max Chart, (SEMLE-max) ad suggests utilizatio of pealized maximum likelihood estimators (MLE) istead of traditioal MLEs because of ubiasedess ad miimum variace amog ubiased estimators. The ew chart, the Pealized SEMLEmax chart is costructed usig similar methodology, ad simulated data is used to compare average ru legths of the proposed chart to those obtaied by the origial chart. INTRODUCTION I the work published by Mukherjee, McCracke, ad Chakraborti (2015), ample justificatio was give to the creatio ad practical implemetatio of o-ormally distributed statistical process cotrol charts. Specifically, the authors focus was upo moitorig time-to-evet processes, as frequetly modeled by the shifted expoetial distributio: f(x) = 1 θ e 1 θ (x η), x > η, < η < I the publicatio, the authors itroduced three cotrol charts explicitly desiged for such processes. The cotrol chart of iterest i this study is their proposed, Shifted Expoetial Maximum Likelihood Estimator Max Chart, or SEMLE-max. I this chartig scheme, the maximum likelihood estimators (MLE) for the scale ad locatio parameters are used to build two plottig statistics based upo the stadard ormal distributio. I order to avoid issues arisig from joit moitorig of two processes, they authors suggest usig a sigle plottig statistic for each respective sample take, which is the absolute value of the maximum of the two respective statistics. The drawback to usig the traditioal MLEs i this circumstace is that they are biased estimators of the parameters. Sarha (1954) derived two ew estimators for the scale ad locatio parameters of the shifted expoetial distributio, usig a least squares techique, which are ubiased ad uiformly miimum variace for ubiased estimators. Zheg (2013) showed these estimators ca be derived from a pealized likelihood fuctio. Thus, i this study we recommed utilizig Sarhas estimators as opposed to the traditioal MLEs for the desirable qualities aforemetioed. We will use Mote Carlo simulatio methods i PROC IML to compare the i-cotrol ad out-of-cotrol average ru legths (IC-ARL & OOC-ARL) betwee the origial chart, ad the proposed modified chart. ESTIMATORS Similar to the derivatio of the traditioal MLEs for the scale ad locatio paramters ( θ ML = X X 1:, η ML = X 1: ), Zheg (2013) itroduced a pealty to the likelihood fuctio of the radom sample. The purpose of the pealty is with specific respect to the overestimatio of η, as P[η ML = X 1: > η] = 1. Thus, 1
2 the differece betwee the miimum order statistic ad η is defied to be the likelihood fuctio s pealty, ad this differece scales the likelihood fuctio, as follows: (X 1: η)l(θ, η) = (X 1: η) 1 θ e 1 θ (x i η) Usig the traditioal MLE derivatio techique, the pealized MLEs for η ad θ are: Ad: l(l(θ, η)) θ l(l(θ, η)) η i=1 = θ + 1 θ 2 (x i η ) = 0 i=1 θ = 1 (x i η ) i=1 1 = X 1: η + θ = 0 η = X 1: θ Pluggig η ito the equatio for θ, the estimators the become: These estimators are ubiased: η = X 1: X 1 θ = 1 (X X 1: ) E[η ] = E [ X 1: X 1 ] = 1 1 (E[X 1: ] E[X ]) = 1 1 ( (θ + η) η θ) = η E[θ ] = E [ 1 (X X 1: )] = 1 (E[X ] E[X 1: ]) = ( 1) (θ + η θ ) = θ Sarha (1954) also proved these estimators to be uiformly miimum variace amog the ubiased estimators. It should also be oted that θ is distributed as a gamma radom variable, with parameters θ/( 1) ad κ = 1, early idetically to the origial MLE, which was show i Mukherjee et all (2015). A closed form for the distributio of η could ot be derived, ad thus, its desity was estimated usig simulatio. CONSTRUCTION OF THE CONTROL CHART As oted i Mukherjee et al (2015), there are a array of issues associated with joitly moitorig two or more ukow parameters. Oe solutio proposed i their work was calculatig two idetically distributed statistics, ad desigatig the plottig statistic for a give sample as the maximum betwee the two. These two statistics are stadardized (as they are ot idetically distributed i their origial form) usig the iverse cumulative desity fuctio (CDF) of the stadard ormal distributio, i.e.: B 1 = Φ 1 (F(x)) B 2 = Φ 1 (G(x)) 2
3 Where Φ 1 deotes the iverse CDF of the stadard ormal distributio, F(x) deotes the CDF of η, ad G(x) deotes the CDF of θ. For a give sample of the process beig moitored, both B 1 ad B 2 will be calculated ad the maximum of their absolute value, say M i, becomes the plottig statistic. Mathematically, M i = max( B 1, B 2 ) Sice oly positive values are yielded by the plottig statistic, oly a upper cotrol limit (UCL) is ecessary to specify. For a desired i-cotrol average ru legth (IC-ARL), called ARL 0, Mukherjee et al (2015) showed the UCL ca be determied by: UCL = Φ 1 (0.5 (1 + 1 (ARL 0 ) 1 )) Thus, for a desired IC-ARL of 500, the UCL to be used is Note, this IC-ARL will be used i comparig out-of-cotrol ARL (OOC-ARL) performace of the modified chart to the chart origially proposed i Mukherjee et al (2015). Therefore, for a give sample, the process is deemed OOC whe M i > UCL. Of ote, a OOC sigal ca assist i idicatig which of the two parameters has shifted away from its specified IC value. Whe M i = B 1 > UCL, this implies η may have shifted, ad coversely, whe M i = B 2 > UCL, this implies θ may have shifted. However, whe both B 1 ad B 2 are greater tha the UCL, it becomes difficult to determie which of the two values (or possibly both) is sigalig the OOC without physical ivestigatio. SIMULTATION METHODS & RESULTS Give the guidelies for the cotrol charts costructio as outlied i the previous sectio, the code for the simulatio was writte i similar steps. First, a radom sample from the OOC shifted expoetial distributio had to be simulated. This was performed usig the iverse probability itegral trasform. This trasformatio utilizes the kowledge of a radom variables cumulative desity fuctio (CDF) beig distributed uiformly betwee 0 & 1 (i.e., Y = F(x) UNIF(0,1)). Values from the uiform distributio are evaluated i a give radom variables iverse CDF. I this case, the shifted expoetial distributios CDF was set equal to Y ad solved for X as give by: Solvig for X: Y = 1 e 1 θ (x η) X = η θ l(1 Y) Thus, values radomly sampled from UNIF(0,1) are iput for Y, ad the resultig X values are distributed as the shifted expoetial for a give value of η ad θ. Usig two DO loops i the DATA step, 10,000 samples, each of sample size = 5, were simulated, ad the estimators for the respective parameters from the OOC distributio were obtaied usig PROC MEANS. At this poit, the plottig statistic for θ, B 2, was obtaied usig the PROBCHI ad QUANTILE fuctios i the DATA step. Here, the probability associated with a give observed value of θ (usig PROBCHI) was iputted ito the iverse CDF of the stadard ormal distributio (i.e., QUANTILE( NORMAL ) ), ad the absolute value of the resultig z-score is the value B 2 for a give sample. Seeig as the desity of η was estimated usig simulatio by ecessity, a hadful of extra procedures were implemeted i order to obtai its plottig statistic, B 1. Usig the same procedure described above, a radom sample from the IC shifted expoetial distributio was obtaied. This radom sample acts as the desity for η. Usig PROC IML, the OOC value of η from each respective sample was iteratively compared to the simulated IC distributio i order to obtai a pseudo-probability. Z-scores for this vector of pseudo-probabilities was obtaied i order to yield the plottig statistic B 1 i a early idetical maer as described for B 2. Now that values for B 1 ad B 2 for each of the respective 10,000 samples are obtaied, they ca be plotted agaist the specified UCL of 3.29 i order to evaluate the effectiveess of these ew estimators. This process was performed iteratively usig PROC IML. Whe a sample plotted above the UCL, thus sigalig a OOC process, the sample umber was 4 recorded, ad the mea differece betwee each 3
4 observed OOC poit was take to be the OOC-ARL for a give pair of OOC values of θ ad η. The etire aforemetioed process was performed 10 times for each pair of OOC values of θ ad η, ad mea of those results were take to be the OOC-ARL used for compariso agaist the Mukherjee et al (2015) chart s performace. The procedure described above was performed for 20 paris of OOC values of η ad θ, where the specified IC values were η = 0 ad θ = 1, respectively, ad the sample size take was = 5. The results of this study are give i the below table, deoted by ARL S, ad the results obtaied i the origial study are deoted ARL C. Table 1: Compariso of Simulated OOC-ARL η θ ARL S ARL C CONCLUSION As show i Table 1, the utilizatio of the pealized MLEs give improved performace over that of the stadard estimators for the dowward shifts i θ coupled with icreases i η. Both chartig schemes, as the locatio parameter shifts upward away from its omial IC value, are highly effective i quickly detectig these large shifts, which is a similar coclusio draw i Mukherjee et al (2015). However, cosiderig the small scale, oe may assume the aalyzed OOC shifts are small, but the magitude of the shifts is quite large. While the proposed modificatios to the estimators yield improvemets i some cases, such as the istace whe OOC value of θ = 0.75 ad η = 0.5, a ))C-ARL of , while a improvemet over , is still far too slow for a 25% shift i θ ad a icreased η. 4
5 A big cotributor to this relatively slow performace is the costructio of the chart itself. While creatively thought out ad mathematically elegat, this techique leds itself oly to detectio of very extreme values, as evideced by the quick detectio rate at large OOC η. However, cosider the case of θ. Sice it ca be trasformed ito a χ 8 2 radom variable i the case whe = 5, it ca easily be determied what the z-score would be for a particular quatile of the χ 8 2 distributio. Say θ is perfectly estimated to be 0.50, a 50% reductio i the IC value of θ. The probability associated with this estimate is about 0.24, which is associated with a z-score of about -0.70, which certaily would ot sigal a shift i θ. I fact, settig the UCL to 3.29 suggests that i order for a OOC sigal to be detected, the lower-tailed probability associated with our estimators would have to be approximately Thus, as show by the results i Table 1, uless η has shifted substatially from 0, it is highly improbable to detect a OOC shift. Cosequetly, it would ot be recommeded to implemet this scheme with such a high UCL uless specific circumstaces warrat it. While a smaller UCL (for example, 2.78) would sacrifice the IC performace of the chartig scheme i both istaces, it would make up the differece by beig far more sesitive to shifts i either parameter. 5
6 REFERENCES Mukherjee, A., McCracke, A., & Chakraborti, S. (2015). Cotrol charts for simultaeous moitorig of parameters of a shifted expoetial distributio. Joural of Quality Techology, 47(2), Sarha, A.E. (1954). Estimatio of the Mea ad Stadard Deviatio by Order Statistics. The Aals of Mathematical Statistics, Vol. 25, Zheg, M. (2013). Pealized Maximum Likelihood Estimatio of Two-parameter Expoetial Distributios (Doctoral dissertatio, Master Thesis, Uiversity of Miesota). 6
7 CONTACT INFORMATION Your commets ad questios are valued ad ecouraged. Cotact the authors at: Austi Brow & Bryce Whitehead Uiversity of Norther Colorado & SAS ad all other SAS Istitute Ic. product or service ames are registered trademarks or trademarks of SAS Istitute Ic. i the USA ad other coutries. idicates USA registratio. Other brad ad product ames are trademarks of their respective compaies. 7
8 8
Topic 9: Sampling Distributions of Estimators
Topic 9: Samplig Distributios of Estimators Course 003, 2018 Page 0 Samplig distributios of estimators Sice our estimators are statistics (particular fuctios of radom variables), their distributio ca be
More informationTopic 9: Sampling Distributions of Estimators
Topic 9: Samplig Distributios of Estimators Course 003, 2016 Page 0 Samplig distributios of estimators Sice our estimators are statistics (particular fuctios of radom variables), their distributio ca be
More informationTopic 9: Sampling Distributions of Estimators
Topic 9: Samplig Distributios of Estimators Course 003, 2018 Page 0 Samplig distributios of estimators Sice our estimators are statistics (particular fuctios of radom variables), their distributio ca be
More informationSince X n /n P p, we know that X n (n. Xn (n X n ) Using the asymptotic result above to obtain an approximation for fixed n, we obtain
Assigmet 9 Exercise 5.5 Let X biomial, p, where p 0, 1 is ukow. Obtai cofidece itervals for p i two differet ways: a Sice X / p d N0, p1 p], the variace of the limitig distributio depeds oly o p. Use the
More informationR. van Zyl 1, A.J. van der Merwe 2. Quintiles International, University of the Free State
Bayesia Cotrol Charts for the Two-parameter Expoetial Distributio if the Locatio Parameter Ca Take o Ay Value Betwee Mius Iity ad Plus Iity R. va Zyl, A.J. va der Merwe 2 Quitiles Iteratioal, ruaavz@gmail.com
More information1 Introduction to reducing variance in Monte Carlo simulations
Copyright c 010 by Karl Sigma 1 Itroductio to reducig variace i Mote Carlo simulatios 11 Review of cofidece itervals for estimatig a mea I statistics, we estimate a ukow mea µ = E(X) of a distributio by
More informationSimulation. Two Rule For Inverting A Distribution Function
Simulatio Two Rule For Ivertig A Distributio Fuctio Rule 1. If F(x) = u is costat o a iterval [x 1, x 2 ), the the uiform value u is mapped oto x 2 through the iversio process. Rule 2. If there is a jump
More informationStatistical Inference (Chapter 10) Statistical inference = learn about a population based on the information provided by a sample.
Statistical Iferece (Chapter 10) Statistical iferece = lear about a populatio based o the iformatio provided by a sample. Populatio: The set of all values of a radom variable X of iterest. Characterized
More informationAAEC/ECON 5126 FINAL EXAM: SOLUTIONS
AAEC/ECON 5126 FINAL EXAM: SOLUTIONS SPRING 2015 / INSTRUCTOR: KLAUS MOELTNER This exam is ope-book, ope-otes, but please work strictly o your ow. Please make sure your ame is o every sheet you re hadig
More informationA new distribution-free quantile estimator
Biometrika (1982), 69, 3, pp. 635-40 Prited i Great Britai 635 A ew distributio-free quatile estimator BY FRANK E. HARRELL Cliical Biostatistics, Duke Uiversity Medical Ceter, Durham, North Carolia, U.S.A.
More informationLecture Note 8 Point Estimators and Point Estimation Methods. MIT Spring 2006 Herman Bennett
Lecture Note 8 Poit Estimators ad Poit Estimatio Methods MIT 14.30 Sprig 2006 Herma Beett Give a parameter with ukow value, the goal of poit estimatio is to use a sample to compute a umber that represets
More informationComparing Two Populations. Topic 15 - Two Sample Inference I. Comparing Two Means. Comparing Two Pop Means. Background Reading
Topic 15 - Two Sample Iferece I STAT 511 Professor Bruce Craig Comparig Two Populatios Research ofte ivolves the compariso of two or more samples from differet populatios Graphical summaries provide visual
More informationCONTROL CHARTS FOR THE LOGNORMAL DISTRIBUTION
CONTROL CHARTS FOR THE LOGNORMAL DISTRIBUTION Petros Maravelakis, Joh Paaretos ad Stelios Psarakis Departmet of Statistics Athes Uiversity of Ecoomics ad Busiess 76 Patisio St., 4 34, Athes, GREECE. Itroductio
More informationRandom Variables, Sampling and Estimation
Chapter 1 Radom Variables, Samplig ad Estimatio 1.1 Itroductio This chapter will cover the most importat basic statistical theory you eed i order to uderstad the ecoometric material that will be comig
More informationDepartment of Mathematics
Departmet of Mathematics Ma 3/103 KC Border Itroductio to Probability ad Statistics Witer 2017 Lecture 19: Estimatio II Relevat textbook passages: Larse Marx [1]: Sectios 5.2 5.7 19.1 The method of momets
More informationGoodness-Of-Fit For The Generalized Exponential Distribution. Abstract
Goodess-Of-Fit For The Geeralized Expoetial Distributio By Amal S. Hassa stitute of Statistical Studies & Research Cairo Uiversity Abstract Recetly a ew distributio called geeralized expoetial or expoetiated
More informationBayesian Control Charts for the Two-parameter Exponential Distribution
Bayesia Cotrol Charts for the Two-parameter Expoetial Distributio R. va Zyl, A.J. va der Merwe 2 Quitiles Iteratioal, ruaavz@gmail.com 2 Uiversity of the Free State Abstract By usig data that are the mileages
More informationStat 421-SP2012 Interval Estimation Section
Stat 41-SP01 Iterval Estimatio Sectio 11.1-11. We ow uderstad (Chapter 10) how to fid poit estimators of a ukow parameter. o However, a poit estimate does ot provide ay iformatio about the ucertaity (possible
More informationProblem Set 4 Due Oct, 12
EE226: Radom Processes i Systems Lecturer: Jea C. Walrad Problem Set 4 Due Oct, 12 Fall 06 GSI: Assae Gueye This problem set essetially reviews detectio theory ad hypothesis testig ad some basic otios
More informationEstimation of Gumbel Parameters under Ranked Set Sampling
Joural of Moder Applied Statistical Methods Volume 13 Issue 2 Article 11-2014 Estimatio of Gumbel Parameters uder Raked Set Samplig Omar M. Yousef Al Balqa' Applied Uiversity, Zarqa, Jorda, abuyaza_o@yahoo.com
More informationECO 312 Fall 2013 Chris Sims LIKELIHOOD, POSTERIORS, DIAGNOSING NON-NORMALITY
ECO 312 Fall 2013 Chris Sims LIKELIHOOD, POSTERIORS, DIAGNOSING NON-NORMALITY (1) A distributio that allows asymmetry differet probabilities for egative ad positive outliers is the asymmetric double expoetial,
More informationResampling Methods. X (1/2), i.e., Pr (X i m) = 1/2. We order the data: X (1) X (2) X (n). Define the sample median: ( n.
Jauary 1, 2019 Resamplig Methods Motivatio We have so may estimators with the property θ θ d N 0, σ 2 We ca also write θ a N θ, σ 2 /, where a meas approximately distributed as Oce we have a cosistet estimator
More information1.010 Uncertainty in Engineering Fall 2008
MIT OpeCourseWare http://ocw.mit.edu.00 Ucertaity i Egieerig Fall 2008 For iformatio about citig these materials or our Terms of Use, visit: http://ocw.mit.edu.terms. .00 - Brief Notes # 9 Poit ad Iterval
More informationOutput Analysis (2, Chapters 10 &11 Law)
B. Maddah ENMG 6 Simulatio Output Aalysis (, Chapters 10 &11 Law) Comparig alterative system cofiguratio Sice the output of a simulatio is radom, the comparig differet systems via simulatio should be doe
More informationData Analysis and Statistical Methods Statistics 651
Data Aalysis ad Statistical Methods Statistics 651 http://www.stat.tamu.edu/~suhasii/teachig.html Suhasii Subba Rao Review of testig: Example The admistrator of a ursig home wats to do a time ad motio
More informationMonte Carlo method and application to random processes
Mote Carlo method ad applicatio to radom processes Lecture 3: Variace reductio techiques (8/3/2017) 1 Lecturer: Eresto Mordecki, Facultad de Ciecias, Uiversidad de la República, Motevideo, Uruguay Graduate
More informationChapter 6 Sampling Distributions
Chapter 6 Samplig Distributios 1 I most experimets, we have more tha oe measuremet for ay give variable, each measuremet beig associated with oe radomly selected a member of a populatio. Hece we eed to
More informationConfidence Interval for Standard Deviation of Normal Distribution with Known Coefficients of Variation
Cofidece Iterval for tadard Deviatio of Normal Distributio with Kow Coefficiets of Variatio uparat Niwitpog Departmet of Applied tatistics, Faculty of Applied ciece Kig Mogkut s Uiversity of Techology
More informationMaximum likelihood estimation from record-breaking data for the generalized Pareto distribution
METRON - Iteratioal Joural of Statistics 004, vol. LXII,. 3, pp. 377-389 NAGI S. ABD-EL-HAKIM KHALAF S. SULTAN Maximum likelihood estimatio from record-breakig data for the geeralized Pareto distributio
More informationEcon 325/327 Notes on Sample Mean, Sample Proportion, Central Limit Theorem, Chi-square Distribution, Student s t distribution 1.
Eco 325/327 Notes o Sample Mea, Sample Proportio, Cetral Limit Theorem, Chi-square Distributio, Studet s t distributio 1 Sample Mea By Hiro Kasahara We cosider a radom sample from a populatio. Defiitio
More informationA proposed discrete distribution for the statistical modeling of
It. Statistical Ist.: Proc. 58th World Statistical Cogress, 0, Dubli (Sessio CPS047) p.5059 A proposed discrete distributio for the statistical modelig of Likert data Kidd, Marti Cetre for Statistical
More informationPower Comparison of Some Goodness-of-fit Tests
Florida Iteratioal Uiversity FIU Digital Commos FIU Electroic Theses ad Dissertatios Uiversity Graduate School 7-6-2016 Power Compariso of Some Goodess-of-fit Tests Tiayi Liu tliu019@fiu.edu DOI: 10.25148/etd.FIDC000750
More informationSoo King Lim Figure 1: Figure 2: Figure 3: Figure 4: Figure 5: Figure 6: Figure 7:
0 Multivariate Cotrol Chart 3 Multivariate Normal Distributio 5 Estimatio of the Mea ad Covariace Matrix 6 Hotellig s Cotrol Chart 6 Hotellig s Square 8 Average Value of k Subgroups 0 Example 3 3 Value
More informationMATH 320: Probability and Statistics 9. Estimation and Testing of Parameters. Readings: Pruim, Chapter 4
MATH 30: Probability ad Statistics 9. Estimatio ad Testig of Parameters Estimatio ad Testig of Parameters We have bee dealig situatios i which we have full kowledge of the distributio of a radom variable.
More informationExpectation and Variance of a random variable
Chapter 11 Expectatio ad Variace of a radom variable The aim of this lecture is to defie ad itroduce mathematical Expectatio ad variace of a fuctio of discrete & cotiuous radom variables ad the distributio
More information1036: Probability & Statistics
036: Probability & Statistics Lecture 0 Oe- ad Two-Sample Tests of Hypotheses 0- Statistical Hypotheses Decisio based o experimetal evidece whether Coffee drikig icreases the risk of cacer i humas. A perso
More informationEstimation for Complete Data
Estimatio for Complete Data complete data: there is o loss of iformatio durig study. complete idividual complete data= grouped data A complete idividual data is the oe i which the complete iformatio of
More informationMathematical Statistics - MS
Paper Specific Istructios. The examiatio is of hours duratio. There are a total of 60 questios carryig 00 marks. The etire paper is divided ito three sectios, A, B ad C. All sectios are compulsory. Questios
More informationc. Explain the basic Newsvendor model. Why is it useful for SC models? e. What additional research do you believe will be helpful in this area?
1. Research Methodology a. What is meat by the supply chai (SC) coordiatio problem ad does it apply to all types of SC s? Does the Bullwhip effect relate to all types of SC s? Also does it relate to SC
More informationClases 7-8: Métodos de reducción de varianza en Monte Carlo *
Clases 7-8: Métodos de reducció de variaza e Mote Carlo * 9 de septiembre de 27 Ídice. Variace reductio 2. Atithetic variates 2 2.. Example: Uiform radom variables................ 3 2.2. Example: Tail
More informationDiscrete Mathematics for CS Spring 2008 David Wagner Note 22
CS 70 Discrete Mathematics for CS Sprig 2008 David Wager Note 22 I.I.D. Radom Variables Estimatig the bias of a coi Questio: We wat to estimate the proportio p of Democrats i the US populatio, by takig
More informationEXAMINATIONS OF THE ROYAL STATISTICAL SOCIETY
EXAMINATIONS OF THE ROYAL STATISTICAL SOCIETY GRADUATE DIPLOMA, 016 MODULE : Statistical Iferece Time allowed: Three hours Cadidates should aswer FIVE questios. All questios carry equal marks. The umber
More informationChapter 6. Sampling and Estimation
Samplig ad Estimatio - 34 Chapter 6. Samplig ad Estimatio 6.. Itroductio Frequetly the egieer is uable to completely characterize the etire populatio. She/he must be satisfied with examiig some subset
More informationStat 319 Theory of Statistics (2) Exercises
Kig Saud Uiversity College of Sciece Statistics ad Operatios Research Departmet Stat 39 Theory of Statistics () Exercises Refereces:. Itroductio to Mathematical Statistics, Sixth Editio, by R. Hogg, J.
More informationLecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting
Lecture 6 Chi Square Distributio (χ ) ad Least Squares Fittig Chi Square Distributio (χ ) Suppose: We have a set of measuremets {x 1, x, x }. We kow the true value of each x i (x t1, x t, x t ). We would
More informationExtreme Value Charts and Analysis of Means (ANOM) Based on the Log Logistic Distribution
Joural of Moder Applied Statistical Methods Volume 11 Issue Article 0 11-1-01 Extreme Value Charts ad Aalysis of Meas (ANOM) Based o the Log Logistic istributio B. Sriivasa Rao R.V.R & J.C. College of
More information1 of 7 7/16/2009 6:06 AM Virtual Laboratories > 6. Radom Samples > 1 2 3 4 5 6 7 6. Order Statistics Defiitios Suppose agai that we have a basic radom experimet, ad that X is a real-valued radom variable
More informationENGI 4421 Probability and Statistics Faculty of Engineering and Applied Science Problem Set 1 Solutions Descriptive Statistics. None at all!
ENGI 44 Probability ad Statistics Faculty of Egieerig ad Applied Sciece Problem Set Solutios Descriptive Statistics. If, i the set of values {,, 3, 4, 5, 6, 7 } a error causes the value 5 to be replaced
More informationLecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting
Lecture 6 Chi Square Distributio (χ ) ad Least Squares Fittig Chi Square Distributio (χ ) Suppose: We have a set of measuremets {x 1, x, x }. We kow the true value of each x i (x t1, x t, x t ). We would
More informationx = Pr ( X (n) βx ) =
Exercise 93 / page 45 The desity of a variable X i i 1 is fx α α a For α kow let say equal to α α > fx α α x α Pr X i x < x < Usig a Pivotal Quatity: x α 1 < x < α > x α 1 ad We solve i a similar way as
More informationIf, for instance, we were required to test whether the population mean μ could be equal to a certain value μ
STATISTICAL INFERENCE INTRODUCTION Statistical iferece is that brach of Statistics i which oe typically makes a statemet about a populatio based upo the results of a sample. I oesample testig, we essetially
More informationModified Lilliefors Test
Joural of Moder Applied Statistical Methods Volume 14 Issue 1 Article 9 5-1-2015 Modified Lilliefors Test Achut Adhikari Uiversity of Norther Colorado, adhi2939@gmail.com Jay Schaffer Uiversity of Norther
More informationCHAPTER 4 BIVARIATE DISTRIBUTION EXTENSION
CHAPTER 4 BIVARIATE DISTRIBUTION EXTENSION 4. Itroductio Numerous bivariate discrete distributios have bee defied ad studied (see Mardia, 97 ad Kocherlakota ad Kocherlakota, 99) based o various methods
More informationLet us give one more example of MLE. Example 3. The uniform distribution U[0, θ] on the interval [0, θ] has p.d.f.
Lecture 5 Let us give oe more example of MLE. Example 3. The uiform distributio U[0, ] o the iterval [0, ] has p.d.f. { 1 f(x =, 0 x, 0, otherwise The likelihood fuctio ϕ( = f(x i = 1 I(X 1,..., X [0,
More informationSample Size Determination (Two or More Samples)
Sample Sie Determiatio (Two or More Samples) STATGRAPHICS Rev. 963 Summary... Data Iput... Aalysis Summary... 5 Power Curve... 5 Calculatios... 6 Summary This procedure determies a suitable sample sie
More information6 Sample Size Calculations
6 Sample Size Calculatios Oe of the major resposibilities of a cliical trial statisticia is to aid the ivestigators i determiig the sample size required to coduct a study The most commo procedure for determiig
More informationLecture 7: Properties of Random Samples
Lecture 7: Properties of Radom Samples 1 Cotiued From Last Class Theorem 1.1. Let X 1, X,...X be a radom sample from a populatio with mea µ ad variace σ
More informationStatistical inference: example 1. Inferential Statistics
Statistical iferece: example 1 Iferetial Statistics POPULATION SAMPLE A clothig store chai regularly buys from a supplier large quatities of a certai piece of clothig. Each item ca be classified either
More informationDirection: This test is worth 250 points. You are required to complete this test within 50 minutes.
Term Test October 3, 003 Name Math 56 Studet Number Directio: This test is worth 50 poits. You are required to complete this test withi 50 miutes. I order to receive full credit, aswer each problem completely
More informationReview Questions, Chapters 8, 9. f(y) = 0, elsewhere. F (y) = f Y(1) = n ( e y/θ) n 1 1 θ e y/θ = n θ e yn
Stat 366 Lab 2 Solutios (September 2, 2006) page TA: Yury Petracheko, CAB 484, yuryp@ualberta.ca, http://www.ualberta.ca/ yuryp/ Review Questios, Chapters 8, 9 8.5 Suppose that Y, Y 2,..., Y deote a radom
More informationConfidence Intervals
Cofidece Itervals Berli Che Deartmet of Comuter Sciece & Iformatio Egieerig Natioal Taiwa Normal Uiversity Referece: 1. W. Navidi. Statistics for Egieerig ad Scietists. Chater 5 & Teachig Material Itroductio
More information1 Inferential Methods for Correlation and Regression Analysis
1 Iferetial Methods for Correlatio ad Regressio Aalysis I the chapter o Correlatio ad Regressio Aalysis tools for describig bivariate cotiuous data were itroduced. The sample Pearso Correlatio Coefficiet
More informationMOMENT-METHOD ESTIMATION BASED ON CENSORED SAMPLE
Vol. 8 o. Joural of Systems Sciece ad Complexity Apr., 5 MOMET-METHOD ESTIMATIO BASED O CESORED SAMPLE I Zhogxi Departmet of Mathematics, East Chia Uiversity of Sciece ad Techology, Shaghai 37, Chia. Email:
More informationFACULTY OF MATHEMATICAL STUDIES MATHEMATICS FOR PART I ENGINEERING. Lectures
FACULTY OF MATHEMATICAL STUDIES MATHEMATICS FOR PART I ENGINEERING Lectures MODULE 5 STATISTICS II. Mea ad stadard error of sample data. Biomial distributio. Normal distributio 4. Samplig 5. Cofidece itervals
More informationHypothesis Testing. Evaluation of Performance of Learned h. Issues. Trade-off Between Bias and Variance
Hypothesis Testig Empirically evaluatig accuracy of hypotheses: importat activity i ML. Three questios: Give observed accuracy over a sample set, how well does this estimate apply over additioal samples?
More informationOverview. p 2. Chapter 9. Pooled Estimate of. q = 1 p. Notation for Two Proportions. Inferences about Two Proportions. Assumptions
Chapter 9 Slide Ifereces from Two Samples 9- Overview 9- Ifereces about Two Proportios 9- Ifereces about Two Meas: Idepedet Samples 9-4 Ifereces about Matched Pairs 9-5 Comparig Variatio i Two Samples
More informationTHE SYSTEMATIC AND THE RANDOM. ERRORS - DUE TO ELEMENT TOLERANCES OF ELECTRICAL NETWORKS
R775 Philips Res. Repts 26,414-423, 1971' THE SYSTEMATIC AND THE RANDOM. ERRORS - DUE TO ELEMENT TOLERANCES OF ELECTRICAL NETWORKS by H. W. HANNEMAN Abstract Usig the law of propagatio of errors, approximated
More informationMinimax Estimation of the Parameter of Maxwell Distribution Under Different Loss Functions
America Joural of heoretical ad Applied Statistics 6; 5(4): -7 http://www.sciecepublishiggroup.com/j/ajtas doi:.648/j.ajtas.654.6 ISSN: 6-8999 (Prit); ISSN: 6-96 (Olie) Miimax Estimatio of the Parameter
More informationThis exam contains 19 pages (including this cover page) and 10 questions. A Formulae sheet is provided with the exam.
Probability ad Statistics FS 07 Secod Sessio Exam 09.0.08 Time Limit: 80 Miutes Name: Studet ID: This exam cotais 9 pages (icludig this cover page) ad 0 questios. A Formulae sheet is provided with the
More informationGoodness-Of-Fit For The Generalized Exponential Distribution. Abstract
Goodess-Of-Fit For The Geeralized Expoetial Distributio By Amal S. Hassa stitute of Statistical Studies & Research Cairo Uiversity Abstract Recetly a ew distributio called geeralized expoetial or expoetiated
More informationStatistical Fundamentals and Control Charts
Statistical Fudametals ad Cotrol Charts 1. Statistical Process Cotrol Basics Chace causes of variatio uavoidable causes of variatios Assigable causes of variatio large variatios related to machies, materials,
More information7-1. Chapter 4. Part I. Sampling Distributions and Confidence Intervals
7-1 Chapter 4 Part I. Samplig Distributios ad Cofidece Itervals 1 7- Sectio 1. Samplig Distributio 7-3 Usig Statistics Statistical Iferece: Predict ad forecast values of populatio parameters... Test hypotheses
More informationUnbiased Estimation. February 7-12, 2008
Ubiased Estimatio February 7-2, 2008 We begi with a sample X = (X,..., X ) of radom variables chose accordig to oe of a family of probabilities P θ where θ is elemet from the parameter space Θ. For radom
More informationThere is no straightforward approach for choosing the warmup period l.
B. Maddah INDE 504 Discrete-Evet Simulatio Output Aalysis () Statistical Aalysis for Steady-State Parameters I a otermiatig simulatio, the iterest is i estimatig the log ru steady state measures of performace.
More informationMONTE CARLO VARIANCE REDUCTION METHODS
MONTE CARLO VARIANCE REDUCTION METHODS M. Ragheb /4/3. INTRODUCTION The questio arises of whether oe ca reduce the variace associated with the samplig of a radom variable? Ideed we ca, but we eed to be
More informationConfidence Intervals For P(X less than Y) In The Exponential Case With Common Location Parameter
Joural of Moder Applied Statistical Methods Volume Issue Article 7 --3 Cofidece Itervals For P(X less tha Y I he Expoetial Case With Commo Locatio Parameter Ayma Baklizi Yarmouk Uiversity, Irbid, Jorda,
More informationUsing SAS to Evaluate Integrals and Reverse Functions in Power and Sample Size Calculations
Usig SAS to Evaluate Itegrals ad everse Fuctios i Power ad Sample Size Calculatios Xigshu Zhu, Merck &Co., Ic. Blue Bell, PA 94 Shupig Zhag, Merck &Co., Ic. Blue Bell, PA 94 ABSTACT We have recetly created
More information6. Sufficient, Complete, and Ancillary Statistics
Sufficiet, Complete ad Acillary Statistics http://www.math.uah.edu/stat/poit/sufficiet.xhtml 1 of 7 7/16/2009 6:13 AM Virtual Laboratories > 7. Poit Estimatio > 1 2 3 4 5 6 6. Sufficiet, Complete, ad Acillary
More informationDirection: This test is worth 150 points. You are required to complete this test within 55 minutes.
Term Test 3 (Part A) November 1, 004 Name Math 6 Studet Number Directio: This test is worth 10 poits. You are required to complete this test withi miutes. I order to receive full credit, aswer each problem
More information1 Models for Matched Pairs
1 Models for Matched Pairs Matched pairs occur whe we aalyse samples such that for each measuremet i oe of the samples there is a measuremet i the other sample that directly relates to the measuremet i
More informationESTIMATION AND PREDICTION BASED ON K-RECORD VALUES FROM NORMAL DISTRIBUTION
STATISTICA, ao LXXIII,. 4, 013 ESTIMATION AND PREDICTION BASED ON K-RECORD VALUES FROM NORMAL DISTRIBUTION Maoj Chacko Departmet of Statistics, Uiversity of Kerala, Trivadrum- 695581, Kerala, Idia M. Shy
More informationENGI 4421 Confidence Intervals (Two Samples) Page 12-01
ENGI 44 Cofidece Itervals (Two Samples) Page -0 Two Sample Cofidece Iterval for a Differece i Populatio Meas [Navidi sectios 5.4-5.7; Devore chapter 9] From the cetral limit theorem, we kow that, for sufficietly
More informationStatistical Theory MT 2008 Problems 1: Solution sketches
Statistical Theory MT 008 Problems : Solutio sketches. Which of the followig desities are withi a expoetial family? Explai your reasoig. a) Let 0 < θ < ad put fx, θ) = θ)θ x ; x = 0,,,... b) c) where α
More informationPSYCHOLOGICAL RESEARCH (PYC 304-C) Lecture 9
Hypothesis testig PSYCHOLOGICAL RESEARCH (PYC 34-C Lecture 9 Statistical iferece is that brach of Statistics i which oe typically makes a statemet about a populatio based upo the results of a sample. I
More informationCHAPTER 10 INFINITE SEQUENCES AND SERIES
CHAPTER 10 INFINITE SEQUENCES AND SERIES 10.1 Sequeces 10.2 Ifiite Series 10.3 The Itegral Tests 10.4 Compariso Tests 10.5 The Ratio ad Root Tests 10.6 Alteratig Series: Absolute ad Coditioal Covergece
More informationEcon 325 Notes on Point Estimator and Confidence Interval 1 By Hiro Kasahara
Poit Estimator Eco 325 Notes o Poit Estimator ad Cofidece Iterval 1 By Hiro Kasahara Parameter, Estimator, ad Estimate The ormal probability desity fuctio is fully characterized by two costats: populatio
More informationDouble Stage Shrinkage Estimator of Two Parameters. Generalized Exponential Distribution
Iteratioal Mathematical Forum, Vol., 3, o. 3, 3-53 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/.9/imf.3.335 Double Stage Shrikage Estimator of Two Parameters Geeralized Expoetial Distributio Alaa M.
More informationTrimmed Mean as an Adaptive Robust Estimator of a Location Parameter for Weibull Distribution
World Academy of Sciece Egieerig ad echology Iteratioal Joural of Mathematical ad Computatioal Scieces Vol: No:6 008 rimmed Mea as a Adaptive Robust Estimator of a Locatio Parameter for Weibull Distributio
More informationControl Charts for Mean for Non-Normally Correlated Data
Joural of Moder Applied Statistical Methods Volume 16 Issue 1 Article 5 5-1-017 Cotrol Charts for Mea for No-Normally Correlated Data J. R. Sigh Vikram Uiversity, Ujjai, Idia Ab Latif Dar School of Studies
More informationChapter 12 - Quality Cotrol Example: The process of llig 12 ouce cas of Dr. Pepper is beig moitored. The compay does ot wat to uderll the cas. Hece, a target llig rate of 12.1-12.5 ouces was established.
More informationElement sampling: Part 2
Chapter 4 Elemet samplig: Part 2 4.1 Itroductio We ow cosider uequal probability samplig desigs which is very popular i practice. I the uequal probability samplig, we ca improve the efficiecy of the resultig
More informationCEE 522 Autumn Uncertainty Concepts for Geotechnical Engineering
CEE 5 Autum 005 Ucertaity Cocepts for Geotechical Egieerig Basic Termiology Set A set is a collectio of (mutually exclusive) objects or evets. The sample space is the (collectively exhaustive) collectio
More informationOutput Analysis and Run-Length Control
IEOR E4703: Mote Carlo Simulatio Columbia Uiversity c 2017 by Marti Haugh Output Aalysis ad Ru-Legth Cotrol I these otes we describe how the Cetral Limit Theorem ca be used to costruct approximate (1 α%
More informationSome Properties of the Exact and Score Methods for Binomial Proportion and Sample Size Calculation
Some Properties of the Exact ad Score Methods for Biomial Proportio ad Sample Size Calculatio K. KRISHNAMOORTHY AND JIE PENG Departmet of Mathematics, Uiversity of Louisiaa at Lafayette Lafayette, LA 70504-1010,
More information( µ /σ)ζ/(ζ+1) µ /σ ( µ /σ)ζ/(ζ 1)
A eective CI for the mea with samples of size 1 ad Melaie Wall James Boe ad Richard Tweedie 1 Abstract It is couterituitive that with a sample of oly oe value from a ormal distributio oe ca costruct a
More informationSTAT 350 Handout 19 Sampling Distribution, Central Limit Theorem (6.6)
STAT 350 Hadout 9 Samplig Distributio, Cetral Limit Theorem (6.6) A radom sample is a sequece of radom variables X, X 2,, X that are idepedet ad idetically distributed. o This property is ofte abbreviated
More informationAn Introduction to Randomized Algorithms
A Itroductio to Radomized Algorithms The focus of this lecture is to study a radomized algorithm for quick sort, aalyze it usig probabilistic recurrece relatios, ad also provide more geeral tools for aalysis
More informationStatistical Theory MT 2009 Problems 1: Solution sketches
Statistical Theory MT 009 Problems : Solutio sketches. Which of the followig desities are withi a expoetial family? Explai your reasoig. (a) Let 0 < θ < ad put f(x, θ) = ( θ)θ x ; x = 0,,,... (b) (c) where
More informationChapter 6 Principles of Data Reduction
Chapter 6 for BST 695: Special Topics i Statistical Theory. Kui Zhag, 0 Chapter 6 Priciples of Data Reductio Sectio 6. Itroductio Goal: To summarize or reduce the data X, X,, X to get iformatio about a
More informationNon-life insurance mathematics. Nils F. Haavardsson, University of Oslo and DNB Skadeforsikring
No-life isurace mathematics Nils F. Haavardsso, Uiversity of Oslo ad DNB Skadeforsikrig Solvecy Fiacial cotrol of liabilities uder early worst-case scearios Target: the reserve which is the upper percetile
More information