Jackknife Empirical Likelihood Methods for Risk Measures and Related Quantities

Size: px
Start display at page:

Download "Jackknife Empirical Likelihood Methods for Risk Measures and Related Quantities"

Transcription

1 Jackkife Empirical Likelihood Methods for Risk Measures ad Related Quatities Liag Peg, Yogcheg Qi 2, Ruodu Wag ad Jigpig Yag 3 Abstract The quatificatio of risks is of importace i isurace. I this paper, we employ the jackkife empirical likelihood methods to costruct cofidece itervals for the risk measures ad some related quatities studied by Joes ad Zitikis (23). A simulatio study shows the advatages of the ew methods over the ormal approximatio method ad aive bootstrap method. Keywords: Cofidece iterval, jackkife empirical likelihood, risk measure. Itroductio I life isurace ad fiace, quatifyig risk is a very importat task for pricig a isurace product or maagig a fiacial portfolio. Geerally speakig, a risk measure is costructed to be a mappig from a set of risks to some real umbers. Some well-kow risk measures iclude coheret risk measures (Yaari (987), Artzer (999)), distortio risk measures, Wag s premium priciple ad proportioal hazards trasform risk measures (Wag, Youg ad Pajer (997); Wag (995), (996), (998); Wirch ad Hardy (999); Necir ad Meraghi (29)). For a risk variable X with distributio fuctio F, Joes ad Zitikis (23) defied a large class of risk measures associated with X as R(F ) = F (t)ψ(t)dt, () where F deotes the geeralized iverse fuctio of F, ad ψ is a oegative fuctio chose for showig the objective opiio about the risk loadig. Differet choices of ψ result i differet risk measures. For example, the proportioal hazards trasform risk measure has ψ(t) = r( t) r School of Mathematics, Georgia Istitute of Techology, Atlata, GA , USA. peg@math.gatech.edu, ruodu.wag@math.gatech.edu 2 Departmet of Mathematics ad Statistics, Uiversity of Miesota Duluth, 7 Uiversity Drive, Duluth, MN yqi@d.um.edu 3 LMAM, Departmet of Fiacial Mathematics, Pekig Uiversity, Beijig 87, Chia. yagjp@math.pku.edu.c

2 ad Wag s premium priciple has ψ(t) = g ( t), where g is a icreasig covex fuctio with derivatives over [, ]; see Joes ad Zitikis (23) for details. Other choices of the fuctio ψ ca be foud i Joes ad Zitikis (27). Joes ad Zitikis (23) also itroduced a related quatity to illustrate the right-tail, left-tail ad two-sided deviatios, which is defied as r(f ) = R(F ) E(X). (2) I this paper, we will focus o the statistical iferece of the risk measure ad its related quatity defied i () ad (2), respectively. How to ifer R(F ) ad r(f ) plays a importat role i the applicatios of risk measures. Recetly, Joes ad Zitikis (23) proposed oparametric estimatio by replacig F ad E(X) by the sample quatile fuctio ad sample mea respectively, ad derived the asymptotic ormality. Therefore, cofidece itervals for R(F ) ad r(f ) ca be costructed via estimatig the asymptotic variace. For comparig two risk measures, we refer to Joes ad Zitikis (25). Joes ad Zitikis (27) ivestigated the oparametric estimatio of the parameter associated with distortio-based risk measures. order to costruct cofidece itervals for R(F ) ad r(f ) without estimatig the asymptotic variace, we ivestigate the possibility of applyig empirical likelihood method i this paper so as to improve the iferece. Empirical likelihood method is a oparametric likelihood approach for statistical iferece, which has bee show to be powerful i iterval estimatio ad hypothesis testig. We refer to Owe (2) for a overview o the method. However, it is kow that empirical likelihood methods are ot effective i dealig with o-liear fuctioals. Recetly, a so-called jackkife empirical likelihood method was proposed by Jig, Yua ad Zhou (29) to deal with oliear fuctioals. The key idea is to formulate a jackkife sample based o estimatig the oliear fuctioal ad the apply the empirical likelihood method for a mea to the jackkife sample. Sice the risk measure R(F ) ad its related quatity r(f ) are o-liear fuctioals, we propose to employ the jackkife empirical likelihood method to obtai iterval estimatio for these two quatities. The paper is orgaized as follows. I Sectio 2, the methodologies ad mai results are preseted. A simulatio study is give i Sectio 3. All proofs are put i Sectio 4. Some coclusios are draw i Sectio 5. I 2

3 2 Methodologies ad mai results Put Ψ(t) = t The the risk measure defied i () ca be writte as R = R(F ) = ψ(s)ds. (Ψ() Ψ(F (t))dt. Let X,, X be idepedet real-valued radom variables with cotiuous distributio fuctio F (x) ad E(X ). Defie the empirical distributio fuctio as F (x) = (X j x). The Joes ad Zitikis (23) proposed to estimate R(F ) ad r(f ) by ˆR = respectively, ad showed that j= (Ψ() Ψ(F (t)))dt, ad ˆr = (Ψ() Ψ(F (t)))dt j= X, j { ˆR R} d N(, σ 2 ) ad {ˆr r(f )} d N(, σ 2 2) (3) uder some regularity coditios, where σ 2 = Q F (Ψ, Ψ), σ 2 2 = µ 2 ( QF (Ψ, Ψ) 2r(F )Q F (Ψ, ) + (r(f )) 2 Q F (, ) ) (4) ad Q F (a, b) = (F (x y) F (x)f (y))a(f (x))b(f (y))dxdy. Based o (3), cofidece itervals for R(F ) ad r(f ) ca be obtaied via estimatig σ 2 ad σ2 2. A alterative way to costruct cofidece itervals is to employ the empirical likelihood method. Sice the risk measure R is o-liear, a commo techique is to liearize the fuctioal by itroducig some lik variables before applyig the profile empirical likelihood method; see the study for variace, ROC curve (Claeskes et al. (23)) ad copulas (Che, Peg ad Zhao (29)). Ufortuately it remais ukow o how to liearize R by itroducig some lik variables. Here we propose to apply the jackkife empirical likelihood method i Jig, Yua ad Zhou (29). This procedure is easy to implemet ad is described as follows. 3

4 Defie F,i = j=,j i (X j x) ad ˆR,i = (Ψ() Ψ(F,i(t)))dt for i =,,. The the jackkife sample is defied as Y i = ˆR ( ) ˆR,i, i =,,. Now we apply the empirical likelihood method to the above jackkife sample. That is, we defie the jackkife empirical likelihood fuctio for θ = R(F ) as L (θ) = sup{ p i : p i, for i =,, ; p i = ; p i Y i = θ}. By Lagrage multiplier techique, we have p i = {+λ(y i θ)} ad 2 log L (θ) = 2 log{+ λ(y i θ)}, where λ = λ(θ) satisfies Y i θ =. (5) + λ(y i θ) The followig theorem shows that the Wilks theorem holds for the proposed jackkife empirical likelihood method. Theorem. Assume that ψ(x) cx α ( x) β, ψ (x) exists ad ψ (x) cx α 2 ( x) β 2 for all < x < ad some costats α > /2, β > /2 ad c >. Further assume E( X i γ ) < for some γ such that γ > /(α /2) ad γ > /(β /2). The we have 2 log L (R ) d χ 2 as, where R deotes the true value of R ad χ 2 freedom. deotes a chi-square distributio with oe degree of Remark. Some well-kow risk measures, such as proportioal hazards trasform risk measure, Wag s right-tail deviatio ad Wag s left-tail deviatio satisfy the assumptios of Theorem ; see Joes ad Zitikis (23). Based o the above theorem, a cofidece iterval for R with level b ca be obtaied as I R b = {R : 2 log L (R) χ 2,b }, where χ 2,b is the b-th quatile of χ2. 4

5 Next we cosider the related quatity r(f ) = R(F )/µ where µ = E(X ). Alteratively, we cosider the quatity R θµ with θ = r(f ). The oe ca estimate this quatity by ˆR θ As before, we defie the jackkife sample as X i = ˆR θ x df (x) = ˆR θ F (x) dx. ( ) ( ) ˆR θ x df (x) ( ) ˆR,i θ x df,i (x) = Y i θx i for i =,,, where Y i s are defied as above. So the jackkife empirical likelihood fuctio for θ = r(f ) is defied as L 2 (θ) = sup{ p i : p i, for i =,, ; p i = ; p i (Y i θx i ) = }. The followig theorem shows that the Wilks theorem holds for the proposed jackkife empirical likelihood method for r(f ). Theorem 2. Assume the coditios of Theorem hold. Further assume EX, ad E(X 2) <. The 2 log L 2 (r ) d χ 2 as, where r deotes the true value of r(f ). Based o the above theorem, a cofidece iterval for r with level b ca be obtaied as I r b = {r : 2 log L 2(r) χ 2,b }. 3 Simulatio study I this sectio we examie the fiite sample behavior of the proposed jackkife empirical likelihood methods i terms of coverage accuracy ad compare with the ormal approximatio method ad the aive bootstrap method. We focus o the proportioal hazards trasform risk measure with ψ(s) = a( s) a. Sice Pareto distributio, log-ormal distributio, Weibull distributio ad Gamma distributio are widely used i fittig the losses data i isurace (Klugma, Pajer ad Willmot (28)), our simulatio study is based o these four distributios. We draw, radom samples of sizes = 3 ad from the followig distributios: 5

6 . Pareto distributio F (x; θ) = x θ for x ; 2. Log-ormal distributio F 2 (x; θ, θ 2 ) = Φ((log x θ )/θ 2 ) for x >, where Φ(x) deotes the stadard ormal distributio; 3. Weibull distributio F 3 (x; θ, θ 2 ) = exp{ (x/θ 2 ) θ } for x > ; 4. Gamma distributio F 4 (x; θ, θ 2 ) = x θ θ 2 Γ(θ ) sθ exp{ θ 2 s}ds for x >. For calculatig the proposed jackkife empirical likelihood itervals for both R(F ) ad r(f ) (JELCI), we employ the proportioal hazards trasform risk measure with ψ(s) = a( s) a, a =.55 ad.85, ad use the R package emplik. For calculatig the cofidece itervals for R(F ) based o the ormal approximatio method (NACI), we use the variace estimatio i Joes ad Zitikis (23). For computig the aive bootstrap cofidece itervals for r(f ) (NBCI), we draw, bootstrap samples with replacemet from each radom sample X,, X. Empirical coverage probabilities are reported i Tables ad 2 for these three cofidece itervals with levels.9,.95 ad.99. From these two tables, we coclude that the proposed jackkife empirical likelihood methods give much more accurate coverage probabilities tha the other two methods. 4 Proofs Throughout we put U i = F (X i ) for i =,,, G (t) = (U i t) ad G,i = ( ) j=,j i (U j t) for i =,,. Sice F is cotiuous, U,, U are idepedet ad uiformly distributed over (, ). Without loss of geerality we also assume o ties i U,, U, ad let U, < < U, deote the order statistics of U,, U. We also use C to deote a geeric costat which may be differet i differet places. First we list some facts which will be employed i the proofs. We assume β α throughout sice proofs for the case of β > α are exactly the same. Therefore we have ψ(x) cx β ( x) β ad ψ (x) cx β 2 ( x) β 2 for all < x <. Sice E X γ < with γ + β < 2, we have which implies P ( X > x) = o(x γ ) as x, (6) (F (x)) β +δ ( F (x)) β +δ dx 2 + C 6 x (β +δ)γ dx < (7)

7 wheever δ ( γ + β, 2 ), ad It follows from the give coditios o ψ that max X j = max F (U j ) = o p ( /γ ). (8) j j Ψ( ) = O( β ) ad Ψ( ) Ψ( ) = O( β ). (9) Lemma. Uder the coditios of Theorem, we have (Y i R ) d N(, σ), 2 () where σ 2 is give i (4). Proof. Write Y i = ( ) = ( ) = ( ) = ( ) + ( ) + ( ) {Ψ(F,i (t)) Ψ(F (t))}dt + ˆR {Ψ(G,i (t)) Ψ(G (t))}df (t) + ˆR {Ψ(G,i (t)) Ψ(G (t))}(u, t < U, )df (t) + ˆR {Ψ(G,i (t)) Ψ(G (t))}(u, t < U,2 )df (t) {Ψ(G,i (t)) Ψ(G (t))}(u,2 t < U, )df (t) {Ψ(G,i (t)) Ψ(G (t))}(u, t < U, )df (t) + ˆR = ( ) {Ψ(G,i (t)) Ψ(G (t))}(u, t < U,2 )df (t) }{{} Z i, + ( ) ψ(g (t)){g,i (t) G (t)}(u,2 t < U, )df (t) }{{} Z i,2 + ψ (ξ,i (t)){g,i (t) G (t)} 2 (U,2 t < U, )df (t) 2 }{{} Z i,3 + ( ) {Ψ(G,i (t)) Ψ(G (t))}(u, t < U, )df (t) + ˆR } {{ } Z i,4 = Z i, + Z i,2 + Z i,3 + Z i,4 + ˆR, 7

8 where for some θ i (t) [, ]. ξ,i (t) = G (t) + θ i (t){g,i (t) G (t)} = G (t) + θ i(t) {G (t) (U i t)} Whe U, t < U,2, we have G (t) = ad G,i(t) = if U i = U, else. () Hece, it follows from (8) ad (9) that Z i, = ( ) {Ψ() Ψ( )}(U, t < U,2 )df (t) + ( ) 2 {Ψ( ) Ψ( )}(U, t < U,2 )df (t) = ( )Ψ( ){F (U,2 ) F (U, )} + ( ) 2 {Ψ( ) Ψ( )}{F (U,2 ) F (U, )} = O(( ) β )o p ( /γ ) + O(( ) 2 β o p ( /γ ) (2) = o p ( /2 β+/γ ) = o p ( ) sice 2 β + γ <. Similarly, we ca show that Sice {G,i(t) G (t)} =, we have Whe t U,2, we have Z i,4 = o p ( ). (3) Z i,2 =. (4) ( ) (U i t) G (t) /( ) 2/ = 2( ), i.e., ξ,i (t) G (t){ 2( ) } uiformly i t U,2. I the same maer, we ca show that the equatio ξ,i (t) ( G (t)){ 8 2( ) }

9 holds uiformly i t < U,. Hece, for large eough, (ξ,i (t), ξ,i (t)) 3 (G (t), G (t)) uiformly for U,2 t < U, ad i. (5) Note that G (t) G (t) sup = O p () ad sup = O p () (6) U,2 t U, t U,2 t U, t (see Page 44 of Shorack ad Weller (986)). It follows from (5) ad (6) that Z i,3 = O p ( ) t β 2 ( t) β 2 {G,i (t) G (t)} 2 (U,2 t < U, )df (t), which coupled with (7) yields Z i,3 = O p ( t β 2 ( t) β 2 ) {G,i (t) G (t)} 2 (U,2 t < U, )df (t) ) = O p ( t β 2 ( t) β 2 ( ) 2 G (t){ G (t)}(u,2 t < U, )df (t) ( ) = O p t β ( t) β (U,2 t < U, )df (t) ( ) = O p t β ( t) β df (t) ( ) = O p δ t β +δ ( t) β +δ df (t) ( ) = O p δ (F (x)) β +δ ( F (x)) β +δ dx = O p ( δ ) (7) for ay δ ( γ + β, 2 ). By Joes ad Zitikis (23), we have { ˆR R} d N(, σ 2 ). (8) Hece, the lemma follows from (2), (4), (7), (3) ad (8). Lemma 2. Uder the coditios of Theorem, we have (Y i R) 2 p σ 2 as. 9

10 Proof. We use the same otatios Z i,j as i the proof of Lemma. The, it follows from () ad (9) that Zi, 2 = ( )2 + = O( = o p (). ( )3 {Ψ() Ψ( )}2 (U, t, t 2 < U,2 )df (t )df (t 2 ) {Ψ( ) Ψ( )}2 (U, t, t 2 < U,2 )df (t )df (t 2 ) ( )2 2β )o p ( 2/γ ( )3 ) + O( 2 2β )o p ( 2/γ ) (9) Similarly, It is easy to check that Zi,2 2 ( )2 = = =2 ψ(g (t ))ψ(g (t 2 )) Zi,4 2 = o p (). (2) {G,i (t ) G (t )}{G,i (t 2 ) G (t 2 )} (U,2 t, t 2 < U, )df (t 2 )df (t ) t ψ(g (t ))ψ(g (t 2 )){G (t t 2 ) G (t )G (t 2 )}(U,2 t, t 2 < U, )df (t 2 )df (t ) ψ(g (t ))ψ(g (t 2 ))G (t 2 ){ G (t )}(U, t, t 2 < U, )df (t 2 )df (t ). } {{ } I By (6), we have ( ) sup ψ(g (t ))ψ(g (t 2 ))G (t 2 ){ G (t )} = O p t β ( t ) β t β 2 ( t 2 ) β t 2 ( t ). U,2 t,t 2 <U, Similar to the proof of (7), we ca show that = t x <. t β ( t ) β t β 2 ( t 2 ) β t 2 ( t )df (t 2 )df (t ) F (x) β ( F (x)) β F (y) β ( F (y)) β F (y)( F (x))dydx By the Gliveko-Catelli theorem, sup <t< G (t) t almost surely. It the follows from the domiated covergece theorem that I p t ψ(t )ψ(t 2 )t 2 ( t )df (t 2 )df (t ).

11 Hece Note that = Z 2 i,2 p {G,i (t ) G (t )} 2 {G,i (t 2 ) G (t 2 )} 2 { G (t ) (U i t ) } 2 { G (t 2 ) (U i t 2 ) ψ(t )ψ(t 2 ){t t 2 t t 2 }df (t 2 )df (t ) = σ 2. (2) = ( ) 4 { 3G2 (t )G 2 (t 2 ) + G 2 (t )G (t 2 ) + G (t )G 2 (t 2 ) + 4G (t )G (t 2 )G (t t 2 ) 2G (t )G (t t 2 ) 2G (t 2 )G (t t 2 ) + G (t t 2 )} = ( ) 4 {3G (t )G (t 2 ) (G (t t 2 ) G (t )G (t 2 )) G }{{} (t ) (G (t t 2 ) G (t )G (t 2 )) }{{} II G (t 2 ) (G (t t 2 ) G (t )G (t 2 )) + ( G }{{} (t ))( G (t 2 ))G (t t 2 ) }{{} II 3 = ( ) 4 {II II 2 II 3 + II 4 }. It follows from (6) that } 2 II 2 II 4 } G (t t 2 ) G (t )G (t 2 ) G (t t 2 )( G (t t 2 ) sup = sup = O p (), U,2 t,t 2 U, t t 2 t t 2 U,2 t,t 2 U, t t 2 ( t t 2 ) which, coupled with (5) ad (6), yields that Z 2 i,3 (( ) 2 =O p t β 2 ( t ) β 2 t β 2 2 ( t 2 ) β 2 4 {G,i (t ) G (t )} 2 {G,i (t 2 ) G (t 2 )} 2 (U,2 t, t 2 < U, )df (t 2 )df (t ) ) =O p ( 2 t β 2 ( t ) β 2 t β 2 2 ( t 2 ) β 2 {II II 2 II 3 + II 4 }(U,2 t, t 2 < U, )df (t 2 )df (t ) ).

12 From the above equatio we ca get that Zi,3 2 U, U, ( =O p 2 t β 2 ( t ) β 2 t β 2 2 ( t 2 ) β 2 t t 2 (t t 2 t t 2 )df (t 2 )df ) (t ) U,2 U,2 }{{} III ( U, U, + O p 2 t β 2 ( t ) β 2 t β 2 2 ( t 2 ) β 2 t (t t 2 t t 2 )df (t 2 )df ) (t ) U,2 U,2 }{{} ( U, U, + O p 2 t β 2 ( t ) β 2 t β 2 2 ( t 2 ) β 2 t 2 (t t 2 t t 2 )df (t 2 )df ) (t ) U,2 U,2 }{{} ( U, + O p 2 U,2 U, U,2 III 2 III 3 t β 2 ( t ) β 2 t β 2 2 ( t 2 ) β 2 ( t )( t 2 )(t t 2 )df (t 2 )df ) (t ) } {{ } III 4 =O p (III ) + O p (III 2 ) + O p (III 3 ) + O p (III 4 ). It is easy to check from (7) that for every δ ( γ + β, 2 ) III 2 + III 3 =2 2 U, U,2 4 2 U, = 2 U,2 =O( 2 U,2 δ ) t U,2 t β 2 ( t ) β 2 t β 2 2 ( t 2 ) β 2 (t + t 2 )t 2 ( t )df (t 2 )df (t ) t β ( t ) β t U,2 t β 2 ( t 2 ) β 2 df (t 2 )df (t ) t β ( t ) β O(U,2 δ + ( t ) δ )df (t ) t β ( t ) β df (t ) + O( 2 ) t β ( t ) β 2 δ df (t ) =O( 2 U,2 δ δ )(U,2 + ( U, ) δ ) + O( 2 )(U,2 δ + ( U, ) 2δ ) =O p ( 2+2δ + +2δ ) =o p (). Similarly, we ca show that III = o p () ad III 4 = o p (). Hece, Zi,3 2 = o p (). (22) Sice ˆR p R, we have ( ˆR R) 2 = o p (). (23) 2

13 It follows from (9), (2), (22) ad (23) that {Z i, + Z i,3 + Z i,4 + ˆR R} 2 = O( Note that Z i,2 {Z i, + Z i,3 + Z i,4 + ˆR R} Therefore, the lemma follows from (9) (25). {Zi, 2 + Zi,3 2 + Zi,4 2 + ( ˆR R) 2 }) = o p (). (24) = o p (). Z 2 i,2 {Z i, + Z i,3 + Z i,4 + ˆR R} 2 Proof of Theorem. First we observe by usig (7) that for ay δ ( γ + β, 2 ) (25) max Z i,2 ψ(g (t))(u,2 t < U, )df (t) i ( ) U, = O p t β ( t) β df (t) U,2 (26) = O p (U δ,2 + ( U, ) δ ) = o p ( /2 ). Similarly we ca show that max Z i,j = o p ( /2 ) for j =, 3, 4. i Hece, max i Y i = o p ( /2 ). By the stadard argumets i the empirical likelihood method (see Chapter of Owe (2)), it follows from Lemmas ad 2 that 2 log L (R) = { (Y i R)} 2 (Y i R) 2 + o p () d χ 2 (). I order to prove Theorem 2, we eed the followig lemmas. Lemma 3. Uder the coditios of Theorem 2, we have ( ) ˆR R(F ) d X i N(, σ 2 ) as, µ where σ 2 = + 2 R(F ) µ ψ(t )ψ(t 2 )(t t 2 t t 2 )df (t )df (t 2 ) + R2 (F ) µ 2 E(X µ) 2 ψ(t )(t t 2 t t 2 )df (t )df (t 2 ). 3

14 Proof. It is kow that there exists a Browia bridge W such that (G (t) t) W (t) sup t t δ ( t) δ = o p () (27) for ay δ (, /2) (see Chapter 4 of Csorgo ad Horvath (993)). The we have { ˆR R(F ) µ X i } = d {Ψ(t) Ψ(G (t))}df (t) + R(F ) {t G (t)}df (t) µ ψ(t)w (t)df (t) R(F ) µ Lemma 4. Uder the coditios of Theorem 2, we have W (t)df (t). ( Y i R(F ) ) µ X d i N(, σ 2 ) as. Proof. It ca be show i a way similar to the proof of Lemma. Lemma 5. Uder the coditios of Theorem 2, we have ( Y i R(F ) ) 2 p µ X i σ 2 as. Proof. It ca be proved i a similar way to the proof of Lemma 2. Proof of Theorem 2. This ca be doe i a way similar to the proof of Theorem. 5 Coclusios This paper proposes jackkife empirical likelihood methods to costruct cofidece itervals for the risk measures ad some related quatities studied by Joes ad Zitikis (23). Ulike the ormal approximatio method, the ew methods do ot eed to estimate the asymptotic variace explicitly ad are easy to implemet by employig the R package emplik. A simulatio study shows that the jackkife empirical likelihood cofidece itervals are more accurate tha the ormal approximatio based cofidece itervals. 4

15 Ackowledgmets Peg s research was supported by NSA grat H , Qi s research was supported by NSA grat H , ad Yag s research was partly supported by the Natioal Basic Research Program (973 Program) of Chia (27CB8495) ad Natioal Natural Sciece Foudatio of Chia (Grats No. 878). Refereces [] P. Artzer (999). Applicatio of coheret risk measures to capital requiremets i isurace. North America Actuarial Joural 3, 29. [2] J. Che, L. Peg ad Y. Zhao (29). Empirical likelihood based cofidece itervals for copulas. J. Multivariate Aal., [3] G. Claeskes, B. Jig, L. Peg ad W. Zhou (23). Empirical likelihood cofidece regios for compariso distributios ad ROC curves. The Caadia Joural of Statistics. 3(2), [4] M. Csorgo ad L. Horvath (993). Weighted Approximatios i Probability ad Statistics. Wiley, Chichester. [5] B. Jig, J. Yua ad W. Zhou (29). Jackkife empirical likelihood. Joural of America Statistical Associatio 4, [6] B.L. Joes ad R. Zitikis (23). Empirical estimatio of risk measures ad related quatities. North America Actuarial Joural 7, [7] B.L. Joes ad R. Zitikis (25). Testig for the order of risk measures: a applicatio of L- statistics i a actuarial sciece. Metro LXIII, [8] B.L. Joes ad R. Zitikis (27). Risk measures, distortio parameters, ad their empirical estimatio. Isurace: Mathematics ad Ecoomics 4, [9] S.A. Klugma, H.H. Pajer ad G.E. Willmot (28). Loss Distributios: From Data to Decisios. Wiley. 5

16 [] A. Necir ad D. Meraghi (29). Empirical estimatio of the proportioal hazard premium for heavy-tailed claim amouts. Isurace: Mathematics ad Ecoomics 45, [] A. Owe (2). Empirical Likelihood. Chapma & Hall/CRC. [2] G.R. Shorack ad J.A. Weller (986). Empirical Processes with Applicatios to Statistics. New York: Wiley. [3] S. Wag (995). Isurace pricig ad icreased limits ratemakig by proportioal hazards trasforms. Isurace: Mathematics ad Ecoomics 7, [4] S. Wag (996). Premium calculatio by trasformig the layer premium desity. ASTIN Bulleti 26, [5] S. Wag (998). A actuarial idex of the right-tail risk. North America Actuarial Joural 2, 88. [6] S. Wag, V.R. Youg ad H.H. Pajer (997). Axiomatic characterizatio of isurace prices. Isurace: Mathematics ad Ecoomics 2, [7] J.L. Wirch ad M.R. Hardy (999). A sythesis of risk measures for capital adequacy. Isurace: Mathematics ad Ecoomics 25, [8] M.E. Yaari (987). The dual theory of choice uder risk. Ecoometrica 55,

17 Table : Coverage probabilities for R(F ) are reported for the itervals based o the proposed jackkife empirical likelihood method (JELCI) ad the ormal approximatio method (NACI). (, a, F ) JELCI NACI JELCI NACI JELCI NACI level.9 level.9 level.95 level.95 level.99 level.99 (3,.55, F (; 4)) (3,.85, F (; 4)) (,.55, F (; 4)) (,.85, F (; 4)) (3,.55, F 2 (;, )) (3,.85, F 2 (;, )) (,.55, F 2 (;, )) (,.85, F 2 (;, )) (3,.55, F 3 (; 4, )) (3,.85, F 3 (; 4, )) (,.55, F 3 (; 4, )) (,.85, F 3 (; 4, )) (3,.55, F 4 (; 4, )) (3,.85, F 4 (; 4, )) (,.55, F 4 (; 4, )) (,.85, F 4 (; 4, ))

18 Table 2: Coverage probabilities for r(f ) are reported for the itervals based o the proposed jackkife empirical likelihood method (JELCI) ad the aive bootstrap method (NBCI). (, a, F ) JELCI NBCI JELCI NBCI JELCI NBCI level.9 level.9 level.95 level.95 level.99 level.99 (3,.55, F (; 4)) (3,.85, F (; 4)) (,.55, F (; 4)) (,.85, F (; 4)) (3,.55, F 2 (;, )) (3,.85, F 2 (;, )) (,.55, F 2 (;, )) (,.85, F 2 (;, )) (3,.55, F 3 (; 4, )) (3,.85, F 3 (; 4, )) (,.55, F 3 (; 4, )) (,.85, F 3 (; 4, )) (3,.55, F 4 (; 4, )) (3,.85, F 4 (; 4, )) (,.55, F 4 (; 4, )) (,.85, F 4 (; 4, ))

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

Comparison of Minimum Initial Capital with Investment and Non-investment Discrete Time Surplus Processes

Comparison of Minimum Initial Capital with Investment and Non-investment Discrete Time Surplus Processes The 22 d Aual Meetig i Mathematics (AMM 207) Departmet of Mathematics, Faculty of Sciece Chiag Mai Uiversity, Chiag Mai, Thailad Compariso of Miimum Iitial Capital with Ivestmet ad -ivestmet Discrete Time

More information

POWER COMPARISON OF EMPIRICAL LIKELIHOOD RATIO TESTS: SMALL SAMPLE PROPERTIES THROUGH MONTE CARLO STUDIES*

POWER COMPARISON OF EMPIRICAL LIKELIHOOD RATIO TESTS: SMALL SAMPLE PROPERTIES THROUGH MONTE CARLO STUDIES* Kobe Uiversity Ecoomic Review 50(2004) 3 POWER COMPARISON OF EMPIRICAL LIKELIHOOD RATIO TESTS: SMALL SAMPLE PROPERTIES THROUGH MONTE CARLO STUDIES* By HISASHI TANIZAKI There are various kids of oparametric

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

UNIFORM INTERVAL ESTIMATION FOR AN AR(1) PROCESS WITH AR ERRORS

UNIFORM INTERVAL ESTIMATION FOR AN AR(1) PROCESS WITH AR ERRORS Statistica Siica 6 (06), 9-36 doi:http://dx.doi.org/0.5705/ss.04.5 UNIFORM INTERVAL ESTIMATION FOR AN AR() PROCESS WITH AR ERRORS Joatha Hill, Deyua Li ad Liag Peg Uiversity of North Carolia, Fuda Uiversity

More information

Self-normalized deviation inequalities with application to t-statistic

Self-normalized deviation inequalities with application to t-statistic Self-ormalized deviatio iequalities with applicatio to t-statistic Xiequa Fa Ceter for Applied Mathematics, Tiaji Uiversity, 30007 Tiaji, Chia Abstract Let ξ i i 1 be a sequece of idepedet ad symmetric

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

THE DATA-BASED CHOICE OF BANDWIDTH FOR KERNEL QUANTILE ESTIMATOR OF VAR

THE DATA-BASED CHOICE OF BANDWIDTH FOR KERNEL QUANTILE ESTIMATOR OF VAR Iteratioal Joural of Iovative Maagemet, Iformatio & Productio ISME Iteratioal c2013 ISSN 2185-5439 Volume 4, Number 1, Jue 2013 PP. 17-24 THE DATA-BASED CHOICE OF BANDWIDTH FOR KERNEL QUANTILE ESTIMATOR

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

32 estimating the cumulative distribution function

32 estimating the cumulative distribution function 32 estimatig the cumulative distributio fuctio 4.6 types of cofidece itervals/bads Let F be a class of distributio fuctios F ad let θ be some quatity of iterest, such as the mea of F or the whole fuctio

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

Simulation. Two Rule For Inverting A Distribution Function

Simulation. Two Rule For Inverting A Distribution Function Simulatio Two Rule For Ivertig A Distributio Fuctio Rule 1. If F(x) = u is costat o a iterval [x 1, x 2 ), the the uiform value u is mapped oto x 2 through the iversio process. Rule 2. If there is a jump

More information

On an Application of Bayesian Estimation

On an Application of Bayesian Estimation O a Applicatio of ayesia Estimatio KIYOHARU TANAKA School of Sciece ad Egieerig, Kiki Uiversity, Kowakae, Higashi-Osaka, JAPAN Email: ktaaka@ifokidaiacjp EVGENIY GRECHNIKOV Departmet of Mathematics, auma

More information

Common Large/Small Sample Tests 1/55

Common Large/Small Sample Tests 1/55 Commo Large/Small Sample Tests 1/55 Test of Hypothesis for the Mea (σ Kow) Covert sample result ( x) to a z value Hypothesis Tests for µ Cosider the test H :μ = μ H 1 :μ > μ σ Kow (Assume the populatio

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

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

ON BARTLETT CORRECTABILITY OF EMPIRICAL LIKELIHOOD IN GENERALIZED POWER DIVERGENCE FAMILY. Lorenzo Camponovo and Taisuke Otsu.

ON BARTLETT CORRECTABILITY OF EMPIRICAL LIKELIHOOD IN GENERALIZED POWER DIVERGENCE FAMILY. Lorenzo Camponovo and Taisuke Otsu. ON BARTLETT CORRECTABILITY OF EMPIRICAL LIKELIHOOD IN GENERALIZED POWER DIVERGENCE FAMILY By Lorezo Campoovo ad Taisuke Otsu October 011 COWLES FOUNDATION DISCUSSION PAPER NO. 185 COWLES FOUNDATION FOR

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

Precise Rates in Complete Moment Convergence for Negatively Associated Sequences

Precise Rates in Complete Moment Convergence for Negatively Associated Sequences Commuicatios of the Korea Statistical Society 29, Vol. 16, No. 5, 841 849 Precise Rates i Complete Momet Covergece for Negatively Associated Sequeces Dae-Hee Ryu 1,a a Departmet of Computer Sciece, ChugWoo

More information

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

A note on self-normalized Dickey-Fuller test for unit root in autoregressive time series with GARCH errors

A note on self-normalized Dickey-Fuller test for unit root in autoregressive time series with GARCH errors Appl. Math. J. Chiese Uiv. 008, 3(): 97-0 A ote o self-ormalized Dickey-Fuller test for uit root i autoregressive time series with GARCH errors YANG Xiao-rog ZHANG Li-xi Abstract. I this article, the uit

More information

STATISTICAL INFERENCE

STATISTICAL INFERENCE STATISTICAL INFERENCE POPULATION AND SAMPLE Populatio = all elemets of iterest Characterized by a distributio F with some parameter θ Sample = the data X 1,..., X, selected subset of the populatio = sample

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

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

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

Convergence of random variables. (telegram style notes) P.J.C. Spreij

Convergence of random variables. (telegram style notes) P.J.C. Spreij Covergece of radom variables (telegram style otes).j.c. Spreij this versio: September 6, 2005 Itroductio As we kow, radom variables are by defiitio measurable fuctios o some uderlyig measurable space

More information

Problem Set 4 Due Oct, 12

Problem Set 4 Due Oct, 12 EE226: Radom Processes i Systems Lecturer: Jea C. Walrad Problem Set 4 Due Oct, 12 Fall 06 GSI: Assae Gueye This problem set essetially reviews detectio theory ad hypothesis testig ad some basic otios

More 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

A RANK STATISTIC FOR NON-PARAMETRIC K-SAMPLE AND CHANGE POINT PROBLEMS

A RANK STATISTIC FOR NON-PARAMETRIC K-SAMPLE AND CHANGE POINT PROBLEMS J. Japa Statist. Soc. Vol. 41 No. 1 2011 67 73 A RANK STATISTIC FOR NON-PARAMETRIC K-SAMPLE AND CHANGE POINT PROBLEMS Yoichi Nishiyama* We cosider k-sample ad chage poit problems for idepedet data i a

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

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

Kolmogorov-Smirnov type Tests for Local Gaussianity in High-Frequency Data

Kolmogorov-Smirnov type Tests for Local Gaussianity in High-Frequency Data Proceedigs 59th ISI World Statistics Cogress, 5-30 August 013, Hog Kog (Sessio STS046) p.09 Kolmogorov-Smirov type Tests for Local Gaussiaity i High-Frequecy Data George Tauche, Duke Uiversity Viktor Todorov,

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

Approximating the ruin probability of finite-time surplus process with Adaptive Moving Total Exponential Least Square

Approximating the ruin probability of finite-time surplus process with Adaptive Moving Total Exponential Least Square WSEAS TRANSACTONS o BUSNESS ad ECONOMCS S. Khotama, S. Boothiem, W. Klogdee Approimatig the rui probability of fiite-time surplus process with Adaptive Movig Total Epoetial Least Square S. KHOTAMA, S.

More information

Stat 319 Theory of Statistics (2) Exercises

Stat 319 Theory of Statistics (2) Exercises Kig Saud Uiversity College of Sciece Statistics ad Operatios Research Departmet Stat 39 Theory of Statistics () Exercises Refereces:. Itroductio to Mathematical Statistics, Sixth Editio, by R. Hogg, J.

More information

MOST PEOPLE WOULD RATHER LIVE WITH A PROBLEM THEY CAN'T SOLVE, THAN ACCEPT A SOLUTION THEY CAN'T UNDERSTAND.

MOST PEOPLE WOULD RATHER LIVE WITH A PROBLEM THEY CAN'T SOLVE, THAN ACCEPT A SOLUTION THEY CAN'T UNDERSTAND. XI-1 (1074) MOST PEOPLE WOULD RATHER LIVE WITH A PROBLEM THEY CAN'T SOLVE, THAN ACCEPT A SOLUTION THEY CAN'T UNDERSTAND. R. E. D. WOOLSEY AND H. S. SWANSON XI-2 (1075) STATISTICAL DECISION MAKING Advaced

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

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

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

Overview. p 2. Chapter 9. Pooled Estimate of. q = 1 p. Notation for Two Proportions. Inferences about Two Proportions. Assumptions

Overview. p 2. Chapter 9. Pooled Estimate of. q = 1 p. Notation for Two Proportions. Inferences about Two Proportions. Assumptions Chapter 9 Slide Ifereces from Two Samples 9- Overview 9- Ifereces about Two Proportios 9- Ifereces about Two Meas: Idepedet Samples 9-4 Ifereces about Matched Pairs 9-5 Comparig Variatio i Two Samples

More information

Introducing a Novel Bivariate Generalized Skew-Symmetric Normal Distribution

Introducing a Novel Bivariate Generalized Skew-Symmetric Normal Distribution Joural of mathematics ad computer Sciece 7 (03) 66-7 Article history: Received April 03 Accepted May 03 Available olie Jue 03 Itroducig a Novel Bivariate Geeralized Skew-Symmetric Normal Distributio Behrouz

More information

Berry-Esseen bounds for self-normalized martingales

Berry-Esseen bounds for self-normalized martingales Berry-Essee bouds for self-ormalized martigales Xiequa Fa a, Qi-Ma Shao b a Ceter for Applied Mathematics, Tiaji Uiversity, Tiaji 30007, Chia b Departmet of Statistics, The Chiese Uiversity of Hog Kog,

More information

EFFECTIVE WLLN, SLLN, AND CLT IN STATISTICAL MODELS

EFFECTIVE WLLN, SLLN, AND CLT IN STATISTICAL MODELS EFFECTIVE WLLN, SLLN, AND CLT IN STATISTICAL MODELS Ryszard Zieliński Ist Math Polish Acad Sc POBox 21, 00-956 Warszawa 10, Polad e-mail: rziel@impagovpl ABSTRACT Weak laws of large umbers (W LLN), strog

More information

ESTIMATING OF THE PROPORTIONAL HAZARD PREMIUM FOR HEAVY-TAILED CLAIM AMOUNTS WITH THE POT METHOD

ESTIMATING OF THE PROPORTIONAL HAZARD PREMIUM FOR HEAVY-TAILED CLAIM AMOUNTS WITH THE POT METHOD REVSTAT Statistical Joural Volume 1, Number 3, November 212, 335 349 ESTIMATING OF THE PROPORTIONAL HAZARD PREMIUM FOR HEAVY-TAILED CLAIM AMOUNTS WITH THE POT METHOD Author: Rassoul Abdelaziz GEE Laboratory,

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

Chapter 3. Strong convergence. 3.1 Definition of almost sure convergence

Chapter 3. Strong convergence. 3.1 Definition of almost sure convergence Chapter 3 Strog covergece As poited out i the Chapter 2, there are multiple ways to defie the otio of covergece of a sequece of radom variables. That chapter defied covergece i probability, covergece i

More information

Bull. Korean Math. Soc. 36 (1999), No. 3, pp. 451{457 THE STRONG CONSISTENCY OF NONLINEAR REGRESSION QUANTILES ESTIMATORS Seung Hoe Choi and Hae Kyung

Bull. Korean Math. Soc. 36 (1999), No. 3, pp. 451{457 THE STRONG CONSISTENCY OF NONLINEAR REGRESSION QUANTILES ESTIMATORS Seung Hoe Choi and Hae Kyung Bull. Korea Math. Soc. 36 (999), No. 3, pp. 45{457 THE STRONG CONSISTENCY OF NONLINEAR REGRESSION QUANTILES ESTIMATORS Abstract. This paper provides suciet coditios which esure the strog cosistecy of regressio

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

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

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

17. Joint distributions of extreme order statistics Lehmann 5.1; Ferguson 15

17. Joint distributions of extreme order statistics Lehmann 5.1; Ferguson 15 17. Joit distributios of extreme order statistics Lehma 5.1; Ferguso 15 I Example 10., we derived the asymptotic distributio of the maximum from a radom sample from a uiform distributio. We did this usig

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

Inferential Statistics. Inference Process. Inferential Statistics and Probability a Holistic Approach. Inference Process.

Inferential Statistics. Inference Process. Inferential Statistics and Probability a Holistic Approach. Inference Process. Iferetial Statistics ad Probability a Holistic Approach Iferece Process Chapter 8 Poit Estimatio ad Cofidece Itervals This Course Material by Maurice Geraghty is licesed uder a Creative Commos Attributio-ShareAlike

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

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall 2013 Lecture 2 9/9/2013. Large Deviations for i.i.d. Random Variables

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall 2013 Lecture 2 9/9/2013. Large Deviations for i.i.d. Random Variables MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall 2013 Lecture 2 9/9/2013 Large Deviatios for i.i.d. Radom Variables Cotet. Cheroff boud usig expoetial momet geeratig fuctios. Properties of a momet

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

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

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

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

Stat 200 -Testing Summary Page 1

Stat 200 -Testing Summary Page 1 Stat 00 -Testig Summary Page 1 Mathematicias are like Frechme; whatever you say to them, they traslate it ito their ow laguage ad forthwith it is somethig etirely differet Goethe 1 Large Sample Cofidece

More information

Math 152. Rumbos Fall Solutions to Review Problems for Exam #2. Number of Heads Frequency

Math 152. Rumbos Fall Solutions to Review Problems for Exam #2. Number of Heads Frequency Math 152. Rumbos Fall 2009 1 Solutios to Review Problems for Exam #2 1. I the book Experimetatio ad Measuremet, by W. J. Youde ad published by the by the Natioal Sciece Teachers Associatio i 1962, the

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

Lecture Notes 15 Hypothesis Testing (Chapter 10)

Lecture Notes 15 Hypothesis Testing (Chapter 10) 1 Itroductio Lecture Notes 15 Hypothesis Testig Chapter 10) Let X 1,..., X p θ x). Suppose we we wat to kow if θ = θ 0 or ot, where θ 0 is a specific value of θ. For example, if we are flippig a coi, we

More information

A Note on the Kolmogorov-Feller Weak Law of Large Numbers

A Note on the Kolmogorov-Feller Weak Law of Large Numbers Joural of Mathematical Research with Applicatios Mar., 015, Vol. 35, No., pp. 3 8 DOI:10.3770/j.iss:095-651.015.0.013 Http://jmre.dlut.edu.c A Note o the Kolmogorov-Feller Weak Law of Large Numbers Yachu

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

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

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

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

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

FACULTY OF MATHEMATICAL STUDIES MATHEMATICS FOR PART I ENGINEERING. Lectures

FACULTY OF MATHEMATICAL STUDIES MATHEMATICS FOR PART I ENGINEERING. Lectures FACULTY OF MATHEMATICAL STUDIES MATHEMATICS FOR PART I ENGINEERING Lectures MODULE 5 STATISTICS II. Mea ad stadard error of sample data. Biomial distributio. Normal distributio 4. Samplig 5. Cofidece itervals

More 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

Research Article On the Strong Laws for Weighted Sums of ρ -Mixing Random Variables

Research Article On the Strong Laws for Weighted Sums of ρ -Mixing Random Variables Hidawi Publishig Corporatio Joural of Iequalities ad Applicatios Volume 2011, Article ID 157816, 8 pages doi:10.1155/2011/157816 Research Article O the Strog Laws for Weighted Sums of ρ -Mixig Radom Variables

More information

Chapter 6 Sampling Distributions

Chapter 6 Sampling Distributions Chapter 6 Samplig Distributios 1 I most experimets, we have more tha oe measuremet for ay give variable, each measuremet beig associated with oe radomly selected a member of a populatio. Hece we eed to

More information

Generalized Law of the Iterated Logarithm and Its Convergence Rate

Generalized Law of the Iterated Logarithm and Its Convergence Rate Stochastic Aalysis ad Applicatios, 25: 89 03, 2007 Copyright Taylor & Fracis Group, LLC ISSN 0736-2994 prit/532-9356 olie DOI: 0.080/073629906005997 Geeralized Law of the Iterated Logarithm ad Its Covergece

More information

Kernel density estimator

Kernel density estimator Jauary, 07 NONPARAMETRIC ERNEL DENSITY ESTIMATION I this lecture, we discuss kerel estimatio of probability desity fuctios PDF Noparametric desity estimatio is oe of the cetral problems i statistics I

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

1.010 Uncertainty in Engineering Fall 2008

1.010 Uncertainty in Engineering Fall 2008 MIT OpeCourseWare http://ocw.mit.edu.00 Ucertaity i Egieerig Fall 2008 For iformatio about citig these materials or our Terms of Use, visit: http://ocw.mit.edu.terms. .00 - Brief Notes # 9 Poit ad Iterval

More information

Journal of Multivariate Analysis. Superefficient estimation of the marginals by exploiting knowledge on the copula

Journal of Multivariate Analysis. Superefficient estimation of the marginals by exploiting knowledge on the copula Joural of Multivariate Aalysis 102 (2011) 1315 1319 Cotets lists available at ScieceDirect Joural of Multivariate Aalysis joural homepage: www.elsevier.com/locate/jmva Superefficiet estimatio of the margials

More information

DERIVING GENERALIZED MEANS AS LEAST SQUARES AND MAXIMUM LIKELIHOOD ESTIMATES

DERIVING GENERALIZED MEANS AS LEAST SQUARES AND MAXIMUM LIKELIHOOD ESTIMATES DERIVING GENERALIZED MEANS AS LEAST SQUARES AND MAXIMUM LIKELIHOOD ESTIMATES Roger L. Berger Statistics Departmet Box 8203 North Carolia State Uiversity Raleigh, NC 27695-8203 ad George Casella 1 Biometrics

More information

7-1. Chapter 4. Part I. Sampling Distributions and Confidence Intervals

7-1. Chapter 4. Part I. Sampling Distributions and Confidence Intervals 7-1 Chapter 4 Part I. Samplig Distributios ad Cofidece Itervals 1 7- Sectio 1. Samplig Distributio 7-3 Usig Statistics Statistical Iferece: Predict ad forecast values of populatio parameters... Test hypotheses

More 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

An almost sure invariance principle for trimmed sums of random vectors

An almost sure invariance principle for trimmed sums of random vectors Proc. Idia Acad. Sci. Math. Sci. Vol. 20, No. 5, November 200, pp. 6 68. Idia Academy of Scieces A almost sure ivariace priciple for trimmed sums of radom vectors KE-ANG FU School of Statistics ad Mathematics,

More information

This is an introductory course in Analysis of Variance and Design of Experiments.

This is an introductory course in Analysis of Variance and Design of Experiments. 1 Notes for M 384E, Wedesday, Jauary 21, 2009 (Please ote: I will ot pass out hard-copy class otes i future classes. If there are writte class otes, they will be posted o the web by the ight before class

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

Lecture 6 Simple alternatives and the Neyman-Pearson lemma

Lecture 6 Simple alternatives and the Neyman-Pearson lemma STATS 00: Itroductio to Statistical Iferece Autum 06 Lecture 6 Simple alteratives ad the Neyma-Pearso lemma Last lecture, we discussed a umber of ways to costruct test statistics for testig a simple ull

More information

Minimax Estimation of the Parameter of Maxwell Distribution Under Different Loss Functions

Minimax Estimation of the Parameter of Maxwell Distribution Under Different Loss Functions America Joural of heoretical ad Applied Statistics 6; 5(4): -7 http://www.sciecepublishiggroup.com/j/ajtas doi:.648/j.ajtas.654.6 ISSN: 6-8999 (Prit); ISSN: 6-96 (Olie) Miimax Estimatio of the Parameter

More information

This chapter focuses on two experimental designs that are crucial to comparative studies: (1) independent samples and (2) matched pair samples.

This chapter focuses on two experimental designs that are crucial to comparative studies: (1) independent samples and (2) matched pair samples. Chapter 9 & : Comparig Two Treatmets: This chapter focuses o two eperimetal desigs that are crucial to comparative studies: () idepedet samples ad () matched pair samples Idepedet Radom amples from Two

More information

Chapter 13: Tests of Hypothesis Section 13.1 Introduction

Chapter 13: Tests of Hypothesis Section 13.1 Introduction Chapter 13: Tests of Hypothesis Sectio 13.1 Itroductio RECAP: Chapter 1 discussed the Likelihood Ratio Method as a geeral approach to fid good test procedures. Testig for the Normal Mea Example, discussed

More information

Section 9.2. Tests About a Population Proportion 12/17/2014. Carrying Out a Significance Test H A N T. Parameters & Hypothesis

Section 9.2. Tests About a Population Proportion 12/17/2014. Carrying Out a Significance Test H A N T. Parameters & Hypothesis Sectio 9.2 Tests About a Populatio Proportio P H A N T O M S Parameters Hypothesis Assess Coditios Name the Test Test Statistic (Calculate) Obtai P value Make a decisio State coclusio Sectio 9.2 Tests

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

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

( θ. sup θ Θ f X (x θ) = L. sup Pr (Λ (X) < c) = α. x : Λ (x) = sup θ H 0. sup θ Θ f X (x θ) = ) < c. NH : θ 1 = θ 2 against AH : θ 1 θ 2

( θ. sup θ Θ f X (x θ) = L. sup Pr (Λ (X) < c) = α. x : Λ (x) = sup θ H 0. sup θ Θ f X (x θ) = ) < c. NH : θ 1 = θ 2 against AH : θ 1 θ 2 82 CHAPTER 4. MAXIMUM IKEIHOOD ESTIMATION Defiitio: et X be a radom sample with joit p.m/d.f. f X x θ. The geeralised likelihood ratio test g.l.r.t. of the NH : θ H 0 agaist the alterative AH : θ H 1,

More information

NEW FAST CONVERGENT SEQUENCES OF EULER-MASCHERONI TYPE

NEW FAST CONVERGENT SEQUENCES OF EULER-MASCHERONI TYPE UPB Sci Bull, Series A, Vol 79, Iss, 207 ISSN 22-7027 NEW FAST CONVERGENT SEQUENCES OF EULER-MASCHERONI TYPE Gabriel Bercu We itroduce two ew sequeces of Euler-Mascheroi type which have fast covergece

More information

Introduction to Extreme Value Theory Laurens de Haan, ISM Japan, Erasmus University Rotterdam, NL University of Lisbon, PT

Introduction to Extreme Value Theory Laurens de Haan, ISM Japan, Erasmus University Rotterdam, NL University of Lisbon, PT Itroductio to Extreme Value Theory Laures de Haa, ISM Japa, 202 Itroductio to Extreme Value Theory Laures de Haa Erasmus Uiversity Rotterdam, NL Uiversity of Lisbo, PT Itroductio to Extreme Value Theory

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

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

Regression with an Evaporating Logarithmic Trend

Regression with an Evaporating Logarithmic Trend Regressio with a Evaporatig Logarithmic Tred Peter C. B. Phillips Cowles Foudatio, Yale Uiversity, Uiversity of Aucklad & Uiversity of York ad Yixiao Su Departmet of Ecoomics Yale Uiversity October 5,

More information

of the matrix is =-85, so it is not positive definite. Thus, the first

of the matrix is =-85, so it is not positive definite. Thus, the first BOSTON COLLEGE Departmet of Ecoomics EC771: Ecoometrics Sprig 4 Prof. Baum, Ms. Uysal Solutio Key for Problem Set 1 1. Are the followig quadratic forms positive for all values of x? (a) y = x 1 8x 1 x

More information

6 Sample Size Calculations

6 Sample Size Calculations 6 Sample Size Calculatios Oe of the major resposibilities of a cliical trial statisticia is to aid the ivestigators i determiig the sample size required to coduct a study The most commo procedure for determiig

More information