25.03.2013 Views

Normal approximation to the hypergeometric distribution in ...

Normal approximation to the hypergeometric distribution in ...

Normal approximation to the hypergeometric distribution in ...

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Note that for all k, n satisfy<strong>in</strong>g |ak,n|,<br />

S.N. Lahiri et al. / Journal of Statistical Plann<strong>in</strong>g and Inference 137 (2007) 3570 –3590 3583<br />

(1 − f)<br />

Np − k Np − (np + (1 − f)npq)>np > 0<br />

2<br />

and<br />

(1 − f)<br />

Nq − (n − k)>nq > 0.<br />

2<br />

Hence, by <strong>the</strong> error bound <strong>in</strong> Stirl<strong>in</strong>g’s <strong>approximation</strong>, for all k, n with |ak,n| and 6(np ∧ nq)1,<br />

∗ 1<br />

<br />

12Np + 1 −<br />

1<br />

12(Np − k) +<br />

1<br />

12Nq + 1 −<br />

1<br />

12(Nq − (n − k)) +<br />

1 1<br />

−<br />

12(N − n) + 1 12N<br />

−<br />

−<br />

12k + 1<br />

(12Np + 1)(12(Np − k)) −<br />

12(n − k) + 1<br />

(12Nq + 1)(12(Nq − n + k))<br />

1<br />

6Np(1 − )(1 − f) −<br />

1<br />

6Nq(1 − )(1 − f)<br />

f<br />

=−<br />

6npq(1 − )(1 − f) ,<br />

∗ <br />

<br />

1 1<br />

f<br />

0 + 0 +<br />

− <br />

12(N − n) + 1 12N 6npq(1 − )(1 − f) ,<br />

∗∗ 1<br />

12n −<br />

1<br />

12k + 1 −<br />

1<br />

12(n − k) + 1 0,<br />

∗∗ 1 1<br />

−<br />

12n + 1 12k −<br />

1<br />

n<br />

1<br />

−<br />

−<br />

12(n − k) 12k(n − k) 6npq(1 − ) .<br />

Hence, <strong>the</strong> lemma follows from (4.21) and <strong>the</strong> above <strong>in</strong>equalities. <br />

Lemma 2. Let g : R −→ [0, ∞) be such that g is ↑ on (−∞,a)andgis↓ on (a, ∞) for some a ∈ R. Then, for any<br />

k ∈ N, b ∈ R and h ∈ (0, ∞),<br />

k<br />

g(b + ih)<br />

i=o<br />

b+hk<br />

where g(x0) = max{g(b + ih) : i = 0, 1,...,k}.<br />

Proof. For ba, by mono<strong>to</strong>nicity,<br />

h<br />

k<br />

g(b + ih)hg(b) +<br />

i=0<br />

b<br />

g(x)dx + 2hg(x0), (4.22)<br />

b+hk<br />

b<br />

g(x)dx.<br />

For b

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!