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

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

Research Note REGRESSION ANALYSIS IN MARKOV CHAIN * A. Y. ALAMUTI AND M. R. MESHKANI **

Estimation of the large covariance matrix with two-step monotone missing data

Introduction to Probability and Statistics

An Improved Generalized Estimation Procedure of Current Population Mean in Two-Occasion Successive Sampling

Estimating function analysis for a class of Tweedie regression models

Published: 14 October 2013

Ratio Estimators in Simple Random Sampling Using Information on Auxiliary Attribute

A Comparison between Biased and Unbiased Estimators in Ordinary Least Squares Regression

ECON 4130 Supplementary Exercises 1-4

A New Asymmetric Interaction Ridge (AIR) Regression Method

Covariance Matrix Estimation for Reinforcement Learning

On Wrapping of Exponentiated Inverted Weibull Distribution

AN OPTIMAL CONTROL CHART FOR NON-NORMAL PROCESSES

On split sample and randomized confidence intervals for binomial proportions

System Reliability Estimation and Confidence Regions from Subsystem and Full System Tests

arxiv: v1 [physics.data-an] 26 Oct 2012

On-Line Appendix. Matching on the Estimated Propensity Score (Abadie and Imbens, 2015)

On the Estimation Of Population Mean Under Systematic Sampling Using Auxiliary Attributes

4. Score normalization technical details We now discuss the technical details of the score normalization method.

DIFFERENTIAL evolution (DE) [3] has become a popular

Combining Logistic Regression with Kriging for Mapping the Risk of Occurrence of Unexploded Ordnance (UXO)

Using the Divergence Information Criterion for the Determination of the Order of an Autoregressive Process

Performance of lag length selection criteria in three different situations

A MIXED CONTROL CHART ADAPTED TO THE TRUNCATED LIFE TEST BASED ON THE WEIBULL DISTRIBUTION

On the capacity of the general trapdoor channel with feedback

DEPARTMENT OF ECONOMICS ISSN DISCUSSION PAPER 20/07 TWO NEW EXPONENTIAL FAMILIES OF LORENZ CURVES

Towards understanding the Lorenz curve using the Uniform distribution. Chris J. Stephens. Newcastle City Council, Newcastle upon Tyne, UK

Approximating min-max k-clustering

ASYMPTOTIC RESULTS OF A HIGH DIMENSIONAL MANOVA TEST AND POWER COMPARISON WHEN THE DIMENSION IS LARGE COMPARED TO THE SAMPLE SIZE

Hotelling s Two- Sample T 2

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

Lecture 23 Maximum Likelihood Estimation and Bayesian Inference

Asymptotically Optimal Simulation Allocation under Dependent Sampling

Lower Confidence Bound for Process-Yield Index S pk with Autocorrelated Process Data

#A64 INTEGERS 18 (2018) APPLYING MODULAR ARITHMETIC TO DIOPHANTINE EQUATIONS

arxiv: v3 [physics.data-an] 23 May 2011

Ž. Ž. Ž. 2 QUADRATIC AND INVERSE REGRESSIONS FOR WISHART DISTRIBUTIONS 1

Robust Beamforming via Matrix Completion

Lecture: Condorcet s Theorem

Probability Estimates for Multi-class Classification by Pairwise Coupling

Slash Distributions and Applications

Statics and dynamics: some elementary concepts

Notes on Instrumental Variables Methods

One-way ANOVA Inference for one-way ANOVA

LECTURE 7 NOTES. x n. d x if. E [g(x n )] E [g(x)]

Supplement to the paper Accurate distributions of Mallows C p and its unbiased modifications with applications to shrinkage estimation

STA 250: Statistics. Notes 7. Bayesian Approach to Statistics. Book chapters: 7.2

E-companion to A risk- and ambiguity-averse extension of the max-min newsvendor order formula

MODELING THE RELIABILITY OF C4ISR SYSTEMS HARDWARE/SOFTWARE COMPONENTS USING AN IMPROVED MARKOV MODEL

Tests for Two Proportions in a Stratified Design (Cochran/Mantel-Haenszel Test)

Supplementary Materials for Robust Estimation of the False Discovery Rate

Deriving Indicator Direct and Cross Variograms from a Normal Scores Variogram Model (bigaus-full) David F. Machuca Mory and Clayton V.

General Linear Model Introduction, Classes of Linear models and Estimation

An Ant Colony Optimization Approach to the Probabilistic Traveling Salesman Problem

Chapter 7 Sampling and Sampling Distributions. Introduction. Selecting a Sample. Introduction. Sampling from a Finite Population

ON THE LEAST SIGNIFICANT p ADIC DIGITS OF CERTAIN LUCAS NUMBERS

Elementary Analysis in Q p

Bayesian Spatially Varying Coefficient Models in the Presence of Collinearity

Discrete alpha-skew-laplace distribution

1. INTRODUCTION. Fn 2 = F j F j+1 (1.1)

1 Probability Spaces and Random Variables

Combinatorics of topmost discs of multi-peg Tower of Hanoi problem

Damage Identification from Power Spectrum Density Transmissibility

Developing A Deterioration Probabilistic Model for Rail Wear

Generalized optimal sub-pattern assignment metric

A Simple Weight Decay Can Improve. Abstract. It has been observed in numerical simulations that a weight decay can improve

arxiv: v2 [stat.me] 3 Nov 2014

CHAPTER 5 STATISTICAL INFERENCE. 1.0 Hypothesis Testing. 2.0 Decision Errors. 3.0 How a Hypothesis is Tested. 4.0 Test for Goodness of Fit

How to Estimate Expected Shortfall When Probabilities Are Known with Interval or Fuzzy Uncertainty

Information collection on a graph

Estimation and Detection

For q 0; 1; : : : ; `? 1, we have m 0; 1; : : : ; q? 1. The set fh j(x) : j 0; 1; ; : : : ; `? 1g forms a basis for the tness functions dened on the i

MATHEMATICAL MODELLING OF THE WIRELESS COMMUNICATION NETWORK

Estimating Time-Series Models

Pairwise active appearance model and its application to echocardiography tracking

Robustness of classifiers to uniform l p and Gaussian noise Supplementary material

Recursive Estimation of the Preisach Density function for a Smart Actuator

Estimation of Sensitive Characteristic using Two-Phase Sampling for Regression Type Estimator

PRELIMINARY SELECTION OF BASIC PARAMETERS OF DIFFERENT TORSIONAL VIBRATION DAMPERS INTENDED FOR USE IN MEDIUM-SPEED DIESEL ENGINES

Re-entry Protocols for Seismically Active Mines Using Statistical Analysis of Aftershock Sequences

Journal of Mathematical Analysis and Applications

On Introducing Asymmetry into Circular Distributions

On Line Parameter Estimation of Electric Systems using the Bacterial Foraging Algorithm

Improved Bounds on Bell Numbers and on Moments of Sums of Random Variables

The Poisson Regression Model

Some Theoretical Properties and Parameter Estimation for the Two-Sided Length Biased Inverse Gaussian Distribution

Session 12 : Monopropellant Thrusters

CERIAS Tech Report The period of the Bell numbers modulo a prime by Peter Montgomery, Sangil Nahm, Samuel Wagstaff Jr Center for Education

Preconditioning techniques for Newton s method for the incompressible Navier Stokes equations

Exponential Ratio Type Estimators of Population Mean under Non-Response

Efficient algorithms for the smallest enclosing ball problem

Information collection on a graph

Maxisets for μ-thresholding rules

VIBRATION ANALYSIS OF BEAMS WITH MULTIPLE CONSTRAINED LAYER DAMPING PATCHES

LINEAR SYSTEMS WITH POLYNOMIAL UNCERTAINTY STRUCTURE: STABILITY MARGINS AND CONTROL

arxiv: v1 [math-ph] 29 Apr 2016

THE search for nearest neighbors (NNs) is a crucial task in

Topic: Lower Bounds on Randomized Algorithms Date: September 22, 2004 Scribe: Srinath Sridhar

Modeling Residual-Geometric Flow Sampling

Transcription:

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 by (Schrodinger, 1915). Stemming from the earlier work by (Tweedie, 1947), the Inverse Gaussian (IG) distribution has received considerable attention in the statistical literature during the last two decades. Over the years it has found alications in different areas ranging from reliability and life testing to meteorology, biology, economics, medicine, market surveys and remote sensing. See (Seshadri, 1999) for an extensive list of alications (Sen and Khattree, 000). The review aer by (Folks and Chhikara, 1978) has resented various interesting roerties of the IG distribution. There are various alternative forms of a IG distribution available in the literature. Out of these we choose to work with the most familiar one having robability density function (.d.f). 1 λ λ ( x µ ) f( x, µλ, ) = ex 3 π x µ x, x > 0, λ > 0, µ > 0 (1) Here, the arameters µ and λ are called location and shae arameters resectively. The mean and the variance of the IG distribution corresonding to (1) 3 are resectively given by E( X ) = µ and V( X ) = µ λ. In cases, where the distribution has arisen from the inverse Gaussian rocess, one may be interested in 1 estimating the recirocal of the mean ie.. µ. Let x1, x,..., x n be a random samle of size n from the Inverse Gaussian distribution (1). It is well known that the maximum likelihood estimators for µ and λ are ˆ µ x n i = x = and ˆ n λ = i = 1 n V n 1 1, where V = i = 1 xi x.

0 The samle mean x follows Inverse Gaussian distribution with arameters µ and nλ. The samle mean x is unbiased for µ where as ˆλ is a biased estimator λ. The statistics ( xv, ) are jointly comlete sufficient statistics for ( µ, λ ) if both µ and λ are unknown. When λ is known, the uniformly minimum variance unbiased estimator (UMVUE) of the recirocal mean (1/ µ ) is ˆ 1 1 θ1 = x nλ () with the variance ˆ 1 1 Var( θ1) = + nλ µ nλ (3) We note that in many ractical situations the value of λ is not known, in such a case the maximum likelihood estimator (mle), uniformly minimum variance unbiased estimator (UMVUE) and minimum mean squared error (MMSE) estimators of (1/ µ ) are resectively given by and ˆ 1 θ ml =, (MLE) (4) x ˆ 1 V θ UMVUE =, (UMVUE) (5) x n( n 1) ˆ 1 V θ MMSE =, (MMSE) (6) x n( n+ 1) It is ointed out in (Sen and Khattree, 000) that the estimators ˆml θ ˆMMSE θ are articular members of the following class of estimators { V ; 0}, ˆUMVUE θ and ˆ 1 θ = α α (7) x The class ˆ θ is a convex subsace of the real line. But it does suffer from certain undesirable features. For instance, this class does not ensure the non-negativity of the estimators and unless α = 0, with ositive robability, any estimator includ-

Estimation of the recirocal of the mean of the Inverse Gaussian distribution etc. 03 ing the UMVUE and MMSE in ˆ θ can take negative values and hence be out of the arameter sace. Following (Lehmann, 1983,. 114) one may comment on ˆUMVUE θ as follows ˆUMVUE θ can take negative values although the estimate is known to be non-negative. Excet when n is small, the robability of such values is not large, but when they do occur they cause an embarrassment. This difficulty can be avoided by relacing it by zero whenever it leads to a negative value, the resulting estimator of course, will no larger be unbiased. This roblem has been further discussed by (Pandey and Malik, 1989). It is to be mentioned that the estimator ˆ θ 1 in () can be used only when the value of the arameter λ is known. The value of λ is not known in most of the ractical situations. However, in many ractical situations the exerimenter has some rior estimate regarding the value of the arameter, either due to his exerience or due to his acquaintance with the behavior of the system. Thus the exerimenter may have evidence that the value of λ is in the neighborhood of λ 0, a known value. We call λ 0 the exerimenter s rior guess see (Pandey and Malik, 1988). The study of the estimators based on rior oint estimate (or guess value) λ 0 revealed that these are better (in terms of mean squared error) than the usual estimators when the guess value is in vicinity of true value. This roerty necessitated the use of reliminary test of hyothesis to decide whether guess value is in vicinity of true value or not. The intention behind reliminary test was that if guess value is in vicinity of the true value, the estimator based on rior oint estimate or guess value λ 0 should be used otherwise usual estimators. In this aer, we define the estimator of the recirocal of the mean (1/ µ ) by incororating the additional information λ 0 of λ. The mean square error criterion will be used to judge the merits of the suggested estimator. Numerical illustrations are given in the suort of resent study. The relevant references in this connection are (Thomson, 1968), (Singh and Saxena, 001), (Singh and Shukla, 000, 003), (Tweedie, 1945, 56, 57a, 57b), (Travadi and Ratani, 1990), (Jani, 1991), (Kourouklis,1994), (Iwase and Seto, 1983, 85), (Iwase, 1987), (Korwar, 1980), (Srivastava et al. 1980), and (Singh and Pandit, 007, 008).. SUGGESTED ESTIMATOR AND ITS PROPERTIES If samle size is sufficiently large the maximum likelihood estimators (MLE) are consistent and efficient in terms of mean squared error. If samle size is small, there is relatively little information about the arameter available from the samle and if there is any rior estimate for the arameter, the shrinkage method can be useful, see (Pandey, 1983). Let the rior oint estimate λ 0 of the arameter λ be available. Then we define the following estimator for (1/ µ ) as

04 ˆ 1 θ x δ 1 n 0 = + ˆ δ W δ0, (8) where δ ˆ ˆ 0 = 1/ λ0, δ = 1/ λ = V/( n 1) is the UMVE of δ = 1/ λ ; W is a constant such that MSE of ˆ θ is least and is a non-zero real number. Taking exectation both sides of (8) we have ˆ 1 1 1 E( V ) E( θ) E δ0 Wδ = + 0 x n ( n 1) (9) We know for I.G. distribution (1) that 1 1 µ E = 1+ x µ nλ (10) and n+ 1 Γ E( V ) = λ n 1 Γ (11) Using (10) and (11) in (9) we have n+ 1 ˆ 1 1 1 Γ 1 E( θ) = + δ0 W + δ0 µ nλ n n 1 λ Γ ( n 1) Thus we get the bias of ˆ θ as n+ 1 ˆ 1 1 Γ 1 B( θ ) 0 W 0 nλ n δ δ = + n 1 λ Γ (1) The mean squared error of ˆ θ is given by

Estimation of the recirocal of the mean of the Inverse Gaussian distribution etc. 05 n+ 4 1 Γ 1 (1 ) δ MSE = MSE η W η. + + x n ( n 1) n 1 Γ n+ 1 Γ + Wη ( n 1) n 1 Γ n+ 1 Γ δ 1 η Wη + n ( n 1) n 1 Γ, (13) 1 where δ = λ, δ0 λ η = δ = λ and 0 1 1 3 δ 1 3δ = + = + n MSE x n λµ λ n µ n (14) The MSE at (13) is minimized for δ0 δ δ0 W = W( ), (15) δ δ where n+ 1 Γ n 1 W( ) = n+ 4 1 Γ Since δ in (15) is unknown, therefore, relacing δ by its estimate ˆ δ, we get an estimate of W as ˆ δ0 δ δ 0 Ŵ = W( ) ˆ δ δ 0 (16)

06 Thus relacing W by Ŵ in (8), we get a class of shrinkage estimators of (1/ ) µ as ˆ θ 1 1 [ ( )( ˆ )] ( ) = 0 W 0 x n δ + δ δ 1 1 1 1 1 = + W( ) x n λo λ λ 0 (17) As ointed out earlier regarding UMVUE ˆUMVUE θ, we note that the resulting family of shrinkage estimators ˆ θ ( ) can also take negative values for smaller values of n. However, this difficulty can be avoided by relacing it by zero whenever it yields negative values. The bias of ˆ θ ( ) is given by (17) ˆ δ B( θ( )) = [{1 W( )}{1 η}] (18) n and the MSE of ˆ θ ( ) is given by ˆ 1 δ ( W( )) MSE( θ( )) = MSE + ( η 1) (1 W( )) 1 + x n ( n 1) (19) which is less than 1 MSE x if 1 G < η < 1+ where G 1 1 1 ( W( )) G1 = 1 (1 W( )) ( n 1) Thus we have the following roosition. Proosition.1 The class of shrinkage estimators ˆ θ ( ) has smaller relative mean squared error (RMSE) than that of ˆ 1 θ ml = for η (1 G1,1 + G1). x The MSE exression at (19) can be re-exressed as

Estimation of the recirocal of the mean of the Inverse Gaussian distribution etc. 07 ˆ ˆ δ (1 + W( )) MSE( θ( ) ) = MSE( θumvue) + ( η 1) (1 W( )) (1 W( )), n ( n 1) (0) or MSE( θ ) MSE( θ ) ( η 1) (1 ( )) n ( n 1) ( n+ 1) ˆ ˆ δ ( W( )) ( ) = MMSE + W + where ˆ 1 MSE( θumvue ) = + nµλ nn ( 1) λ and ˆ 1 ( n + ) MSE( θmmse ) = + nµλ n ( n+ 1) λ, (1) () (3) Now we state the following roositions which can be easily roved from (0) and (1). Proosition. The class of shrinkage estimators ˆ θ ( ) has smaller relative mean squared error (RMSE) than that of UMVUE ˆ 1 V θumvue =, for η (1 G, 1 + G ). x n( n 1) where G 1 (1 ( W( )) ) = ( n 1)(1 W( )). Proosition.3 The class of shrinkage estimators ˆ θ ( ) has smaller RMSE than that of MMSE estimator ˆ 1 V θmmse = x n( n+ 1) for η (1 G3, 1 + G3),

08 where G 1 ( W( )) = (1 W( )) ( n+ 1) ( n 1) 3 1. 3. SPECIAL CASE In this section we discuss the roerties of an estimator ˆ θ (1) (say), which is a articular member of the roosed class of shrinkage estimators ˆ θ ( ) defined at (17). For =1, the estimator θ ( ) reduces to the estimator ˆ θ (1), thus ˆ ˆ 1 1 n 1 ˆ θ(1) = δ0 + ( δ δ0 ) x n n+ 1 ˆ δ 0 = θmmse nn ( + 1) (4) Putting =1 in (18) and (19) we get the bias and MSE of ˆ θ (1) resectively as ˆ (1 η ) δ B( θ(1) ) = nn ( + 1) and (5) ˆ ˆ 4δ MSE( θ(1) ) = MSE( θmmse ) + η( η ), (6) n ( n+ 1) where MSE( ˆ θ MMSE ) is given by (3). In order to comare the bias of ˆMMSE θ with B( ˆ θ (1)), we write the bias of ˆMMSE θ as ˆ δ B( θ MMSE ) = nn ( + 1) (7) It follows from (5) and (7) that B( ˆ θ ) < B( ˆ θ ) if (1) MMSE 0< η < (8)

Estimation of the recirocal of the mean of the Inverse Gaussian distribution etc. 09 Further from (6) we note that MSE( ˆ ) < MSE( ˆ ) if θ(1) θ MMSE 0< η < (9) Thus we state the following theorem. Theorem 3.1 The estimator ˆ θ (1) is less biased as well more efficient than MMSE estimator ˆMMSE θ iff 0< η <. Remark 3.1 It is to be noted that the estimator ˆMMSE θ is the minimum mean squared error estimator so it has less MSE than that of usual estimator ˆ θ ml = 1/ x and UMVUE ˆUMVUE θ. Thus it is interesting to note from (9) that the estimator ˆ θ (1) is more efficient than the estimators ˆ θ, ˆ ml θ UMVUE and ˆMMSE θ under the condition 0 < η <. 4. NUMERICAL ILLUSTRATION AND CONCLUSIONS To get tangible idea about the erformance of the roosed class of shrinkage estimators ˆ θ ( ) over usual estimator ˆ 1 θ ml =, UMVUE x ˆ 1 V θumvue = and MMSE estimator ˆ 1 V θ MMSE =, we x n( n 1) x n( n+ 1) have comuted the ercent relative efficiencies (PRE s) of ˆ θ ( ) with resect to ˆ θ, ˆ ml θ UMVUE and ˆMMSE θ resectively using the following formulae: ˆ ˆ A PRE( θ( ), θ ml ) = 100 ( A+ A ) ˆ ˆ PRE( θ( ), θ UMVUE) = 100 ( A+ A1 ) ˆ ˆ 1 na na 3 PRE( θ( ), θ MMSE) = 100 ( A+ A1 ) (30) (31) (3) n where A = 3 + C, (33)

10 ( W( )) A1 = ( η 1) (1 W( )) + 1, (34) ( n 1) A A 1 = + C n 1, (35) 1 ( n + ) = +, (36) C nn ( + 1) 3 µ C = λ 1, is the coefficient of variation. The PRE s have been comuted for different values of = ± 1, ±. n = 5, 10, 15; η = 0.5(0.5)1.75 ; and C = 1,5. The comuted values are dislayed in tables 4.1(a). 4.(a) and 4.3(a). We have also comuted the range of η for different values of, n, η and C; and resented in tables 4.1(b), 4.(b) and 4.3(b). We note that ( n 7)( n 9) W ( ) = ( n 1) W ( 1) = ( n 5) ( n 1) W (1) = ( n 1) ( n + 1) ( n 1) W () = ( n+ 5)( n+ 3) It is observed from tables 4.1(a), 4.(a) and 4.3(a) that the ercent relative efficiencies of ˆ θ ( ), =±, ± 1; with resect to ˆml θ, ˆUMVUE θ and ˆMMSE θ resectively (.. i e PRE( ˆ θ, ˆ θ ), PRE( ˆ θ, ˆ θ ) and PRE( ˆ θ, ˆ θ )): ( ) ml ( ) UMVUE ( ) MMSE (i) attain their maximum at η = 1 (i.e. when λ coincide with λ 0 ), (ii) decrease as n increases, (iii) increase as the value of coefficient of variation (C) increases, (iv) is more than 100% for the largest range of dominance of η when =1, (i.e. ˆ θ (1)), but the gain in efficiency is smaller comared to other estimators ˆ θ ( ),

Estimation of the recirocal of the mean of the Inverse Gaussian distribution etc. 11 ˆ θ( 1) and ˆ θ (). However, the estimator ˆ θ () seems to be the good intermediate choice between the estimators ˆ θ( ), ˆ θ( 1) and ˆ θ (1). Comaring the results of tables 4.1(a), 4.(a) and 4.3(a) it is seen that the gain in efficiency by using ˆ θ ( ) over ˆml θ is the largest followed by ˆUMVUE θ and then ˆMMSE θ. Also the roosed family of shrinkage estimators ˆ ( ) θ is more efficient than ˆml θ for widest range of η followed by ˆUMVUE θ and ˆMMSE θ see tables 4.1(b), 4.(b) and 4.3(b). It is exected too. Thus we conclude that the suggested family of shrinkage estimators ˆ θ ( ) is to be recommended for its use in ractice when (i) the guessed value λ 0 moves in the vicinity of the true value λ. (ii) samle size n is small (i.e. in the situations where the samling is costly) and (iii) the oulation is heterogeneous. TABLE 4.1(a) Percent relative efficiency of ˆ θ ( ) over ˆml θ C P 1 η n 5 10 15 5 10 15 0.5 110.11 103.8 103.87 105.79 106.74 105.15 0.50 111.30 106.8 104.95 110.35 107.8 105.30 1 0.75 11.04 107.81 105.61 113.7 107.61 105.40 1.00 11.8 108.33 105.83 114.9 107.7 105.43 1.5 11.04 107.81 105.61 113.7 107.61 105.40 1.50 111.30 106.8 104.95 110.35 107.8 105.30 1.75 110.11 103.8 103.87 105.79 106.74 105.15 0.5 19.78 116.36 1.90 115.84 131.80 13.41 0.50 134.03 19.17 130.86 130.61 135.03 133.67 0.75 136.7 138.31 136.15 141.44 137.05 134.43 5 1.00 137.63 141.65 138.01 145.45 137.73 134.69 1.5 136.7 138.31 136.15 141.44 137.05 134.43 1.50 134.03 19.17 130.86 130.61 135.03 133.67 1.75 19.79 116.36 1.90 115.84 131.80 13.41 1 0.5 109.8 106.84 105.15 108.40 106.9 104.90 0.50 110.35 106.93 105.18 111.4 107.3 105.30 0.75 110.66 106.99 105.0 113.31 107.80 105.54 1 1.00 110.77 107.01 105.1 113.96 107.99 105.6 1.5 110.64 106.99 105.0 113.31 107.80 105.54 1.50 110.35 106.93 105.18 111.4 107.3 105.30 1.75 109.8 106.84 105.15 108.40 106.9 104.90 0.5 18.79 13.43 13.44 14.03 19.5 130.48 0.50 130.61 13.97 13.68 134.45 134.7 133.6 0.75 131.73 133.9 13.8 141.59 138.3 135.57 5 1.00 13.11 133.40 13.87 144.14 139.44 136.4 1.5 131.73 133.9 13.8 141.59 138.3 135.58 1.50 130.61 13.97 13.68 134.45 134.7 133.6 1.75 18.79 13.43 13.44 14.03 19.5 130.48

1 TABLE 4.(a) Percent relative efficiency of ˆ θ ( ) over ˆUMVUE θ C P 1 1 5 1 5 η n 5 10 15 5 10 15 0.5 103.3 97.61 98.9 99.17 100.35 100.10 0.50 104.35 99.10 99.95 103.40 100.86 100.30 0.75 105.03 101.36 100.60 106.0 101.17 100.40 1.00 105.6 101.85 100.80 107.10 101.7 100.40 1.5 105.03 101.36 100.60 106.0 101.17 100.40 1.50 104.35 99.10 99.95 103.40 100.86 100.30 1.75 103.3 97.61 98.9 99.17 100.35 100.10 0.5 109.51 89.74 93.64 97.74 101.65 100.90 0.50 113.09 99.6 99.61 110.0 104.14 101.80 0.75 115.35 106.67 103.70 119.30 105.70 10.40 1.00 116.13 109.5 105.10 1.70 106. 10.60 1.5 115.35 106.67 103.70 119.30 105.70 10.40 1.50 113.09 99.6 99.69 110.0 104.14 101.80 1.75 109.51 89.74 93.64 97.74 101.65 100.90 1 0.5 10.96 100.45 100.10 101.60 99.93 99.91 0.50 103.45 100.54 100.0 104.50 100.81 100.30 0.75 103.75 100.59 100.0 106.0 101.35 100.50 1.00 103.85 100.61 100.0 106.80 101.53 100.60 1.5 103.75 100.59 100.0 106.0 101.35 100.50 1.50 103.45 100.54 100.0 104.50 100.81 100.30 1.75 10.96 100.45 100.10 101.60 99.93 99.91 0.5 108.66 10.14 100.90 104.70 99.68 99.41 0.50 110.0 10.55 101.10 113.40 103.90 101.80 0.75 111.15 10.80 101.0 119.50 106.61 103.30 1.00 111.47 10.88 101.0 11.60 107.54 103.80 1.5 111.15 10.80 101.0 119.50 106.61 103.30 1.50 110.0 10.55 101.10 113.40 103.90 101.80 1.75 108.66 10.14 100.90 104.70 99.68 99.41

Estimation of the recirocal of the mean of the Inverse Gaussian distribution etc. 13 TABLE 4.3(a) Percent relative efficiency of ˆ θ ( ) over ˆMMSE θ C P 1 1 5 1 5 η n 5 10 15 5 10 15 0.5 100.93 97.8 98.8 96.98 100.0 100.04 0.50 10.03 99.59 99.85 101.15 100.53 100.18 0.75 10.70 101.0 100.47 103.83 100.83 100.7 1.00 10.9 101.51 100.68 104.76 100.94 100.30 1.5 10.70 101.0 100.47 103.83 100.83 100.7 1.50 10.03 99.59 99.85 101.15 100.53 100.18 1.75 100.93 97.8 98.8 96.97 100.0 100.04 0.5 10.75 88.36 93.0 91.70 100.08 100.3 0.50 106.11 98.09 99.05 103.40 10.54 101.18 0.75 108.3 105.03 103.06 111.97 104.07 101.75 1.00 108.96 107.56 104.46 115.15 104.59 101.95 1.5 108.3 105.03 103.06 111.97 104.07 101.75 1.50 106.1 98.09 99.05 103.40 10.54 101.18 1.75 10.75 88.36 93.03 91.70 100.08 100.3 1 0.5 100.67 100.11 100.04 99.37 99.60 99.80 0.50 101.15 100.0 100.07 10.13 100.48 100.18 0.75 101.44 100.5 100.09 103.87 101.01 100.41 1.00 101.54 100.7 100.09 104.46 101.19 100.48 1.5 101.44 100.5 100.09 103.87 101.01 100.41 1.50 101.15 100.0 100.07 10.16 100.48 100.18 1.75 100.67 100.1 100.04 99.37 99.60 99.80 0.5 101.96 100.56 100.5 98.19 98.15 98.76 0.50 103.40 100.97 100.43 106.44 10.30 101.14 0.75 104.9 101.1 100.54 11.09 104.96 10.6 1.00 104.59 101.30 100.58 114.11 105.88 103.13 1.5 104.9 101.1 100.54 11.09 104.96 10.6 1.50 103.40 100.96 100.43 106.44 10.30 101.14 1.75 101.96 100.56 100.5 98.19 98.15 98.76

14 n Estimator TABLE 4.1(b) Range of η for different values of n 5 10 15 ˆ θ( ) (0.0000,.8708) (0.0000,.0383) (0.0000,.3186) ˆ θ( 1) (0.0000,.0000) (0.0000, 3.1715) (0.0000, 4.3700) ˆ θ (1) (0.0000, 3.6458) (0.0000, 6.0745) (0.0000, 8.5498) ˆ θ () (0.0000,.374) (0.0000,.6774) (0.0000, 3.1481) n Estimator TABLE 4.(b) Range of η for different values of n 5 10 15 ˆ θ( ) (0.0000,.1067) (0.3006, 1.6994) (0.3033, 1.6966) ˆ θ( 1) (0.1591, 1.8409) (0.0609, 1.9391) (0.3780, 1.96) ˆ θ (1) (0.0000,.574) (0.0000,.09) (0.0000,.099) ˆ θ () (0.0693, 1.9306) (0.1436, 1.8564) (0.1658, 1.834) n Estimator TABLE 4.3(b) Range of η for different values of n 5 10 15 ˆ θ( ) (0.0446,1.9554) (0.3348,1.6651) (0.378,1.67) ˆ θ( 1) (0.40,1.7598) (0.199,1.8701) (0.1061,1.8939) ˆ θ (1) (0.0000,.0000) (0.0000,.0000) (0.0000,.0000) ˆ θ () (0.1635,1.8365) (0.1951,1.8049) (0.056,1.7944) ACKNOWLEDGEMENT The authors are thankful to the executive editor Professor ssa. Marilena Pillati and an anonymous referee for their constructive suggestions regarding imrovement of the aer. School of Studies in Statistics Vikram University, Ujjain, India HOUSILA P. SINGH SHEETAL PANDIT

Estimation of the recirocal of the mean of the Inverse Gaussian distribution etc. 15 REFERENCES J.L. FOLKS, R.S. CHHIKARA, (1978), The Inverse Gaussian distribution and its alication, a review, Journal of Royal Statistical Society, B, 40,. 63-89. R.V. HOGG, (1974), Adotive robust rocedures, A artial review and some suggestions for future alications and theory, Journal of the American Statistical Association, 69,. 909-93. KOSEI IWASE, (1987), UMVU estimation for the Inverse Gaussian distribution I ( µ, c µ ) with known C, Communication Statistics, 16, 5,. 1315-130. KOSEI IWASE, NORAKI SETO, (1983), UMVUE s for the Inverse Gaussian distribution, Journal of the American Statistical Association, 78,. 660-663. KOSEI IWASE, NORAKI SETO, (1985), UMVU estimators of the mode and limits of an interval for the Inverse Gaussian distribution, Communication in Statistics, Theory and Methods, 14,. 1151-1161. P.N. JANI, (1991), A class of shrinkage estimators for the scale arameter of the exonential distribution, IEEE Transaction on Reliability, 40,. 60-70. R.M. KORWAR, (1980), On the UMVUE s of the variation and its recirocal of an Inverse Gaussian distribution, Journal of the American Statistical Association, 75,. 734-735. S. KOUROUKLIS, (1994), Estimation in the two-arameter exonential distribution with rior information, IEEE Transaction on Reliability, 43, 3,. 446-450. E.L. LEHMANN, (1983), Theory of oint estimation, Wiley, New York,. 114. B.N. PANDEY, (1983), Shrinkage estimation of the exonential scale arameter, IEEE Transaction on Reliability, 3,,. 03-05. B.N. PANDEY, H.J. MALIK, (1988), Some imroved estimators for a measure of disersion of an Inverse Gaussian distribution, Communication in Statistics, Theory and Methods, 17, 11,. 3935-3949. B.N. PANDEY, H.J. MALIK, (1989), Estimation of the mean and the recirocal of the mean of the Inverse Gaussian distribution, Communication Statistics, Simulation, 18(3),. 1187-101. E. SCHRONDINGER, (1915), Zur theories der fall-und Steigvursuche on Tailchenmit Brownsher Bewegung, Physikaliche Zeitschrift, 16,. 89-95. SEN, R. KHATTREE, (000), Revisiting the roblem of estimating the Inverse Gaussian arameters, IAPQR Transactions, 5,,. 63-79. R. SESHADRI, (1999), The Inverse Gaussian distribution, Statistical Theory and Alication, Sringer Verlag, New York. H.P. SINGH, S.K. SHUKLA, (000), Estimation in the two-arameter Weibull distribution with rior information, IAPQR. Transactions, 5,,. 107-117. H.P. SINGH, S. SAXENA, (001), Imroved estimation in one-arameter Exonential distribution with rior information, Gujarat Statistical Review, 8, 1-,. 5-35. H.P. SINGH, S.K. SHUKLA, (003), A family of shrinkage estimators of the square of mean in normal distribution, Statistical Paers, Verlag, 44,. 433-44. H.P. SINGH, S. PANDIT, (007), Imroved estimation if Inverse Gaussian shae arameter and measure of disersion with rior information, Journal of statistical Theory and Practice, 1,,. 05-6. H.P. SINGH, S. PANDIT, (008), Estimation of shae arameter and measure of disersion of Inverse Gaussian distribution using rior information, Bulletin of statistics and economics,, S08,. 81-100. V.K. SRIVASTAVA, T.D. DWIVEDI, S. BHATNAGAR, (1980), Estimation of the square of mean in normal oulation, Statistica, XL, 4,. 455-466. J.R. THOMPSON, (1968), Some shrinkage technique for estimating the mean, Journal of the American Statistical Association, 63,. 113-13.

16 R.J. TRAVEDI, R.T. RATANI, (1990), On estimation of reliability function for Inverse Gaussian distribution with known coefficient of variation, IAPQR, Transactions, 5,,. 9-37. M.C.K. TWEEDIE, (1945), Inverse Statistical Variate, Nature,. 155-453. M.C.K. TWEEDIE, (1947), Functions of statistical variate with given means, with secial inference to lalacian distributions, Proceedings of the Cambridge Philosohical Society, 43,. 41-49. M.C.K. TWEEDIE, (1956), Some statistical roerties of Inverse Gaussian distribution. Virginia Journal Science, 7,. 160-165. M.C.K. TWEEDIE, (1957 a), Statistical roerties of Inverse Gaussian distribution, I, Annals of Mathematical Statistics, 8,. 36-377. M.C.K. TWEEDIE, (1957 b), Statistical roerties of Inverse Gaussian distribution, II, Annals of Mathematical Statistics, 8,. 696-705. SUMMARY Estimation of the recirocal of the mean of the Inverse Gaussian distribution with rior information This aer considers the roblem of estimating the recirocal of the mean (1/ µ ) of the Inverse Gaussian distribution when a rior estimate or guessed value λ 0 of the shae arameter λ is available. We have roosed a class of estimators ˆ θ ( ), say, for (1/ µ ) with its mean squared error formula. Realistic conditions are obtained in which the estimator ˆ θ ( ) is better than usual estimator, uniformly minimum variance unbiased estimator (UMVUE) and the minimum mean squared error estimator (MMSE). Numerical illustrations are given in suort of the resent study.