the Presence of k Outliers
|
|
- Frederick West
- 6 years ago
- Views:
Transcription
1 RESEARCH ARTICLE OPEN ACCESS On the Estimation of the Presence of k Outliers for Weibull Distribution in Amal S. Hassan 1, Elsayed A. Elsherpieny 2, and Rania M. Shalaby 3 1 (Department of Mathematical Statistics, Institute of Statistical Studies& Research, Cairo University, Orman, Giza, Egypt) 2 (Department of Mathematical Statistics, Institute of Statistical Studies& Research, Cairo University, Orman, Giza, Egypt) 3 (The Higher Institute of Managerial Science, Culture and Science City, 6 th of October, Egypt) ABSTRACT This paper deals with the estimation problem of reliability where X, Y and Z are independently Weibull distribution with common known shape parameter but with different scale parameters, in presence of k outliers in strength X. The moment, maximum likelihood and mixture estimators of R are derived. An extensive computer simulation is used to compare the performance of the proposed estimators using MathCAD (14). Simulation study showed that the mixture estimators are better and easier than the maximum likelihood and moment estimators. Keywords-Maximum likelihood estimator, mixture estimator, moment estimator, outliers and Weibull distribution. I. Introduction In the context of reliability the stress-strength model describes the life of a component which has a random strength X and is subjected to random stress Y. This problem arises in the classical stress strength reliability where one is interested in estimating the proportion of the times the random strength X of a component exceeds the random stress Y to which the component is subjected. If X Y, then either the component fails or the system that uses the component may malfunction, and there is no failure when Y<X. The stress-strength models of the types P(Y<X), P(Y<X<Z) have extensive applications in various subareas of engineering, psychology, genetics, clinical trials and so on. [See,Kotz[1]]. The germ of this idea was introduced by Birnbaum [2] and developed by Birnbaum and McCarty [3]. An important particular case is estimation of R= P(Y<X<Z) which represents the situation where the strength X should not only be greater than stress Y but also be smaller than stress Z. For example, many devices cannot function at high temperatures; neither can do at very low ones. Similarly, person's blood pressure has two limits systolic and diastolic and his/her blood pressure should lie within these limits. Chandra and Owen [4] constructed maximum likelihood estimators (MLEs) and uniform minimum unbiased estimators (UMVUEs) for R= P(Y<X< Z). Singh [5] presented the minimum variance unbiased, maximum likelihood and empirical estimators of R= P(Y<X<Z), where X, Y and Z are mutually independent random variables and follow the normal distribution. Dutta and Sriwastav[6] dealt with the estimation of R when X, Y and Z are exponentially distributed. Ivshin[7] investigated the MLE and UMVUE of R when X, Y and Z are either uniform or exponential random variables with unknown location parameters. Wang et al. [8] make statistical inference for P(X<Y<Z) via two methods, the nonparametric normal approximation and the jackknife empirical likelihood, since the usual empirical likelihood method for U-statistics is too complicated to apply. The results of the simulation studies indicate that these two methods work promisingly compared to other existing methods. Some classical and real data sets were analyzed using these two proposed methods. The main aim of this article is to focus on the estimation of R= P(Y<X<Z), under the assumption that, X, Y and Z are independent. The stresses Y and Z have Weibull distribution with common known shape parameter β and scale parameters α and θ respectively. While the strength X has Weibull distribution with known shape parameter β and scale parameter γ in presence of k outliers. Maximum likelihood estimator, moment estimator (ME) and mixture estimator (Mix) are obtained. Monte Carlo simulation is performed for comparing different methods of estimation. The rest of the paper is organized as follows. In Section 2, the estimation of R= P(Y<X< Z) will be derived. The moment estimator of R derived in Section 3. Section 4 discussed MLE of R. The mixture estimator of R is obtained in Section 5. Monte Carlo simulation results are laid out in Section 6. Finally, conclusionsare presented in Section P a g e
2 II. Estimation of This Section deals with estimate the reliability of where Z, Y and X have Weibull distribution in the presence of k outliers in the strength X, such that X, Y and Z are independent. is the strength of independent observations such that k of them are distributed Weibully with scale parameter shape parameter and has the following probability density function (pdf), where is the cdf of Z at, is the cdf of Y at and is the survival function of Z at. Then, the reliability R in the presence of k outliers is given by while the remaining random variables are distributed Weibully with scale parameter known shape parameter and has the following pdf, According to Dixit [9] and Dixit and Nasiri [10], the joint distribution of in the presence of k outliers, can be expressed as The corresponding cumulative distribution function (cdf) Now to compute the estimate of R, the estimate of the parameters and must be obtained firstly. Three well known methods, namely, maximum likelihood, moments and mixture will be used to obtain the estimate of the parameters and will be discussed in details the following Sections. III. Moment Estimator of R In this Section, the moment estimator of R denoted by will be obtained. The moment estimators and of the unknown parameters and will be obtained by equating the population moments with the corresponding sample moments. The population means ofrandom stresses Y and Z are given by where is the stress of independent observations of Weibull distribution with scale parameter known shape parameter and has the following pdf, In addition, let is the stress of independent observations of random Z with known shape parameter and scale parameter and has the following pdf, According to Singh [5], the reliability takes the following form In addition the first and second population moments of strength X are given by, Suppose that be a random sample of size and be a random sample of size drawn from Weibull distribution with known shape parameter and scale parameter respectively. Then the means of these samples are given by be a random sample of size drawn from Weibull distribution with known shape parameter and scale parameter, then the first and second sample moments are given by,. By equating the samples moments with the corresponding population moments, then 1728 P a g e
3 and The moment estimator of and denoted by can be obtained from (8), respectively as, The moment estimator of denoted by can be obtained from (9) as follow IV. Maximum Likelihood Estimator of R This Section deals with MLE of reliability when X, Y and Z are independent Weibull distribution with parameters and respectively and the shape parameter is known. To compute the MLE for, firstly the MLEs, and must be obtained. The MLEs, and of the parameters and are the values which maximize the likelihood function. To obtain the maximum Likelihood estimator of, let be a random sample of size drawn from Weibull distribution with parameters the likelihood function of observed sample is given by of Thus it must be first obtain the moment estimator denoted by, by substitute (12) in (10), therefore The log-likelihood function denoted by is given by, and The maximum likelihood estimate of say, is obtained by setting the first partial derivatives of (15) to zero Therefore, if is non-negative then the roots are real. Therefore the moment estimator of denoted by, takes the following form, Hence the moment estimator can be obtained by substitute (13)in (12). Finally, the moment estimator of, denoted by is obtained by substitute and in (7), therefore takes the following form, By the similar way, to obtain the maximum Likelihood estimator of, let be a random sample of size drawn from Weibull distribution with parameters, then takes the following form To obtain the maximum Likelihood estimator and of, let be a random sample of size drawn from Weibull distribution with presence of k outliers with parameters the likelihood function of observed sample can be written as 1729 P a g e
4 where, The first partial derivatives of the log likelihood (18) with respect to are given, respectively by, The mixture estimator of, denoted by obtained by substitute, and in (7). is MLE s of denoted by and are solution to the system of equations obtained by setting the partial derivatives of the logarithm of likelihood function (19) and (20) to be zero. Obviously, it is not easy to obtain a closed form solution to this system of equations. Therefore, an iterative method must be applied to solve this equation numerically to estimate. The MLE of, denoted by is obtained by substitute in (7). V. Mixture Estimator of R To avoid the difficulty of complicated in the system of likelihood equations, the mixture estimator of R denoted by will be obtained. Following Read [11], the mixture estimator of and denoted by, and will be derived by mixing between moment estimates and MLEs which are obtained previously in Sections 3 and 4. The mixture estimator of can be obtained from moment estimator as follows The mixture estimator of can be obtained from likelihood estimators by substitute as follows VI. Numerical Illustration In this Section, an extensive numerical investigation will be carried out to compare the performance of the different estimators for different sample sizes and parameter values for Weibull distribution in the presence of k outliers. The investigated properties are biases and mean square errors (MSEs).All the computation is performed via MathCAD (14) statistical package. The algorithm for the parameter estimation can be summarized in the following steps: Step (1): Generate 1000 random samples, and from Weibull distribution with the sample sizes; (15,15,15),(20,20,20),(25,25,25), (15,15,20), (15,15,25), (15,15,20), (15,25,20), (15,25,25), (20,15,15), (20,15,20), (20,15,25), (20,25,20), (20,25,25), (25,15,15), (25,15,20), (25,15,25), (25,20,20), (25,20,25). Step (2): The parameter values are selected as Step (3): The moment estimator of and are obtained by solving (11). ME of is obtained numerically from (12). The ME of is obtained by substitute in (13). Once the estimate of these estimator are computed, then will be obtainedusing (7). Step (4): The MLE of and are obtained from (16) and (17). The nonlinear equations (19) and (20) of the MLEs will be solved iteratively using NetwonRaphson method. The MLE of will be obtained by substitute the MLEs of in (7). Step (5): The mixture estimator of and are computed by using (21). NetwonRaphson method is used for solving (22) and (23), to obtain the 1730 P a g e
5 estimators of he mixture estimator of will be obtained by substitute in (7). Step (6):The performance of the estimators can be evaluated through some measures of accuracy which are Biases and MSEs of Simulation results are summarized in Tables 1-4.From these Tables, the following conclusions can be observed on the properties of estimated parameters from R. 1. The estimated value of increases as the value of outliers, increases. 2. The estimated values of based on moment method is the smallest value, on the other hand the estimated values of based on mixture method is the highest one. 3. The biases of are the smallest relative to the biases of and. 4. Comparing the MSEs of all estimators, the mixture estimators perform the best estimator. 5. The biases and MSEs of based on different estimators increase as the value of outliers increases in almost all cases expect for some few cases. VII. Conclusions This article considered the problem of estimating the reliability for the Weibull distribution with presence of k outliers in strength X. Assuming that, X,Y and Z are independent with common known shape parameter but with different scale parameters. The moment, maximum likelihood and mixture estimators of are derived. Performance of estimators is usually evaluated through their biases and MSEs. Comparison study revealed that the mixture estimator works the best with respect to biases and MSEs, so the researcher strongly feels that mixture estimator is better and easy to calculate than the maximum likelihood and moment estimators. In general, the mixture method for estimating of the Weibull distribution in the presence of k outliers issuggested to be used. Table 1: Estimates of R, Biases and MSE s of the point estimates from Weibull Distribution, when (15,15,15) (20,20,20) (25,25,25) (15,15,20) (15,15,25) (15,20,20) (15,25,20) (15,25,25) (20,15,15) (20,15,20) (20,15,25) (20,25,20) (20,25,25) (25,15,15) (25,15,20) (25,15,25) (25,20,20) (25,20,25) P a g e
6 Table 2: Estimates of R, Biases and MSE s of the point estimates from Weibull Distribution, when (15,15,15) (20,20,20) (25,25,25) (15,15,20) (15,15,25) (15,20,20) (15,25,20) (15,25,25) (20,15,15) (20,15,20) (20,15,25) (20,25,20) (20,25,25) (25,15,15) (25,15,20) (25,15,25) (25,20,20) (25,20,25) Table 3: Estimates of R, Biases and MSE s of the point estimates from Weibull Distribution, when (15,15,15) (20,20,20) (25,25,25) (15,15,20) (15,15,25) (15,20,20) (15,25,20) (15,25,25) (20,15,15) (20,15,20) (20,15,25) (20,25,20) (20,25,25) (25,15,15) (25,15,20) (25,15,25) (25,20,20) (25,20,25) P a g e
7 Table 4: Estimates of R, Biases and MSE s of the point estimates from Weibull Distribution, when (15,15,15) (20,20,20) (25,25,25) (15,15,20) (15,15,25) (15,20,20) (15,25,20) (15,25,25) (20,15,15) (20,15,20) (20,15,25) (20,25,20) (20,25,25) (25,15,15) (25,15,20) (25,15,25) (25,20,20) (25,20,25) References [1] S.Kotz, Y. Lumelskii, and M. Pensky, The stress- strength model and its generalizations, Theory and Applications (World Scientific Publishing Co.,2003). [2] Z. W.Birnbaum, On a use of the Mann- Whitney statistic, proceedings of the third barkeley symposium on mathematical statistics and probability, I, 1956, [3] Z. W.Birnbaum, and R. C. McCarty,A distribution free upper confidence bound for Pr (Y<X), based on independent samples of X and Y, Annals of Mathematical Statistics, 29, 1958, [4] S. Chandraand D. B.Owen, On estimating the reliability of acomponent subject to several different stresses (strengths). Naval Research Logistics Quarterly,22,1975, [5] N. Singh, On the estimation of. Communication in Statistics Theory & Methods, 9,1980, [6] K. Dutta,and G. L.Sriwastav, An n-standby system with P(X < Y < Z). IAPQR Transaction, 12, 1986, [7] V. V.Ivshin, O`n the estimation of the probabilities of a doublelinear inequality in the case of uniform and two-parameter exponential distributions, Journal of Mathematical Science, 88, 1998, [8] Z.Wang, W.Xiping, and P.Guangming, Nonparametric statistical inference for P(X <Y <Z). The Indian Journal of Statistics, 75- A (1),2013, [9] U. J.Dixit, Estimation of parameters of the gamma distribution in the presence of outliers. Communications in statistics theory and methods, 18,1989, [10] U. J. Dixit, and P.F.Nasiri, Estimation of parameters of the exponential distribution in the presence of outliers generated from uniform distribution. Metron, 49(3-4),2001, [11] R. R.Read, Representation of certain covariance matrices with application to asymptotic efficiency. Journal of the American Statistical Association, 76, 1981, P a g e
Moments of the Reliability, R = P(Y<X), As a Random Variable
International Journal of Computational Engineering Research Vol, 03 Issue, 8 Moments of the Reliability, R = P(Y
More informationPoint and Interval Estimation for Gaussian Distribution, Based on Progressively Type-II Censored Samples
90 IEEE TRANSACTIONS ON RELIABILITY, VOL. 52, NO. 1, MARCH 2003 Point and Interval Estimation for Gaussian Distribution, Based on Progressively Type-II Censored Samples N. Balakrishnan, N. Kannan, C. T.
More informationEstimation of Stress-Strength Reliability for Kumaraswamy Exponential Distribution Based on Upper Record Values
International Journal of Contemporary Mathematical Sciences Vol. 12, 2017, no. 2, 59-71 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ijcms.2017.7210 Estimation of Stress-Strength Reliability for
More informationESTIMATOR IN BURR XII DISTRIBUTION
Journal of Reliability and Statistical Studies; ISSN (Print): 0974-804, (Online): 9-5666 Vol. 0, Issue (07): 7-6 ON THE VARIANCE OF P ( Y < X) ESTIMATOR IN BURR XII DISTRIBUTION M. Khorashadizadeh*, S.
More informationStep-Stress Models and Associated Inference
Department of Mathematics & Statistics Indian Institute of Technology Kanpur August 19, 2014 Outline Accelerated Life Test 1 Accelerated Life Test 2 3 4 5 6 7 Outline Accelerated Life Test 1 Accelerated
More informationTwo hours. To be supplied by the Examinations Office: Mathematical Formula Tables THE UNIVERSITY OF MANCHESTER. 21 June :45 11:45
Two hours MATH20802 To be supplied by the Examinations Office: Mathematical Formula Tables THE UNIVERSITY OF MANCHESTER STATISTICAL METHODS 21 June 2010 9:45 11:45 Answer any FOUR of the questions. University-approved
More informationEstimation of Parameters of the Weibull Distribution Based on Progressively Censored Data
International Mathematical Forum, 2, 2007, no. 41, 2031-2043 Estimation of Parameters of the Weibull Distribution Based on Progressively Censored Data K. S. Sultan 1 Department of Statistics Operations
More informationMathematical statistics
October 1 st, 2018 Lecture 11: Sufficient statistic Where are we? Week 1 Week 2 Week 4 Week 7 Week 10 Week 14 Probability reviews Chapter 6: Statistics and Sampling Distributions Chapter 7: Point Estimation
More informationSome Theoretical Properties and Parameter Estimation for the Two-Sided Length Biased Inverse Gaussian Distribution
Journal of Probability and Statistical Science 14(), 11-4, Aug 016 Some Theoretical Properties and Parameter Estimation for the Two-Sided Length Biased Inverse Gaussian Distribution Teerawat Simmachan
More informationOptimal Unbiased Estimates of P {X < Y } for Some Families of Distributions
Metodološki zvezki, Vol., No., 04, -9 Optimal Unbiased Estimates of P {X < } for Some Families of Distributions Marko Obradović, Milan Jovanović, Bojana Milošević Abstract In reliability theory, one of
More informationinferences on stress-strength reliability from lindley distributions
inferences on stress-strength reliability from lindley distributions D.K. Al-Mutairi, M.E. Ghitany & Debasis Kundu Abstract This paper deals with the estimation of the stress-strength parameter R = P (Y
More informationINVERTED KUMARASWAMY DISTRIBUTION: PROPERTIES AND ESTIMATION
Pak. J. Statist. 2017 Vol. 33(1), 37-61 INVERTED KUMARASWAMY DISTRIBUTION: PROPERTIES AND ESTIMATION A. M. Abd AL-Fattah, A.A. EL-Helbawy G.R. AL-Dayian Statistics Department, Faculty of Commerce, AL-Azhar
More informationAMCS243/CS243/EE243 Probability and Statistics. Fall Final Exam: Sunday Dec. 8, 3:00pm- 5:50pm VERSION A
AMCS243/CS243/EE243 Probability and Statistics Fall 2013 Final Exam: Sunday Dec. 8, 3:00pm- 5:50pm VERSION A *********************************************************** ID: ***********************************************************
More information0, otherwise. U = Y 1 Y 2 Hint: Use either the method of distribution functions or a bivariate transformation. (b) Find E(U).
1. Suppose Y U(0, 2) so that the probability density function (pdf) of Y is 1 2, 0 < y < 2 (a) Find the pdf of U = Y 4 + 1. Make sure to note the support. (c) Suppose Y 1, Y 2,..., Y n is an iid sample
More informationEstimation of Stress-Strength Reliability Using Record Ranked Set Sampling Scheme from the Exponential Distribution
Filomat 9:5 015, 1149 116 DOI 10.98/FIL1505149S Published by Faculty of Sciences and Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat Estimation of Stress-Strength eliability
More informationEstimation of P (X > Y ) for Weibull distribution based on hybrid censored samples
Estimation of P (X > Y ) for Weibull distribution based on hybrid censored samples A. Asgharzadeh a, M. Kazemi a, D. Kundu b a Department of Statistics, Faculty of Mathematical Sciences, University of
More informationReview. December 4 th, Review
December 4 th, 2017 Att. Final exam: Course evaluation Friday, 12/14/2018, 10:30am 12:30pm Gore Hall 115 Overview Week 2 Week 4 Week 7 Week 10 Week 12 Chapter 6: Statistics and Sampling Distributions Chapter
More informationInference on reliability in two-parameter exponential stress strength model
Metrika DOI 10.1007/s00184-006-0074-7 Inference on reliability in two-parameter exponential stress strength model K. Krishnamoorthy Shubhabrata Mukherjee Huizhen Guo Received: 19 January 2005 Springer-Verlag
More informationACTEX CAS EXAM 3 STUDY GUIDE FOR MATHEMATICAL STATISTICS
ACTEX CAS EXAM 3 STUDY GUIDE FOR MATHEMATICAL STATISTICS TABLE OF CONTENTS INTRODUCTORY NOTE NOTES AND PROBLEM SETS Section 1 - Point Estimation 1 Problem Set 1 15 Section 2 - Confidence Intervals and
More informationOptimum Test Plan for 3-Step, Step-Stress Accelerated Life Tests
International Journal of Performability Engineering, Vol., No., January 24, pp.3-4. RAMS Consultants Printed in India Optimum Test Plan for 3-Step, Step-Stress Accelerated Life Tests N. CHANDRA *, MASHROOR
More informationEstimation for generalized half logistic distribution based on records
Journal of the Korean Data & Information Science Society 202, 236, 249 257 http://dx.doi.org/0.7465/jkdi.202.23.6.249 한국데이터정보과학회지 Estimation for generalized half logistic distribution based on records
More informationProblem Selected Scores
Statistics Ph.D. Qualifying Exam: Part II November 20, 2010 Student Name: 1. Answer 8 out of 12 problems. Mark the problems you selected in the following table. Problem 1 2 3 4 5 6 7 8 9 10 11 12 Selected
More informationISI Web of Knowledge (Articles )
ISI Web of Knowledge (Articles 1 -- 18) Record 1 of 18 Title: Estimation and prediction from gamma distribution based on record values Author(s): Sultan, KS; Al-Dayian, GR; Mohammad, HH Source: COMPUTATIONAL
More informationEstimation of the Bivariate Generalized. Lomax Distribution Parameters. Based on Censored Samples
Int. J. Contemp. Math. Sciences, Vol. 9, 2014, no. 6, 257-267 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ijcms.2014.4329 Estimation of the Bivariate Generalized Lomax Distribution Parameters
More informationEstimation in an Exponentiated Half Logistic Distribution under Progressively Type-II Censoring
Communications of the Korean Statistical Society 2011, Vol. 18, No. 5, 657 666 DOI: http://dx.doi.org/10.5351/ckss.2011.18.5.657 Estimation in an Exponentiated Half Logistic Distribution under Progressively
More informationOn the Comparison of Fisher Information of the Weibull and GE Distributions
On the Comparison of Fisher Information of the Weibull and GE Distributions Rameshwar D. Gupta Debasis Kundu Abstract In this paper we consider the Fisher information matrices of the generalized exponential
More informationExponentiated Rayleigh Distribution: A Bayes Study Using MCMC Approach Based on Unified Hybrid Censored Data
Exponentiated Rayleigh Distribution: A Bayes Study Using MCMC Approach Based on Unified Hybrid Censored Data ABSTRACT M. G. M. Ghazal 1, H. M. Hasaballah 2 1 Mathematics Department, Faculty of Science,
More informationThe comparative studies on reliability for Rayleigh models
Journal of the Korean Data & Information Science Society 018, 9, 533 545 http://dx.doi.org/10.7465/jkdi.018.9..533 한국데이터정보과학회지 The comparative studies on reliability for Rayleigh models Ji Eun Oh 1 Joong
More informationBayesian Analysis of Simple Step-stress Model under Weibull Lifetimes
Bayesian Analysis of Simple Step-stress Model under Weibull Lifetimes A. Ganguly 1, D. Kundu 2,3, S. Mitra 2 Abstract Step-stress model is becoming quite popular in recent times for analyzing lifetime
More informationReview Quiz. 1. Prove that in a one-dimensional canonical exponential family, the complete and sufficient statistic achieves the
Review Quiz 1. Prove that in a one-dimensional canonical exponential family, the complete and sufficient statistic achieves the Cramér Rao lower bound (CRLB). That is, if where { } and are scalars, then
More informationTABLE OF CONTENTS CHAPTER 1 COMBINATORIAL PROBABILITY 1
TABLE OF CONTENTS CHAPTER 1 COMBINATORIAL PROBABILITY 1 1.1 The Probability Model...1 1.2 Finite Discrete Models with Equally Likely Outcomes...5 1.2.1 Tree Diagrams...6 1.2.2 The Multiplication Principle...8
More informationHybrid Censoring; An Introduction 2
Hybrid Censoring; An Introduction 2 Debasis Kundu Department of Mathematics & Statistics Indian Institute of Technology Kanpur 23-rd November, 2010 2 This is a joint work with N. Balakrishnan Debasis Kundu
More informationThe Distributions of Sums, Products and Ratios of Inverted Bivariate Beta Distribution 1
Applied Mathematical Sciences, Vol. 2, 28, no. 48, 2377-2391 The Distributions of Sums, Products and Ratios of Inverted Bivariate Beta Distribution 1 A. S. Al-Ruzaiza and Awad El-Gohary 2 Department of
More informationMathematical statistics
October 4 th, 2018 Lecture 12: Information Where are we? Week 1 Week 2 Week 4 Week 7 Week 10 Week 14 Probability reviews Chapter 6: Statistics and Sampling Distributions Chapter 7: Point Estimation Chapter
More informationSimultaneous Prediction Intervals for the (Log)- Location-Scale Family of Distributions
Statistics Preprints Statistics 10-2014 Simultaneous Prediction Intervals for the (Log)- Location-Scale Family of Distributions Yimeng Xie Virginia Tech Yili Hong Virginia Tech Luis A. Escobar Louisiana
More informationEstimation theory. Parametric estimation. Properties of estimators. Minimum variance estimator. Cramer-Rao bound. Maximum likelihood estimators
Estimation theory Parametric estimation Properties of estimators Minimum variance estimator Cramer-Rao bound Maximum likelihood estimators Confidence intervals Bayesian estimation 1 Random Variables Let
More informationABC methods for phase-type distributions with applications in insurance risk problems
ABC methods for phase-type with applications problems Concepcion Ausin, Department of Statistics, Universidad Carlos III de Madrid Joint work with: Pedro Galeano, Universidad Carlos III de Madrid Simon
More informationReview and continuation from last week Properties of MLEs
Review and continuation from last week Properties of MLEs As we have mentioned, MLEs have a nice intuitive property, and as we have seen, they have a certain equivariance property. We will see later that
More informationComparison of Least Square Estimators with Rank Regression Estimators of Weibull Distribution - A Simulation Study
ISSN 168-80 Journal of Statistics Volume 0, 01. pp. 1-10 Comparison of Least Square Estimators with Rank Regression Estimators of Weibull Distribution - A Simulation Study Abstract Chandika Rama Mohan
More informationStatistics - Lecture One. Outline. Charlotte Wickham 1. Basic ideas about estimation
Statistics - Lecture One Charlotte Wickham wickham@stat.berkeley.edu http://www.stat.berkeley.edu/~wickham/ Outline 1. Basic ideas about estimation 2. Method of Moments 3. Maximum Likelihood 4. Confidence
More informationMathematical statistics
October 18 th, 2018 Lecture 16: Midterm review Countdown to mid-term exam: 7 days Week 1 Chapter 1: Probability review Week 2 Week 4 Week 7 Chapter 6: Statistics Chapter 7: Point Estimation Chapter 8:
More informationAnswers and expectations
Answers and expectations For a function f(x) and distribution P(x), the expectation of f with respect to P is The expectation is the average of f, when x is drawn from the probability distribution P E
More informationConstant Stress Partially Accelerated Life Test Design for Inverted Weibull Distribution with Type-I Censoring
Algorithms Research 013, (): 43-49 DOI: 10.593/j.algorithms.01300.0 Constant Stress Partially Accelerated Life Test Design for Mustafa Kamal *, Shazia Zarrin, Arif-Ul-Islam Department of Statistics & Operations
More informationMS&E 226: Small Data. Lecture 11: Maximum likelihood (v2) Ramesh Johari
MS&E 226: Small Data Lecture 11: Maximum likelihood (v2) Ramesh Johari ramesh.johari@stanford.edu 1 / 18 The likelihood function 2 / 18 Estimating the parameter This lecture develops the methodology behind
More informationSTATISTICAL INFERENCE IN ACCELERATED LIFE TESTING WITH GEOMETRIC PROCESS MODEL. A Thesis. Presented to the. Faculty of. San Diego State University
STATISTICAL INFERENCE IN ACCELERATED LIFE TESTING WITH GEOMETRIC PROCESS MODEL A Thesis Presented to the Faculty of San Diego State University In Partial Fulfillment of the Requirements for the Degree
More informationInferences on stress-strength reliability from weighted Lindley distributions
Inferences on stress-strength reliability from weighted Lindley distributions D.K. Al-Mutairi, M.E. Ghitany & Debasis Kundu Abstract This paper deals with the estimation of the stress-strength parameter
More informationStatistical Inference Using Progressively Type-II Censored Data with Random Scheme
International Mathematical Forum, 3, 28, no. 35, 1713-1725 Statistical Inference Using Progressively Type-II Censored Data with Random Scheme Ammar M. Sarhan 1 and A. Abuammoh Department of Statistics
More informationChapter 9 Regression with a Binary Dependent Variable. Multiple Choice. 1) The binary dependent variable model is an example of a
Chapter 9 Regression with a Binary Dependent Variable Multiple Choice ) The binary dependent variable model is an example of a a. regression model, which has as a regressor, among others, a binary variable.
More informationI I FINAL, 01 Jun 8.4 to 31 May TITLE AND SUBTITLE 5 * _- N, '. ', -;
R AD-A237 850 E........ I N 11111IIIII U 1 1I!til II II... 1. AGENCY USE ONLY Leave 'VanK) I2. REPORT DATE 3 REPORT TYPE AND " - - I I FINAL, 01 Jun 8.4 to 31 May 88 4. TITLE AND SUBTITLE 5 * _- N, '.
More informationINFERENCE FOR MULTIPLE LINEAR REGRESSION MODEL WITH EXTENDED SKEW NORMAL ERRORS
Pak. J. Statist. 2016 Vol. 32(2), 81-96 INFERENCE FOR MULTIPLE LINEAR REGRESSION MODEL WITH EXTENDED SKEW NORMAL ERRORS A.A. Alhamide 1, K. Ibrahim 1 M.T. Alodat 2 1 Statistics Program, School of Mathematical
More informationBayesian and Non Bayesian Estimations for. Birnbaum-Saunders Distribution under Accelerated. Life Testing Based oncensoring sampling
Applied Mathematical Sciences, Vol. 7, 2013, no. 66, 3255-3269 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2013.34232 Bayesian and Non Bayesian Estimations for Birnbaum-Saunders Distribution
More informationAnalysis of Type-II Progressively Hybrid Censored Data
Analysis of Type-II Progressively Hybrid Censored Data Debasis Kundu & Avijit Joarder Abstract The mixture of Type-I and Type-II censoring schemes, called the hybrid censoring scheme is quite common in
More informationBTRY 4090: Spring 2009 Theory of Statistics
BTRY 4090: Spring 2009 Theory of Statistics Guozhang Wang September 25, 2010 1 Review of Probability We begin with a real example of using probability to solve computationally intensive (or infeasible)
More informationA comparison of inverse transform and composition methods of data simulation from the Lindley distribution
Communications for Statistical Applications and Methods 2016, Vol. 23, No. 6, 517 529 http://dx.doi.org/10.5351/csam.2016.23.6.517 Print ISSN 2287-7843 / Online ISSN 2383-4757 A comparison of inverse transform
More informationCOPYRIGHTED MATERIAL CONTENTS. Preface Preface to the First Edition
Preface Preface to the First Edition xi xiii 1 Basic Probability Theory 1 1.1 Introduction 1 1.2 Sample Spaces and Events 3 1.3 The Axioms of Probability 7 1.4 Finite Sample Spaces and Combinatorics 15
More informationMonte Carlo Studies. The response in a Monte Carlo study is a random variable.
Monte Carlo Studies The response in a Monte Carlo study is a random variable. The response in a Monte Carlo study has a variance that comes from the variance of the stochastic elements in the data-generating
More informationBayesian Estimation for the Generalized Logistic Distribution Type-II Censored Accelerated Life Testing
Int. J. Contemp. Math. Sciences, Vol. 8, 2013, no. 20, 969-986 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ijcms.2013.39111 Bayesian Estimation for the Generalized Logistic Distribution Type-II
More information1. Point Estimators, Review
AMS571 Prof. Wei Zhu 1. Point Estimators, Review Example 1. Let be a random sample from. Please find a good point estimator for Solutions. There are the typical estimators for and. Both are unbiased estimators.
More informationDistribution Fitting (Censored Data)
Distribution Fitting (Censored Data) Summary... 1 Data Input... 2 Analysis Summary... 3 Analysis Options... 4 Goodness-of-Fit Tests... 6 Frequency Histogram... 8 Comparison of Alternative Distributions...
More informationMultivariate Survival Analysis
Multivariate Survival Analysis Previously we have assumed that either (X i, δ i ) or (X i, δ i, Z i ), i = 1,..., n, are i.i.d.. This may not always be the case. Multivariate survival data can arise in
More informationPENALIZED LIKELIHOOD PARAMETER ESTIMATION FOR ADDITIVE HAZARD MODELS WITH INTERVAL CENSORED DATA
PENALIZED LIKELIHOOD PARAMETER ESTIMATION FOR ADDITIVE HAZARD MODELS WITH INTERVAL CENSORED DATA Kasun Rathnayake ; A/Prof Jun Ma Department of Statistics Faculty of Science and Engineering Macquarie University
More informationA General Overview of Parametric Estimation and Inference Techniques.
A General Overview of Parametric Estimation and Inference Techniques. Moulinath Banerjee University of Michigan September 11, 2012 The object of statistical inference is to glean information about an underlying
More informationClosed-Form Estimators for the Gamma Distribution Derived from Likelihood Equations
The American Statistician ISSN: 3-135 (Print) 1537-2731 (Online) Journal homepage: http://www.tandfonline.com/loi/utas2 Closed-Form Estimators for the Gamma Distribution Derived from Likelihood Equations
More informationSTATISTICS 4, S4 (4769) A2
(4769) A2 Objectives To provide students with the opportunity to explore ideas in more advanced statistics to a greater depth. Assessment Examination (72 marks) 1 hour 30 minutes There are four options
More informationTwo hours. To be supplied by the Examinations Office: Mathematical Formula Tables and Statistical Tables THE UNIVERSITY OF MANCHESTER.
Two hours MATH38181 To be supplied by the Examinations Office: Mathematical Formula Tables and Statistical Tables THE UNIVERSITY OF MANCHESTER EXTREME VALUES AND FINANCIAL RISK Examiner: Answer any FOUR
More informationPractice Problems Section Problems
Practice Problems Section 4-4-3 4-4 4-5 4-6 4-7 4-8 4-10 Supplemental Problems 4-1 to 4-9 4-13, 14, 15, 17, 19, 0 4-3, 34, 36, 38 4-47, 49, 5, 54, 55 4-59, 60, 63 4-66, 68, 69, 70, 74 4-79, 81, 84 4-85,
More informationStat 5101 Lecture Notes
Stat 5101 Lecture Notes Charles J. Geyer Copyright 1998, 1999, 2000, 2001 by Charles J. Geyer May 7, 2001 ii Stat 5101 (Geyer) Course Notes Contents 1 Random Variables and Change of Variables 1 1.1 Random
More informationEstimation for Mean and Standard Deviation of Normal Distribution under Type II Censoring
Communications for Statistical Applications and Methods 2014, Vol. 21, No. 6, 529 538 DOI: http://dx.doi.org/10.5351/csam.2014.21.6.529 Print ISSN 2287-7843 / Online ISSN 2383-4757 Estimation for Mean
More informationDistributions of Functions of Random Variables. 5.1 Functions of One Random Variable
Distributions of Functions of Random Variables 5.1 Functions of One Random Variable 5.2 Transformations of Two Random Variables 5.3 Several Random Variables 5.4 The Moment-Generating Function Technique
More informationEstimation of the Exponential Distribution based on Multiply Progressive Type II Censored Sample
Communications of the Korean Statistical Society 2012 Vol. 19 No. 5 697 70 DOI: http://dx.doi.org/10.5351/ckss.2012.19.5.697 Estimation of the Exponential Distribution based on Multiply Progressive Type
More informationThis midterm covers Chapters 6 and 7 in WMS (and the notes). The following problems are stratified by chapter.
This midterm covers Chapters 6 and 7 in WMS (and the notes). The following problems are stratified by chapter. Chapter 6 Problems 1. Suppose that Y U(0, 2) so that the probability density function (pdf)
More informationThe Kumaraswamy-Burr Type III Distribution: Properties and Estimation
British Journal of Mathematics & Computer Science 14(2): 1-21, 2016, Article no.bjmcs.19958 ISSN: 2231-0851 SCIENCEDOMAIN international www.sciencedomain.org The Kumaraswamy-Burr Type III Distribution:
More informationECE 275A Homework 7 Solutions
ECE 275A Homework 7 Solutions Solutions 1. For the same specification as in Homework Problem 6.11 we want to determine an estimator for θ using the Method of Moments (MOM). In general, the MOM estimator
More informationParameters Estimation for a Linear Exponential Distribution Based on Grouped Data
International Mathematical Forum, 3, 2008, no. 33, 1643-1654 Parameters Estimation for a Linear Exponential Distribution Based on Grouped Data A. Al-khedhairi Department of Statistics and O.R. Faculty
More informationSurvival Analysis. Lu Tian and Richard Olshen Stanford University
1 Survival Analysis Lu Tian and Richard Olshen Stanford University 2 Survival Time/ Failure Time/Event Time We will introduce various statistical methods for analyzing survival outcomes What is the survival
More informationStock Sampling with Interval-Censored Elapsed Duration: A Monte Carlo Analysis
Stock Sampling with Interval-Censored Elapsed Duration: A Monte Carlo Analysis Michael P. Babington and Javier Cano-Urbina August 31, 2018 Abstract Duration data obtained from a given stock of individuals
More informationFall 2017 STAT 532 Homework Peter Hoff. 1. Let P be a probability measure on a collection of sets A.
1. Let P be a probability measure on a collection of sets A. (a) For each n N, let H n be a set in A such that H n H n+1. Show that P (H n ) monotonically converges to P ( k=1 H k) as n. (b) For each n
More informationIntroduction to Maximum Likelihood Estimation
Introduction to Maximum Likelihood Estimation Eric Zivot July 26, 2012 The Likelihood Function Let 1 be an iid sample with pdf ( ; ) where is a ( 1) vector of parameters that characterize ( ; ) Example:
More informationMaximum Likelihood Estimation. only training data is available to design a classifier
Introduction to Pattern Recognition [ Part 5 ] Mahdi Vasighi Introduction Bayesian Decision Theory shows that we could design an optimal classifier if we knew: P( i ) : priors p(x i ) : class-conditional
More informationProbability and Estimation. Alan Moses
Probability and Estimation Alan Moses Random variables and probability A random variable is like a variable in algebra (e.g., y=e x ), but where at least part of the variability is taken to be stochastic.
More informationParameter estimation and forecasting. Cristiano Porciani AIfA, Uni-Bonn
Parameter estimation and forecasting Cristiano Porciani AIfA, Uni-Bonn Questions? C. Porciani Estimation & forecasting 2 Temperature fluctuations Variance at multipole l (angle ~180o/l) C. Porciani Estimation
More informationThree hours. To be supplied by the Examinations Office: Mathematical Formula Tables and Statistical Tables THE UNIVERSITY OF MANCHESTER.
Three hours To be supplied by the Examinations Office: Mathematical Formula Tables and Statistical Tables THE UNIVERSITY OF MANCHESTER EXTREME VALUES AND FINANCIAL RISK Examiner: Answer QUESTION 1, QUESTION
More informationImproving Efficiency of Inferences in Randomized Clinical Trials Using Auxiliary Covariates
Improving Efficiency of Inferences in Randomized Clinical Trials Using Auxiliary Covariates Anastasios (Butch) Tsiatis Department of Statistics North Carolina State University http://www.stat.ncsu.edu/
More informationMATH4427 Notebook 2 Fall Semester 2017/2018
MATH4427 Notebook 2 Fall Semester 2017/2018 prepared by Professor Jenny Baglivo c Copyright 2009-2018 by Jenny A. Baglivo. All Rights Reserved. 2 MATH4427 Notebook 2 3 2.1 Definitions and Examples...................................
More informationSTA 2201/442 Assignment 2
STA 2201/442 Assignment 2 1. This is about how to simulate from a continuous univariate distribution. Let the random variable X have a continuous distribution with density f X (x) and cumulative distribution
More informationON THE FAILURE RATE ESTIMATION OF THE INVERSE GAUSSIAN DISTRIBUTION
ON THE FAILURE RATE ESTIMATION OF THE INVERSE GAUSSIAN DISTRIBUTION ZHENLINYANGandRONNIET.C.LEE Department of Statistics and Applied Probability, National University of Singapore, 3 Science Drive 2, Singapore
More informationOn Sarhan-Balakrishnan Bivariate Distribution
J. Stat. Appl. Pro. 1, No. 3, 163-17 (212) 163 Journal of Statistics Applications & Probability An International Journal c 212 NSP On Sarhan-Balakrishnan Bivariate Distribution D. Kundu 1, A. Sarhan 2
More informationMax. Likelihood Estimation. Outline. Econometrics II. Ricardo Mora. Notes. Notes
Maximum Likelihood Estimation Econometrics II Department of Economics Universidad Carlos III de Madrid Máster Universitario en Desarrollo y Crecimiento Económico Outline 1 3 4 General Approaches to Parameter
More informationDouble Gamma Principal Components Analysis
Applied Mathematical Sciences, Vol. 12, 2018, no. 11, 523-533 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ams.2018.8455 Double Gamma Principal Components Analysis Ameerah O. Bahashwan, Zakiah
More informationRobust Monte Carlo Methods for Sequential Planning and Decision Making
Robust Monte Carlo Methods for Sequential Planning and Decision Making Sue Zheng, Jason Pacheco, & John Fisher Sensing, Learning, & Inference Group Computer Science & Artificial Intelligence Laboratory
More informationLogistic-Modified Weibull Distribution and Parameter Estimation
International Journal of Contemporary Mathematical Sciences Vol. 13, 2018, no. 1, 11-23 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ijcms.2018.71135 Logistic-Modified Weibull Distribution and
More informationDefine characteristic function. State its properties. State and prove inversion theorem.
ASSIGNMENT - 1, MAY 013. Paper I PROBABILITY AND DISTRIBUTION THEORY (DMSTT 01) 1. (a) Give the Kolmogorov definition of probability. State and prove Borel cantelli lemma. Define : (i) distribution function
More informationStatistical Inference on Constant Stress Accelerated Life Tests Under Generalized Gamma Lifetime Distributions
Int. Statistical Inst.: Proc. 58th World Statistical Congress, 2011, Dublin (Session CPS040) p.4828 Statistical Inference on Constant Stress Accelerated Life Tests Under Generalized Gamma Lifetime Distributions
More informationHybrid Censoring Scheme: An Introduction
Department of Mathematics & Statistics Indian Institute of Technology Kanpur August 19, 2014 Outline 1 2 3 4 5 Outline 1 2 3 4 5 What is? Lifetime data analysis is used to analyze data in which the time
More informationQuality Control Using Inferential Statistics In Weibull Based Reliability Analyses S. F. Duffy 1 and A. Parikh 2
Quality Control Using Inferential Statistics In Weibull Based Reliability Analyses S. F. Duffy 1 and A. Parikh 2 1 Cleveland State University 2 N & R Engineering www.inl.gov ASTM Symposium on Graphite
More informationA New Class of Positively Quadrant Dependent Bivariate Distributions with Pareto
International Mathematical Forum, 2, 27, no. 26, 1259-1273 A New Class of Positively Quadrant Dependent Bivariate Distributions with Pareto A. S. Al-Ruzaiza and Awad El-Gohary 1 Department of Statistics
More informationReliability Engineering I
Happiness is taking the reliability final exam. Reliability Engineering I ENM/MSC 565 Review for the Final Exam Vital Statistics What R&M concepts covered in the course When Monday April 29 from 4:30 6:00
More informationParametric Models. Dr. Shuang LIANG. School of Software Engineering TongJi University Fall, 2012
Parametric Models Dr. Shuang LIANG School of Software Engineering TongJi University Fall, 2012 Today s Topics Maximum Likelihood Estimation Bayesian Density Estimation Today s Topics Maximum Likelihood
More informationTHE WEIBULL GENERALIZED FLEXIBLE WEIBULL EXTENSION DISTRIBUTION
Journal of Data Science 14(2016), 453-478 THE WEIBULL GENERALIZED FLEXIBLE WEIBULL EXTENSION DISTRIBUTION Abdelfattah Mustafa, Beih S. El-Desouky, Shamsan AL-Garash Department of Mathematics, Faculty of
More information