Modied moment estimation for the two-parameter Birnbaum Saunders distribution

Size: px
Start display at page:

Download "Modied moment estimation for the two-parameter Birnbaum Saunders distribution"

Transcription

1 Computatioal Statistics & Data Aalysis 43 (23) Modied momet estimatio for the two-parameter Birbaum Sauders distributio H.K.T. Ng a, D. Kudu b, N. Balakrisha a; a Departmet of Mathematics ad Statistics, McMaster Uiversity, Hamilto, Ot., Caada L8S 4K1 b Departmet of Mathematics, Idia Istitute of Techology Kapur, Kapur 2816, Idia Received 1 March 22; received i revised form 1 July 22 Abstract The maximum likelihood estimators ad a modicatio of the momet estimators of a twoparameter Birbaum Sauders distributio are discussed. A simple bias-reductio method is proposed to reduce the bias of the maximum likelihood estimators ad the modied momet estimators. The jackkife techique is also used to reduce the bias of these estimators. Mote Carlo simulatio is used to compare the performace of all these estimators. The probability coverages of codece itervals based o iferetial quatities associated with all these estimators are evaluated usig Mote Carlo simulatios for small, moderate ad large sample sizes. Two illustrative examples ad some cocludig remarks are ally preseted. c 22 Elsevier Sciece B.V. All rights reserved. Keywords: Birbaum Sauders distributio; Jackkife estimate; Bias-reduced estimator; Mote Carlo simulatios; Probability coverage; Codece iterval 1. Itroductio The two-parameter Birbaum Sauders distributio was origially proposed (Birbaum ad Sauders, 1969a) as a failure time distributio for fatigue failure caused uder cyclic loadig. It was also assumed that the failure is due to the developmet ad growth of a domiat crack. A more geeral derivatio was provided by Desmod (1985) based o a biological model. Desmod (1985) also stregtheed the physical justicatio for the use of this distributio by relaxig the assumptios made by Birbaum ad Sauders (1969a). Desmod (1986) ivestigated the relatioship betwee the Birbaum Sauders distributio ad the iverse Gaussia distributio. Some Correspodig author. address: bala@mcmail.cis.mcmaster.ca (N. Balakrisha) /3/$ - see frot matter c 22 Elsevier Sciece B.V. All rights reserved. PII: S (2)254-2

2 284 H.K.T. Ng et al. / Computatioal Statistics & Data Aalysis 43 (23) recet works o Birbaum Sauders distributio ca be foud i Chag ad Tag (1993, 1994), Dupuis ad Mills (1998) ad Rieck (1995, 1999), ad a review of these developmets ca be foud i Johso et al. (1995). The cumulative distributio fuctio (CDF) of a two-parameter Birbaum Sauders radom variable T ca be writte as F T (t; ; )= [ 1 { ( t ) 1=2 ( ) }] 1=2 ; t ; ; ; (1) t where ( ) is the stadard ormal CDF. The parameters ad are the shape ad the scale parameters, respectively. It is kow that the desity fuctio of the Birbaum Sauders distributio is uimodal ad although the hazard rate is ot a icreasig fuctio of t but the average hazard rate is early a o-decreasig fuctio of t (Ma et al., 1974, p. 155). The maximum likelihood estimators (MLEs) were discussed origially by Birbaum ad Sauders (1969b) ad their asymptotic distributios were obtaied by Egelhardt et al. (1981). Although the MLEs have several optimal properties, oe still eeds to solve a o-liear equatio i to obtai the solutio; for this purpose, Birbaum ad Sauders (1969b) suggested some iterative schemes to solve the required o-liear equatio. Also, the exact distributio of the MLEs are ot available. Therefore, for costructig codece itervals of the ukow parameters ad, the asymptotic distributios of the MLEs eed to be used. However, it is ot kow how these asymptotic codece itervals behave i the case of small sample sizes. Moreover, the covetioal momet estimators also have a diculty i that they may ot always exist ad eve if they do, they may ot be uique. For this reaso, modied momet estimators (MMEs) for ad are rst proposed. The MMEs are very easy to compute ad they have explicit expressios i terms of the sample observatios. Ulike the momet estimators, MMEs always exist uiquely. The asymptotic distributios of the MMEs are derived which are the used to costruct co- dece itervals for the ukow parameters. The performace of all these estimators is evaluated through simulatios. Eve though the MLEs ad MMEs are asymptotically ubiased, these simulatio results reveal that they are highly biased i case of small sample sizes. A simple bias correctio techique is proposed which performs quite well eve for small sample sizes. Jackkife procedure is also used to reduce the bias of the MLEs ad MMEs ad is show to work very well i this case; but, this procedure becomes computatioally quite ivolved i case of large sample sizes. 2. The Birbaum--Sauders distributio The CDF of a two-parameter Birbaum Sauders radom variable T is give by (1) ad the correspodig probability desity fuctio (PDF) is [ ( ) 1=2 ( ) ] 3=2 [ 1 f T (t; ; )= 2 + exp 1 ( t 2 t t )] t 2 ; t ;; : (2)

3 or H.K.T. Ng et al. / Computatioal Statistics & Data Aalysis 43 (23) Cosider the followig mootoe trasformatio: [ (T X = 1 ) 1=2 ( ) ] 1=2 T 2 T = (1+2X 2 +2X (1 + X 2 ) 1=2 ); the, from (1), we have X to be a ormally distributed with mea zero ad variace Usig the above trasformatio, the expected value, variace, ad coeciets of skewess ad kurtosis ca be easily obtaied as E(T)=( ); (3) Var(T )=() 2 ( ); (4) 1 (T )= 162 ( ) (5 2 +4) 3 ; (5) 2 (T )=3+ 62 ( ) (5 2 +4) 2 : (6) Moreover, if T has a Birbaum Sauders distributio with parameters ad, the T 1 also has a Birbaum Sauders distributio with the correspodig parameters ad 1, respectively (Birbaum ad Sauders, 1969a). Therefore, we also readily have ad E(T 1 )= 1 ( ) (7) Var(T 1 )= 2 2 ( ): (8) 3. Maximum likelihood estimators Let {t 1 ;t 2 ;:::;t } be a radom sample of size from the Birbaum Sauders distributio with the PDF as give i (2). The sample arithmetic ad harmoic meas are deed by s = 1 t i ; r = [ 1 t 1 i ] 1 : Let us further dee the harmoic mea fuctio K by [ ] 1 1 K(x)= (x + t i ) 1 for x ; so that r K().

4 286 H.K.T. Ng et al. / Computatioal Statistics & Data Aalysis 43 (23) The MLE of (deoted by ˆ) ca be obtaied as the uique positive root of the equatio 2 [2r + K()] + r[s + K()]=: (9) Oce ˆ is obtaied as a solutio of (9), the MLE of (deoted by ˆ) ca be obtaied explicitly as [ s ˆ = ˆ + ˆ ] 1=2 r 2 : Sice (9) is a o-liear equatio i, oe eeds to use a iterative procedure to solve for ˆ. Birbaum ad Sauders (1969b) proposed two iterative procedures (oe simple ad oe complicated) to compute ˆ, but oted that the simple oe works very well for small ( 1 2 ) but may ot work at all for large ( 2). The complicated oe also does ot work i certai rage of the sample space. I this paper, Newto Raphso method is used for the computatio of the MLEs. Egelhardt et al. (1981) showed that the asymptotic joit distributio of ˆ ad ˆ is bivariate ormal ad is give by ( ) ( ) ˆ N ; ˆ where I()= [: I()] {[1 + g(x)] 1 :5} 2 d(x); ; (1) ) 1=2 g(y)=1+ (1+ y2 2 + y y2 : 4 It is iterestig to observe that ˆ ad ˆ are asymptotically idepedet of each other. The asymptotic codece iterval of ca be easily obtaied from (1). Moreover, the asymptotic codece iterval of, for a give, ca also be obtaied from (1). 4. Modied momet estimators For the usual momet estimators i a two-parameter case, the rst ad secod populatio momets are equated with the correspodig sample momets. I this case, the sample mea ad the sample variace ca be equated to the right-had sides of (3) ad (4), respectively, ad the correspodig momet estimators of ad ca the be obtaied as solutios of ad to these equatios. It ca be easily see from these equatios that if the sample coeciet of variatio is greater tha 5, the the momet estimators do ot exist. If the sample coeciet of variatio is less tha 5, the momet estimators exist; however, the momet estimator of may ot be uique. Istead of usig (3) ad (4), we propose to use (3) ad (7) ad equate them with the correspodig sample estimates to obtai the MMEs. I this case, we have the

5 H.K.T. Ng et al. / Computatioal Statistics & Data Aalysis 43 (23) followig two momet equatios: s = ( ) ; (11) r 1 = 1 ( ): (12) Solvig Eqs. (11) ad (12) for ad, we obtai the MMEs for ad deoted by ad as { [ ( s ) ]} 1=2 1=2 = 2 1 ; r =(sr) 1=2 : The asymptotic joit distributio of ad is bivariate ormal ad is give by ( ) ( ) 2 2 N ; ( 1+ 3 ) : (13) 4 2 () 2 ( ) 2 The proof of this result is preseted i the appedix. Note that the MMEs ad are also asymptotically idepedet of each other, just as i the case of the MLEs. 5. Bias-reduced estimators Based o the results of a extesive Mote Carlo simulatio study, we observed that the MLEs ad the MMEs performed very similarly i terms of both bias ad mea square error, especially for small values of. Upo ispectig the patter of the bias of the MLEs ad MMEs, we observed that Bias(ˆ) Bias( ) ; Bias( ˆ) Bias( ) 2 4 : The, by employig a stadard bias reductio method, we ca simply costruct almost ubiased maximum likelihood estimators (UMLEs, deoted by ˆ ad ˆ ) ad almost ubiased modied momet estimators (UMMEs, deoted by ad )of ad. These bias-reduced estimators are give by ˆ = ( ) 1 ( ) = ; 1 = ) 1 ) 1 ˆ; ˆ = (1+ 2 ˆ (1+ ˆ 2 ˆ; (14) 4 4 ( ) 1 ) 1 (1+ 2 : (15) 4

6 288 H.K.T. Ng et al. / Computatioal Statistics & Data Aalysis 43 (23) From the distributioal results preseted i (1), we readily have the asymptotic joit distributio of ˆ ad ˆ to be bivariate ormal ad is give by ( ˆ ) ( ) 2 2( 1) N ; 2 ; (16) ˆ 16 2 (4+ 2 ) 2 [: I()] similarly, from (13), the asymptotic joit distributio of ad is bivariate ormal ad is give by ( ) ( ) 2 2( 1) 2 N ; 1+ 3 ) : (17) 4 2 ( 16() 2 (4+ 2 ) 2 ( ) 2 6. Jackkife estimators Jackkig is based o sequetially deletig oe sample poit t i ad recomputig the MLEs ad MMEs from the reduced sample of size 1. We remove the poit t j from the data set, ad the recompute r ad s ad also the fuctio K as s (j) = 1 1 i j 1 r (j) = 1 t i = s t j 1 ; i j 1 K (j) (x)= 1 ti 1 1 = (x + t i ) 1 i j r t 1 j 1 ; 1 = K(x) (x + t j) 1 : 1 The, we obtai ˆ (j) as the uique positive root of the equatio 2 [2r (j) + K (j) ()] + r (j) [s (j) + K (j) ()]= ad [ s (j) ˆ (j) = + ˆ ] 1=2 (j) 2 ; ˆ (j) r (j) Similarly, we d { [ (s(j) ) ]} 1=2 1=2 (j) = 2 1 ; r (j) (j) =(s (j) r (j) ) 1=2 :

7 H.K.T. Ng et al. / Computatioal Statistics & Data Aalysis 43 (23) Let us ow dee ˆ ( ) = 1 ˆ (j) ; j=1 ˆ( ) = 1 ˆ (j) ; j=1 ( ) = 1 (j) ; j=1 ( ) = 1 (j) : j=1 The, the bias-corrected jackkife maximum likelihood estimates (JMLEs) of ad (see, for example, Efro, 1982) are give by ˆ J = ˆ ( 1) ˆ ( ) ; ˆ J = ˆ ( 1) ˆ ( ) ; similarly, the bias-corrected jackkife modied momet estimates (JMMEs) of ad are give by J = ( 1) ( ) ; J = ( 1) ( ) : 7. Mote Carlo simulatio results I order to compare the performace of all the above estimators, we performed a simulatio study for dieret sample sizes ad for dieret parameter values. We took the sample size as =5; 1; 2; 5; 1, ad the shape parameter as =:1; :25; :5; 1:; 2:. Sice is the scale parameter, was kept xed at 1., without loss of ay geerality. Of course, the values of bias ad stadard deviatio of estimates of simply eed to be multiplied by (if it is dieret from 1) while the bias ad stadard deviatio of estimates of will ot be aected. All the results were based o 1, Mote Carlo rus. We computed the MLE, MME, UMLE, UMME, JMLE ad JMME for each ru, ad the computed the average estimates ad the stadard deviatios over the 1, rus for all these estimators. The results so obtaied are reported i Tables 1 ad 2. We also computed the 9% ad 95% probability coverages of codece itervals based o iferetial quatities associated with all these estimators usig the asymptotic distributios give earlier. Specically, the 1(1 )% codece itervals for ad based o the MLEs ad UMLEs are give by [ ( 1 ( ) ] 1 z=2 z1 =2 ˆ +1) ; ˆ +1 ; 2 2 ( ) 1 ( ) 1 z =2 ˆ h1 (ˆ) +1 ; ˆ z 1 =2 h1 (ˆ) +1 ;

8 29 H.K.T. Ng et al. / Computatioal Statistics & Data Aalysis 43 (23) Table 1 Meas of estimates based o Mote Carlo simulatio ( = 1:) Estimate of Estimate of MLE MME UMLE UMME JMLE JMME MLE MME UMLE UMME JMLE JMME

9 H.K.T. Ng et al. / Computatioal Statistics & Data Aalysis 43 (23) Table 2 Stadard deviatios of estimates based o Mote Carlo simulatio ( = 1:) Estimate of Estimate of MLE MME UMLE UMME JMLE JMME MLE MME UMLE UMME JMLE JMME

10 292 H.K.T. Ng et al. / Computatioal Statistics & Data Aalysis 43 (23) ad [ ( ) 1 ( ) 1 ] ˆ z =2 2 ( 1) +1 ; ˆ z 1 =2 2 ( 1) +1 ; [ ( ) 1 ( ) ] 1 ˆ 4z =2 h 1 (ˆ) (4 +ˆ 2 ) +1 ; ˆ 4z 1 =2 h 1 (ˆ) (4 +ˆ 2 ) +1 ; where h 1 (x)=:25 + x 2 + I(x) ad z p is the 1pth percetile of the stadard ormal distributio. Similarly, the 1(1 )% codece itervals for ad based o the MMEs ad UMMEs are give by [ ( z=2 2 +1) 1 ; ( z1 = ) ] 1 ; ( ) 1 ( ) 1 z =2 h2 ( ) +1 ; z 1 =2 h2 ( ) +1 ; ad [ ( ) 1 ( ) 1 ] z =2 2 ( 1) +1 ; z 1 =2 2 ( 1) +1 ; [ ( ) 1 ( ) ] 1 4z =2 h 2 ( ) (4 + 2 ) +1 ; 4z 1 =2 h 2 ( ) (4 + 2 ) +1 ; where h 2 (x)= x2 ( x2 ) : 2 These results are reported i Tables 3 ad 4, respectively. From the simulatio results, it is clear that the performace of the MLEs ad MMEs are almost idetical for dieret sample sizes ad if the shape parameter is ot too large. The average estimates of the MLEs ad MMEs ad their stadard deviatios coicide upto four decimal places if is less tha.5. It is also evidet from these results that the MLEs ad the MMEs are both highly biased if is small ad is large. The bias reductio method works very well i this case for both the parameters eve for small samples; but as expected, it icreases the correspodig stadard deviatio of the estimators. The performace of the JMLE ad JMME are almost idetical at least for small values of, ad they perform better tha UMLE ad UMME i terms of bias. However, JMLE ad JMME possess larger stadard deviatios tha UMLE ad UMME. I additio, UMLE ad UMME also are computatioally simpler while the jackkifed estimators will demad cosiderably high computatioal time i case of large sample sizes. The asymptotic codece itervals do ot work very well whe the sample size is very small as the coverage probabilities are much lower tha the correspodig omial

11 H.K.T. Ng et al. / Computatioal Statistics & Data Aalysis 43 (23) Table 3 Probability coverages of 9% codece itervals based o Mote Carlo simulatio ( = 1:) Probability coverages for Probability coverages for MLE MME UMLE UMME MLE MME UMLE UMME levels. But for sample sizes 2 or more, the performaces are quite satisfactory for codece itervals for both ad. The bias reductio techique deitely helps to improve the coverage probabilities i both cases to a certai extet. 8. Illustrative examples Practical applicatio of the above estimators is illustrated here with two examples with oe ivolvig a large sample ad the other with a small sample. Example 1. The data set is give by Birbaum ad Sauders (1969b) o the fatigue life of 661-T6 alumium coupos cut parallel to the directio of rollig ad oscillated at 18 cycles per secod (cps). The data set cosists of 11 observatios with maximum stress per cycle 31; psi. The data are preseted i Table 5.

12 294 H.K.T. Ng et al. / Computatioal Statistics & Data Aalysis 43 (23) Table 4 Probability coverages of 95% codece itervals based o Mote Carlo simulatio ( = 1:) Probability coverages for Probability coverages for MLE MME UMLE UMME MLE MME UMLE UMME Table 5 Fatigue lifetime data preseted by Birbaum ad Sauders (1969b) I summary, we have i this case =11, s=133:73267, ad r =129: For this example, the poit estimates of ad obtaied by all the methods are summarized i Table 6.

13 H.K.T. Ng et al. / Computatioal Statistics & Data Aalysis 43 (23) Table 6 Poit estimates of ad for Example 1 Estimator MLE MME UMLE UMME JMLE JMME Table 7 Stadard deviatios of estimates ad iterval estimates of ad for Example 1 Estimator SD 9% CI 95% CI SD 9% CI 95% CI MLE.12 (.1527,.1927) (.1497,.1976) ( , ) ( , ) MME.12 (.1527,.1927) (.1497,.1976) ( , ) ( , ) UMLE.122 (.1541,.1949) (.1511,.1999) ( , ) ( , ) UMME.122 (.1541,.1949) (.1511,.1999) ( , ) ( , ) Table 8 Fatigue lifetime data preseted by McCool (1974) From Eqs. (1), (13), (16) ad (17), the asymptotic variace of the estimators ca be readily obtaied, ad also the codece itervals for ad based o the MLEs, MMEs, UMLEs ad UMMEs ca be readily costructed usig the asymptotic ormality. The results so obtaied are preseted i Table 7. Example 2. This example is from McCool (1974) o the fatigue life i hours of 1 bearigs of a certai type. These data were used as a illustrative example for the three-parameter Weibull distributio by Cohe et al. (1984). The data are preseted i Table 8. I this case, we d = 1; s = 22:48 ad r = 23:8853. The poit estimates ad iterval estimates obtaied from this data are summarized i Tables 9 ad 1, respectively. It is hearteig to observe i both these examples that the MMEs are very early the same as the MLEs ad also the correspodig codece itervals. Also, while the ubiased estimators correct for the bias, they do result i a larger stadard deviatio.

14 296 H.K.T. Ng et al. / Computatioal Statistics & Data Aalysis 43 (23) Table 9 Poit estimates of ad for Example 2 Estimator MLE MME UMLE UMME JMLE JMME Table 1 Stadard deviatios of estimates ad iterval estimates of ad for Example 2 Estimator SD 9% CI 95% CI SD 9% CI 95% CI MLE.632 (.265,.4468) (.1964,.529) ( , ) ( , ) MME.632 (.265,.4468) (.1964,.529) ( , ) ( , ) UMLE.78 (.2228,.538) (.2111,.6118) ( , ) ( , ) UMME.78 (.2228,.538) (.2111,.6118) ( , ) ( , ) 9. Cocludig remarks Sice the maximum likelihood estimators ad the modied momet estimators behave very similarly i almost all cases cosidered, we recommed the use of the modied momet estimators because they are explicit estimators ad are very easy to compute. Although jackkife estimators also work very well, they caot be recommeded for large sample sizes. If the sample size is very small (say, less tha 1), the bias corrected modied momet estimators or bias corrected jackkife estimators should be used as i this case the origial estimators ca be highly biased. The asymptotic codece itervals behave very well for large sample sizes (at least 2), but ot for small sample sizes. I the latter case, oe may rely o simulated percetage poits rather tha o the asymptotic ormality. Appedix. Derivatio of the asymptotic distributio of the modied momet estimators Let T; T 1 ;:::;T be idepedet ad idetically distributed Birbaum Sauders radom variables with desity fuctio as i (2). Let us dee the radom variables [ ] 1 S = 1 1 T i ad R = Ti 1 : From the strog law of large umbers, it is kow that S ad R 1 coverge almost surely to E(T ) ad E(T 1 ), respectively. Also from the cetral limit theorem, we

15 H.K.T. Ng et al. / Computatioal Statistics & Data Aalysis 43 (23) readily observe that S ad R 1 are asymptotically ormally distributed; furthermore, ay liear combiatio of the form as + br 1 = 1 [at i + bti 1 ] is also asymptotically ormally distributed for all a ad b. Therefore, (S R 1 ) is asymptotically distributed as bivariate ormal. Specially, we have ( ) [( ) ] S E(T ) N ; ; R 1 E(T 1 ) where ( ) = ; = Var(T )=() 2 ( ); 12 = 21 = Cov(T; T 1 )=E(1) E(T)E(T 1 )=1 ( ) 2 ; 22 = Var(T 1 )= 2 2 ( ): We ow eed to d the asymptotic joit distributio of ( ) [ ] f1 (S; R) = ; f 2 (S; R) where { [ ( ) 1=2 1=2 x f 1 (x; y)= 2 1]} y ad f 2 (x; y)=(xy) 1=2 : By usig the Taylor series expasio, we obtai ( ) [( ) ] N ; ; @y = ( () ( ) ) : T x=e(t );y=e(t 1 )

16 298 H.K.T. Ng et al. / Computatioal Statistics & Data Aalysis 43 (23) Of course, we are usig here the followig 1 = x=e(t );y=e(t 1 ) = x=e(t );y=e(t 1 ) 2 ; x=e(t );y=e(t 1 ) x=e(t );y=e(t 1 ) = ; = : Refereces Birbaum, Z.W., Sauders, S.C., 1969a. A ew family of life distributio. J. Appl. Probab. 6, Birbaum, Z.W., Sauders, S.C., 1969b. Estimatio for a family of life distributios with applicatios to fatigue. J. Appl. Probab. 6, Chag, D.S., Tag, L.C., Reliability bouds ad critical time for the Birbaum Sauders distributio. IEEE Tras. Reliab. 42, Chag, D.S., Tag, L.C., Percetile bouds ad tolerace limits for the Birbaum Sauders distributio. Comm. Statist. Theory Methods 23, Cohe, A.C., Whitte, B.J., Dig, Y., Modied momet estimatio for the three-parameter Weibull distributio. J. Qual. Techol. 16, Desmod, A.F., Stochastic models of failure i radom eviromets. Caad. J. Statist. 13, Desmod, A.F., O the relatioship betwee two fatigue-life models. IEEE Tras. Reliab. 35, Dupuis, D.J., Mills, J.E., Robust estimatio of the Birbaum Sauders distributio. IEEE Tras. Reliab. 47, Efro, B., The Jackkife, the Bootstrap ad Other Resamplig Plas. Society for Idustrial ad Applied Mathematics, Philadelphia. Egelhardt, M., Bai, L.J., Wright, F.T., Ifereces o the parameters of the Birbaum Sauders fatigue life distributio based o maximum likelihood estimatio. Techometrics 23, Johso, N.L., Kotz, S., Balakrisha, N., Cotiuous Uivariate Distributios Vol. 2, 2d Editio. Wiley, New York. Ma, N.R., Schafer, R.E., Sigpurwalla, N.D., Methods for Statistical Aalysis of Reliability ad Life Data. Wiley, New York. McCool, J.I., Iferetial techiques for Weibull populatios. Aerospace Research Laboratories Report ARL TR Wright Patterso Air Force Base, Dayto, OH. Rieck, J.R., Parametric estimatio for the Birbaum Sauders distributio based o symmetrically cesored samples. Comm. Statist. Theory Methods 24, Rieck, J.R., A momet-geeratig fuctio with applicatio to the Birbaum Sauders distributio. Comm. Statist. Theory Methods 28,

MOMENT-METHOD ESTIMATION BASED ON CENSORED SAMPLE

MOMENT-METHOD ESTIMATION BASED ON CENSORED SAMPLE Vol. 8 o. Joural of Systems Sciece ad Complexity Apr., 5 MOMET-METHOD ESTIMATIO BASED O CESORED SAMPLE I Zhogxi Departmet of Mathematics, East Chia Uiversity of Sciece ad Techology, Shaghai 37, Chia. Email:

More information

Bootstrap Intervals of the Parameters of Lognormal Distribution Using Power Rule Model and Accelerated Life Tests

Bootstrap Intervals of the Parameters of Lognormal Distribution Using Power Rule Model and Accelerated Life Tests Joural of Moder Applied Statistical Methods Volume 5 Issue Article --5 Bootstrap Itervals of the Parameters of Logormal Distributio Usig Power Rule Model ad Accelerated Life Tests Mohammed Al-Ha Ebrahem

More information

Resampling Methods. X (1/2), i.e., Pr (X i m) = 1/2. We order the data: X (1) X (2) X (n). Define the sample median: ( n.

Resampling Methods. X (1/2), i.e., Pr (X i m) = 1/2. We order the data: X (1) X (2) X (n). Define the sample median: ( n. Jauary 1, 2019 Resamplig Methods Motivatio We have so may estimators with the property θ θ d N 0, σ 2 We ca also write θ a N θ, σ 2 /, where a meas approximately distributed as Oce we have a cosistet estimator

More information

Confidence Interval for Standard Deviation of Normal Distribution with Known Coefficients of Variation

Confidence Interval for Standard Deviation of Normal Distribution with Known Coefficients of Variation Cofidece Iterval for tadard Deviatio of Normal Distributio with Kow Coefficiets of Variatio uparat Niwitpog Departmet of Applied tatistics, Faculty of Applied ciece Kig Mogkut s Uiversity of Techology

More information

R. van Zyl 1, A.J. van der Merwe 2. Quintiles International, University of the Free State

R. van Zyl 1, A.J. van der Merwe 2. Quintiles International, University of the Free State Bayesia Cotrol Charts for the Two-parameter Expoetial Distributio if the Locatio Parameter Ca Take o Ay Value Betwee Mius Iity ad Plus Iity R. va Zyl, A.J. va der Merwe 2 Quitiles Iteratioal, ruaavz@gmail.com

More information

Maximum likelihood estimation from record-breaking data for the generalized Pareto distribution

Maximum likelihood estimation from record-breaking data for the generalized Pareto distribution METRON - Iteratioal Joural of Statistics 004, vol. LXII,. 3, pp. 377-389 NAGI S. ABD-EL-HAKIM KHALAF S. SULTAN Maximum likelihood estimatio from record-breakig data for the geeralized Pareto distributio

More information

Confidence interval for the two-parameter exponentiated Gumbel distribution based on record values

Confidence interval for the two-parameter exponentiated Gumbel distribution based on record values Iteratioal Joural of Applied Operatioal Research Vol. 4 No. 1 pp. 61-68 Witer 2014 Joural homepage: www.ijorlu.ir Cofidece iterval for the two-parameter expoetiated Gumbel distributio based o record values

More information

Topic 9: Sampling Distributions of Estimators

Topic 9: Sampling Distributions of Estimators Topic 9: Samplig Distributios of Estimators Course 003, 2016 Page 0 Samplig distributios of estimators Sice our estimators are statistics (particular fuctios of radom variables), their distributio ca be

More information

Goodness-Of-Fit For The Generalized Exponential Distribution. Abstract

Goodness-Of-Fit For The Generalized Exponential Distribution. Abstract Goodess-Of-Fit For The Geeralized Expoetial Distributio By Amal S. Hassa stitute of Statistical Studies & Research Cairo Uiversity Abstract Recetly a ew distributio called geeralized expoetial or expoetiated

More information

Lecture 7: Properties of Random Samples

Lecture 7: Properties of Random Samples Lecture 7: Properties of Radom Samples 1 Cotiued From Last Class Theorem 1.1. Let X 1, X,...X be a radom sample from a populatio with mea µ ad variace σ

More information

Statistical Inference (Chapter 10) Statistical inference = learn about a population based on the information provided by a sample.

Statistical Inference (Chapter 10) Statistical inference = learn about a population based on the information provided by a sample. Statistical Iferece (Chapter 10) Statistical iferece = lear about a populatio based o the iformatio provided by a sample. Populatio: The set of all values of a radom variable X of iterest. Characterized

More information

Approximate Confidence Interval for the Reciprocal of a Normal Mean with a Known Coefficient of Variation

Approximate Confidence Interval for the Reciprocal of a Normal Mean with a Known Coefficient of Variation Metodološki zvezki, Vol. 13, No., 016, 117-130 Approximate Cofidece Iterval for the Reciprocal of a Normal Mea with a Kow Coefficiet of Variatio Wararit Paichkitkosolkul 1 Abstract A approximate cofidece

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

1 Inferential Methods for Correlation and Regression Analysis

1 Inferential Methods for Correlation and Regression Analysis 1 Iferetial Methods for Correlatio ad Regressio Aalysis I the chapter o Correlatio ad Regressio Aalysis tools for describig bivariate cotiuous data were itroduced. The sample Pearso Correlatio Coefficiet

More information

Mathematical Modeling of Optimum 3 Step Stress Accelerated Life Testing for Generalized Pareto Distribution

Mathematical Modeling of Optimum 3 Step Stress Accelerated Life Testing for Generalized Pareto Distribution America Joural of Theoretical ad Applied Statistics 05; 4(: 6-69 Published olie May 8, 05 (http://www.sciecepublishiggroup.com/j/ajtas doi: 0.648/j.ajtas.05040. ISSN: 6-8999 (Prit; ISSN: 6-9006 (Olie Mathematical

More information

A Note on Box-Cox Quantile Regression Estimation of the Parameters of the Generalized Pareto Distribution

A Note on Box-Cox Quantile Regression Estimation of the Parameters of the Generalized Pareto Distribution A Note o Box-Cox Quatile Regressio Estimatio of the Parameters of the Geeralized Pareto Distributio JM va Zyl Abstract: Makig use of the quatile equatio, Box-Cox regressio ad Laplace distributed disturbaces,

More information

Topic 9: Sampling Distributions of Estimators

Topic 9: Sampling Distributions of Estimators Topic 9: Samplig Distributios of Estimators Course 003, 2018 Page 0 Samplig distributios of estimators Sice our estimators are statistics (particular fuctios of radom variables), their distributio ca be

More information

Department of Mathematics

Department of Mathematics Departmet of Mathematics Ma 3/103 KC Border Itroductio to Probability ad Statistics Witer 2017 Lecture 19: Estimatio II Relevat textbook passages: Larse Marx [1]: Sectios 5.2 5.7 19.1 The method of momets

More information

Bayesian Control Charts for the Two-parameter Exponential Distribution

Bayesian Control Charts for the Two-parameter Exponential Distribution Bayesia Cotrol Charts for the Two-parameter Expoetial Distributio R. va Zyl, A.J. va der Merwe 2 Quitiles Iteratioal, ruaavz@gmail.com 2 Uiversity of the Free State Abstract By usig data that are the mileages

More information

Access to the published version may require journal subscription. Published with permission from: Elsevier.

Access to the published version may require journal subscription. Published with permission from: Elsevier. This is a author produced versio of a paper published i Statistics ad Probability Letters. This paper has bee peer-reviewed, it does ot iclude the joural pagiatio. Citatio for the published paper: Forkma,

More information

A goodness-of-fit test based on the empirical characteristic function and a comparison of tests for normality

A goodness-of-fit test based on the empirical characteristic function and a comparison of tests for normality A goodess-of-fit test based o the empirical characteristic fuctio ad a compariso of tests for ormality J. Marti va Zyl Departmet of Mathematical Statistics ad Actuarial Sciece, Uiversity of the Free State,

More information

MATH 320: Probability and Statistics 9. Estimation and Testing of Parameters. Readings: Pruim, Chapter 4

MATH 320: Probability and Statistics 9. Estimation and Testing of Parameters. Readings: Pruim, Chapter 4 MATH 30: Probability ad Statistics 9. Estimatio ad Testig of Parameters Estimatio ad Testig of Parameters We have bee dealig situatios i which we have full kowledge of the distributio of a radom variable.

More information

ECONOMETRIC THEORY. MODULE XIII Lecture - 34 Asymptotic Theory and Stochastic Regressors

ECONOMETRIC THEORY. MODULE XIII Lecture - 34 Asymptotic Theory and Stochastic Regressors ECONOMETRIC THEORY MODULE XIII Lecture - 34 Asymptotic Theory ad Stochastic Regressors Dr. Shalabh Departmet of Mathematics ad Statistics Idia Istitute of Techology Kapur Asymptotic theory The asymptotic

More information

Lecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting

Lecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting Lecture 6 Chi Square Distributio (χ ) ad Least Squares Fittig Chi Square Distributio (χ ) Suppose: We have a set of measuremets {x 1, x, x }. We kow the true value of each x i (x t1, x t, x t ). We would

More information

Estimation of Gumbel Parameters under Ranked Set Sampling

Estimation of Gumbel Parameters under Ranked Set Sampling Joural of Moder Applied Statistical Methods Volume 13 Issue 2 Article 11-2014 Estimatio of Gumbel Parameters uder Raked Set Samplig Omar M. Yousef Al Balqa' Applied Uiversity, Zarqa, Jorda, abuyaza_o@yahoo.com

More information

Lecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting

Lecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting Lecture 6 Chi Square Distributio (χ ) ad Least Squares Fittig Chi Square Distributio (χ ) Suppose: We have a set of measuremets {x 1, x, x }. We kow the true value of each x i (x t1, x t, x t ). We would

More information

Topic 9: Sampling Distributions of Estimators

Topic 9: Sampling Distributions of Estimators Topic 9: Samplig Distributios of Estimators Course 003, 2018 Page 0 Samplig distributios of estimators Sice our estimators are statistics (particular fuctios of radom variables), their distributio ca be

More information

ENGI 4421 Confidence Intervals (Two Samples) Page 12-01

ENGI 4421 Confidence Intervals (Two Samples) Page 12-01 ENGI 44 Cofidece Itervals (Two Samples) Page -0 Two Sample Cofidece Iterval for a Differece i Populatio Meas [Navidi sectios 5.4-5.7; Devore chapter 9] From the cetral limit theorem, we kow that, for sufficietly

More information

Estimation for Complete Data

Estimation for Complete Data Estimatio for Complete Data complete data: there is o loss of iformatio durig study. complete idividual complete data= grouped data A complete idividual data is the oe i which the complete iformatio of

More information

Lecture 2: Monte Carlo Simulation

Lecture 2: Monte Carlo Simulation STAT/Q SCI 43: Itroductio to Resamplig ethods Sprig 27 Istructor: Ye-Chi Che Lecture 2: ote Carlo Simulatio 2 ote Carlo Itegratio Assume we wat to evaluate the followig itegratio: e x3 dx What ca we do?

More information

Lecture 19: Convergence

Lecture 19: Convergence Lecture 19: Covergece Asymptotic approach I statistical aalysis or iferece, a key to the success of fidig a good procedure is beig able to fid some momets ad/or distributios of various statistics. I may

More information

x = Pr ( X (n) βx ) =

x = Pr ( X (n) βx ) = Exercise 93 / page 45 The desity of a variable X i i 1 is fx α α a For α kow let say equal to α α > fx α α x α Pr X i x < x < Usig a Pivotal Quatity: x α 1 < x < α > x α 1 ad We solve i a similar way as

More information

Stat 421-SP2012 Interval Estimation Section

Stat 421-SP2012 Interval Estimation Section Stat 41-SP01 Iterval Estimatio Sectio 11.1-11. We ow uderstad (Chapter 10) how to fid poit estimators of a ukow parameter. o However, a poit estimate does ot provide ay iformatio about the ucertaity (possible

More information

1 Introduction to reducing variance in Monte Carlo simulations

1 Introduction to reducing variance in Monte Carlo simulations Copyright c 010 by Karl Sigma 1 Itroductio to reducig variace i Mote Carlo simulatios 11 Review of cofidece itervals for estimatig a mea I statistics, we estimate a ukow mea µ = E(X) of a distributio by

More information

It should be unbiased, or approximately unbiased. Variance of the variance estimator should be small. That is, the variance estimator is stable.

It should be unbiased, or approximately unbiased. Variance of the variance estimator should be small. That is, the variance estimator is stable. Chapter 10 Variace Estimatio 10.1 Itroductio Variace estimatio is a importat practical problem i survey samplig. Variace estimates are used i two purposes. Oe is the aalytic purpose such as costructig

More information

POWER AKASH DISTRIBUTION AND ITS APPLICATION

POWER AKASH DISTRIBUTION AND ITS APPLICATION POWER AKASH DISTRIBUTION AND ITS APPLICATION Rama SHANKER PhD, Uiversity Professor, Departmet of Statistics, College of Sciece, Eritrea Istitute of Techology, Asmara, Eritrea E-mail: shakerrama009@gmail.com

More information

EECS564 Estimation, Filtering, and Detection Hwk 2 Solns. Winter p θ (z) = (2θz + 1 θ), 0 z 1

EECS564 Estimation, Filtering, and Detection Hwk 2 Solns. Winter p θ (z) = (2θz + 1 θ), 0 z 1 EECS564 Estimatio, Filterig, ad Detectio Hwk 2 Sols. Witer 25 4. Let Z be a sigle observatio havig desity fuctio where. p (z) = (2z + ), z (a) Assumig that is a oradom parameter, fid ad plot the maximum

More information

Limit distributions for products of sums

Limit distributions for products of sums Statistics & Probability Letters 62 (23) 93 Limit distributios for products of sums Yogcheg Qi Departmet of Mathematics ad Statistics, Uiversity of Miesota-Duluth, Campus Ceter 4, 7 Uiversity Drive, Duluth,

More information

Expectation and Variance of a random variable

Expectation and Variance of a random variable Chapter 11 Expectatio ad Variace of a radom variable The aim of this lecture is to defie ad itroduce mathematical Expectatio ad variace of a fuctio of discrete & cotiuous radom variables ad the distributio

More information

Basis for simulation techniques

Basis for simulation techniques Basis for simulatio techiques M. Veeraraghava, March 7, 004 Estimatio is based o a collectio of experimetal outcomes, x, x,, x, where each experimetal outcome is a value of a radom variable. x i. Defiitios

More information

The standard deviation of the mean

The standard deviation of the mean Physics 6C Fall 20 The stadard deviatio of the mea These otes provide some clarificatio o the distictio betwee the stadard deviatio ad the stadard deviatio of the mea.. The sample mea ad variace Cosider

More information

Direction: This test is worth 250 points. You are required to complete this test within 50 minutes.

Direction: This test is worth 250 points. You are required to complete this test within 50 minutes. Term Test October 3, 003 Name Math 56 Studet Number Directio: This test is worth 50 poits. You are required to complete this test withi 50 miutes. I order to receive full credit, aswer each problem completely

More information

Since X n /n P p, we know that X n (n. Xn (n X n ) Using the asymptotic result above to obtain an approximation for fixed n, we obtain

Since X n /n P p, we know that X n (n. Xn (n X n ) Using the asymptotic result above to obtain an approximation for fixed n, we obtain Assigmet 9 Exercise 5.5 Let X biomial, p, where p 0, 1 is ukow. Obtai cofidece itervals for p i two differet ways: a Sice X / p d N0, p1 p], the variace of the limitig distributio depeds oly o p. Use the

More information

A new distribution-free quantile estimator

A new distribution-free quantile estimator Biometrika (1982), 69, 3, pp. 635-40 Prited i Great Britai 635 A ew distributio-free quatile estimator BY FRANK E. HARRELL Cliical Biostatistics, Duke Uiversity Medical Ceter, Durham, North Carolia, U.S.A.

More information

Goodness-Of-Fit For The Generalized Exponential Distribution. Abstract

Goodness-Of-Fit For The Generalized Exponential Distribution. Abstract Goodess-Of-Fit For The Geeralized Expoetial Distributio By Amal S. Hassa stitute of Statistical Studies & Research Cairo Uiversity Abstract Recetly a ew distributio called geeralized expoetial or expoetiated

More information

ANOTHER WEIGHTED WEIBULL DISTRIBUTION FROM AZZALINI S FAMILY

ANOTHER WEIGHTED WEIBULL DISTRIBUTION FROM AZZALINI S FAMILY ANOTHER WEIGHTED WEIBULL DISTRIBUTION FROM AZZALINI S FAMILY Sulema Nasiru, MSc. Departmet of Statistics, Faculty of Mathematical Scieces, Uiversity for Developmet Studies, Navrogo, Upper East Regio, Ghaa,

More information

A Generalized Gamma-Weibull Distribution: Model, Properties and Applications

A Generalized Gamma-Weibull Distribution: Model, Properties and Applications Marquette Uiversity e-publicatios@marquette Mathematics, Statistics ad Computer Sciece Faculty Research ad Publicatios Mathematics, Statistics ad Computer Sciece, Departmet of --06 A Geeralized Gamma-Weibull

More information

CEE 522 Autumn Uncertainty Concepts for Geotechnical Engineering

CEE 522 Autumn Uncertainty Concepts for Geotechnical Engineering CEE 5 Autum 005 Ucertaity Cocepts for Geotechical Egieerig Basic Termiology Set A set is a collectio of (mutually exclusive) objects or evets. The sample space is the (collectively exhaustive) collectio

More information

G. R. Pasha Department of Statistics Bahauddin Zakariya University Multan, Pakistan

G. R. Pasha Department of Statistics Bahauddin Zakariya University Multan, Pakistan Deviatio of the Variaces of Classical Estimators ad Negative Iteger Momet Estimator from Miimum Variace Boud with Referece to Maxwell Distributio G. R. Pasha Departmet of Statistics Bahauddi Zakariya Uiversity

More information

Probability and statistics: basic terms

Probability and statistics: basic terms Probability ad statistics: basic terms M. Veeraraghava August 203 A radom variable is a rule that assigs a umerical value to each possible outcome of a experimet. Outcomes of a experimet form the sample

More information

7.1 Convergence of sequences of random variables

7.1 Convergence of sequences of random variables Chapter 7 Limit Theorems Throughout this sectio we will assume a probability space (, F, P), i which is defied a ifiite sequece of radom variables (X ) ad a radom variable X. The fact that for every ifiite

More information

Random Variables, Sampling and Estimation

Random Variables, Sampling and Estimation Chapter 1 Radom Variables, Samplig ad Estimatio 1.1 Itroductio This chapter will cover the most importat basic statistical theory you eed i order to uderstad the ecoometric material that will be comig

More information

Record Values from T-X Family of. Pareto-Exponential Distribution with. Properties and Simulations

Record Values from T-X Family of. Pareto-Exponential Distribution with. Properties and Simulations Applied Mathematical Scieces, Vol. 3, 209, o., 33-44 HIKARI Ltd, www.m-hikari.com https://doi.org/0.2988/ams.209.879 Record Values from T-X Family of Pareto-Epoetial Distributio with Properties ad Simulatios

More information

Discrete Orthogonal Moment Features Using Chebyshev Polynomials

Discrete Orthogonal Moment Features Using Chebyshev Polynomials Discrete Orthogoal Momet Features Usig Chebyshev Polyomials R. Mukuda, 1 S.H.Og ad P.A. Lee 3 1 Faculty of Iformatio Sciece ad Techology, Multimedia Uiversity 75450 Malacca, Malaysia. Istitute of Mathematical

More information

MATH/STAT 352: Lecture 15

MATH/STAT 352: Lecture 15 MATH/STAT 352: Lecture 15 Sectios 5.2 ad 5.3. Large sample CI for a proportio ad small sample CI for a mea. 1 5.2: Cofidece Iterval for a Proportio Estimatig proportio of successes i a biomial experimet

More information

Parameter, Statistic and Random Samples

Parameter, Statistic and Random Samples Parameter, Statistic ad Radom Samples A parameter is a umber that describes the populatio. It is a fixed umber, but i practice we do ot kow its value. A statistic is a fuctio of the sample data, i.e.,

More information

Econ 325/327 Notes on Sample Mean, Sample Proportion, Central Limit Theorem, Chi-square Distribution, Student s t distribution 1.

Econ 325/327 Notes on Sample Mean, Sample Proportion, Central Limit Theorem, Chi-square Distribution, Student s t distribution 1. Eco 325/327 Notes o Sample Mea, Sample Proportio, Cetral Limit Theorem, Chi-square Distributio, Studet s t distributio 1 Sample Mea By Hiro Kasahara We cosider a radom sample from a populatio. Defiitio

More information

Direction: This test is worth 150 points. You are required to complete this test within 55 minutes.

Direction: This test is worth 150 points. You are required to complete this test within 55 minutes. Term Test 3 (Part A) November 1, 004 Name Math 6 Studet Number Directio: This test is worth 10 poits. You are required to complete this test withi miutes. I order to receive full credit, aswer each problem

More information

A LARGER SAMPLE SIZE IS NOT ALWAYS BETTER!!!

A LARGER SAMPLE SIZE IS NOT ALWAYS BETTER!!! A LARGER SAMLE SIZE IS NOT ALWAYS BETTER!!! Nagaraj K. Neerchal Departmet of Mathematics ad Statistics Uiversity of Marylad Baltimore Couty, Baltimore, MD 2250 Herbert Lacayo ad Barry D. Nussbaum Uited

More information

Properties and Hypothesis Testing

Properties and Hypothesis Testing Chapter 3 Properties ad Hypothesis Testig 3.1 Types of data The regressio techiques developed i previous chapters ca be applied to three differet kids of data. 1. Cross-sectioal data. 2. Time series data.

More information

Investigating the Significance of a Correlation Coefficient using Jackknife Estimates

Investigating the Significance of a Correlation Coefficient using Jackknife Estimates Iteratioal Joural of Scieces: Basic ad Applied Research (IJSBAR) ISSN 2307-4531 (Prit & Olie) http://gssrr.org/idex.php?joural=jouralofbasicadapplied ---------------------------------------------------------------------------------------------------------------------------

More information

Lecture 33: Bootstrap

Lecture 33: Bootstrap Lecture 33: ootstrap Motivatio To evaluate ad compare differet estimators, we eed cosistet estimators of variaces or asymptotic variaces of estimators. This is also importat for hypothesis testig ad cofidece

More information

Confidence Intervals For P(X less than Y) In The Exponential Case With Common Location Parameter

Confidence Intervals For P(X less than Y) In The Exponential Case With Common Location Parameter Joural of Moder Applied Statistical Methods Volume Issue Article 7 --3 Cofidece Itervals For P(X less tha Y I he Expoetial Case With Commo Locatio Parameter Ayma Baklizi Yarmouk Uiversity, Irbid, Jorda,

More information

Statistics 511 Additional Materials

Statistics 511 Additional Materials Cofidece Itervals o mu Statistics 511 Additioal Materials This topic officially moves us from probability to statistics. We begi to discuss makig ifereces about the populatio. Oe way to differetiate probability

More information

Chapter 2 The Monte Carlo Method

Chapter 2 The Monte Carlo Method Chapter 2 The Mote Carlo Method The Mote Carlo Method stads for a broad class of computatioal algorithms that rely o radom sampligs. It is ofte used i physical ad mathematical problems ad is most useful

More information

Using the IML Procedure to Examine the Efficacy of a New Control Charting Technique

Using the IML Procedure to Examine the Efficacy of a New Control Charting Technique Paper 2894-2018 Usig the IML Procedure to Examie the Efficacy of a New Cotrol Chartig Techique Austi Brow, M.S., Uiversity of Norther Colorado; Bryce Whitehead, M.S., Uiversity of Norther Colorado ABSTRACT

More information

Some Properties of the Exact and Score Methods for Binomial Proportion and Sample Size Calculation

Some Properties of the Exact and Score Methods for Binomial Proportion and Sample Size Calculation Some Properties of the Exact ad Score Methods for Biomial Proportio ad Sample Size Calculatio K. KRISHNAMOORTHY AND JIE PENG Departmet of Mathematics, Uiversity of Louisiaa at Lafayette Lafayette, LA 70504-1010,

More information

Chapter 6 Principles of Data Reduction

Chapter 6 Principles of Data Reduction Chapter 6 for BST 695: Special Topics i Statistical Theory. Kui Zhag, 0 Chapter 6 Priciples of Data Reductio Sectio 6. Itroductio Goal: To summarize or reduce the data X, X,, X to get iformatio about a

More information

EXAMINATIONS OF THE ROYAL STATISTICAL SOCIETY

EXAMINATIONS OF THE ROYAL STATISTICAL SOCIETY EXAMINATIONS OF THE ROYAL STATISTICAL SOCIETY GRADUATE DIPLOMA, 016 MODULE : Statistical Iferece Time allowed: Three hours Cadidates should aswer FIVE questios. All questios carry equal marks. The umber

More information

Control Charts for Mean for Non-Normally Correlated Data

Control Charts for Mean for Non-Normally Correlated Data Joural of Moder Applied Statistical Methods Volume 16 Issue 1 Article 5 5-1-017 Cotrol Charts for Mea for No-Normally Correlated Data J. R. Sigh Vikram Uiversity, Ujjai, Idia Ab Latif Dar School of Studies

More information

Economics 241B Relation to Method of Moments and Maximum Likelihood OLSE as a Maximum Likelihood Estimator

Economics 241B Relation to Method of Moments and Maximum Likelihood OLSE as a Maximum Likelihood Estimator Ecoomics 24B Relatio to Method of Momets ad Maximum Likelihood OLSE as a Maximum Likelihood Estimator Uder Assumptio 5 we have speci ed the distributio of the error, so we ca estimate the model parameters

More information

Distribution of Random Samples & Limit theorems

Distribution of Random Samples & Limit theorems STAT/MATH 395 A - PROBABILITY II UW Witer Quarter 2017 Néhémy Lim Distributio of Radom Samples & Limit theorems 1 Distributio of i.i.d. Samples Motivatig example. Assume that the goal of a study is to

More information

Statistical inference: example 1. Inferential Statistics

Statistical inference: example 1. Inferential Statistics Statistical iferece: example 1 Iferetial Statistics POPULATION SAMPLE A clothig store chai regularly buys from a supplier large quatities of a certai piece of clothig. Each item ca be classified either

More information

On Marshall-Olkin Extended Weibull Distribution

On Marshall-Olkin Extended Weibull Distribution Joural of Statistical Theory ad Applicatios, Vol. 6, No. March 27) 7 O Marshall-Olki Exteded Weibull Distributio Haa Haj Ahmad Departmet of Mathematics, Uiversity of Hail Hail, KSA haaahm@yahoo.com Omar

More information

Output Analysis and Run-Length Control

Output Analysis and Run-Length Control IEOR E4703: Mote Carlo Simulatio Columbia Uiversity c 2017 by Marti Haugh Output Aalysis ad Ru-Legth Cotrol I these otes we describe how the Cetral Limit Theorem ca be used to costruct approximate (1 α%

More information

LECTURE 8: ASYMPTOTICS I

LECTURE 8: ASYMPTOTICS I LECTURE 8: ASYMPTOTICS I We are iterested i the properties of estimators as. Cosider a sequece of radom variables {, X 1}. N. M. Kiefer, Corell Uiversity, Ecoomics 60 1 Defiitio: (Weak covergece) A sequece

More information

Exponential Families and Bayesian Inference

Exponential Families and Bayesian Inference Computer Visio Expoetial Families ad Bayesia Iferece Lecture Expoetial Families A expoetial family of distributios is a d-parameter family f(x; havig the followig form: f(x; = h(xe g(t T (x B(, (. where

More information

Lecture 3. Properties of Summary Statistics: Sampling Distribution

Lecture 3. Properties of Summary Statistics: Sampling Distribution Lecture 3 Properties of Summary Statistics: Samplig Distributio Mai Theme How ca we use math to justify that our umerical summaries from the sample are good summaries of the populatio? Lecture Summary

More information

Statistical Intervals for a Single Sample

Statistical Intervals for a Single Sample 3/5/06 Applied Statistics ad Probability for Egieers Sixth Editio Douglas C. Motgomery George C. Ruger Chapter 8 Statistical Itervals for a Sigle Sample 8 CHAPTER OUTLINE 8- Cofidece Iterval o the Mea

More information

The estimation of the m parameter of the Nakagami distribution

The estimation of the m parameter of the Nakagami distribution The estimatio of the m parameter of the Nakagami distributio Li-Fei Huag Mig Chua Uiversity Dept. of Applied Stat. ad Ifo. Sciece 5 Teh-Mig Rd., Gwei-Sha Taoyua City Taiwa lhuag@mail.mcu.edu.tw Je-Je Li

More information

Study the bias (due to the nite dimensional approximation) and variance of the estimators

Study the bias (due to the nite dimensional approximation) and variance of the estimators 2 Series Methods 2. Geeral Approach A model has parameters (; ) where is ite-dimesioal ad is oparametric. (Sometimes, there is o :) We will focus o regressio. The fuctio is approximated by a series a ite

More information

DS 100: Principles and Techniques of Data Science Date: April 13, Discussion #10

DS 100: Principles and Techniques of Data Science Date: April 13, Discussion #10 DS 00: Priciples ad Techiques of Data Sciece Date: April 3, 208 Name: Hypothesis Testig Discussio #0. Defie these terms below as they relate to hypothesis testig. a) Data Geeratio Model: Solutio: A set

More information

CS434a/541a: Pattern Recognition Prof. Olga Veksler. Lecture 5

CS434a/541a: Pattern Recognition Prof. Olga Veksler. Lecture 5 CS434a/54a: Patter Recogitio Prof. Olga Veksler Lecture 5 Today Itroductio to parameter estimatio Two methods for parameter estimatio Maimum Likelihood Estimatio Bayesia Estimatio Itroducto Bayesia Decisio

More information

A statistical method to determine sample size to estimate characteristic value of soil parameters

A statistical method to determine sample size to estimate characteristic value of soil parameters A statistical method to determie sample size to estimate characteristic value of soil parameters Y. Hojo, B. Setiawa 2 ad M. Suzuki 3 Abstract Sample size is a importat factor to be cosidered i determiig

More information

4. Partial Sums and the Central Limit Theorem

4. Partial Sums and the Central Limit Theorem 1 of 10 7/16/2009 6:05 AM Virtual Laboratories > 6. Radom Samples > 1 2 3 4 5 6 7 4. Partial Sums ad the Cetral Limit Theorem The cetral limit theorem ad the law of large umbers are the two fudametal theorems

More information

CONTROL CHARTS FOR THE LOGNORMAL DISTRIBUTION

CONTROL CHARTS FOR THE LOGNORMAL DISTRIBUTION CONTROL CHARTS FOR THE LOGNORMAL DISTRIBUTION Petros Maravelakis, Joh Paaretos ad Stelios Psarakis Departmet of Statistics Athes Uiversity of Ecoomics ad Busiess 76 Patisio St., 4 34, Athes, GREECE. Itroductio

More information

Parameter, Statistic and Random Samples

Parameter, Statistic and Random Samples Parameter, Statistic ad Radom Samples A parameter is a umber that describes the populatio. It is a fixed umber, but i practice we do ot kow its value. A statistic is a fuctio of the sample data, i.e.,

More information

Chapter 8: Estimating with Confidence

Chapter 8: Estimating with Confidence Chapter 8: Estimatig with Cofidece Sectio 8.2 The Practice of Statistics, 4 th editio For AP* STARNES, YATES, MOORE Chapter 8 Estimatig with Cofidece 8.1 Cofidece Itervals: The Basics 8.2 8.3 Estimatig

More information

6.3 Testing Series With Positive Terms

6.3 Testing Series With Positive Terms 6.3. TESTING SERIES WITH POSITIVE TERMS 307 6.3 Testig Series With Positive Terms 6.3. Review of what is kow up to ow I theory, testig a series a i for covergece amouts to fidig the i= sequece of partial

More information

Statistical Inference Based on Extremum Estimators

Statistical Inference Based on Extremum Estimators T. Rotheberg Fall, 2007 Statistical Iferece Based o Extremum Estimators Itroductio Suppose 0, the true value of a p-dimesioal parameter, is kow to lie i some subset S R p : Ofte we choose to estimate 0

More information

PREDICTION INTERVALS FOR FUTURE SAMPLE MEAN FROM INVERSE GAUSSIAN DISTRIBUTION

PREDICTION INTERVALS FOR FUTURE SAMPLE MEAN FROM INVERSE GAUSSIAN DISTRIBUTION Qatar Uiv. Sci. J. (1991), 11: 19-26 PREDICTION INTERVALS FOR FUTURE SAMPLE MEAN FROM INVERSE GAUSSIAN DISTRIBUTION By MUHAMMAD S. ABU-SALIH ad RAFIQ K. AL-BAITAT Departmet of Statistics, Yarmouk Uiversity,

More information

The Sampling Distribution of the Maximum. Likelihood Estimators for the Parameters of. Beta-Binomial Distribution

The Sampling Distribution of the Maximum. Likelihood Estimators for the Parameters of. Beta-Binomial Distribution Iteratioal Mathematical Forum, Vol. 8, 2013, o. 26, 1263-1277 HIKARI Ltd, www.m-hikari.com http://d.doi.org/10.12988/imf.2013.3475 The Samplig Distributio of the Maimum Likelihood Estimators for the Parameters

More information

Big Picture. 5. Data, Estimates, and Models: quantifying the accuracy of estimates.

Big Picture. 5. Data, Estimates, and Models: quantifying the accuracy of estimates. 5. Data, Estimates, ad Models: quatifyig the accuracy of estimates. 5. Estimatig a Normal Mea 5.2 The Distributio of the Normal Sample Mea 5.3 Normal data, cofidece iterval for, kow 5.4 Normal data, cofidece

More information

Asymptotic distribution of products of sums of independent random variables

Asymptotic distribution of products of sums of independent random variables Proc. Idia Acad. Sci. Math. Sci. Vol. 3, No., May 03, pp. 83 9. c Idia Academy of Scieces Asymptotic distributio of products of sums of idepedet radom variables YANLING WANG, SUXIA YAO ad HONGXIA DU ollege

More information

TMA4245 Statistics. Corrected 30 May and 4 June Norwegian University of Science and Technology Department of Mathematical Sciences.

TMA4245 Statistics. Corrected 30 May and 4 June Norwegian University of Science and Technology Department of Mathematical Sciences. Norwegia Uiversity of Sciece ad Techology Departmet of Mathematical Scieces Corrected 3 May ad 4 Jue Solutios TMA445 Statistics Saturday 6 May 9: 3: Problem Sow desity a The probability is.9.5 6x x dx

More information

Bayesian and E- Bayesian Method of Estimation of Parameter of Rayleigh Distribution- A Bayesian Approach under Linex Loss Function

Bayesian and E- Bayesian Method of Estimation of Parameter of Rayleigh Distribution- A Bayesian Approach under Linex Loss Function Iteratioal Joural of Statistics ad Systems ISSN 973-2675 Volume 12, Number 4 (217), pp. 791-796 Research Idia Publicatios http://www.ripublicatio.com Bayesia ad E- Bayesia Method of Estimatio of Parameter

More information

Simple Linear Regression

Simple Linear Regression Simple Liear Regressio 1. Model ad Parameter Estimatio (a) Suppose our data cosist of a collectio of pairs (x i, y i ), where x i is a observed value of variable X ad y i is the correspodig observatio

More information

A NEW METHOD FOR CONSTRUCTING APPROXIMATE CONFIDENCE INTERVALS FOR M-ESTU1ATES. Dennis D. Boos

A NEW METHOD FOR CONSTRUCTING APPROXIMATE CONFIDENCE INTERVALS FOR M-ESTU1ATES. Dennis D. Boos .- A NEW METHOD FOR CONSTRUCTING APPROXIMATE CONFIDENCE INTERVALS FOR M-ESTU1ATES by Deis D. Boos Departmet of Statistics North Carolia State Uiversity Istitute of Statistics Mimeo Series #1198 September,

More information

Department of Civil Engineering-I.I.T. Delhi CEL 899: Environmental Risk Assessment HW5 Solution

Department of Civil Engineering-I.I.T. Delhi CEL 899: Environmental Risk Assessment HW5 Solution Departmet of Civil Egieerig-I.I.T. Delhi CEL 899: Evirometal Risk Assessmet HW5 Solutio Note: Assume missig data (if ay) ad metio the same. Q. Suppose X has a ormal distributio defied as N (mea=5, variace=

More information

Module 1 Fundamentals in statistics

Module 1 Fundamentals in statistics Normal Distributio Repeated observatios that differ because of experimetal error ofte vary about some cetral value i a roughly symmetrical distributio i which small deviatios occur much more frequetly

More information