Section 6.1

Performance Evaluations in Large Samples

[39/0%]
Show
LaTeX
Problem 6.1.1

Let X1,…,XnX_{1}, \ldots , X_{n} be iid with E[Xi]=ξ\mathbb {E}\left[X_{i}\right]=\xi.

?
(a)

If the Xi sX_{i} \mathrm{~ s} have a finite fourth moment, establish (1.3)

(b)

For kk a positive integer, show that E[(Xˉ−ξ)2k−1]\mathbb {E}\left[(\bar{X}-\xi )^{2 k-1}\right] and E[(Xˉ−ξ)2k]\mathbb {E}\left[(\bar{X}-\xi )^{2 k}\right], if they exist, are both O(1/nk)O\left(1 / n^{k}\right).

Problem 6.1.2

For fixed nn, describe the relative error in Example 1.3 as a function of pp.

?
Problem 6.1.3

Theorem 1.1 states: let X1,…,XnX_1, \ldots , X_n be iid with E[X1]=ξ\mathbb {E}\left[X_1\right]=\xi, Var⁡[X1]=σ2\operatorname {Var}\left[X_1\right]=\sigma^2, and finite fourth moment, and suppose hh is a function of a real variable whose first four derivatives h′,h′′,h′′′,h(iv)h', h'', h''', h^{(iv)} exist for all x∈Ix \in I, where II is an interval with P(X1∈I)=1\mathbb {P}\left(X_1 \in I\right)=1. Suppose further that ∣h(iv)(x)∣≤M\left|h^{(iv)}(x)\right| \leq M for all x∈Ix \in I, for some M<∞M<\infty. Then,

E[h(Xˉ)]=h(ξ)+σ22nh′′(ξ)+Rn, \mathbb {E}\left[h(\bar{X})\right]=h(\xi )+\frac{\sigma ^{2}}{2 n} h''(\xi )+R_{n},

and if, in addition, the fourth derivative of h2h^2 is also bounded,

Var⁡[h(Xˉ)]=σ2n[h′(ξ)]2+Rn, \operatorname {Var}\left[h(\bar{X})\right]=\frac{\sigma ^{2}}{n}\left[h'(\xi )\right]^{2}+R_{n},

where the remainder RnR_n in both cases is O(1/n2)O(1/n^2).

Theorem 1.5 states: suppose that the assumptions of Theorem 1.1 hold, and that cnc_n is a sequence of constants satisfying cn=1+an+O ⁣(1n2)c_{n}=1+\dfrac {a}{n}+O\! \left(\dfrac {1}{n^{2}}\right). Then, the variance of δn(Xˉ)=h(cnXˉ)\delta_{n}(\bar{X})=h\left(c_{n} \bar{X}\right) satisfies

Var⁡[δn(Xˉ)]=σ2n[h′(ξ)]2+O ⁣(1n2). \operatorname {Var}\left[\delta _{n}(\bar{X})\right]=\frac{\sigma ^{2}}{n}\left[h'(\xi )\right]^{2}+O\! \left(\frac{1}{n^{2}}\right).

Prove Theorem 1.5.

?
Problem 6.1.4

Let X1,…,XnX_{1}, \ldots , X_{n} be iid as N(ξ,σ2),σ2N\left(\xi , \sigma^{2}\right), \sigma^{2} known, and let g(ξ)=ξr,r=2,3,4g(\xi )=\xi^{r}, r=2,3,4. Determine, up to terms of order 1/n1 / n,

?
(a)

the variance of the UMVU estimator of g(ξ)g(\xi );

(b)

the bias of the MLE of g(ξ)g(\xi ).

Problem 6.1.5

Let X1,…,XnX_{1}, \ldots , X_{n} be iid as N(ξ,σ2),ξN\left(\xi , \sigma^{2}\right), \xi known. For even rr, determine the variance of the UMVU estimator (2.2.4) of σr\sigma^{r} up to terms of order rr.

?
Problem 6.1.6

Solve the preceding problem for the case that ξ\xi is unknown.

?
Problem 6.1.7

For estimating pmp^{m} in Example 3.3.1, determine, up to order 1/n1 / n,

?
(a)

the variance of the UMVU estimator (2.3.2);

(b)

the bias of the MLE.

Problem 6.1.8

Solve the preceding problem if pmp^{m} is replaced by the estimand of Problem 2.3.3.

?
Problem 6.1.9

Let X1,…,XnX_{1}, \ldots , X_{n} be iid as Poisson P(θ)P(\theta ).

?
(a)

Determine the UMVU estimator of P(Xi=0)=e−θ\mathbb {P}\left(X_{i}=0\right)=e^{-\theta }.

(b)

Calculate the variance of the estimator of (a) up to terms of order 1/n1 / n.

Problem 6.1.10

Solve part (b) of the preceding problem for the estimator (2.3.22).

?
Problem 6.1.11

Under the assumptions of Problem 1.1, show that E[∣Xˉ−ξ∣2k−1]=O(n−k+1/2)\mathbb {E}\left[\left|\bar{X}-\xi \right|^{2 k-1}\right]=O\left(n^{-k+1 / 2}\right).

?
Problem 6.1.12

Obtain a variant of Theorem 1.1, which requires existence and boundedness of only h′′′h^{\prime \prime \prime } instead of h(iv)h^{(i v)}, but where RnR_{n} is only O(n−3/2)O\left(n^{-3 / 2}\right).

?
Problem 6.1.13

To see that Theorem 1.1 is not necessarily valid without boundedness of the fourth (or some higher) derivative, suppose that the XX 's are distributed as N(ξ,σ2)N\left(\xi , \sigma^{2}\right) and let h(X)=ex4h(X)=e^{x^{4}}. Then, all moments of the XX 's and all derivatives of hh exist.

?
(a)

Show that the expectation of h(Xˉ)h(\bar{X}) does not exist for any nn, and hence that E[(n[h(Xˉ)−h(ξ)])2]=∞\mathbb {E}\left[\left(\sqrt{n}[h(\bar{X}) -h(\xi )]\right)^{2}\right]=\infty for all values of nn.

(b)

On the other hand, show that n[h(Xˉ)−h(ξ)]\sqrt{n}[h(\bar{X})-h(\xi )] has a normal limit distribution with finite variance, and determine that variance.

Problem 6.1.14

Let X1,…,XnX_{1}, \ldots , X_{n} be iid from the exponential distribution with density (1/θ)e−x/θ,x>0(1 / \theta ) e^{-x / \theta }, x>0, and θ>0\theta >0.

?
(a)

Use Theorem 1.1 to find approximations to E[Xˉ]\mathbb {E}\left[\sqrt{\bar{X}}\right] and Var⁡[Xˉ]\operatorname {Var}\left[\sqrt{\bar{X}}\right].

(b)

Verify the exact calculation

Var⁡[Xˉ]=[1−1n(Γ(n+1/2)Γ(n))2]θ \operatorname {Var}\left[\sqrt{\bar{X}}\right]=\left[1-\frac{1}{n}\left(\frac{\Gamma (n+1 / 2)}{\Gamma (n)}\right)^{2}\right] \theta

and show that lim⁡n→∞nVar⁡[Xˉ]=θ/4\lim_{n \rightarrow \infty } n \operatorname {Var}\left[\sqrt{\bar{X}}\right]=\theta / 4.

(c)

Reconcile the results in parts (a) and (b). Explain why, even though Theorem 1.1 did not apply, it gave the correct answer.

(d)

Show that a similar conclusion holds for h(x)=1/xh(x)=1 / x.

Problem 6.1.15

Let X1,…,XnX_{1}, \ldots , X_{n} be iid according to U(0,θ)U(0, \theta ). Determine the variance of the UMVU estimator of θk\theta^{k}, where kk is an integer, k>−nk>-n.

?
Problem 6.1.16

Under the assumptions of Problem 1.15, find the MLE of θk\theta^{k} and compare its expected squared error with the variance of the UMVU estimator.

?
Problem 6.1.17

Let X1,…,XnX_{1}, \ldots , X_{n} be iid according to U(0,θ)U(0, \theta ), let T=max⁡(X1,…,Xn)T=\max \left(X_{1}, \ldots , X_{n}\right), and let hh be a function satisfying the conditions of Theorem 1.1. Show that

E[h(T)]=h(θ)−θnh′(θ)+1n2[θh′(θ)+θ2h′′(θ)]+O(1n3) \mathbb {E}\left[h(T)\right]=h(\theta )-\frac{\theta }{n} h^{\prime }(\theta )+\frac{1}{n^{2}}\left[\theta h^{\prime }(\theta )+\theta ^{2} h^{\prime \prime }(\theta )\right]+O\left(\frac{1}{n^{3}}\right)

and

Var⁡[h(T)]=θ2n2[h′(θ)]2+O(1n3). \operatorname {Var}\left[h(T)\right]=\frac{\theta ^{2}}{n^{2}}\left[h^{\prime }(\theta )\right]^{2}+O\left(\frac{1}{n^{3}}\right).
?
Problem 6.1.18

Apply the results of Problem 1.17 to obtain approximate answers to Problems 1.15 and 1.16, and compare the answers with the exact solutions.

?
Problem 6.1.19

If the XX 's are as in Theorem 1.1 and if the first five derivatives of hh exist and the fifth derivative is bounded, show that

E[h(Xˉ)]=h(ξ)+12h′′σ2n+124n2[4h′′′μ3+3h(iv)σ4]+O(n−5/2) \mathbb {E}\left[h(\bar{X})\right]=h(\xi )+\frac{1}{2} h^{\prime \prime } \frac{\sigma ^{2}}{n}+\frac{1}{24 n^{2}}\left[4 h^{\prime \prime \prime } \mu _{3}+3 h^{(i v)} \sigma ^{4}\right]+O\left(n^{-5 / 2}\right)

and if the fifth derivative of h2h^{2} is also bounded

Var⁡[h(Xˉ)]=(h′2)σ2n+1n2[h′h′′μ3+(h′h′′′+12h′′2)σ4]+O(n−5/2) \operatorname {Var}\left[h(\bar{X})\right]=\left(h^{\prime 2}\right) \frac{\sigma ^{2}}{n}+\frac{1}{n^{2}}\left[h^{\prime } h^{\prime \prime } \mu _{3}+\left(h^{\prime } h^{\prime \prime \prime }+\frac{1}{2} h^{\prime \prime 2}\right) \sigma ^{4}\right]+O\left(n^{-5 / 2}\right)

where μ3=E[(X−ξ)3]\mu_{3}=\mathbb {E}\left[(X-\xi )^{3}\right].

?
Problem 6.1.20

Under the assumptions of the preceding problem, carry the calculation of the variance (1.16) to terms of order 1/n21 / n^{2}, and compare the result with that of the preceding problem.

?
Problem 6.1.21

Carry the calculation of Problem 1.4 to terms of order 1/n21 / n^{2}.

?
Problem 6.1.22

For the estimands of Problem 1.4, calculate the expected squared error of the MLE to terms of order 1/n21 / n^{2}, and compare it with the variance calculated in Problem 1.21.

?
Problem 6.1.23

Calculate the variance (1.18) to terms of order 1/n21 / n^{2} and compare it with the expected squared error of the MLE carried to the same order.

?
Problem 6.1.24

Find the variance of the estimator (2.3.17) up to terms of the order 1/n31 / n^{3}.

?
Problem 6.1.25

For the situation of Example 1.12, show that the UMVU estimator δ1n\delta_{1 n} is the biascorrected MLE, where the MLE is δ3n\delta_{3 n}.

?
Problem 6.1.26

For the estimators of Example 1.13:

?
(a)

Calculate their exact variances.

(b)

Use the result of part (a) to verify (1.27).

Problem 6.1.27
?
(a)

Under the assumptions of Theorem 1.5, if all fourth moments of the XiνX_{i \nu } are finite, show that E[Xˉi−ξi](Xˉj−ξj)=σij/n\mathbb {E}\left[\bar{X}_{i}-\xi_{i}\right]\left(\bar{X}_{j}-\xi_{j}\right)=\sigma_{i j} / n and that all third and fourth moments E[Xˉi−ξi](Xˉj−ξj)(Xˉk−ξk)\mathbb {E}\left[\bar{X}_{i}-\xi_{i}\right]\left(\bar{X}_{j}-\xi_{j}\right)\left(\bar{X}_{k}-\xi_{k}\right), and so on are of the order 1/n21 / n^{2}.

(b)

If, in addition, all derivatives of hh of total order ≤4\leq 4 exist and those of order 4 are uniformly bounded, then

E[h(Xˉ1,…,Xˉs)]=h(ξ1,…,ξs)+12n∑i=1s∑j=1sσij∂2h(ξ1,…,ξs)∂ξi∂ξj+Rn, \mathbb {E}\left[h\left(\bar{X}_{1}, \ldots , \bar{X}_{s}\right)\right]=h\left(\xi _{1}, \ldots , \xi _{s}\right)+\frac{1}{2 n} \sum _{i=1}^{s} \sum _{j=1}^{s} \sigma _{i j} \frac{\partial ^{2} h\left(\xi _{1}, \ldots , \xi _{s}\right)}{\partial \xi _{i} \partial \xi _{j}}+R_{n},

and if the derivatives of h2h^{2} of order 4 are also bounded,

Var⁡[h(Xˉ1,…,Xˉs)]=1nΣΣσij∂h∂ξi∂h∂ξj+Rn \operatorname {Var}\left[h\left(\bar{X}_{1}, \ldots , \bar{X}_{s}\right)\right]=\frac{1}{n} \Sigma \Sigma \sigma _{i j} \frac{\partial h}{\partial \xi _{i}} \frac{\partial h}{\partial \xi _{j}}+R_{n}

where the remainder RnR_{n} in both cases is O(1/n2)O\left(1 / n^{2}\right).

Problem 6.1.28

On the basis of a sample from N(ξ,σ2)N\left(\xi , \sigma^{2}\right), let Pn(ξ,σ)P_{n}(\xi , \sigma ) be the probability that the UMVU estimator Xˉ2−σ2/n\bar{X}^{2}-\sigma^{2} / n of ξ2\xi^{2} ( σ\sigma known) is negative.

?
(a)

Show that Pn(ξ,σ)P_{n}(\xi , \sigma ) is a decreasing function of n∣ξ∣/σ\sqrt{n}\left|\xi \right| / \sigma.

(b)

Show that Pn(ξ,σ)→0P_{n}(\xi , \sigma ) \rightarrow 0 as n→∞n \rightarrow \infty for any fixed ξ≠0\xi \neq 0 and σ\sigma.

(c)

Determine the value of Pn(0,σ)P_{n}(0, \sigma ).

Problem 6.1.29

Use the tt-distribution to find the value of Pn(0,σ)P_{n}(0, \sigma ) in the preceding problem for the UMVU estimator of ξ2\xi^{2} when σ\sigma is unknown for representative values of nn.

?
Problem 6.1.30

Theorem 1.9 states: if a sequence of random variables TnT_{n} tends to ξ\xi in probability and if hh is continuous at ξ\xi, then h(Tn)→Ph(ξ)h(T_{n}) \xrightarrow {P} h(\xi ). To show that P(∣h(Tn)−h(ξ)∣<a)→1\mathbb {P}\left(\left|h(T_{n})-h(\xi )\right|<a\right) \rightarrow 1, it is enough to notice that by the continuity of hh, the difference ∣h(Tn)−h(ξ)∣\left|h(T_{n})-h(\xi )\right| will be arbitrarily small if TnT_{n} is close to ξ\xi, and that for any aa, the probability that ∣Tn−ξ∣<a\left|T_{n}-\xi \right|<a tends to 11 as n→∞n \rightarrow \infty.

Fill in the details of the proof of Theorem 1.9. (See also Problem 1.8.8.)

?
Problem 6.1.31

In Example 8.13 with θ=0\theta =0, show that δ2n\delta_{2 n} is not exactly distributed as σ2(χ12−1)/n\sigma^{2}\left(\chi_{1}^{2}-1\right) / n.

?
Problem 6.1.32

In Example 8.13, let δ4n=max⁡(0,Xˉ2−σ2/n)\delta_{4 n}=\max \left(0, \bar{X}^{2}-\sigma^{2} / n\right), which is an improvement over δ1n\delta_{1 n}.

?
(a)

Show that n(δ4n−θ2)\sqrt{n}\left(\delta_{4 n}-\theta^{2}\right) has the same limit distribution as n(δ1n−θ2)\sqrt{n}\left(\delta_{1 n}-\theta^{2}\right) when θ≠0\theta \neq 0.

(b)

Describe the limit distribution of nδ4nn \delta_{4 n} when θ=0\theta =0.

Problem 6.1.33

Let XX have the binomial distribution b(p,n)b(p, n), and let g(p)=pqg(p)=p q. The UMVU estimator of g(p)g(p) is δ=X(n−X)/n(n−1)\delta =X(n-X) / n(n-1). Determine the limit distribution of n(δ−pq)\sqrt{n}(\delta - p q) and n(δ−pq)n(\delta -p q) when g′(p)≠0g^{\prime }(p) \neq 0 and g′(p)=0g^{\prime }(p)=0, respectively.

?
Problem 6.1.34

Let X1,…,XnX_{1}, \ldots , X_{n} be iid as N(ξ,1)N(\xi , 1). Determine the limit behavior of the distribution of the UMVU estimator of p=P(∣Xi∣≤u)p=\mathbb {P}\left(\left|X_{i}\right| \leq u\right).

?
Problem 6.1.35

Determine the limit behavior of the estimator (2.3.22) as n→∞n \rightarrow \infty.

?
Problem 6.1.36

Let X1,…,XnX_{1}, \ldots , X_{n} be iid with distribution PθP_{\theta }, and suppose δn\delta_{n} is UMVU for estimating g(θ)g(\theta ) on the basis of X1,…,XnX_{1}, \ldots , X_{n}. If there exists n0n_{0} and an unbiased estimator δ0(X1,…,Xn0)\delta_{0}\left(X_{1}, \ldots , X_{n_{0}}\right) which has finite variance for all θ\theta, then δn\delta_{n} is consistent for g(θ)g(\theta ).

?
Problem 6.1.37

Let YnY_{n} be distributed as N(0,1)N(0,1) with probability πn\pi_{n} and as N(0,τn2)N\left(0, \tau_{n}^{2}\right) with probability 1−πn1-\pi_{n}. If τn→∞\tau_{n} \rightarrow \infty and πn→π\pi_{n} \rightarrow \pi, determine for what values of π\pi the sequence {Yn}\left\{ Y_{n}\right\} does and does not have a limit distribution.

?
Problem 6.1.38
?
(a)

In Problem 1.37, determine to what values Var⁡[Yn]\operatorname {Var}\left[Y_{n}\right] can tend as n→∞n \rightarrow \infty if πn→1\pi_{n} \rightarrow 1 and τn→∞\tau_{n} \rightarrow \infty but otherwise both are arbitrary.

(b)

Use (a) to show that the limit of the variance need not agree with the variance of the limit distribution.

Problem 6.1.39

Let bm,n,m,n=1,2,…b_{m, n}, m, n=1,2, \ldots, be a double sequence of real numbers, which for each fixed mm is nondecreasing in nn. Show that lim⁡n→∞lim⁡m→∞bm,n=lim⁡m,n→∞inf⁡bm,n\lim_{n \rightarrow \infty } \lim_{m \rightarrow \infty } b_{m, n}=\lim_{m, n \rightarrow \infty } \inf b_{m, n} and lim⁡m→∞lim⁡n→∞bm,n=lim⁡m,n→∞sup⁡bm,n\lim_{m \rightarrow \infty } \lim_{n \rightarrow \infty } b_{m, n}=\lim_{m, n \rightarrow \infty } \sup b_{m, n} provided the indicated limits exist (they may be infinite) and where lim⁡inf⁡bm,n\lim \inf b_{m, n} and lim⁡sup⁡bm,n\lim \sup b_{m, n} denote, respectively, the smallest and the largest limit points attainable by a sequence bmk,nk,k=1,2,…b_{m_{k}, n_{k}}, k=1,2, \ldots, with mk→∞m_{k} \rightarrow \infty and nk→∞n_{k} \rightarrow \infty.

?