Single Stage Shrinkage Estimator for the Shape Parameter of the Pareto Distribution

Size: px
Start display at page:

Download "Single Stage Shrinkage Estimator for the Shape Parameter of the Pareto Distribution"

Transcription

1 International Journal of Mathematics Trends and Technology IJMTT - Volume 35 Number 3- July 20 Single Stage Shrinkage Estimator for the Shape Parameter of the Pareto Distribution Abbas Najim Salman*, Adel Abdulkadhim Hussein** Department of Mathematics-Ibn-Al-Haitham College of Education - University of Baghdad Abbasnajim@yahoocom* adilabed57@yahoocom** Abstract: In this paper, we employee single stage shrinkage estimator for estimate the shape parameter of the Pareto distribution when the scale parameter is known The proposed estimator is shown to have a smaller mean squared error in a region around when comparison with the usual and existing estimators Keywords: Pareto Distribution, Shrinkage Estimation, Prior Guess, Pre Test Region, Bias, Mean Squared Error and Relative Efficiency Introduction The Model "Pareto distribution is a power law probability distribution that coincides with social, scientific, geophysical and many other types of observable phenomena It is a continuous distribution bounded on the lower side and it has two parameters, shape and mode and it is a highly skewed distribution It is a decreasing function and it has a finite value at the minimum value It is a heavy tailed distribution meaning that a random variable following a Pareto distribution can have extreme values The mode parameter for Pareto distribution sets the position of the left edge of the probability density function The only outcomes that can be observed from this distribution are greater than or equal to the value of the mode parameter Changes in the mode simply shift the boundary to the left or right Pareto distribution named after the Italian economist and sociologist Vilfredo Pareto -923 It was proposed first by Pareto in 97 at the university of Lausanne as a model for the distribution of income The pareto model is very often used as a basis of Excess of loss quotations as it gives a pretty good description of the random behavior of large losses, and often used to model the distribution of income ; [], [7],[],[2] and [] Consider a random sample t, t2,, tn from the Pareto distribution with the following pdf ISSN: Page 5

2 International Journal of Mathematics Trends and Technology IJMTT - Volume 35 Number 3- July 20 <0 =0 >0 0 ft t Figure Shown pdf of Pareto distribution In conventional notation, we write t par,, where and are the shape and scale parameter respectively Sometimes we may have a prior guess value due past experiences or from an acquaintance with similar situation of the parameter to be estimated, if this value is very close to the true value, the shrinkage technique is useful to obtain an improved estimator According to Thompson [3], such prior estimate may arise for any one of a number of reasons, eg, we are estimating and i we believe is close to the true value of ; or ii We fear that may be near the true value of ; ie; A bad thing can happen when you do not we use it if = For such cases, this prior guess value may be utilized to improve the estimation procedure with usual estimator MLE via shrinkage estimator which including shrinkage weight factor as follows, see[5] = 2 The authors in [9],[3],[5],[] and others suggested shrinkage estimator 2 for estimating the parameters of different distribution when a guess prior value is available They showed that these estimators perform better than usual estimator in the term of mean square error when the guess value is close to the true value Consequently, to make sure whether the prior guess value is approximately or close to the true value or not, we may test H0: = vs HA: and a preliminary test shrunken estimator of significance was employed as follows: ISSN: / / 2 3 Page 57

3 International Journal of Mathematics Trends and Technology IJMTT - Volume 35 Number 3- July 20 Using pre-assigned level of significance, the critical region R for such test is pretest region using specific test statistic, where, i =, 2 is shrinkage weight factors such that 0 and is the usual estimator MLE for example Preliminary test shrinkage estimators 3 have been considered in different contexts by [],[2],[3], [0], [], [5] The aim of this paper is to estimate the shape parameter of pareto distribution using mentioned preliminary test shrinkage estimator 3 and study its behavior when we derived the expressions of its Bias, Mean Square Error and Relative Efficiency and study their performance Numerical results and conclusions for these expressions were made to show the effective of the proposed estimator as well as some comparisons with the usual and existing estimators were made 2 Usual Estimation of the Shape Parameter Let t, t2,, tn be identically independent pareto distributed random variables when the scale parameter is known, then the maximum likelihood estimator of as below: See[2] It follows easily that Ln t will be exponentially distributed with mean / y= will be Gamma distributed with parameters n and ie; y Gn, Therefore, when E Then, we get and var and the mean square error is 2 Single Stage Shrunken Estimator Recall the estimator 3 when Where R is the pre test region for test H0: ISSN: = vs HA: 3 using test statistic Page 5

4 International Journal of Mathematics Trends and Technology IJMTT - Volume 35 Number 3- July 20, ie; - For simplest we denote to R as below R=[a,b], where and are the lower and upper 00 /2 percentile of chi square distribution with degree of freedom 2n The Expression for Bias of the proposed estimator is defined as = [ ] where R is the complement region of R in real space and f is a pdf of which has the following form: [ ] 5 Then the Bias will be *, - +, we denote to the Bias Ratio as B which is defined as: 7 The expression of the Mean Squared Error MSE of, *,, The Efficiency of ISSN: is, - relative to the, - denoted by REff + R is defined as below: Page 59

5 International Journal of Mathematics Trends and Technology IJMTT - Volume 35 Number 3- July 20 9 see [],[2],[3] 3 Numerical Results and Discussion The computations of the statistical indicators Relative Efficiency [REff ] and Bias Ratio [B ] expressions were used for the considered testimators These computations were performed for the constants = 00,005,0, n =,,,0,2 and = and k=exp-n Some of these computations are displayed in attached table for some samples of these constants The observation mentioned in the table leads to the following results: i The Relative Efficiency [REff ] of are adversely proportional with small value of especially when =, ie = 00 yield highest efficiency iithe Relative Efficiency [REff ] of has maximum value when = =, for each n,, and decreasing otherwise This feature shown the important usefulness of prior knowledge which given higher effects of proposed estimator as well as the important role of shrinkage technique iiibias ratio [B ] of increases when increases ivbias ratio [B ] of are reasonably small when = for each n,, and increases otherwise This property shown that the proposed estimator is very close to unbiased property especially when = vthe Relative efficiency [REff ] of increases function with increases value of n for =, and vice versa otherwise for each vithe Effective Interval [the value of that makes REff greater than one] using proposed estimator is [075,25] Here the pretest criterion is very important for guarantee that prior information is very closely to the actual value and prevent it far away from it, which get optimal effect of the considered estimator to obtain high efficiency viithe considered estimator is better than the classical estimator especially when, which is given the effective of and important weight of prior knowledge as well as the increment of efficiency may be reach to tens time viiithe proposed estimator has smaller MSE than some existing estimators introduced by authors, see for examples [] Conclusions From the above discussions it is obvious that by using guess point value one can improve the usual estimator It can be noted that if the guess point is very close to the true value of the parameter ie; =, the proposed estimators perform better than the usual estimator If one has no confidence in the guessed value then proposed preliminary test shrinkage estimators can be suggested We can safely use the proposed estimators for small sample size at usual level of significance and moderate value of shrinkage weight factors ISSN: Page 0

6 International Journal of Mathematics Trends and Technology IJMTT - Volume 35 Number 3- July 20 References Al-Hemyari,ZA, Khurshid, A and Al-Joboori, A N,2009, On Thompson Type Estimators for the Mean of Normal Distribution, REVISTA INVESTIGACION OPERACIONAL J Vol 30, No2,09- Al-Joboori, AN, et al 20, Single and Double Stage Shrinkage Estimators for the normal Mean with the Variance Cases, International Journal of statistic, Vol3,2,PP, 27-3 Al-Joboori, AN, Khalef, BA and Hayder SK, 20, Improved of Stein-Type Estimator of the Variance of Normal Distribution via Shrinkage Estimation Technique, International Journal of Inventive Engineering and Science IJIES, Vol2,, PP-7 Al-Joboori, AN,AA Hussein and MAMohammd, 20,On Significance Testimator in Pareto Distribution Via Shrinkage Technique, Journal of College of E Casella, G, 2002, Shrinkage, In: Encyclopedia of Envirnmetrics,, EdsAHELShaarawi and wwpiegorsch, John Wiley :New York Cirillo R, "The Economics of Vilfredo Pareto", Frank Case and Company Limited, London, 979 Malik, HJ, 970, Estimation of the Parameter of the Pareto Distribution, Metrika, Vol5, pp 2-32 Ord, JK, 975, Statistical Models for Personal income Distributions in Statistical 5, Dordrecht, Reidel Mehta,JS and Srinivasan,R, 97, "Estimation of the Mean by Shrinkage to a Point", J Amer Statist Assoc,, p-90 Saleh, AKE, 200, Theory of Preliminary Test and Stein-Type Estimators with Application, Wiley and Sons, NewYork Rytgaard, M, 990, Estimation in the Pareto Distribution, Astin Bulletin, Vol20, No2, pp20-2 Rashid, H A and Al-Gazi, N A A, 20, Bayes Estimators for the Shape Parameter of Pareto Distribution under Generalized Square Errors Loss Function, Mathematical Theory and Modeling, Vol,No pp Thompson, JR, 9, Some Shrinkage Techniques for Estimating the Mean, J Amer Statist Assoc, Vol3, pp3-22 Prakash, G, Singh, DC and Singh, RD, 200, Some Test Estimator for the Scale Parameter of Classical Pareto Distribution, Journal of Statistical Research, Vol02, pp-5 Singh, DC, Prakash, G and Singh, P, 2007, Shrinkage Testimators of the Shape Parameter of Pareto Distribution Using LINEX Loss Function, Communication in Statistics-Theory and Methods, Vol3, pp7-753 Singh, DC, Singh, P and Singh, PR, 99, Shrunken Estimators for the Scale Parameter of Classical Pareto Distribution, Microelectron Reliability, Vol33, pp3539 ISSN: Page

7 International Journal of Mathematics Trends and Technology IJMTT - Volume 35 Number 3- July 20 Table : Shown Bais Ratio B and Relative Efficiency REff of estimator Δ n REff B - REff 07 B - REff 0232 B - REff B REff 0222 B REff B - REff 07 B - REff B - REff B - REff 0222 B REff B - REff 07 B - REff 0232 B - REff B - REff 0222 B ISSN: e+09 57e e e-0 599e e e e e e-02 77e+05 75e e+020 5e-0 52e e-0 955e e e e e e e+0 735e e e e e e e Page 2

Two Different Shrinkage Estimator Classes for the Shape Parameter of Classical Pareto Distribution

Two Different Shrinkage Estimator Classes for the Shape Parameter of Classical Pareto Distribution Two Different Shrinkage Estimator Classes for the Shape Parameter of Classical Pareto Distribution Meral EBEGIL and Senay OZDEMIR Abstract In this study, biased estimators for the shape parameter of a

More information

SHRINKAGE ESTIMATOR AND TESTIMATORS FOR SHAPE PARAMETER OF CLASSICAL PARETO DISTRIBUTION

SHRINKAGE ESTIMATOR AND TESTIMATORS FOR SHAPE PARAMETER OF CLASSICAL PARETO DISTRIBUTION Journal of Scientific Research Vol. 55, 0 : 8-07 Banaras Hindu University, Varanasi ISSN : 0447-9483 SHRINKAGE ESTIMATOR AND TESTIMATORS FOR SHAPE PARAMETER OF CLASSICAL PARETO DISTRIBUTION Ganesh Singh,

More information

An Efficient Single Stage Shrinkage Estimator for the Scale parameter of Inverted Gamma Distribution

An Efficient Single Stage Shrinkage Estimator for the Scale parameter of Inverted Gamma Distribution مجلة إبن الهيثم للعلوم الصرفة و التطبيقية المجلد )30 ) العدد )1( عام 2012 Ibn Al-Haitham J. for Pure & Appl. Sci. Vol.30 (1) 2017 An Efficient Single Stage Shrinkage Estimator for the Scale parameter of

More information

ISSN: ISO 9001:2008 Certified International Journal of Engineering and Innovative Technology (IJEIT) Volume 4, Issue 6, December 2014

ISSN: ISO 9001:2008 Certified International Journal of Engineering and Innovative Technology (IJEIT) Volume 4, Issue 6, December 2014 Bayesian Estimation for the Reliability Function of Pareto Type I Distribution under Generalized Square Error Loss Function Dr. Huda A. Rasheed, Najam A. Aleawy Al-Gazi Abstract The main objective of this

More information

Journal of Statistical Research 2007, Vol. 41, No. 1, pp Bangladesh

Journal of Statistical Research 2007, Vol. 41, No. 1, pp Bangladesh Journal of Statistical Research 007, Vol. 4, No., pp. 5 Bangladesh ISSN 056-4 X ESTIMATION OF AUTOREGRESSIVE COEFFICIENT IN AN ARMA(, ) MODEL WITH VAGUE INFORMATION ON THE MA COMPONENT M. Ould Haye School

More information

Terminology Suppose we have N observations {x(n)} N 1. Estimators as Random Variables. {x(n)} N 1

Terminology Suppose we have N observations {x(n)} N 1. Estimators as Random Variables. {x(n)} N 1 Estimation Theory Overview Properties Bias, Variance, and Mean Square Error Cramér-Rao lower bound Maximum likelihood Consistency Confidence intervals Properties of the mean estimator Properties of the

More information

New Statistical Methods That Improve on MLE and GLM Including for Reserve Modeling GARY G VENTER

New Statistical Methods That Improve on MLE and GLM Including for Reserve Modeling GARY G VENTER New Statistical Methods That Improve on MLE and GLM Including for Reserve Modeling GARY G VENTER MLE Going the Way of the Buggy Whip Used to be gold standard of statistical estimation Minimum variance

More information

Estimators as Random Variables

Estimators as Random Variables Estimation Theory Overview Properties Bias, Variance, and Mean Square Error Cramér-Rao lower bound Maimum likelihood Consistency Confidence intervals Properties of the mean estimator Introduction Up until

More information

Intro to probability concepts

Intro to probability concepts October 31, 2017 Serge Lang lecture This year s Serge Lang Undergraduate Lecture will be given by Keith Devlin of our main athletic rival. The title is When the precision of mathematics meets the messiness

More information

Probability Distributions Columns (a) through (d)

Probability Distributions Columns (a) through (d) Discrete Probability Distributions Columns (a) through (d) Probability Mass Distribution Description Notes Notation or Density Function --------------------(PMF or PDF)-------------------- (a) (b) (c)

More information

Ordinary and Bayesian Shrinkage Estimation. Mohammad Qabaha

Ordinary and Bayesian Shrinkage Estimation. Mohammad Qabaha An - Najah Univ. J. Res. (N. Sc.) Vol., 007 Ordinary and ayesian Shrinkage Estimation التقدير باستخدام طريقتي التقلص العادية والتقلص لبييز Mohammad Qabaha محمد قبها Department of Mathematics, Faculty of

More information

ACTEX CAS EXAM 3 STUDY GUIDE FOR MATHEMATICAL STATISTICS

ACTEX 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 information

Bayesian Predictive Modeling for Exponential-Pareto Composite Distribution

Bayesian Predictive Modeling for Exponential-Pareto Composite Distribution Bayesian Predictive Modeling for Exponential-Pareto Composite Distribution ABSTRCT Composite distributions have well-known applications in the insurance industry. In this article a composite Exponential-Pareto

More information

INVERTED KUMARASWAMY DISTRIBUTION: PROPERTIES AND ESTIMATION

INVERTED 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 information

Robustness of location estimators under t- distributions: a literature review

Robustness of location estimators under t- distributions: a literature review IOP Conference Series: Earth and Environmental Science PAPER OPEN ACCESS Robustness of location estimators under t- distributions: a literature review o cite this article: C Sumarni et al 07 IOP Conf.

More information

Part 4: Multi-parameter and normal models

Part 4: Multi-parameter and normal models Part 4: Multi-parameter and normal models 1 The normal model Perhaps the most useful (or utilized) probability model for data analysis is the normal distribution There are several reasons for this, e.g.,

More information

ON THE FAILURE RATE ESTIMATION OF THE INVERSE GAUSSIAN DISTRIBUTION

ON 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 information

INTRODUCTION TO ANALYSIS OF VARIANCE

INTRODUCTION TO ANALYSIS OF VARIANCE CHAPTER 22 INTRODUCTION TO ANALYSIS OF VARIANCE Chapter 18 on inferences about population means illustrated two hypothesis testing situations: for one population mean and for the difference between two

More information

Statistics Ph.D. Qualifying Exam: Part I October 18, 2003

Statistics Ph.D. Qualifying Exam: Part I October 18, 2003 Statistics Ph.D. Qualifying Exam: Part I October 18, 2003 Student Name: 1. Answer 8 out of 12 problems. Mark the problems you selected in the following table. 1 2 3 4 5 6 7 8 9 10 11 12 2. Write your answer

More information

Complex Systems Methods 11. Power laws an indicator of complexity?

Complex Systems Methods 11. Power laws an indicator of complexity? Complex Systems Methods 11. Power laws an indicator of complexity? Eckehard Olbrich e.olbrich@gmx.de http://personal-homepages.mis.mpg.de/olbrich/complex systems.html Potsdam WS 2007/08 Olbrich (Leipzig)

More information

Statistical Models with Uncertain Error Parameters (G. Cowan, arxiv: )

Statistical Models with Uncertain Error Parameters (G. Cowan, arxiv: ) Statistical Models with Uncertain Error Parameters (G. Cowan, arxiv:1809.05778) Workshop on Advanced Statistics for Physics Discovery aspd.stat.unipd.it Department of Statistical Sciences, University of

More information

Parameters Estimation for a Linear Exponential Distribution Based on Grouped Data

Parameters 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 information

Sampling Distributions of Statistics Corresponds to Chapter 5 of Tamhane and Dunlop

Sampling Distributions of Statistics Corresponds to Chapter 5 of Tamhane and Dunlop Sampling Distributions of Statistics Corresponds to Chapter 5 of Tamhane and Dunlop Slides prepared by Elizabeth Newton (MIT), with some slides by Jacqueline Telford (Johns Hopkins University) 1 Sampling

More information

Best linear unbiased and invariant reconstructors for the past records

Best linear unbiased and invariant reconstructors for the past records Best linear unbiased and invariant reconstructors for the past records B. Khatib and Jafar Ahmadi Department of Statistics, Ordered and Spatial Data Center of Excellence, Ferdowsi University of Mashhad,

More information

p y (1 p) 1 y, y = 0, 1 p Y (y p) = 0, otherwise.

p y (1 p) 1 y, y = 0, 1 p Y (y p) = 0, otherwise. 1. Suppose Y 1, Y 2,..., Y n is an iid sample from a Bernoulli(p) population distribution, where 0 < p < 1 is unknown. The population pmf is p y (1 p) 1 y, y = 0, 1 p Y (y p) = (a) Prove that Y is the

More information

Economics 573 Problem Set 4 Fall 2002 Due: 20 September

Economics 573 Problem Set 4 Fall 2002 Due: 20 September Economics 573 Problem Set 4 Fall 2002 Due: 20 September 1. Ten students selected at random have the following "final averages" in physics and economics. Students 1 2 3 4 5 6 7 8 9 10 Physics 66 70 50 80

More information

STA 2201/442 Assignment 2

STA 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 information

Economics 520. Lecture Note 19: Hypothesis Testing via the Neyman-Pearson Lemma CB 8.1,

Economics 520. Lecture Note 19: Hypothesis Testing via the Neyman-Pearson Lemma CB 8.1, Economics 520 Lecture Note 9: Hypothesis Testing via the Neyman-Pearson Lemma CB 8., 8.3.-8.3.3 Uniformly Most Powerful Tests and the Neyman-Pearson Lemma Let s return to the hypothesis testing problem

More information

Comparing Different Estimators of Reliability Function for Proposed Probability Distribution

Comparing Different Estimators of Reliability Function for Proposed Probability Distribution American Journal of Mathematics and Statistics 21, (2): 84-94 DOI: 192/j.ajms.212. Comparing Different Estimators of Reliability Function for Proposed Probability Distribution Dhwyia S. Hassun 1, Nathie

More information

Reconstruction of Order Statistics in Exponential Distribution

Reconstruction of Order Statistics in Exponential Distribution JIRSS 200) Vol. 9, No., pp 2-40 Reconstruction of Order Statistics in Exponential Distribution M. Razmkhah, B. Khatib, Jafar Ahmadi Department of Statistics, Ordered and Spatial Data Center of Excellence,

More information

ECON 4130 Supplementary Exercises 1-4

ECON 4130 Supplementary Exercises 1-4 HG Set. 0 ECON 430 Sulementary Exercises - 4 Exercise Quantiles (ercentiles). Let X be a continuous random variable (rv.) with df f( x ) and cdf F( x ). For 0< < we define -th quantile (or 00-th ercentile),

More information

Week 1 Quantitative Analysis of Financial Markets Distributions A

Week 1 Quantitative Analysis of Financial Markets Distributions A Week 1 Quantitative Analysis of Financial Markets Distributions A Christopher Ting http://www.mysmu.edu/faculty/christophert/ Christopher Ting : christopherting@smu.edu.sg : 6828 0364 : LKCSB 5036 October

More information

Asymptotic distribution of the sample average value-at-risk

Asymptotic distribution of the sample average value-at-risk Asymptotic distribution of the sample average value-at-risk Stoyan V. Stoyanov Svetlozar T. Rachev September 3, 7 Abstract In this paper, we prove a result for the asymptotic distribution of the sample

More information

Inference on reliability in two-parameter exponential stress strength model

Inference 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 information

Inference and Regression

Inference and Regression Name Inference and Regression Final Examination, 2016 Department of IOMS This course and this examination are governed by the Stern Honor Code. Instructions Please write your name at the top of this page.

More information

Moments of the Reliability, R = P(Y<X), As a Random Variable

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 information

Testing Hypothesis. Maura Mezzetti. Department of Economics and Finance Università Tor Vergata

Testing Hypothesis. Maura Mezzetti. Department of Economics and Finance Università Tor Vergata Maura Department of Economics and Finance Università Tor Vergata Hypothesis Testing Outline It is a mistake to confound strangeness with mystery Sherlock Holmes A Study in Scarlet Outline 1 The Power Function

More information

Module 22: Bayesian Methods Lecture 9 A: Default prior selection

Module 22: Bayesian Methods Lecture 9 A: Default prior selection Module 22: Bayesian Methods Lecture 9 A: Default prior selection Peter Hoff Departments of Statistics and Biostatistics University of Washington Outline Jeffreys prior Unit information priors Empirical

More information

High-dimensional regression

High-dimensional regression High-dimensional regression Advanced Methods for Data Analysis 36-402/36-608) Spring 2014 1 Back to linear regression 1.1 Shortcomings Suppose that we are given outcome measurements y 1,... y n R, and

More information

Unobservable Parameter. Observed Random Sample. Calculate Posterior. Choosing Prior. Conjugate prior. population proportion, p prior:

Unobservable Parameter. Observed Random Sample. Calculate Posterior. Choosing Prior. Conjugate prior. population proportion, p prior: Pi Priors Unobservable Parameter population proportion, p prior: π ( p) Conjugate prior π ( p) ~ Beta( a, b) same PDF family exponential family only Posterior π ( p y) ~ Beta( a + y, b + n y) Observed

More information

Chapter 14 Stein-Rule Estimation

Chapter 14 Stein-Rule Estimation Chapter 14 Stein-Rule Estimation The ordinary least squares estimation of regression coefficients in linear regression model provides the estimators having minimum variance in the class of linear and unbiased

More information

3 Continuous Random Variables

3 Continuous Random Variables Jinguo Lian Math437 Notes January 15, 016 3 Continuous Random Variables Remember that discrete random variables can take only a countable number of possible values. On the other hand, a continuous random

More information

Estimation of Parameters

Estimation of Parameters CHAPTER Probability, Statistics, and Reliability for Engineers and Scientists FUNDAMENTALS OF STATISTICAL ANALYSIS Second Edition A. J. Clark School of Engineering Department of Civil and Environmental

More information

WU Weiterbildung. Linear Mixed Models

WU Weiterbildung. Linear Mixed Models Linear Mixed Effects Models WU Weiterbildung SLIDE 1 Outline 1 Estimation: ML vs. REML 2 Special Models On Two Levels Mixed ANOVA Or Random ANOVA Random Intercept Model Random Coefficients Model Intercept-and-Slopes-as-Outcomes

More information

Testing and Model Selection

Testing and Model Selection Testing and Model Selection This is another digression on general statistics: see PE App C.8.4. The EViews output for least squares, probit and logit includes some statistics relevant to testing hypotheses

More information

Keywords: Estimator, Bias, Mean-squared error, normality, generalized Pareto distribution

Keywords: Estimator, Bias, Mean-squared error, normality, generalized Pareto distribution Testing approxiate norality of an estiator using the estiated MSE and bias with an application to the shape paraeter of the generalized Pareto distribution J. Martin van Zyl Abstract In this work the norality

More information

Two hours. To be supplied by the Examinations Office: Mathematical Formula Tables THE UNIVERSITY OF MANCHESTER. 21 June :45 11:45

Two 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 information

Pareto and Laplace Distributions

Pareto and Laplace Distributions Pareto and Laplace Distributions In this chapter we consider the Pareto distribution and Laplace distribution (double exponential distribution) 6.1. The Pareto distribution The Pareto distribution is a

More information

Introduction to Econometrics. Review of Probability & Statistics

Introduction to Econometrics. Review of Probability & Statistics 1 Introduction to Econometrics Review of Probability & Statistics Peerapat Wongchaiwat, Ph.D. wongchaiwat@hotmail.com Introduction 2 What is Econometrics? Econometrics consists of the application of mathematical

More information

ESTIMATION OF THE RECIPROCAL OF THE MEAN OF THE INVERSE GAUSSIAN DISTRIBUTION WITH PRIOR INFORMATION

ESTIMATION OF THE RECIPROCAL OF THE MEAN OF THE INVERSE GAUSSIAN DISTRIBUTION WITH PRIOR INFORMATION STATISTICA, anno LXVIII, n., 008 ESTIMATION OF THE RECIPROCAL OF THE MEAN OF THE INVERSE GAUSSIAN DISTRIBUTION WITH PRIOR INFORMATION 1. INTRODUCTION The Inverse Gaussian distribution was first introduced

More information

Optimum Test Plan for 3-Step, Step-Stress Accelerated Life Tests

Optimum 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 information

Chapter 5 Confidence Intervals

Chapter 5 Confidence Intervals Chapter 5 Confidence Intervals Confidence Intervals about a Population Mean, σ, Known Abbas Motamedi Tennessee Tech University A point estimate: a single number, calculated from a set of data, that is

More information

A comparison of inverse transform and composition methods of data simulation from the Lindley distribution

A 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 information

Chapter 6. Estimation of Confidence Intervals for Nodal Maximum Power Consumption per Customer

Chapter 6. Estimation of Confidence Intervals for Nodal Maximum Power Consumption per Customer Chapter 6 Estimation of Confidence Intervals for Nodal Maximum Power Consumption per Customer The aim of this chapter is to calculate confidence intervals for the maximum power consumption per customer

More information

Estimation of parametric functions in Downton s bivariate exponential distribution

Estimation of parametric functions in Downton s bivariate exponential distribution Estimation of parametric functions in Downton s bivariate exponential distribution George Iliopoulos Department of Mathematics University of the Aegean 83200 Karlovasi, Samos, Greece e-mail: geh@aegean.gr

More information

CSE 103 Homework 8: Solutions November 30, var(x) = np(1 p) = P r( X ) 0.95 P r( X ) 0.

CSE 103 Homework 8: Solutions November 30, var(x) = np(1 p) = P r( X ) 0.95 P r( X ) 0. () () a. X is a binomial distribution with n = 000, p = /6 b. The expected value, variance, and standard deviation of X is: E(X) = np = 000 = 000 6 var(x) = np( p) = 000 5 6 666 stdev(x) = np( p) = 000

More information

2008 Winton. Statistical Testing of RNGs

2008 Winton. Statistical Testing of RNGs 1 Statistical Testing of RNGs Criteria for Randomness For a sequence of numbers to be considered a sequence of randomly acquired numbers, it must have two basic statistical properties: Uniformly distributed

More information

Bayesian Regression Linear and Logistic Regression

Bayesian Regression Linear and Logistic Regression When we want more than point estimates Bayesian Regression Linear and Logistic Regression Nicole Beckage Ordinary Least Squares Regression and Lasso Regression return only point estimates But what if we

More information

Exam C Solutions Spring 2005

Exam C Solutions Spring 2005 Exam C Solutions Spring 005 Question # The CDF is F( x) = 4 ( + x) Observation (x) F(x) compare to: Maximum difference 0. 0.58 0, 0. 0.58 0.7 0.880 0., 0.4 0.680 0.9 0.93 0.4, 0.6 0.53. 0.949 0.6, 0.8

More information

EXAMINERS REPORT & SOLUTIONS STATISTICS 1 (MATH 11400) May-June 2009

EXAMINERS REPORT & SOLUTIONS STATISTICS 1 (MATH 11400) May-June 2009 EAMINERS REPORT & SOLUTIONS STATISTICS (MATH 400) May-June 2009 Examiners Report A. Most plots were well done. Some candidates muddled hinges and quartiles and gave the wrong one. Generally candidates

More information

Chapter 8 Heteroskedasticity

Chapter 8 Heteroskedasticity Chapter 8 Walter R. Paczkowski Rutgers University Page 1 Chapter Contents 8.1 The Nature of 8. Detecting 8.3 -Consistent Standard Errors 8.4 Generalized Least Squares: Known Form of Variance 8.5 Generalized

More information

MS&E 226: Small Data

MS&E 226: Small Data MS&E 226: Small Data Lecture 12: Frequentist properties of estimators (v4) Ramesh Johari ramesh.johari@stanford.edu 1 / 39 Frequentist inference 2 / 39 Thinking like a frequentist Suppose that for some

More information

Introduction to Maximum Likelihood Estimation

Introduction 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 information

Math/Stats 425, Sec. 1, Fall 04: Introduction to Probability. Final Exam: Solutions

Math/Stats 425, Sec. 1, Fall 04: Introduction to Probability. Final Exam: Solutions Math/Stats 45, Sec., Fall 4: Introduction to Probability Final Exam: Solutions. In a game, a contestant is shown two identical envelopes containing money. The contestant does not know how much money is

More information

ISI Web of Knowledge (Articles )

ISI 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 information

Do not copy, post, or distribute

Do not copy, post, or distribute 14 CORRELATION ANALYSIS AND LINEAR REGRESSION Assessing the Covariability of Two Quantitative Properties 14.0 LEARNING OBJECTIVES In this chapter, we discuss two related techniques for assessing a possible

More information

Week 2: Review of probability and statistics

Week 2: Review of probability and statistics Week 2: Review of probability and statistics Marcelo Coca Perraillon University of Colorado Anschutz Medical Campus Health Services Research Methods I HSMP 7607 2017 c 2017 PERRAILLON ALL RIGHTS RESERVED

More information

Probability. Lecture Notes. Adolfo J. Rumbos

Probability. Lecture Notes. Adolfo J. Rumbos Probability Lecture Notes Adolfo J. Rumbos October 20, 204 2 Contents Introduction 5. An example from statistical inference................ 5 2 Probability Spaces 9 2. Sample Spaces and σ fields.....................

More information

Confidence Intervals. Confidence interval for sample mean. Confidence interval for sample mean. Confidence interval for sample mean

Confidence Intervals. Confidence interval for sample mean. Confidence interval for sample mean. Confidence interval for sample mean Confidence Intervals Confidence interval for sample mean The CLT tells us: as the sample size n increases, the sample mean is approximately Normal with mean and standard deviation Thus, we have a standard

More information

GENERALIZED CONFIDENCE INTERVALS FOR THE SCALE PARAMETER OF THE INVERTED EXPONENTIAL DISTRIBUTION

GENERALIZED CONFIDENCE INTERVALS FOR THE SCALE PARAMETER OF THE INVERTED EXPONENTIAL DISTRIBUTION Internation Journ of Latest Research in Science and Technology ISSN (Online):7- Volume, Issue : Page No.-, November-December 0 (speci Issue Paper ) http://www.mnkjourns.com/ijlrst.htm Speci Issue on Internation

More information

Linear model selection and regularization

Linear model selection and regularization Linear model selection and regularization Problems with linear regression with least square 1. Prediction Accuracy: linear regression has low bias but suffer from high variance, especially when n p. It

More information

y ˆ i = ˆ " T u i ( i th fitted value or i th fit)

y ˆ i = ˆ  T u i ( i th fitted value or i th fit) 1 2 INFERENCE FOR MULTIPLE LINEAR REGRESSION Recall Terminology: p predictors x 1, x 2,, x p Some might be indicator variables for categorical variables) k-1 non-constant terms u 1, u 2,, u k-1 Each u

More information

Chapter 9. Non-Parametric Density Function Estimation

Chapter 9. Non-Parametric Density Function Estimation 9-1 Density Estimation Version 1.2 Chapter 9 Non-Parametric Density Function Estimation 9.1. Introduction We have discussed several estimation techniques: method of moments, maximum likelihood, and least

More information

Introduction to Rare Event Simulation

Introduction to Rare Event Simulation Introduction to Rare Event Simulation Brown University: Summer School on Rare Event Simulation Jose Blanchet Columbia University. Department of Statistics, Department of IEOR. Blanchet (Columbia) 1 / 31

More information

Maximum-Likelihood Estimation: Basic Ideas

Maximum-Likelihood Estimation: Basic Ideas Sociology 740 John Fox Lecture Notes Maximum-Likelihood Estimation: Basic Ideas Copyright 2014 by John Fox Maximum-Likelihood Estimation: Basic Ideas 1 I The method of maximum likelihood provides estimators

More information

Improving Reproducibility Probability estimation, preserving RP-testing

Improving Reproducibility Probability estimation, preserving RP-testing Improving Reproducibility Probability estimation, preserving -testing L. De Capitani, D. De Martini Department of Statistics and Quantitative Methods, University of Milano-Bicocca via Bicocca degli Arcimboldi

More information

Estimation for inverse Gaussian Distribution Under First-failure Progressive Hybird Censored Samples

Estimation for inverse Gaussian Distribution Under First-failure Progressive Hybird Censored Samples Filomat 31:18 (217), 5743 5752 https://doi.org/1.2298/fil1718743j Published by Faculty of Sciences and Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat Estimation for

More information

Problem 1 (20) Log-normal. f(x) Cauchy

Problem 1 (20) Log-normal. f(x) Cauchy ORF 245. Rigollet Date: 11/21/2008 Problem 1 (20) f(x) f(x) 0.0 0.1 0.2 0.3 0.4 0.0 0.2 0.4 0.6 0.8 4 2 0 2 4 Normal (with mean -1) 4 2 0 2 4 Negative-exponential x x f(x) f(x) 0.0 0.1 0.2 0.3 0.4 0.5

More information

9. Linear Regression and Correlation

9. Linear Regression and Correlation 9. Linear Regression and Correlation Data: y a quantitative response variable x a quantitative explanatory variable (Chap. 8: Recall that both variables were categorical) For example, y = annual income,

More information

Lecture 5. G. Cowan Lectures on Statistical Data Analysis Lecture 5 page 1

Lecture 5. G. Cowan Lectures on Statistical Data Analysis Lecture 5 page 1 Lecture 5 1 Probability (90 min.) Definition, Bayes theorem, probability densities and their properties, catalogue of pdfs, Monte Carlo 2 Statistical tests (90 min.) general concepts, test statistics,

More information

Statistical Inference Using Progressively Type-II Censored Data with Random Scheme

Statistical 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 information

A CHAIN RATIO EXPONENTIAL TYPE ESTIMATOR IN TWO- PHASE SAMPLING USING AUXILIARY INFORMATION

A CHAIN RATIO EXPONENTIAL TYPE ESTIMATOR IN TWO- PHASE SAMPLING USING AUXILIARY INFORMATION STATISTICA, anno LXXIII, n. 2, 2013 A CHAIN RATIO EXPONENTIAL TYPE ESTIMATOR IN TWO- PHASE SAMPLING USING AUXILIARY INFORMATION Rohini Yadav Department of Applied Mathematics, Indian School of Mines, Dhanbad,

More information

M(t) = 1 t. (1 t), 6 M (0) = 20 P (95. X i 110) i=1

M(t) = 1 t. (1 t), 6 M (0) = 20 P (95. X i 110) i=1 Math 66/566 - Midterm Solutions NOTE: These solutions are for both the 66 and 566 exam. The problems are the same until questions and 5. 1. The moment generating function of a random variable X is M(t)

More information

Institute of Actuaries of India

Institute of Actuaries of India Institute of Actuaries of India Subject CT3 Probability & Mathematical Statistics May 2011 Examinations INDICATIVE SOLUTION Introduction The indicative solution has been written by the Examiners with the

More information

Eco517 Fall 2004 C. Sims MIDTERM EXAM

Eco517 Fall 2004 C. Sims MIDTERM EXAM Eco517 Fall 2004 C. Sims MIDTERM EXAM Answer all four questions. Each is worth 23 points. Do not devote disproportionate time to any one question unless you have answered all the others. (1) We are considering

More information

Binary Logistic Regression

Binary Logistic Regression The coefficients of the multiple regression model are estimated using sample data with k independent variables Estimated (or predicted) value of Y Estimated intercept Estimated slope coefficients Ŷ = b

More information

Statistics - Lecture One. Outline. Charlotte Wickham 1. Basic ideas about estimation

Statistics - 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 information

Financial Econometrics and Quantitative Risk Managenent Return Properties

Financial Econometrics and Quantitative Risk Managenent Return Properties Financial Econometrics and Quantitative Risk Managenent Return Properties Eric Zivot Updated: April 1, 2013 Lecture Outline Course introduction Return definitions Empirical properties of returns Reading

More information

Double Stage Shrinkage Estimator of Two Parameters. Generalized Exponential Distribution

Double 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 information

A Very Brief Summary of Statistical Inference, and Examples

A Very Brief Summary of Statistical Inference, and Examples A Very Brief Summary of Statistical Inference, and Examples Trinity Term 2009 Prof. Gesine Reinert Our standard situation is that we have data x = x 1, x 2,..., x n, which we view as realisations of random

More information

t x 1 e t dt, and simplify the answer when possible (for example, when r is a positive even number). In particular, confirm that EX 4 = 3.

t x 1 e t dt, and simplify the answer when possible (for example, when r is a positive even number). In particular, confirm that EX 4 = 3. Mathematical Statistics: Homework problems General guideline. While working outside the classroom, use any help you want, including people, computer algebra systems, Internet, and solution manuals, but

More information

Applied Statistics and Econometrics

Applied Statistics and Econometrics Applied Statistics and Econometrics Lecture 5 Saul Lach September 2017 Saul Lach () Applied Statistics and Econometrics September 2017 1 / 44 Outline of Lecture 5 Now that we know the sampling distribution

More information

A Study of Five Parameter Type I Generalized Half Logistic Distribution

A Study of Five Parameter Type I Generalized Half Logistic Distribution Pure and Applied Mathematics Journal 2017; 6(6) 177-181 http//www.sciencepublishinggroup.com/j/pamj doi 10.11648/j.pamj.20170606.14 ISSN 2326-9790 (Print); ISSN 2326-9812 (nline) A Study of Five Parameter

More information

ECE531 Lecture 10b: Maximum Likelihood Estimation

ECE531 Lecture 10b: Maximum Likelihood Estimation ECE531 Lecture 10b: Maximum Likelihood Estimation D. Richard Brown III Worcester Polytechnic Institute 05-Apr-2011 Worcester Polytechnic Institute D. Richard Brown III 05-Apr-2011 1 / 23 Introduction So

More information

COMPARISON OF RELATIVE RISK FUNCTIONS OF THE RAYLEIGH DISTRIBUTION UNDER TYPE-II CENSORED SAMPLES: BAYESIAN APPROACH *

COMPARISON OF RELATIVE RISK FUNCTIONS OF THE RAYLEIGH DISTRIBUTION UNDER TYPE-II CENSORED SAMPLES: BAYESIAN APPROACH * Jordan Journal of Mathematics and Statistics JJMS 4,, pp. 6-78 COMPARISON OF RELATIVE RISK FUNCTIONS OF THE RAYLEIGH DISTRIBUTION UNDER TYPE-II CENSORED SAMPLES: BAYESIAN APPROACH * Sanku Dey ABSTRACT:

More information

IEOR 165 Lecture 7 1 Bias-Variance Tradeoff

IEOR 165 Lecture 7 1 Bias-Variance Tradeoff IEOR 165 Lecture 7 Bias-Variance Tradeoff 1 Bias-Variance Tradeoff Consider the case of parametric regression with β R, and suppose we would like to analyze the error of the estimate ˆβ in comparison to

More information

40.2. Interval Estimation for the Variance. Introduction. Prerequisites. Learning Outcomes

40.2. Interval Estimation for the Variance. Introduction. Prerequisites. Learning Outcomes Interval Estimation for the Variance 40.2 Introduction In Section 40.1 we have seen that the sampling distribution of the sample mean, when the data come from a normal distribution (and even, in large

More information

Chapter 5: HYPOTHESIS TESTING

Chapter 5: HYPOTHESIS TESTING MATH411: Applied Statistics Dr. YU, Chi Wai Chapter 5: HYPOTHESIS TESTING 1 WHAT IS HYPOTHESIS TESTING? As its name indicates, it is about a test of hypothesis. To be more precise, we would first translate

More information

Variations. ECE 6540, Lecture 10 Maximum Likelihood Estimation

Variations. ECE 6540, Lecture 10 Maximum Likelihood Estimation Variations ECE 6540, Lecture 10 Last Time BLUE (Best Linear Unbiased Estimator) Formulation Advantages Disadvantages 2 The BLUE A simplification Assume the estimator is a linear system For a single parameter

More information

STAT 536: Genetic Statistics

STAT 536: Genetic Statistics STAT 536: Genetic Statistics Tests for Hardy Weinberg Equilibrium Karin S. Dorman Department of Statistics Iowa State University September 7, 2006 Statistical Hypothesis Testing Identify a hypothesis,

More information