Chapter 2

Unbiasedness

[136/0%]
Show
LaTeX
Problem 2.1.1

Verify

?
(a)

that (1.4) defines a probability distribution and

(b)

condition (1.5).

Problem 2.1.2

In Example 1.5, show that ai∗a_{i}^{*} minimizes (1.6) for i=0,1i=0,1, and simplify the expression for a0∗a_{0}^{*}.

?
Problem 2.1.3

Let XX take on the values -1, 0, 1, 2, 3 with probabilities P(X=−1)=2pq\mathbb {P}\left(X=-1\right)=2 p q and P(X=k)=pkq3−k\mathbb {P}\left(X=k\right)=p^{k} q^{3-k} for k=0,1,2,3k=0,1,2,3.

?
(a)

Check that this is a probability distribution.

(b)

Determine the LMVU estimator at p0p_{0} of (i) pp, and (ii) pqp q, and decide for each whether it is UMVU.

Problem 2.1.4

For a sample of size nn, suppose that the estimator T(x)T(\mathbf{x}) of τ(θ)\tau (\theta ) has expectation

E[T(X)]=τ(θ)+∑k=1∞aknk, \mathbb {E}\left[T(\mathbf{X})\right]=\tau (\theta )+\sum _{k=1}^{\infty } \frac{a_{k}}{n^{k}},

where aka_{k} may depend on θ\theta but not on nn.

?
(a)

Show that the expectation of the jackknife estimator TJT_{J} of (1.3) is

E[TJ(X)]=τ(θ)−a2n2+O(1/n3). \mathbb {E}\left[T_{J}(\mathbf{X})\right]=\tau (\theta )-\frac{a_{2}}{n^{2}}+O\left(1 / n^{3}\right).
(b)

Show that if var T∼c/nT \sim c / n for some constant cc, then var TJ∼c′/nT_{J} \sim c^{\prime } / n for some constant c′c^{\prime }. Thus, the jackknife will reduce bias and not increase variance.

A second-order jackknife can be defined by jackknifing TJT_{J}, and this will result in further bias reduction, but may not maintain a variance of the same order (Robson and Whitlock 1964; see also Thorburn 1976 and Note 8.3).

Problem 2.1.5
?
(a)

Any two random variables XX and YY with finite second moments satisfy the covariance inequality [cov⁡(X,Y)]2≤var⁡(X)⋅var⁡(Y)[\operatorname {cov}(X, Y)]^{2} \leq \operatorname {var}(X) \cdot \operatorname {var}(Y).

(b)

The inequality in part (a) is an equality if and only if there exist constants aa and bb for which P(X=aY+b)=1\mathbb {P}\left(X=a Y+b\right)=1.

Problem 2.1.6

An alternative proof of the Schwarz inequality is obtained by noting that

∫(f+λg)2dP=∫f2dP+2λ∫fgdP+λ2∫g2dP≥0 for all λ, \int (f+\lambda g)^{2} d P=\int f^{2} d P+2 \lambda \int f g d P+\lambda ^{2} \int g^{2} d P \geq 0 \text{ for all } \lambda ,

so that this quadratic in λ\lambda has at most one root.

?
Problem 2.1.7

Suppose XX is distributed on ( 0,1 ) with probability density pθ(x)=(1−θ)+θ/2xp_{\theta }(x)=(1-\theta )+\theta / 2 \sqrt{x} for all 0<x<1,0≤θ≤10<x<1,0 \leq \theta \leq 1. Show that there does not exist an LMVU estimator of θ\theta.

?
Problem 2.1.8

If δ\delta and δ′\delta^{\prime } have finite variance, so does δ′−δ\delta^{\prime }-\delta.

?
Problem 2.1.9

In Example 1.9,

?
(a)

determine all unbiased estimators of zero;

(b)

show that no nonconstant estimator is UMVU.

Problem 2.1.10

If estimators are restricted to the class of linear estimators, characterization of best unbiased estimators is somewhat easier. Although the following is a consequence of Theorem 1.7, it should be established without using that theorem.

Let Xp×1\mathbf{X}_{p \times 1} satisfy E[X]=Bψ\mathbb {E}\left[\mathbf{X}\right]=B \psi and var⁡(X)=I\operatorname {var}(\mathbf{X})=I, where Bp×rB_{p \times r} is known, and ψr×1\psi_{r \times 1} is unknown. A linear estimator is an estimator of the form a′Xa^{\prime } \mathbf{X}, where ap×1a_{p \times 1} is a known vector. We are concerned with the class of estimators

D={δ(x):δ(x)=a′x, for some known vector a}. \mathcal{D}=\left\{ \delta (\mathbf{x}): \delta (\mathbf{x})=a^{\prime } \mathbf{x}, \text{ for some known vector } a\right\} .
?
(a)

For a known vector cc, show that the estimators in D\mathcal{D} that are unbiased estimators of c′ψc^{\prime } \psi satisfy a′B=c′a^{\prime } B=c^{\prime }.

(b)

Let Dc={δ(x):δ(x)=a′x,a′B=c′}\mathcal{D}_{c}=\left\{ \delta (\mathbf{x}): \delta (\mathbf{x})=a^{\prime } \mathbf{x}, a^{\prime } B=c^{\prime }\right\} be the class of linear unbiased estimators of c′ψc^{\prime } \psi. Show that the best linear unbiased estimator (BLUE) of c′ψc^{\prime } \psi, the linear unbiased estimator with minimum variance, is δ∗(x)=a∗′x\delta^{*}(\mathbf{x})=a^{* \prime } \mathbf{x}, where a∗′=a′B(B′B)−1B′a^{* \prime }=a^{\prime } B\left(B^{\prime } B\right)^{-1} B^{\prime } and a∗′B=ca^{* \prime } B=c with variance var⁡(δ∗)=c′c\operatorname {var}\left(\delta^{*}\right)=c^{\prime } c.

(c)

Let D0={δ(x):δ(x)=a′x,a′B=0.}\mathcal{D}_{0}=\left\{ \delta (\mathbf{x}): \delta (\mathbf{x})=a^{\prime } \mathbf{x}, a^{\prime } B=0.\right\} be the class of linear unbiased estimators of zero. Show that if δ∈D0\delta \in \mathcal{D}_{0}, then cov⁡(δ,δ∗)=0\operatorname {cov}\left(\delta , \delta^{*}\right)=0.

(d)

Hence, establish the analog of Theorem 1.7 for linear estimators: Theorem. An estimator δ∗∈Dc\delta^{*} \in \mathcal{D}_{c} satisfies var⁡(δ∗)=min⁡δ∈Dcvar⁡(δ)\operatorname {var}\left(\delta^{*}\right)=\min_{\delta \in \mathcal{D}_{c}} \operatorname {var}(\delta ) if and only if cov⁡(δ∗,U)=0\operatorname {cov}\left(\delta^{*}, U\right)=0, where UU is any estimator in D0\mathcal{D}_{0}.

(e)

Show that the results here can be directly extended to the case of var⁡(X)=Σ\operatorname {var}(\mathbf{X})=\Sigma, where Σp×p\Sigma_{p \times p} is a known matrix, by considering the transformed problem with X∗=Σ1/2X\mathbf{X}^{*}=\Sigma^{1 / 2} \mathbf{X} and B∗=Σ1/2BB^{*}=\Sigma^{1 / 2} B.

Problem 2.1.11

Use Theorem 1.7 to find UMVU estimators of some of the ηθ(di)\eta_{\theta }\left(d_{i}\right) in the doseresponse model (1.6.16), with the restriction (1.6.17) (Messig and Strawderman 1993). Let the classes Δ\Delta and U\mathcal{U} be defined as in Theorem 1.7.

?
(a)

Show that an estimator U∈UU \in \mathcal{U} if and only if U(x1,x2)=a[I(x1=0)−I(x2=0)]U\left(x_{1}, x_{2}\right)=a\left[I\left(x_{1}=0\right)-I\left(x_{2}=0\right)\right] for an arbitrary constant a<∞a<\infty.

(b)

Using part (a) and (1.7), show that an estimator δ\delta is UMVU for its expectation only if it is of the form δ(x1,x2)=aI(0,0)(x1,x2)+bI(0,1),(1,0),(2,0)(x1,x2)+cI(1,1)(x1,x2)+dI(2,1)(x1,x2)\delta \left(x_{1}, x_{2}\right)=a I_{(0,0)}\left(x_{1}, x_{2}\right)+b I_{(0,1),(1,0),(2,0)}\left(x_{1}, x_{2}\right)+c I_{(1,1)}\left(x_{1}, x_{2}\right)+ d I_{(2,1)}\left(x_{1}, x_{2}\right) where a,b,ca, b, c, and dd are arbitrary constants.

(c)

Show that there does not exist a UMVU estimator of ηθ(d1)=1−e−θ\eta_{\theta }\left(d_{1}\right)=1-e^{-\theta }, but the UMVU estimator of ηθ(d2)=1−e−2θ\eta_{\theta }\left(d_{2}\right)=1-e^{-2 \theta } is δ(x1,x2)=1−12[I(x1=0)+I(x2=0)]\delta \left(x_{1}, x_{2}\right)=1-\frac{1}{2}\left[I\left(x_{1}=0\right)+I\left(x_{2}=0\right)\right].

(d)

Show that the LMVU estimator of 1−e−θ1-e^{-\theta } is δ(x1,x2)=x12+12(1+e−θ)[I(x1=0)−I(x2=0)]\delta \left(x_{1}, x_{2}\right)=\frac{x_{1}}{2}+\frac{1}{2\left(1+e^{-\theta }\right)}\left[I\left(x_{1}=0\right)-I\left(x_{2}=0\right)\right].

Problem 2.1.12

Show that if δ(X)\delta (X) is a UMVU estimator of g(θ)g(\theta ), it is the unique UMVU estimator of g(θ)g(\theta ). (Do not assume completeness, but rather use the covariance inequality and the conditions under which it is an equality.)

?
Problem 2.1.13

If δ1\delta_{1} and δ2\delta_{2} are in Δ\Delta and are UMVU estimators of g1(θ)g_{1}(\theta ) and g2(θ)g_{2}(\theta ), respectively, then a1δ1+a2δ2a_{1} \delta_{1}+a_{2} \delta_{2} is also in Δ\Delta and is UMVU for estimating a1g1(θ)+a2g2(θ)a_{1} g_{1}(\theta )+a_{2} g_{2}(\theta ), for any real a1a_{1} and a2a_{2}.

?
Problem 2.1.14

Completeness of TT is not only sufficient but also necessary so that every g(θ)g(\theta ) that can be estimated unbiasedly has only one unbiased estimator that is a function of TT.

?
Problem 2.1.15

Suppose X1,…,XnX_{1}, \ldots , X_{n} are iid Poisson (λ)(\lambda ).

?
(a)

Show that Xˉ\bar{X} is the UMVU estimator for λ\lambda.

(b)

For S2=∑i=1n(Xi−Xˉ)2/(n−1)S^{2}=\sum_{i=1}^{n}\left(X_{i}-\bar{X}\right)^{2} /(n-1), we have that E[S2]=E[Xˉ]=λ\mathbb {E}\left[S^{2}\right]=\mathbb {E}\left[\bar{X}\right]=\lambda. To directly establish that var S2>S^{2}> var Xˉ\bar{X}, prove that E[S2∣Xˉ]=Xˉ\mathbb {E}\left[S^{2} \mid \bar{X}\right]=\bar{X}.

Note: The identity E[S2∣Xˉ]=Xˉ\mathbb {E}\left[S^{2} \mid \bar{X}\right]=\bar{X} shows how completeness can be used in calculating conditional expectations.

Problem 2.1.16
?
(a)

If X1,…,XnX_{1}, \ldots , X_{n} are iid (not necessarily normal) with var⁡(Xi)=σ2<∞\operatorname {var}\left(X_{i}\right)=\sigma^{2}<\infty, show that δ=Σ(Xi−Xˉ)2/(n−1)\delta =\Sigma \left(X_{i}-\bar{X}\right)^{2} /(n-1) is an unbiased estimator of σ2\sigma^{2}.

(b)

If the XiX_{i} take on the values 1 and 0 with probabilities pp and q=1−pq=1-p, the estimator δ\delta of (a) depends only on T=ΣXiT=\Sigma X_{i} and hence is UMVU for estimating σ2=pq\sigma^{2}=p q. Compare this result with that of Example 1.13.

Problem 2.1.17

If TT has the binomial distribution b(p,n)b(p, n) with n>3n>3, use Method 1 to find the UMVU estimator of p3p^{3}.

?
Problem 2.1.18

Let X1,…,XnX_{1}, \ldots , X_{n} be iid according to the Poisson distribution P(λ)P(\lambda ). Use Method 1 to find the UMVU estimator of

?
(a)

λk\lambda^{k} for any positive integer kk and

(b)

e−λe^{-\lambda }.

Problem 2.1.19

Let X1,…,XnX_{1}, \ldots , X_{n} be distributed as in Example 1.14. Use Method 1 to find the UMVU estimator of θk\theta^{k} for any integer k>−nk>-n.

?
Problem 2.1.20

Solve Problem 1.18(b) by Method 2, using the fact that an unbiased estimator of e−λe^{-\lambda } is δ=1\delta =1 if X1=0X_{1}=0, and δ=0\delta =0 otherwise.

?
Problem 2.1.21

In nn Bernoulli trials, let Xi=1X_{i}=1 or 0 as the ii th trial is a success or failure, and let T=ΣXiT=\Sigma X_{i}. Solve Problem 1.17 by Method 2, using the fact that an unbiased estimator of p3p^{3} is δ=1\delta =1 if X1=X2=X3=1X_{1}=X_{2}=X_{3}=1, and δ=0\delta =0 otherwise.

?
Problem 2.1.22

Let XX take on the values 1 and 0 with probability pp and qq, respectively, and assume that 1/4<p<3/41 / 4<p<3 / 4. Consider the problem of estimating pp with loss function L(p,d)=L(p, d)= 1 if ∣d−p∣≥1/4\left|d-p\right| \geq 1 / 4, and 0 otherwise. Let δ∗\delta^{*} be the randomized estimator which is Y0Y_{0} or Y1Y_{1} when X=0X=0 or 1 where Y0Y_{0} and Y1Y_{1} are distributed as U(−1/2,1/2)U(-1 / 2,1 / 2) and U(1/2,3/2)U(1 / 2,3 / 2), respectively.

?
(a)

Show that δ∗\delta^{*} is unbiased.

(b)

Compare the risk function of δ∗\delta^{*} with that of XX.

Problem 2.2.1

If X1,…,XnX_{1}, \ldots , X_{n} are iid as N(ξ,σ2)N\left(\xi , \sigma^{2}\right) with σ2\sigma^{2} known, find the UMVU estimator of

?
(a)

ξ2\xi^{2},

(b)

ξ3\xi^{3}, and

(c)

ξ4\xi^{4}.

Problem 2.2.2

Solve the preceding problem when σ\sigma is unknown.

?
Problem 2.2.3

In Example 2.1 with σ\sigma known, let δ=ΣciXi\delta =\Sigma c_{i} X_{i} be any linear estimator of ξ\xi. If δ\delta is biased, show that its risk E[(δ−ξ)2]\mathbb {E}\left[(\delta -\xi )^{2}\right] is unbounded.

?
Problem 2.2.4

Suppose, as in Example 2.1, that X1,…,XnX_{1}, \ldots , X_{n} are iid as N(ξ,σ2)N\left(\xi , \sigma^{2}\right), with one of the parameters known, and that the estimand is a polynomial in ξ\xi or σ\sigma. Then, the UMVU estimator is a polynomial in Xˉ\bar{X} or S2=∑(Xi−ξ)2S^{2}=\sum \left(X_{i}-\xi \right)^{2}. The variance of any such polynomial can be estimated if one knows the moments E[Xˉk]\mathbb {E}\left[\bar{X}^{k}\right] and E[Sk]\mathbb {E}\left[S^{k}\right] for all k=1,2,…k=1,2, \ldots. To determine E[Xˉk]\mathbb {E}\left[\bar{X}^{k}\right], write Xˉ=Y+ξ\bar{X}=Y+\xi, where YY is distributed as N(0,σ2/n)N\left(0, \sigma^{2} / n\right). Show that

?
(a)
E[Xˉk]=∑r=0k(kr)ξk−rE[Yr] \mathbb {E}\left[\bar{X}^{k}\right]=\sum _{r=0}^{k}\binom {k}{r} \xi ^{k-r} \mathbb {E}\left[Y^{r}\right]

with

E[Yr]={(r−1)(r−3)⋯3⋅1(σ2/n)r/2 when r≥2 is even 0 when r is odd . \mathbb {E}\left[Y^{r}\right]= \begin{cases} (r-1)(r-3) \cdots 3 \cdot 1\left(\sigma ^{2} / n\right)^{r / 2} & \text{ when } r \geq 2 \text{ is even } \\ 0 & \text{ when } r \text{ is odd }.\end{cases}
(b)

As an example, consider the UMVU estimator S2/nS^{2} / n of σ2\sigma^{2}. Show that E[S4]=n(n+2)σ2\mathbb {E}\left[S^{4}\right]= n(n+2) \sigma^{2} and var⁡(S2n)=2σ4n\operatorname {var}\left(\frac{S^{2}}{n}\right)=\frac{2 \sigma^{4}}{n} and that the UMVU estimator of this variance is 2S4/n2(n+2)2 S^{4} / n^{2}(n+2).

Problem 2.2.5

In Example 2.1, when both parameters are unknown, show that the UMVU estimator of ξ2\xi^{2} is given by δ=Xˉ2−S2n(n−1)\delta =\bar{X}^{2}-\frac{S^{2}}{n(n-1)} where now S2=∑(Xi−Xˉ)2S^{2}=\sum \left(X_{i}-\bar{X}\right)^{2}.

?
Problem 2.2.6
?
(a)

Determine the variance of the estimator Problem 2.5.

(b)

Find the UMVU estimator of the variance in part (a).

Problem 2.2.7

If XX is a single observation from N(ξ,σ2)N\left(\xi , \sigma^{2}\right), show that no unbiased estimator δ\delta of σ2\sigma^{2} exists when ξ\xi is unknown.

?
Problem 2.2.8

Let Xi,i=1,…,nX_{i}, i=1, \ldots , n, be independently distributed as N(α+βti,σ2)N\left(\alpha +\beta t_{i}, \sigma^{2}\right) where α,β\alpha , \beta, and σ2\sigma^{2} are unknown, and the tt 's are known constants that are not all equal. Find the UMVU estimators of α\alpha and β\beta.

?
Problem 2.2.9

In Example 2.2 with n=1n=1, the UMVU estimator of pp is the indicator of the event X1≤uX_{1} \leq u whether σ\sigma is known or unknown.

?
Problem 2.2.10

Verify Equation (2.14), the density of (X1−Xˉ)/S\left(X_{1}-\bar{X}\right) / S in normal sampling.

?
Remark.
?

The UMVU estimator in (2.13) is used by Kiefer (1977) as an example of his estimated confidence approach.

Problem 2.2.11

Assuming (2.15) with σ=τ\sigma =\tau, determine the UMVU estimators of σ2\sigma^{2} and (η−ξ)/σ(\eta -\xi ) / \sigma.

?
Problem 2.2.12

Assuming (2.15) with η=ξ\eta =\xi and σ2/τ2=γ\sigma^{2} / \tau^{2}=\gamma, show that when γ\gamma is known:

?
(a)

T′T^{\prime } defined in Example 2.3(iii) is a complete sufficient statistic;

(b)

δγ\delta_{\gamma } is UMVU for ξ\xi.

Problem 2.2.13

Show that in the preceding problem with γ\gamma unknown,

?
(a)

a UMVU estimator of ξ\xi does not exist;

(b)

the estimator ξ^\hat{\xi } is unbiased under the conditions stated in Example 2.3.

Problem 2.2.14

For the model (2.15) find the UMVU estimator of P(X1<Y1)\mathbb {P}\left(X_{1}<Y_{1}\right) when

?
(a)

σ=τ\sigma =\tau and

(b)

σ\sigma and τ\tau are arbitrary.

Problem 2.2.15

If (X1,Y1),…,(Xn,Yn)\left(X_{1}, Y_{1}\right), \ldots ,\left(X_{n}, Y_{n}\right) are iid according to any bivariate distribution with finite second moments, show that SXY/(n−1)S_{X Y} /(n-1) given by (2.17) is an unbiased estimator of cov⁡(Xi,Yi)\operatorname {cov}\left(X_{i}, Y_{i}\right).

?
Problem 2.2.16

In a sample size N=n+k+1N=n+k+1, some of the observations are missing. Assume that (Xi,Yi),i=1,…,n\left(X_{i}, Y_{i}\right), i=1, \ldots , n, are iid according to the bivariate normal distribution (2.16), and that U1,…,UkU_{1}, \ldots , U_{k} and V1,…,VlV_{1}, \ldots , V_{l} are independent N(ξ,σ2)N\left(\xi , \sigma^{2}\right) and N(η,τ2)N\left(\eta , \tau^{2}\right), respectively.

?
(a)

Show that the minimal sufficient statistics are complete when ξ\xi and η\eta are known but not when they are unknown.

(b)

When ξ\xi and η\eta are known, find the UMVU estimators for σ2,τ2\sigma^{2}, \tau^{2}, and ρστ\rho \sigma \tau, and suggest reasonable unbiased estimators for these parameters when ξ\xi and η\eta are unknown.

Problem 2.2.17

For the family (2.22), show that the UMVU estimator of aa when bb is known and the UMVU estimator of bb is known are as stated in Example 2.5.

?
Problem 2.2.18

Show that the estimators (2.23) are UMVU.

?
Problem 2.2.19

For the family (2.22) with b=1b=1, find the UMVU estimator of P(X1≥u)\mathbb {P}\left(X_{1} \geq u\right) and of the density e−(u−a)e^{-(u-a)} of X1X_{1} at uu.

?
Problem 2.2.20

Find the UMVU estimator of P(X1≥u)\mathbb {P}\left(X_{1} \geq u\right) for the family (2.22) when both aa and bb are unknown.

?
Problem 2.2.21

Let X1,…,XmX_{1}, \ldots , X_{m} and Y1,…,YnY_{1}, \ldots , Y_{n} be independently distributed as E(a,b)E(a, b) and E(a′,b′)E(a^{\prime }, b^{\prime }), respectively.

?
(a)

If a,b,a′a, b, a^{\prime }, and b′b^{\prime } are completely unknown, X(1),Y(1),Σ[Xi−X(1)]X_{(1)}, Y_{(1)}, \Sigma \left[X_{i}-X_{(1)}\right], and Σ[Yj−Y(1)]\Sigma \left[Y_{j}-Y_{(1)}\right] jointly are sufficient and complete.

(b)

Find the UMVU estimators of a′−aa^{\prime }-a and b′/bb^{\prime } / b.

Problem 2.2.22

In the preceding problem, suppose that b′=bb^{\prime }=b.

?
(a)

Show that X(1),Y(1)X_{(1)}, Y_{(1)}, and Σ[Xi−X(1)]+Σ[Yj−Y(1)]\Sigma \left[X_{i}-X_{(1)}\right]+\Sigma \left[Y_{j}-Y_{(1)}\right] are sufficient and complete.

(b)

Find the UMVU estimators of bb and (a′−a)/b\left(a^{\prime }-a\right) / b.

Problem 2.2.23

In Problem 2.21, suppose that a′=aa^{\prime }=a.

?
(a)

Show that the complete sufficient statistic of Problem 2.21(a) is still minimal sufficient but no longer complete.

(b)

Show that a UMVU estimator for a′=aa^{\prime }=a does not exist.

(c)

Suggest a reasonable unbiased estimator for a′=aa^{\prime }=a.

Problem 2.2.24

Let X1,…,XnX_{1}, \ldots , X_{n} be iid according to the uniform distribution U(ξ−b,ξ+b)U(\xi -b, \xi +b). If ξ,b\xi , b are both unknown, find the UMVU estimators of ξ,b\xi , b, and ξ/b\xi / b.

?
Problem 2.2.25

Let X1,…,XmX_{1}, \ldots , X_{m} and Y1,…,YnY_{1}, \ldots , Y_{n} be iid as U(0,θ)U(0, \theta ) and U(0,θ′)U\left(0, \theta^{\prime }\right), respectively. If n>1n>1, determine the UMVU estimator of θ/θ′\theta / \theta^{\prime }.

?
Problem 2.2.26

Verify the ML estimators given in (2.24).

?
Problem 2.2.27

In Example 2.6(b), show that

?
(a)

The bias of the ML estimator is 0 when ξ=u\xi =u.

(b)

At ξ=u\xi =u, the ML estimator has smaller expected squared error than the UMVU estimator.

Problem 2.2.28

Verify (2.26).

?
Problem 2.2.29

Under the assumptions of Lemma 2.7, show that:

?
(a)

If bb is replaced by any random variable BB which is independent of XX and not 0 with probability 1, then Rδ(θ)<Rδ∗(θ)R_{\delta }(\theta )<R_{\delta^{*}}(\theta ).

(b)

If squared error is replaced by any loss function of the form L(θ,δ)=ρ(d−θ)L(\theta , \delta )=\rho (d-\theta ) and δ\delta is risk unbiased with respect to LL, then Rδ(θ)<Rδ∗(θ)R_{\delta }(\theta )<R_{\delta^{*}}(\theta ).

Problem 2.3.1
?
(a)

In Example 3.1, show that Σ(Xi−Xˉ)2=T(n−T)/n\Sigma \left(X_{i}-\bar{X}\right)^{2}=T(n-T) / n.

(b)

The variance of T(n−T)/n(n−1)T(n-T) / n(n-1) in Example 3.1 is (pq/n)[(q−p)2+2pq/(n−1)](p q / n)\left[(q-p)^{2}+2 p q /(n-1)\right].

Problem 2.3.2

If TT is distributed as b(p,n)b(p, n), find an unbiased estimator δ(T)\delta (T) of pm(m≤n)p^{m}(m \leq n) by Method 1, that is, using (1.10).

?
Problem 2.3.3
?
(a)

Use the method leading to (3.2) to find the UMVU estimator πk(T)\pi_{k}(T) of P(X1+⋯+Xm=k)=(mk)pkqm−k(m≤n)\mathbb {P}\left(X_{1}+\cdots +X_{m}=k\right)=\binom {m}{k} p^{k} q^{m-k}(m \leq n).

(b)

For fixed tt and varying kk, show that the πk(t)\pi_{k}(t) are the probabilities of a hypergeometric distribution.

Problem 2.3.4

If YY is distributed according to (3.3), use Method 1 of Section 2.1

?
(a)

to show that the UMVU estimator of pr(r<m)p^{r}(r<m) is

δ(y)=(m−r+y−1)(m−r+y−2)…(m−r)(m+y−1)(m+y−2)⋯m, \delta (y)=\frac{(m-r+y-1)(m-r+y-2) \ldots (m-r)}{(m+y-1)(m+y-2) \cdots m},

and hence in particular that the UMVU estimator of 1/p,1/p21 / p, 1 / p^{2} and pp are, respectively, (m+y)/m,(m+y)(m+y+1)/m(m+1)(m+y) / m,(m+y)(m+y+1) / m(m+1), and (m−1)/(m+y−1)(m-1) /(m+y-1);

(b)

to determine the UMVU estimator of var⁡(Y)\operatorname {var}(Y);

(c)

to show how to calculate the UMVU estimator δ\delta of log⁡p\log p.

Problem 2.3.5

Consider the scheme in which binomial sampling is continued until at least aa successes and bb failures have been obtained. Show how to calculate a reasonable estimator of log⁡(p/q)\log (p / q).

?
Problem 2.3.6

If binomial sampling is continued until mm successes have been obtained, let XiX_{i} ( i=1,…,mi= 1, \ldots , m ) be the number of failures between the (i−1)(i-1) st and ii th success.

?
(a)

The XiX_{i} are iid according to the geometric distribution P(Xi=x)=pqx,x=0,1,…\mathbb {P}\left(X_{i}=x\right)=p q^{x}, x= 0,1, \ldots.

(b)

The statistic Y=ΣXiY=\Sigma X_{i} is sufficient for (X1,…,Xm)\left(X_{1}, \ldots , X_{m}\right) and has the distribution (3.3).

Problem 2.3.7

Suppose that binomial sampling is continued until the number of successes equals the number of failures.

?
(a)

This rule is closed if p=1/2p=1 / 2 but not otherwise.

(b)

If p=1/2p=1 / 2 and NN denotes the number of trials required, E[N]=∞\mathbb {E}\left[N\right]=\infty.

Problem 2.3.8

Verify Equation (3.7) with the appropriate definition of N′(x,y)N^{\prime }(x, y)

?
(a)

for the estimation of pp and

(b)

for the estimation of paqbp^{a} q^{b}.

Problem 2.3.9

Consider sequential binomial sampling with the stopping points (0,1)(0,1) and (2,y)(2, y), y=0,1,…y=0,1, \ldots.

?
(a)

Show that this plan is closed and simple.

(b)

Show that (X,Y)(X, Y) is not complete by finding a nontrivial unbiased estimator of zero.

Problem 2.3.10

In Example 3.4(ii),

?
(a)

show that the plan is closed but not simple,

(b)

show that (X,Y)(X, Y) is not complete, and

(c)

evaluate the unbiased estimator (3.7) of pp.

Problem 2.3.11

Curtailed single sampling. Let a,b<na, b<n be three non-negative integers. Continue observation until either aa successes, bb failures, or nn observations have been obtained. Determine the UMVU estimator of pp.

?
Problem 2.3.12

For any sequential binomial sampling plan, the coordinates (X,Y)(X, Y) of the end point of the sample path are minimal sufficient.

?
Problem 2.3.13

Consider any closed sequential binomial sampling plan with a set BB of stopping points, and let B′B^{\prime } be the set B∪{(x0,y0)}B \cup \left\{ \left(x_{0}, y_{0}\right)\right\} where (x0,y0)\left(x_{0}, y_{0}\right) is a point not in BB that has positive probability of being reached under plan BB. Show that the sufficient statistic T=(X,Y)T=(X, Y) is not complete for the sampling plan which has B′B^{\prime } as its set of stopping points.

?
Problem 2.3.14

For any sequential binomial sampling plan under which the point (1,1)(1,1) is reached with positive probability but is not a stopping point, find an unbiased estimator of pqp q depending only on (X,Y)(X, Y). Evaluate this estimator for

?
(a)

taking a sample of fixed size n>2n>2;

(b)

inverse binomial sampling.

Problem 2.3.15

Use (3.3) to determine A(t,n)A(t, n) in (3.11) for the negative binomial distribution with m=nm=n, and evaluate the estimators (3.13) of qrq^{r}, and (3.14).

?
Problem 2.3.16

Consider nn binomial trials with success probability pp, and let rr and ss be two positive integers with r+s<nr+s<n. To the boundary x+y=nx+y=n, add the boundary point (r,s)(r, s), that is, if the number of successes in the first r+sr+s trials is exactly rr, the process is stopped and the remaining n−(r+s)n-(r+s) trials are not performed.

?
(a)

Show that UU is an unbiased estimator of zero if and only if U(k,n−k)=0U(k, n-k)=0 for k=0,1,…,r−1k=0,1, \ldots , r-1 and k=n−s+1,n−s+2,…,nk=n-s+1, n-s+2, \ldots , n, and U(k,n−k)=ckU(r,s)U(k, n-k)=c_{k} U(r, s) for k=r,…,n−sk=r, \ldots , n-s, where the cc 's are given constants ≠0\neq 0.

(b)

Show that δ\delta is the UMVU estimator of its expectation if and only if

δ(k,n−k)=δ(r,s) for k=r,…,n−s. \delta (k, n-k)=\delta (r, s) \text{ for } k=r, \ldots , n-s.
Problem 2.3.17

Generalize the preceding problem to the case that two points (r1,s1)\left(r_{1}, s_{1}\right) and (r2,s2)\left(r_{2}, s_{2}\right) with ri+si<nr_{i}+s_{i}<n are added to the boundary. Assume that these two points are such that all n+1n+1 points x+y=nx+y=n remain boundary points.

?
Problem 2.3.18

If XX has the Poisson distribution P(θ)P(\theta ), show that 1/θ1 / \theta does not have an unbiased estimator.

?
Problem 2.3.19

If X1,…,XnX_{1}, \ldots , X_{n} are iid according to (3.18), the Poisson distribution truncated on the left at 0, find the UMVU estimator of θ\theta when

?
(a)

n=1n=1 and

(b)

n=2n=2.

Problem 2.3.20

Let X1,…,XnX_{1}, \ldots , X_{n} be a sample from the Poisson distribution truncated on the left at 0, i.e., with sample space X={1,2,3,…}\mathcal{X}=\left\{ 1,2,3, \ldots \right\}.

?
(a)

For t=Σxit=\Sigma x_{i}, the UMVU estimator of λ\lambda is (Tate and Goen 1958) λ^=Ct−1nCtnt\hat{\lambda }=\frac{\mathcal{C}_{t-1}^{n}}{\mathcal{C}_{t}^{n}} t where Ctn=(−1)nn!∑k=0∞(nk)(−1)kkt\mathcal{C}_{t}^{n}=\frac{(-1)^{n}}{n!} \sum_{k=0}^{\infty }\binom {n}{k}(-1)^{k} k^{t} is a Stirling number of the second kind.

(b)

An alternate form of the UMVU estimator is λ^=tn(1−Ct−1n−1Ctn)\hat{\lambda }=\frac{t}{n}\left(1-\frac{\mathcal{C}_{t-1}^{n-1}}{\mathcal{C}_{t}^{n}}\right).

(c)

The Cramér-Rao lower bound for the variance of unbiased estimators of λ\lambda is λ(1−e−λ)2/[n(1−e−λ−λe−λ)]\lambda \left(1-e^{-\lambda }\right)^{2} /\left[n\left(1-e^{-\lambda }-\lambda e^{-\lambda }\right)\right], and it is not attained by the UMVU estimator. (It is, however, the asymptotic variance of the ML estimator.)

Problem 2.3.21

Suppose that XX has the Poisson distribution truncated on the right at aa, so that it has the conditional distribution of YY given Y≤aY \leq a, where YY is distributed as P(λ)P(\lambda ). Show that λ\lambda does not have an unbiased estimator.

?
Problem 2.3.22

For the negative binomial distribution truncated at zero, evaluate the estimators (3.13) and (3.14) for m=1,2m=1,2, and 3.

?
Problem 2.3.23

If X1,…,XnX_{1}, \ldots , X_{n} are iid P(λ)P(\lambda ), consider estimation of e−bλe^{-b \lambda }, where bb is known.

?
(a)

Show that δ∗=(1−b/n)t\delta^{*}=(1-b / n)^{t} is the UMVU estimator of e−bλe^{-b \lambda }.

(b)

For b>nb>n, describe the behavior of δ∗\delta^{*}, and suggest why it might not be a reasonable estimator.

(The probability e−bλe^{-b \lambda }, for b>nb>n, is that of an "unobservable" event, in that it can be interpreted as the probability of no occurrence in a time interval of length bb. A number of such situations are described and analyzed in Lehmann (1983), where it is suggested that, in these problems, no reasonable estimator may exist.)

Problem 2.3.24

If X1,…,XnX_{1}, \ldots , X_{n} are iid according to the logarithmic series distribution of Problem 1.5.14, evaluate the estimators (3.13) and (3.14) for n=1,2n=1,2, and 3.

?
Problem 2.3.25

For the multinomial distribution of Example 3.8,

?
(a)

show that p0r0⋯psrsp_{0}^{r_{0}} \cdots p_{s}^{r_{s}} has an unbiased estimator provided r0,…,rsr_{0}, \ldots , r_{s} are nonnegative integers with Σri≤n\Sigma r_{i} \leq n;

(b)

find the totality of functions that can be estimated unbiasedly;

(c)

determine the UMVU estimator of the estimand of (a).

Problem 2.3.26

In Example 3.9 when pij=pi+p+jp_{i j}=p_{i+} p_{+j}, determine the variances of the two unbiased estimators δ0=nij/n\delta_{0}=n_{i j} / n and δ1=ni+n+j/n2\delta_{1}=n_{i+} n_{+j} / n^{2} of pijp_{i j}, and show directly that var⁡(δ0)>var⁡(δ1)\operatorname {var}\left(\delta_{0}\right)>\operatorname {var}\left(\delta_{1}\right) for all n>1n>1.

?
Problem 2.3.27

In Example 3.9, show that independence of AA and BB implies that (n1+,…,nI+)(n_{1+}, \ldots , n_{I^{+}}) and (n+1,…,n+J)\left(n_{+1}, \ldots , n_{+J}\right) are independent with multinomial distributions as stated.

?
Problem 2.3.28

Verify (3.20).

?
Problem 2.3.29

Let X,YX, Y, and gg be such that E[g(X,Y)∣y]\mathbb {E}\left[g(X, Y) \mid y\right] is independent of yy. Then, E[f(Y)g(X,Y)]=E[f(Y)]E[g(X,Y)]\mathbb {E}\left[f(Y) g(X, Y)\right]=\mathbb {E}\left[f(Y)\right] \mathbb {E}\left[g(X, Y)\right], and hence f(Y)f(Y) and g(X,Y)g(X, Y) are uncorrelated, for all ff.

?
Problem 2.3.30

In Example 3.10, show that the estimator δ1\delta_{1} of pijkp_{i j k} is unbiased for the model (3.20).

?
Problem 2.4.1

Let X1,…,XnX_{1}, \ldots , X_{n} be iid with distribution FF.

?
(a)

Characterize the totality of functions f(X1,…,Xn)f\left(X_{1}, \ldots , X_{n}\right) which are unbiased estimators of zero for the class F0\mathcal{F}_{0} of all distributions FF having a density.

(b)

Give one example of a nontrivial unbiased estimator of zero when (i) n=2n=2 and (ii) n=3n=3.

Problem 2.4.2

Let F\mathcal{F} be the class of all univariate distribution functions FF that have a probability density function ff and finite mm th moment.

?
(a)

Let X1,…,XnX_{1}, \ldots , X_{n} be independently distributed with common distribution F∈FF \in \mathcal{F}. For n≥mn \geq m, find the UMVU estimator of ξm\xi^{m} where ξ=ξ(F)=E[Xi]\xi =\xi (F)=\mathbb {E}\left[X_{i}\right].

(b)

Show that for the case that P(Xi=1)=p,P(Xi=0)=q,p+q=1\mathbb {P}\left(X_{i}=1\right)=p, \mathbb {P}\left(X_{i}=0\right)=q, p+q=1, the estimator of (a) reduces to (3.2).

Problem 2.4.3

In the preceding problem, show that 1/var⁡FXi1 / \operatorname {var}_{F} X_{i} does not have an unbiased estimator for any nn.

?
Problem 2.4.4

Let X1,…,XnX_{1}, \ldots , X_{n} be iid with distribution F∈FF \in \mathcal{F} where F\mathcal{F} is the class of all symmetric distributions with a probability density. There exists no UMVU estimator of the center of symmetry θ\theta of FF (if unbiasedness is required only for the distributions FF for which the expectation of the estimator exists).

?
Problem 2.4.5

If X1,…,XmX_{1}, \ldots , X_{m} and Y1,…,YnY_{1}, \ldots , Y_{n} are independently distributed according to FF and G∈F0G \in \mathcal{F}_{0}, defined in Problem 4.1, the order statistics X(1)<⋯<X(m)X_{(1)}<\cdots <X_{(m)} and Y(1)<⋯<Y(n)Y_{(1)}<\cdots < Y_{(n)} are sufficient and complete.

?
Problem 2.4.6

Under the assumptions of the preceding problem, find the UMVU estimator of P(Xi<Yj)\mathbb {P}\left(X_{i}<Y_{j}\right).

?
Problem 2.4.7

Under the assumptions of Problem 4.5, let ξ=E[Xi]\xi =\mathbb {E}\left[X_{i}\right] and η=E[Yj]\eta =\mathbb {E}\left[Y_{j}\right]. Show that ξ2η2\xi^{2} \eta^{2} possesses an unbiased estimator if and only if m≥2m \geq 2 and n≥2n \geq 2.

?
Problem 2.4.8

Let (X1,Y1),…,(Xn,Yn)\left(X_{1}, Y_{1}\right), \ldots ,\left(X_{n}, Y_{n}\right) be iid F∈FF \in \mathcal{F}, where F\mathcal{F} is the family of all distributions with probability density and finite second moments. Show that δ(X,Y)=∑(Xi−Xˉ)(Yi−Yˉ)/(n−1)\delta (X, Y)=\sum \left(X_{i}-\bar{X}\right)\left(Y_{i}-\bar{Y}\right) /(n-1) is UMVU for cov⁡(X,Y)\operatorname {cov}(X, Y).

?
Problem 2.4.9

Under the assumptions of the preceding problem, find the UMVU estimator of

?
(a)

P(Xi≤Yi)\mathbb {P}\left(X_{i} \leq Y_{i}\right);

(b)

P(Xi≤Xj and Yi≤Yj),i≠j\mathbb {P}\left(X_{i} \leq X_{j} \text{ and } Y_{i} \leq Y_{j}\right), i \neq j.

Problem 2.4.10

Let (X1,Y1),…,(Xn,Yn)\left(X_{1}, Y_{1}\right), \ldots ,\left(X_{n}, Y_{n}\right) be iid with F∈FF \in \mathcal{F}, where F\mathcal{F} is the family of all bivariate densities. Show that the sufficient statistic TT, which generalizes the order statistics to the bivariate case, is complete.

?
Problem 2.5.1

Under the assumptions of Problem 1.3, determine for each p1p_{1}, the value LV(p1)L_{V}\left(p_{1}\right) of the LMVU estimator of pp at p1p_{1} and compare the function LV(p),0<p<1L_{V}(p), 0<p<1 with the variance Vp0(p)V_{p_{0}}(p) of the estimator which is LMVU at

?
(a)

p0=1/3p_{0}=1 / 3 and

(b)

p0=1/2p_{0}=1 / 2.

Problem 2.5.2

Determine the conditions under which equality holds in (5.1).

?
Problem 2.5.3

Verify I(θ)I(\theta ) for the distributions of Table 5.1.

?
Problem 2.5.4

If XX is normal with mean zero and standard deviation σ\sigma, determine I(σ)I(\sigma ).

?
Problem 2.5.5

Find I(p)I(p) for the negative binomial distribution.

?
Problem 2.5.6

If XX is distributed as P(λ)P(\lambda ), show that the information it contains about λ\sqrt{\lambda } is independent of λ\lambda.

?
Problem 2.5.7

Verify the following statements, asserted by Basu (1988, Chapter 1), which illustrate the relationship between information, sufficiency, and ancillarity. Suppose that we let I(θ)=Eθ[−∂2/∂θ2log⁡f(x∣θ)]I(\theta )=\mathbb {E}_{\theta }\left[-\partial^{2} / \partial \theta^{2} \log f(x \mid \theta )\right] be the information in XX about θ\theta and let J(θ)=Eθ[−∂2/∂θ2log⁡g(T∣θ)]J(\theta )=\mathbb {E}_{\theta }\left[-\partial^{2} / \partial \theta^{2} \log g(T \mid \theta )\right] be the information about θ\theta contained in a statistic TT, where g(⋅∣θ)g(\cdot \mid \theta ) is the density function of TT. Define λ(θ)=I(θ)−J(θ)\lambda (\theta )=I(\theta )-J(\theta ), a measure of information lost by using TT instead of XX. Under suitable regularity conditions, show that

?
(a)

λ(θ)≥0\lambda (\theta ) \geq 0 for all θ\theta

(b)

λ(θ)=0\lambda (\theta )=0 if and only if TT is sufficient for θ\theta.

(c)

If YY is ancillary but (T,Y)(T, Y) is sufficient, then I(θ)=Eθ[J(θ∣Y)]I(\theta )=\mathbb {E}_{\theta }\left[J(\theta \mid Y)\right], where

J(θ∣y)=Eθ[−∂2∂θ2log⁡h(T∣y,θ)∣ Y=y] J(\theta \mid y)=\mathbb {E}_{\theta }\left[\left.-\frac{\partial ^{2}}{\partial \theta ^{2}} \log h(T \mid y, \theta ) \right\rvert \, Y=y\right]

and h(t∣y,θ)h(t \mid y, \theta ) is the conditional density of TT given Y=yY=y.

(Basu's "regularity conditions" are mainly concerned with interchange of integration and differentiation. Assume any such interchanges are valid.)

Problem 2.5.8

Find a function of θ\theta for which the amount of information is independent of θ\theta :

?
(a)

for the gamma distribution Γ(α,β)\Gamma (\alpha , \beta ) with α\alpha known and with θ=β\theta =\beta;

(b)

for the binomial distribution b(p,n)b(p, n) with θ=p\theta =p.

Problem 2.5.9

For inverse binomial sampling (see Example 3.2):

?
(a)

Show that the best unbiased estimator of pp is given by δ∗(Y)=(m−1)/(Y+m−1)\delta^{*}(Y)=(m-1) /(Y+m-1).

(b)

Show that the information contained in YY about PP is I(p)=m/p2(1−p)I(p)=m / p^{2}(1-p).

(c)

Show that var⁡δ∗>1/I(p)\operatorname {var} \delta^{*}>1 / I(p).

(The estimator δ∗\delta^{*} can be interpreted as the success rate if we ignore the last trial, which we know must be a success.)

Problem 2.5.10

Show that (5.13) can be differentiated by differentiating under the integral sign when pθ(x)p_{\theta }(x) is given by (5.24), for each of the distributions of Table 5.2.

?
Problem 2.5.11

Verify the entries of Table 5.2.

?
Problem 2.5.12

Evaluate (5.25) when ff is the density of Student's tt-distribution with vv degrees of freedom.

?
Problem 2.5.13

For the distribution with density (5.24), show that I(θ)I(\theta ) is independent of θ\theta.

?
Problem 2.5.14

Verify

?
(a)

formula (5.25) and

(b)

formula (5.27).

Problem 2.5.15

For the location tt density, calculate the information inequality bound for unbiased estimators of θ\theta.

?
Problem 2.5.16
?
(a)

For the scale family with density (1/θ)f(x/θ),θ>0(1 / \theta ) f(x / \theta ), \theta >0, the amount of information a single observation XX has about θ\theta is

1θ2∫[yf′(y)f(y)+1]2f(y)dy \frac{1}{\theta ^{2}} \int \left[\frac{y f^{\prime }(y)}{f(y)}+1\right]^{2} f(y) d y
(b)

Show that the information XX contains about ξ=log⁡θ\xi =\log \theta is independent of θ\theta.

(c)

For the Cauchy distribution C(0,θ),I(θ)=1/(2θ2)\mathrm{C}(0, \theta ), I(\theta )=1 /\left(2 \theta^{2}\right).

Problem 2.5.17

If pθ(x)p_{\theta }(x) is given by 1.5.1 with s=1s=1 and T(x)=δ(x)T(x)=\delta (x), show that var⁡[δ(X)]\operatorname {var}[\delta (X)] attains the lower bound (5.31) and is the only estimator to do so.

?
Problem 2.5.18

Show that if a given function g(θ)g(\theta ) has an unbiased estimator, there exists an unbiased estimator δ\delta which for all θ\theta values attains the lower bound (5.1) for some ψ(x,θ)\psi (x, \theta ) satisfying (5.2) if and only if g(θ)g(\theta ) has a UMVU estimator δ0\delta_{0}. (Blyth 1974).

?
Problem 2.5.19

Show that if Eθ[δ]=g(θ)\mathbb {E}_{\theta }\left[\delta \right]=g(\theta ), and var⁡(δ)\operatorname {var}(\delta ) attains the information inequality bound (5.31), then

δ(x)=g(θ)+g′(θ)I(θ)∂∂θpθ(x). \delta (x)=g(\theta )+\frac{g^{\prime }(\theta )}{I(\theta )} \frac{\partial }{\partial \theta } p_{\theta }(x).
?
Problem 2.5.20

If Eθ[δ]=g(θ)\mathbb {E}_{\theta }\left[\delta \right]=g(\theta ), the information inequality lower bound is IB(θ)=[g′(θ)]2/I(θ)I B(\theta )=\left[g^{\prime }(\theta )\right]^{2} / I(\theta ). If θ=h(ξ)\theta =h(\xi ) where hh is differentiable, show that IB(ξ)=IB(θ)I B(\xi )=I B(\theta ).

?
Problem 2.5.21

Let XX be an observation from the normal mixture density

pθ(x)=122π{e−(1/2)(x−θ)2+e−(1/2)(x+θ)2},θ∈Ω, p_{\theta }(x)=\frac{1}{2 \sqrt{2 \pi }}\left\{ e^{-(1 / 2)(x-\theta )^{2}}+e^{-(1 / 2)(x+\theta )^{2}}\right\} , \theta \in \Omega ,

where Ω\Omega is any neighborhood of zero. Thus, the random variable XX is either N(θ,1)N(\theta , 1) or N(−θ,1)N(-\theta , 1), each with probability 1/21 / 2. Show that θ=0\theta =0 is a singular point, that is, if there exists an unbiased estimator of θ\theta it will have infinite variance at θ=0\theta =0.

?
Problem 2.5.22

Let X1,…,XnX_{1}, \ldots , X_{n} be a sample from the Poisson (λ)(\lambda ) distribution truncated on the left at 0, i.e., with sample space X={1,2,3,…}\mathcal{X}=\left\{ 1,2,3, \ldots \right\} (see Problem 3.20). Show that the Cramér-Rao lower bound for the variance of unbiased estimators of λ\lambda is

λ(1−e−λ)2n(1−e−λ−λe−λ) \frac{\lambda \left(1-e^{-\lambda }\right)^{2}}{n\left(1-e^{-\lambda }-\lambda e^{-\lambda }\right)}

and is not attained by the UMVU estimator. (It is, however, the asymptotic variance of the ML estimator.)

?
Problem 2.5.23

Let X1,…,XnX_{1}, \ldots , X_{n} be iid according to a density p(x,θ)p(x, \theta ) which is positive for all xx. Then, the variance of any unbiased estimator δ\delta of θ\theta satisfies

var⁡θ0(δ)≥(θ−θ0)2{∫−∞∞[p(x,θ)]2p(x,θ0)}n−1,θ≠θ0 \operatorname {var}_{\theta _{0}}(\delta ) \geq \frac{\left(\theta -\theta _{0}\right)^{2}}{\left\{ \int _{-\infty }^{\infty } \frac{[p(x, \theta )]^{2}}{p\left(x, \theta _{0}\right)}\right\} ^{n}-1}, \quad \theta \neq \theta _{0}
?
Problem 2.5.24

If X1,…,XnX_{1}, \ldots , X_{n} are iid as N(θ,σ2)N\left(\theta , \sigma^{2}\right) where σ\sigma is known and θ\theta is known to have one of the values 0, ±1, ±2, . . ., the inequality of the preceding problem shows that any unbiased estimator δ\delta of the restricted parameter θ\theta satisfies

var⁡θ0(δ)≥Δ2enΔ2/σ2−1,Δ≠0 \operatorname {var}_{\theta _{0}}(\delta ) \geq \frac{\Delta ^{2}}{e^{n \Delta ^{2} / \sigma ^{2}}-1}, \quad \Delta \neq 0

where Δ=θ−θ0\Delta =\theta -\theta_{0}, and hence sup⁡Δ≠0var⁡θ0(δ)≥1/[en/σ2−1]\sup_{\Delta \neq 0} \operatorname {var}_{\theta_{0}}(\delta ) \geq 1 /\left[e^{n / \sigma^{2}}-1\right].

?
Problem 2.5.25

Under the assumptions of the preceding problem, let Xˉ∗\bar{X}^{*} be the integer closest to Xˉ\bar{X}.

?
(a)

The estimator Xˉ∗\bar{X}^{*} is unbiased for the restricted parameter θ\theta.

(b)

There exist positive constants aa and bb such that for all sufficiently large n,var⁡θ(Xˉ∗)≤ae−bnn, \operatorname {var}_{\theta }\left(\bar{X}^{*}\right) \leq a e^{-b n} for all integers θ\theta.

Problem 2.5.26

Kiefer inequality.

?
(a)

Let XX have density (with respect to μ)p(x,θ)\mu ) p(x, \theta ) which is > 0 for all xx, and let Λ1\Lambda_{1} and Λ2\Lambda_{2} be two distributions on the real line with finite first moments. Then, any unbiased estimator δ\delta of θ\theta satisfies

var⁡(δ)≥[∫ΔdΛ1(Δ)−∫ΔdΛ2(Δ)]2∫ψ2(x,θ)p(x,θ)dμ(x) \operatorname {var}(\delta ) \geq \frac{\left[\int \Delta d \Lambda _{1}(\Delta )-\int \Delta d \Lambda _{2}(\Delta )\right]^{2}}{\int \psi ^{2}(x, \theta ) p(x, \theta ) d \mu (x)}

where

ψ(x,θ)=∫Ωθp(x,θ+Δ)[dΛ1(Δ)−dΛ2(Δ)]p(x,θ) \psi (x, \theta )=\frac{\int _{\Omega _{\theta }} p(x, \theta +\Delta )\left[d \Lambda _{1}(\Delta )-d \Lambda _{2}(\Delta )\right]}{p(x, \theta )}

with Ωθ={Δ:θ+Δ∈Ω}\Omega_{\theta }=\left\{ \Delta : \theta +\Delta \in \Omega \right\}.

(b)

If Λ1\Lambda_{1} and Λ2\Lambda_{2} assign probability 1 to Δ=0\Delta =0 and Δ\Delta, respectively, the inequality reduces to (5.6) with g(θ)=θg(\theta )=\theta. (Kiefer 1952.)

Problem 2.5.27

Verify directly that the following families of densities satisfy (5.38).

?
(a)

The exponential family of (1.5.1),

pη(x)=h(x)eηT(x)−A(η). p_{\eta }(x)=h(x) e^{\eta T(x)-A(\eta )}.
(b)

The location tt family of Example 5.16.

(c)

The logistic density of Table 1.4.1.

Problem 2.5.28

Extend condition (5.38) to vector-valued parameters, and show that it is satisfied by the exponential family (1.5.1) for s>1s>1.

?
Problem 2.5.29

Show that the assumption (5.36(b)) implies (5.38), so Theorem 5.15 is, in fact, a corollary of Theorem 5.10.

?
Problem 2.5.30

Show that (5.38) is satisfied if either of the following is true:

?
(a)

∣∂log⁡pθ/∂θ∣\left|\partial \log p_{\theta } / \partial \theta \right| is bounded.

(b)

[pθ+Δ(x)−pθ(x)]/Δ→∂log⁡pθ/∂θ\left[p_{\theta +\Delta }(x)-p_{\theta }(x)\right] / \Delta \rightarrow \partial \log p_{\theta } / \partial \theta uniformly.

Problem 2.5.31
?
(a)

Show that if (5.38) holds, then the family of densities is strongly differentiable (see Note 8.6).

(b)

Show that weak differentiability is implied by strong differentiability.

Problem 2.5.32

Brown and Gajek (1990) give two different sufficient conditions for (8.2) to hold, which are given below. Show that each implies (8.2). (Note that, in the progression from (a) to (b) the conditions become weaker, thus more widely applicable and harder to check.)

?
(a)

For some B<∞B<\infty,

Eθ0[∂2∂θ2pθ(X)/pθ0(X)]2<B \mathbb {E}_{\theta _{0}}\left[\frac{\partial ^{2}}{\partial \theta ^{2}} p_{\theta }(X) / p_{\theta _{0}}(X)\right]^{2}<B

for all θ\theta in a neighborhood of θ\theta.

(b)

If pt∗(x)=∂/∂θpθ(x)∣θ=tp_{t}^{*}(x)=\partial /\left.\partial \theta p_{\theta }(x)\right|_{\theta =t}, then

lim⁡Δ→0Eθ0[pθ0+Δ∗(X)−pθ0∗(X)pθ0(X)]2=0. \lim _{\Delta \rightarrow 0} \mathbb {E}_{\theta _{0}}\left[\frac{p_{\theta _{0}+\Delta }^{*}(X)-p_{\theta _{0}}^{*}(X)}{p_{\theta _{0}}(X)}\right]^{2}=0.
Problem 2.5.33

Let F\mathcal{F} be the class of all unimodal symmetric densities or, more generally, densities symmetric around zero and satisfying f(x)≤f(0)f(x) \leq f(0) for all xx. Show that

min⁡f∈F∫x2f(x)dx=112, \min _{f \in \mathcal{F}} \int x^{2} f(x) d x=\frac{1}{12},

and that the minimum is attained by the uniform (−12,12)\left(-\frac{1}{2}, \frac{1}{2}\right) distribution. Thus, the uniform distribution has minimum variance among symmetric unimodal distributions. (See Example 4.8.6 for large-sample properties of the scale uniform.)

?
Problem 2.6.1

For any random variables (ψ1,…,ψs)\left(\psi_{1}, \ldots , \psi_{s}\right), show that the matrices ∥Eψiψj∥\left\| E \psi_{i} \psi_{j}\right\| and C=∥cov⁡(ψi,ψj)∥C= \left\| \operatorname {cov}\left(\psi_{i}, \psi_{j}\right)\right\| are positive semidefinite.

?
Problem 2.6.2

In this problem, we establish some facts about eigenvalues and eigenvectors of square matrices. (For a more general treatment, see, for example, Marshall and Olkin 1979, Chapter 20.) We use the facts that a scalar λ>0\lambda >0 is an eigenvalue of the n×nn \times n symmetric matrix AA if there exists an n×1n \times 1 vector pp, the corresponding eigenvector, satisfying Ap=λpA p=\lambda p. If A is nonsingular, there are nn eigenvalues with corresponding linearly independent eigenvectors.

?
(a)

Show that A=P′DλPA=P^{\prime } D_{\lambda } P, where DλD_{\lambda } is a diagonal matrix of eigenvalues of AA and PP is and n×nn \times n matrix whose rows are the corresponding eigenvalues that satisfies P′P=PP′=IP^{\prime } P=P P^{\prime }=I, the identity matrix.

(b)

Show that max⁡xx′Axx′x=\max_{x} \frac{x^{\prime } A x}{x^{\prime } x}= largest eigenvalue of AA.

(c)

If BB is a nonsingular symmetric matrix with eigenvector-eigenvalue representation B=Q′DβQB=Q^{\prime } D_{\beta } Q, then max⁡xx′Axx′Bx=\max_{x} \frac{x^{\prime } A x}{x^{\prime } B x}= largest eigenvalue of A∗A^{*}, where A∗=Dβ−1/2QAQ′Dβ−1/2A^{*}= D_{\beta }^{-1 / 2} Q A Q^{\prime } D_{\beta }^{-1 / 2} and Dβ−1/2D_{\beta }^{-1 / 2} is a diagonal matrix whose elements are the reciprocals of the square roots of the eigenvalues of BB.

(d)

For any square matrices CC and DD, show that the eigenvalues of the matrix CDC D are the same as the eigenvalues of the matrix DCD C, and hence that max⁡xx′Axx′Bx=\max_{x} \frac{x^{\prime } A x}{x^{\prime } B x}= largest eigenvalue of AB−1A B^{-1}.

(e)

If A=aa′A=a a^{\prime }, where aa is a n×1n \times 1 vector ( AA is thus a rank-one matrix), then max⁡xx′aa′xx′Bx=a′B−1a\max_{x} \frac{x^{\prime } a a^{\prime } x}{x^{\prime } B x}= a^{\prime } B^{-1} a.

Problem 2.6.3

An alternate proof of Theorem 6.1 uses the method of Lagrange (or undetermined) multipliers. Show that, for fixed γ\gamma, the maximum value of a′γa^{\prime } \gamma, subject to the constraint that a′Ca=1a^{\prime } C a=1, is obtained by the solutions to

∂∂ai{a′γ−12λ[a′Ca−1]}=0, \frac{\partial }{\partial a_{i}}\left\{ a^{\prime } \gamma -\frac{1}{2} \lambda \left[a^{\prime } C a-1\right]\right\} =0,

where λ\lambda is the undetermined multiplier. (The solution is a=±C−1γ/γ′C−1γa= \pm C^{-1} \gamma / \sqrt{\gamma^{\prime } C^{-1} \gamma }.)

?
Problem 2.6.4

Prove (6.11) under the assumptions of the text.

?
Problem 2.6.5

Verify

?
(a)

the information matrices of Table 6.1 and

(b)

Equations (6.15) and (6.16).

Problem 2.6.6

If p(x)=(1−ε)ϕ(x−ξ)+(ε/τ)ϕ[(x−ξ)/τ]p(x)=(1-\varepsilon ) \phi (x-\xi )+(\varepsilon / \tau ) \phi [(x-\xi ) / \tau ] where ϕ\phi is the standard normal density, find I(ε,ξ,τ)I(\varepsilon , \xi , \tau ).

?
Problem 2.6.7

Verify the expressions (6.20) and (6.21).

?
Problem 2.6.8

Let A=(A11A12A21A22)A=\left(\begin{array}{ll}A_{11} & A_{12} \\ A_{21} & A_{22}\end{array}\right) be a partitioned matrix with A22A_{22} square and nonsingular, and let

B=(I−A12A22−10I). B=\left(\begin{array}{cc} I-A_{12} A_{22}^{-1} 0 & I \end{array}\right).

Show that ∣A∣=∣A11−A12A22−1A21∣⋅∣A22∣\left|A\right|=\left|A_{11}-A_{12} A_{22}^{-1} A_{21}\right| \cdot \left|A_{22}\right|.

?
Problem 2.6.9
?
(a)

Let

A=(ab′bC) A=\left(\begin{array}{ll} a & b^{\prime } \\ b & C \\ \end{array}\right)

where aa is a scalar and bb a column matrix, and suppose that AA is positive definite. Show that ∣A∣≤a∣C∣\left|A\right| \leq a\left|C\right| with equality holding if and only if b=0b=0.

(b)

More generally, if the matrix AA of Problem 6.8 is positive definite, show that ∣A∣≤∣A11∣⋅∣A22∣\left|A\right| \leq \left|A_{11}\right| \cdot \left|A_{22}\right| with equality holding if and only if A12=0A_{12}=0.

Problem 2.6.10
?
(a)

Show that if the matrix AA is nonsingular, then for any vector x,(x′Ax)(x′A−1x)>(x′x)2x,\left(x^{\prime } A x\right)\left(x^{\prime } A^{-1} x\right)> \left(x^{\prime } x\right)^{2}.

(b)

Show that, in the notation of Theorem 6.6 and the following discussion,

[∂∂θiEθ[δ]]2Iii(θ)=(εi′α)2εi′I(θ)εi, \frac{\left[\frac{\partial }{\partial \theta _{i}} \mathbb {E}_{\theta }\left[\delta \right]\right]^{2}}{I_{i i}(\theta )}=\frac{\left(\varepsilon _{i}^{\prime } \alpha \right)^{2}}{\varepsilon _{i}^{\prime } I(\theta ) \varepsilon _{i}},

and if α=(0,…,0,αi,0,…0),α′I(θ)−1α=(εi′α)2εi′I(θ)−1εi\alpha =\left(0, \ldots , 0, \alpha_{i}, 0, \ldots 0\right), \alpha^{\prime } I(\theta )^{-1} \alpha =\left(\varepsilon_{i}^{\prime } \alpha \right)^{2} \varepsilon_{i}^{\prime } I(\theta )^{-1} \varepsilon_{i}, and hence establish (6.25).

Problem 2.6.11

Prove that (6.26) is necessary for equality in (6.25).

?
Problem 2.6.12

Prove the Bhattacharyya inequality (6.29) and show that the condition of equality is as stated.

?