Chapter 5

Estimation in Nonparametric Models

[40/100%]
Show
LaTeX
Problem 5.1

Let p≥1p\geq 1 and Fp\mathcal F_p be the set of cumulative distribution functions on R\mathbb {R} having finite ppth moment. Mallows' distance between FF and GG in Fp\mathcal F_p is defined to be

ϱMp(F,G)=inf⁡(E[∥X−Y∥p])1/p, \varrho _{M_p}(F,G)=\inf \left(\mathbb {E}\left[\left\| X-Y\right\| ^p\right]\right)^{1/p},

where the infimum is taken over all pairs of random variables XX and YY having marginal distributions FF and GG, respectively. Show that ϱMp\varrho_{M_p} is a distance on Fp\mathcal F_p.

?
Remark.
?

Corresponding exercise in Shao (2003): 5.3.

Problem 5.2

Let F1\mathcal F_1 be the collection of cumulative distribution functions on R\mathbb {R} with finite means and ϱM1\varrho_{M_1} be as defined in Exercise 5.1. Show that

?
(i)

ϱM1(F,G)=∫01∣F−1(t)−G−1(t)∣ dt\varrho_{M_1}(F,G)=\int_0^1\left|F^{-1}(t)-G^{-1}(t)\right|\, dt;

(ii)

ϱM1(F,G)=∫−∞∞∣F(x)−G(x)∣ dx\varrho_{M_1}(F,G)=\int_{-\infty }^{\infty }\left|F(x)-G(x)\right|\, dx.

Remark.
?

Recall that Fp\mathcal F_p is the collection of cumulative distribution functions on R\mathbb {R} with finite ppth moment, and ϱMp(F,G)=inf⁡[E[∣X−Y∣p]]1/p\varrho_{M_p}(F,G)=\inf \left[\mathbb {E}\left[\left|X-Y\right|^p\right]\right]^{1/p}, where the infimum ranges over joint distributions with marginals FF and GG with p=1p=1.

Corresponding exercise in Shao (2003): 5.5.

Problem 5.3

Let ϱMp\varrho_{M_p} be the Mallows' distance defined in Exercise 5.1 and {G,G1,G2,…}⊂Fp\left\{ G,G_1,G_2,\ldots \right\} \subset \mathcal F_p. Show that lim⁡nϱMp(Gn,G)=0\lim_n\varrho_{M_p}(G_n,G)=0 if and only if

lim⁡n∫∣x∣p dGn(x)=∫∣x∣p dG(x) \lim _n\int \left|x\right|^p\, dG_n(x)=\int \left|x\right|^p\, dG(x)

and lim⁡nGn(x)=G(x)\lim_nG_n(x)=G(x) for any xx at which GG is continuous.

?
Remark.
?

Recall that Fp\mathcal F_p is the collection of cumulative distribution functions on R\mathbb {R} with finite ppth moment, and ϱMp(F,G)=inf⁡[E[∣X−Y∣p]]1/p\varrho_{M_p}(F,G)=\inf \left[\mathbb {E}\left[\left|X-Y\right|^p\right]\right]^{1/p}, where the infimum ranges over joint distributions with marginals FF and GG (p≥1p\geq 1).

Problem 5.4

Find an example of cumulative distribution functions G,G1,G2,…G,G_1,G_2,\ldots on R\mathbb {R} such that

?
(i)

lim⁡nϱ∞(Gn,G)=0\lim_n\varrho_\infty (G_n,G)=0 but ϱMp(Gn,G)\varrho_{M_p}(G_n,G) does not converge to 0, where ϱMp\varrho_{M_p} is the distance defined in Exercise 5.1 and ϱ∞\varrho_\infty is the sup-norm distance defined as ϱ∞(F,G)=sup⁡x∣F(x)−G(x)∣\varrho_\infty (F,G)=\sup_x\left|F(x)-G(x)\right| for any cumulative distribution functions FF and GG;

(ii)

lim⁡nϱMp(Gn,G)=0\lim_n\varrho_{M_p}(G_n,G)=0 but ϱ∞(Gn,G)\varrho_\infty (G_n,G) does not converge to 0.

Remark.
?

Recall that Fp\mathcal F_p is the collection of cumulative distribution functions on R\mathbb {R} with finite ppth moment, and ϱMp(F,G)=inf⁡[E[∣X−Y∣p]]1/p\varrho_{M_p}(F,G)=\inf \left[\mathbb {E}\left[\left|X-Y\right|^p\right]\right]^{1/p}, where the infimum ranges over joint distributions with marginals FF and GG (p≥1p\geq 1).

Corresponding exercise in Shao (2003): 5.6.

Problem 5.5

Let XX be a random variable having cumulative distribution function FF. Show that

?
(i)

E[∣X∣2]<∞\mathbb {E}\left[\left|X\right|^2\right]<\infty implies ∫{F(t)[1−F(t)]}p/2 dt<∞\int \left\{ F(t)[1-F(t)]\right\}^{p/2}\, dt<\infty for p>1p>1;

(ii)

E[∣X∣2+δ]<∞\mathbb {E}\left[\left|X\right|^{2+\delta }\right]<\infty with some δ>0\delta >0 implies ∫{F(t)[1−F(t)]}1/2 dt<∞\int \left\{ F(t)[1-F(t)]\right\}^{1/2}\, dt<\infty.

Remark.
?

Corresponding exercise in Shao (2003): 5.8.

Problem 5.6

Show that pi=c/np_i=c/n, i=1,…,ni=1,\ldots ,n, λ=−(c/n)n−1\lambda =-(c/n)^{n-1} is a maximum of the function

H(p1,…,pn,λ)=∏i=1npi+λ(∑i=1npi−c) H(p_1,\ldots ,p_n,\lambda )=\prod _{i=1}^n p_i+\lambda \left(\sum _{i=1}^n p_i-c\right)

over pi>0p_i>0, i=1,…,ni=1,\ldots ,n, ∑i=1npi=c\sum_{i=1}^n p_i=c.

?
Remark.
?

This exercise shows that the empirical distribution function (which puts mass n−1n^{-1} to each of nn observations) is a maximum likelihood estimator.

Corresponding exercise in Shao (2003): 5.10.

Problem 5.7

Consider the problem of maximizing ∏i=1npi\prod_{i=1}^n p_i over

pi>0,i=1,…,n,∑i=1npi=1,and∑i=1npiui=0, p_i>0,\quad i=1,\ldots ,n,\quad \sum _{i=1}^n p_i=1,\quad \text{and}\quad \sum _{i=1}^n p_i u_i=0,

where uiu_i's are ss-vectors. Show that the solution is

p^i=1n(1+λτui),i=1,…,n, \widehat{p}_i=\frac1{n(1+\lambda ^\tau u_i)},\quad i=1,\ldots ,n,

and λ∈Rs\lambda \in \mathbb {R}^{s} satisfying

∑i=1nui1+λτui=0. \sum _{i=1}^n\frac{u_i}{1+\lambda ^\tau u_i}=0.
?
Remark.
?

Corresponding exercise in Shao (2003): 5.11.

Problem 5.8

Let δ1,…,δn\delta_1,\ldots ,\delta_n be nn observations from a binary random variable and

ℓ(p1,…,pn+1)=∏i=1npiδi(∑j=i+1n+1pj)1−δi. \ell (p_1,\ldots ,p_{n+1})=\prod _{i=1}^n p_i^{\delta _i}\left(\sum _{j=i+1}^{n+1}p_j\right)^{1-\delta _i}.
?
(i)

Show that maximizing ℓ(p1,…,pn+1)\ell (p_1,\ldots ,p_{n+1}) subject to pi≥0p_i\geq 0, i=1,…,n+1i=1,\ldots ,n+1, ∑i=1n+1pi=1\sum_{i=1}^{n+1}p_i=1 is equivalent to maximizing

∏i=1nqiδi(1−qi)n−i+1−δi, \prod _{i=1}^n q_i^{\delta _i}(1-q_i)^{n-i+1-\delta _i},

where qi=pi/∑j=in+1pjq_i=p_i/\sum_{j=i}^{n+1}p_j, i=1,…,ni=1,\ldots ,n.

(ii)

Show that

p^i=δin−i+1∏j=1i−1(1−δjn−j+1),i=1,…,n,p^n+1=1−∑i=1np^i \widehat{p}_i=\frac{\delta _i}{n-i+1}\prod _{j=1}^{i-1}\left(1-\frac{\delta _j}{n-j+1}\right),\quad i=1,\ldots ,n,\quad \widehat{p}_{n+1}=1-\sum _{i=1}^n\widehat{p}_i

maximizes ℓ(p1,…,pn+1)\ell (p_1,\ldots ,p_{n+1}) subject to pi≥0p_i\geq 0, i=1,…,n+1i=1,\ldots ,n+1, ∑i=1n+1pi=1\sum_{i=1}^{n+1}p_i=1.

(iii)

For any x1≤x2≤⋯≤xnx_1\leq x_2\leq \cdots \leq x_n, show that

F^(t)=∑i=1n+1p^i  1(0,t](xi)=1−∏xi≤t(1−δin−i+1). \widehat{F}(t)=\sum _{i=1}^{n+1}\widehat{p}_i\; \mathbb {1}_{(0,t]}(x_i)=1-\prod _{x_i\leq t}\left(1-\frac{\delta _i}{n-i+1}\right).
(iv)

When δi=1\delta_i=1 for all ii, show that p^i=n−1\widehat{p}_i=n^{-1}, i=1,…,ni=1,\ldots ,n, and p^n+1=0\widehat{p}_{n+1}=0.

Remark.
?

This exercise shows that the well-known Kaplan--Meier product-limit estimator is a maximum likelihood estimator.

Corresponding exercise in Shao (2003): 5.15.

Problem 5.9

Let (X1,…,Xn)(X_1,\ldots ,X_n) be a random sample of random variables with Lebesgue density ff,

Fn(t)=1n∑i=1n  1(−∞,t](Xi),t∈R, F_n(t)=\frac1n\sum _{i=1}^n\; \mathbb {1}_{(-\infty ,t]}(X_i),\quad t\in \mathbb {R},

be the empirical distribution, and

fn(t)=Fn(t+λn)−Fn(t−λn)2λn,t∈R, f_n(t)=\frac{F_n(t+\lambda _n)-F_n(t-\lambda _n)}{2\lambda _n},\quad t\in \mathbb {R},

where {λn}\left\{ \lambda_n\right\} is a sequence of positive constants.

?
(i)

Show that fnf_n is a Lebesgue density on R\mathbb {R}.

(ii)

Suppose that ff is continuously differentiable at tt, lim⁡nλn=0\lim_n\lambda_n=0, and lim⁡nnλn=∞\lim_n n\lambda_n=\infty. Show that the mean squared error of fn(t)f_n(t) as an estimator of f(t)f(t) equals

f(t)2nλn+o(1nλn)+O(λn2) \frac{f(t)}{2n\lambda _n}+o\left(\frac1{n\lambda _n}\right)+O(\lambda _n^2)

as n→∞n\to \infty.

(iii)

Under lim⁡nnλn3=0\lim_n n\lambda_n^3=0 and the conditions in (ii), show that

nλn[fn(t)−f(t)]→dN(0,f(t)/2). \sqrt{n\lambda _n}[f_n(t)-f(t)]\to _d\mathscr {N}\left(0, f(t)/2\right).
(iv)

Suppose that ff is continuous on [a,b][a,b], −∞<a<b<∞-\infty <a<b<\infty, lim⁡nλn=0\lim_n\lambda_n=0, and lim⁡nnλn=∞\lim_n n\lambda_n=\infty. Show that

∫abfn(t) dt→p∫abf(t) dt. \int _a^b f_n(t)\, dt\to _p\int _a^b f(t)\, dt.
Remark.
?

Corresponding exercise in Shao (2003): 5.16.

Problem 5.10

Let (X1,…,Xn)(X_1,\ldots ,X_n) be a random sample of random variables with Lebesgue density ff, ww be a known Lebesgue density on R\mathbb {R}, and

f^(t)=1nλn∑i=1nw(t−Xiλn), \widehat{f}(t)=\frac1{n\lambda _n}\sum _{i=1}^n w\left(\frac{t-X_i}{\lambda _n}\right),

where {λn}\left\{ \lambda_n\right\} is a sequence of positive constants.

?
(i)

Show that f^\widehat{f} is a Lebesgue density on R\mathbb {R}.

(ii)

Show that if λn→0\lambda_n\to 0, nλn→∞n\lambda_n\to \infty, and ff is bounded and continuous at tt and w0=∫−∞∞[w(t)]2 dt<∞w_0=\int_{-\infty }^\infty [w(t)]^2\, dt<\infty, then

nλn{f^(t)−E[f^(t)]}→dN(0,w0f(t)). \sqrt{n\lambda _n}\left\{ \widehat{f}(t)-\mathbb {E}\left[\widehat{f}(t)\right]\right\} \to _d\mathscr {N}\left(0, w_0f(t)\right).
(iii)

Assume that lim⁡nnλn3=0\lim_n n\lambda_n^3=0, ∫−∞∞∣t∣w(t) dt<∞\int_{-\infty }^\infty \left|t\right|w(t)\, dt<\infty, f′f' is bounded and continuous at tt, and the conditions in (ii) hold. Show that

nλn[f^(t)−f(t)]→dN(0,w0f(t)). \sqrt{n\lambda _n}[\widehat{f}(t)-f(t)]\to _d\mathscr {N}\left(0, w_0f(t)\right).
(iv)

Suppose that λn→0\lambda_n\to 0, nλn→∞n\lambda_n\to \infty, ww is bounded, and ff is bounded and continuous on [a,b][a,b], −∞<a<b<∞-\infty <a<b<\infty. Show that ∫abf^(t) dt→p∫abf(t) dt\int_a^b\widehat{f}(t)\, dt\to_p\int_a^b f(t)\, dt.

Remark.
?

Corresponding exercise in Shao (2003): 5.17.

Problem 5.11

Let ℓ(θ,ξ)\ell (\theta ,\xi ) be a likelihood function. Show that a maximum profile likelihood estimator θ^\widehat{\theta } of θ\theta is an MLE if ξ(θ)\xi (\theta ), the maximum of sup⁡ξℓ(θ,ξ)\sup_\xi \ell (\theta ,\xi ) for a fixed θ\theta, does not depend on θ\theta.

?
Remark.
?

A maximum profile likelihood estimator θ^\widehat{\theta } maximizes the profile likelihood function ℓP(θ)=ℓ(θ,ξ(θ))\ell_P(\theta )=\ell (\theta ,\xi (\theta )), where ℓ(θ,ξ(θ))=sup⁡ξℓ(θ,ξ)\ell (\theta ,\xi (\theta ))=\sup_\xi \ell (\theta ,\xi ) for each fixed θ\theta.

Corresponding exercise in Shao (2003): 5.20.

Problem 5.12

Let (X1,…,Xn)(X_1,\ldots ,X_n) be a random sample from N(μ,σ2)\mathscr {N}\left(\mu , \sigma^2\right). Derive the profile likelihood function for μ\mu or σ2\sigma^2. Discuss in each case whether the maximum profile likelihood estimator is the same as the MLE.

?
Remark.
?

Corresponding exercise in Shao (2003): 5.21.

Problem 5.13

Let (X1,…,Xn)(X_1,\ldots ,X_n) be a random sample from a distribution FF and let π(x)=P(δi=1∣Xi=x)\pi (x)=\mathbb {P}\left( \delta_i=1 \mid X_i=x \right), where δi=1\delta_i=1 if XiX_i is observed and δi=0\delta_i=0 if XiX_i is missing. Assume that 0<π=∫π(x) dF(x)<10<\pi =\int \pi (x)\, dF(x)<1.

?
(i)

Let F1(x)=P(Xi≤x∣δi=1)F_1(x)=\mathbb {P}\left( X_i\leq x \mid \delta_i=1 \right). Show that FF and F1F_1 are the same if and only if π(x)≡π\pi (x)\equiv \pi.

(ii)

Let F^\widehat{F} be the empirical distribution putting mass r−1r^{-1} to each observed XiX_i, where rr is the number of observed XiX_i's. Show that F^(x)\widehat{F}(x) is unbiased and consistent for F1(x)F_1(x), x∈Rx\in \mathbb {R}.

(iii)

When π(x)≡π\pi (x)\equiv \pi, show that F^(x)\widehat{F}(x) in part (ii) is unbiased and consistent for F(x)F(x), x∈Rx\in \mathbb {R}. When π(x)\pi (x) is not constant, show that F^(x)\widehat{F}(x) is biased and inconsistent for F(x)F(x) for some x∈Rx\in \mathbb {R}.

Remark.
?

Corresponding exercise in Shao (2003): 5.23.

Problem 5.14

Let F\mathcal F be a collection of distributions on Rd\mathbb {R}^{d}. A functional TT defined on F\mathcal F is Gateaux differentiable at G∈FG\in \mathcal F if and only if there is a linear functional LGL_G on D={c(G1−G2):c∈R,Gj∈F,j=1,2}\mathcal D=\left\{ c(G_1-G_2):c\in \mathbb {R},G_j\in \mathcal F,j=1,2\right\} (i.e., LG(c1Δ1+c2Δ2)=c1LG(Δ1)+c2LG(Δ2)L_G(c_1\Delta_1+c_2\Delta_2)=c_1L_G(\Delta_1)+c_2L_G(\Delta_2) for any Δj∈D\Delta_j\in \mathcal D and cj∈Rc_j\in \mathbb {R}) such that Δ∈D\Delta \in \mathcal D and G+tΔ∈FG+t\Delta \in \mathcal F imply

lim⁡t→0[T(G+tΔ)−T(G)t−LG(Δ)]=0. \lim _{t\to 0}\left[\frac{T(G+t\Delta )-T(G)}t-L_G(\Delta )\right]=0.

Assume that the functional LFL_F is continuous in the sense that ∥Δj−Δ∥∞→0\left\| \Delta_j-\Delta \right\|_{\infty }\to 0 implies LF(Δj)→LF(Δ)L_F(\Delta_j)\to L_F(\Delta ), where D∈DD\in \mathcal D, Dj∈DD_j\in \mathcal D, and ∥D∥∞=sup⁡x∣D(x)∣\left\| D\right\|_{\infty }=\sup_x\left|D(x)\right| for any D∈DD\in \mathcal D is the sup-norm. Show that ϕF(x)=LF(δx−F)\phi_F(x)=L_F(\delta_x-F) is a bounded function of xx, where δx\delta_x is the degenerated distribution at xx.

?
Remark.
?

Corresponding exercise in Shao (2003): 5.25.

Problem 5.15

Suppose that a functional TT is Gateaux differentiable at FF with a bounded and continuous influence function ϕF(x)=LF(δx−F)\phi_F(x)=L_F(\delta_x-F), where δx\delta_x is the degenerated distribution at xx. Show that LFL_F is continuous in the sense described in the previous exercise.

?
Remark.
?

Here D={c(G1−G2):c∈R, G1,G2∈F}\mathcal D=\left\{ c(G_1-G_2):c\in \mathbb {R},\ G_1,G_2\in \mathcal F\right\} and LFL_F is the linear differential satisfying lim⁡t→0[T(F+tΔ)−T(F)]/t=LF(Δ)\lim_{t\to 0}[T(F+t\Delta )-T(F)]/t=L_F(\Delta ) for admissible Δ∈D\Delta \in \mathcal D. Continuity means LF(Dj)→LF(D)L_F(D_j)\to L_F(D) whenever Dj,D∈DD_j,D\in \mathcal D and sup⁡x∣Dj(x)−D(x)∣→0\sup_x\left|D_j(x)-D(x)\right|\to 0.

Corresponding exercise in Shao (2003): 5.26.

Problem 5.16

Let F\mathcal F be the collection of all distributions on R\mathbb {R} and zz be a fixed real number. Define

T(G)=∫G(z−y) dG(y),G∈F. T(G)=\int G(z-y)\, dG(y),\quad G\in \mathcal F.

Obtain the influence function ϕF\phi_F for TT and show that ϕF\phi_F is continuous if and only if FF is continuous.

?
Remark.
?

Corresponding exercise in Shao (2003): 5.29.

Problem 5.17

An L-functional is defined as

T(G)=∫xJ(G(x)) dG(x),G∈F0, T(G)=\int xJ(G(x))\, dG(x),\quad G\in \mathcal F_0,

where F0\mathcal F_0 contains all distributions on R\mathbb {R} for which TT is well defined and J(t)J(t) is a Borel function on [0,1][0,1].

?
(i)

Show that the influence function is

ϕF(x)=−∫−∞∞(δx−F)(y)J(F(y)) dy, \phi _F(x)=-\int _{-\infty }^\infty (\delta _x-F)(y)J(F(y))\, dy,

where δx\delta_x is the degenerated distribution at xx.

(ii)

Show that ∫ϕF(x) dF(x)=0\int \phi_F(x)\, dF(x)=0 and, if JJ is bounded and FF has a finite second moment, then ∫[ϕF(x)]2 dF(x)<∞\int [\phi_F(x)]^2\, dF(x)<\infty.

Remark.
?

Corresponding exercise in Shao (2003): 5.34.

Problem 5.18

Obtain explicit forms of the influence functions for L-functionals (Exercise 5.17) in the following cases and discuss which of them are bounded and continuous.

?
(i)

J≡1J\equiv 1.

(ii)

J(t)=4t−2J(t)=4t-2.

(iii)

J(t)=(β−α)−1  1(α,β)(t)J(t)=(\beta -\alpha )^{-1}\; \mathbb {1}_{(\alpha ,\beta )}(t) for some constants α<β\alpha <\beta.

Remark.
?

An L-functional is T(G)=∫xJ(G(x)) dG(x)T(G)=\int xJ(G(x))\, dG(x) on the distributions for which the integral is well defined, with JJ a Borel function on [0,1][0,1]. Its influence function is ϕF(x)=LF(δx−F)\phi_F(x)=L_F(\delta_x-F), where LFL_F is the Gateaux differential and δx\delta_x is the point mass at xx.

Corresponding exercise in Shao (2003): 5.37.

Problem 5.19

In part (iii) of the previous exercise, show that if FF is continuous at F−1(α)F^{-1}(\alpha ) and F−1(β)F^{-1}(\beta ), then

ϕF(x)={F−1(α)(1−α)−F−1(β)(1−β)β−α−T(F)x<F−1(α),]6pt]x−F−1(α)α−F−1(β)(1−β)β−α−T(F)F−1(α)≤x≤F−1(β),F−1(β)β−F−1(α)αβ−α−T(F)x>F−1(β). \phi _F(x)=\begin{cases} \displaystyle \frac{F^{-1}(\alpha )(1-\alpha )-F^{-1}(\beta )(1-\beta )}{\beta -\alpha }-T(F)& x<F^{-1}(\alpha ),\\ ]6pt] \displaystyle \frac{x-F^{-1}(\alpha )\alpha -F^{-1}(\beta )(1-\beta )}{\beta -\alpha }-T(F)& F^{-1}(\alpha )\leq x\leq F^{-1}(\beta ),\\[6pt] \displaystyle \frac{F^{-1}(\beta )\beta -F^{-1}(\alpha )\alpha }{\beta -\alpha }-T(F)& x>F^{-1}(\beta ). \end{cases}
?
Remark.
?

An L-functional is T(G)=∫xJ(G(x)) dG(x)T(G)=\int xJ(G(x))\, dG(x) on the distributions for which the integral is well defined, with JJ a Borel function on [0,1][0,1]. Its influence function is ϕF(x)=LF(δx−F)\phi_F(x)=L_F(\delta_x-F), where LFL_F is the Gateaux differential and δx\delta_x is the point mass at xx. Here J(t)=(β−α)−1  1(α,β)(t)J(t)=(\beta -\alpha )^{-1}\; \mathbb {1}_{(\alpha ,\beta )}(t), as in Exercise 5.18(iii).

Problem 5.20

Let TT be an L-functional defined in Exercise 5.17.

?
(i)

Show that T(F)=θT(F)=\theta if FF is symmetric about θ\theta, JJ is symmetric about 12\frac12, and ∫01J(t) dt=1\int_0^1J(t)\, dt=1.

(ii)

Assume that

σF2=∫−∞∞∫−∞∞J(F(x))J(F(y))[F(min⁡{x,y})−F(x)F(y)] dx dy \sigma _F^2=\int _{-\infty }^\infty \int _{-\infty }^\infty J(F(x))J(F(y))[F(\min \left\{ x,y\right\} )-F(x)F(y)]\, dx\, dy

is finite. Show that σ2=∫[ϕF(x)]2 dF(x)\sigma^2=\int [\phi_F(x)]^2\, dF(x), where ϕF\phi_F is the influence function of TT.

(iii)

Show that if J≡1J\equiv 1, then σF2\sigma_F^2 in (ii) is equal to the variance of FF.

Remark.
?

An L-functional is T(G)=∫xJ(G(x)) dG(x)T(G)=\int xJ(G(x))\, dG(x) on the distributions for which the integral is well defined, with JJ a Borel function on [0,1][0,1]. Its influence function is ϕF(x)=LF(δx−F)\phi_F(x)=L_F(\delta_x-F), where LFL_F is the Gateaux differential and δx\delta_x is the point mass at xx.

Corresponding exercises in Shao (2003): 5.67, 5.69, 5.74.

Problem 5.21

Let TT be an L-functional given in Exercise 5.18(iii) with β=1−α\beta =1-\alpha and α∈(0,12)\alpha \in (0,\frac12), and let FnF_n be the empirical distribution based on a random sample from a distribution FF.

?
(i)

Let X(1)≤⋯≤X(n)X_{(1)}\leq \cdots \leq X_{(n)} be the order statistics. Show that

T(Fn)=X‾α=1(1−2α)n∑j=mα+1n−mαX(j), T(F_n)=\overline{X}_\alpha =\frac1{(1-2\alpha )n}\sum _{j=m_\alpha +1}^{n-m_\alpha }X_{(j)},

which is called the α\alpha-trimmed sample mean. Here mαm_\alpha is the integer part of nαn\alpha.

(ii)

Assume that FF is continuous at F−1(α)F^{-1}(\alpha ) and F−1(1−α)F^{-1}(1-\alpha ) and is symmetric about θ\theta. Show that

σα2=2(1−2α)2{∫0F0−1(1−α)x2 dF0(x)+α[F0−1(1−α)]2}, \sigma _\alpha ^2=\frac2{(1-2\alpha )^2}\left\{ \int _0^{F_0^{-1}(1-\alpha )}x^2\, dF_0(x)+\alpha [F_0^{-1}(1-\alpha )]^2\right\} ,

where F0(x−θ)=F(x)F_0(x-\theta )=F(x), is equal to the σF2\sigma_F^2 in part (ii) of the previous exercise with J(t)=(1−2α)−1  1(α,1−α)(t)J(t)=(1-2\alpha )^{-1}\; \mathbb {1}_{(\alpha ,1-\alpha )}(t).

(iii)

Show that if F0′(0)F_0'(0) exists and is positive, then lim⁡α→12σα2=1/[2F0′(0)]2\lim_{\alpha \to \frac12}\sigma_\alpha^2=1/[2F_0'(0)]^2.

(iv)

Show that if σ2=∫x2 dF0(x)<∞\sigma^2=\int x^2\, dF_0(x)<\infty, then lim⁡α→0σα2=σ2\lim_{\alpha \to 0}\sigma_\alpha^2=\sigma^2.

Remark.
?

An L-functional is T(G)=∫xJ(G(x)) dG(x)T(G)=\int xJ(G(x))\, dG(x) on the distributions for which the integral is well defined, with JJ a Borel function on [0,1][0,1]. Its influence function is ϕF(x)=LF(δx−F)\phi_F(x)=L_F(\delta_x-F), where LFL_F is the Gateaux differential and δx\delta_x is the point mass at xx. Here J(t)=(1−2α)−1  1(α,1−α)(t)J(t)=(1-2\alpha )^{-1}\; \mathbb {1}_{(\alpha ,1-\alpha )}(t).

Corresponding exercises in Shao (2003): 5.65, 5.72, 5.73.

Problem 5.22

Calculate σF2\sigma_F^2 defined in Exercise 5.20(ii) with J(t)=4t−2J(t)=4t-2 and FF being the double exponential distribution with location parameter μ∈R\mu \in \mathbb {R} and scale parameter 1.

?
Remark.
?

The quantity in Exercise 5.20(ii) is σF2=∫−∞∞∫−∞∞J(F(x))J(F(y))[F(min⁡{x,y})−F(x)F(y)] dx dy\sigma_F^2=\int_{-\infty }^{\infty }\int_{-\infty }^{\infty }J(F(x))J(F(y))[F(\min \left\{ x,y\right\} )-F(x)F(y)]\, dx\, dy.

Corresponding exercise in Shao (2003): 5.75.

Problem 5.23

Let T(G)=G−1(p)T(G)=G^{-1}(p) be the ppth quantile functional. Suppose that FF has a positive derivative F′F' in a neighborhood of θ=F−1(p)\theta =F^{-1}(p). Show that TT is Gateaux differentiable at FF and obtain the influence function.

?
Remark.
?

Corresponding exercise in Shao (2003): 5.59.

Problem 5.24

Let FnF_n be the empirical distribution based on a random sample of size nn from a distribution FF on R\mathbb {R} having Lebesgue density ff. Let φn(t)\varphi_n(t) be the Lebesgue density of ppth sample quantile Fn−1(p)F_n^{-1}(p). Prove that

φn(t)=n(n−1lp−1)[F(t)]lp−1[1−F(t)]n−lpf(t), \varphi _n(t)=n\binom {n-1}{l_p-1}[F(t)]^{l_p-1}[1-F(t)]^{n-l_p}f(t),

where lp=npl_p=np if npnp is an integer and lp=1+l_p=1+ the integer part of npnp if npnp is not an integer, by

?
(i)

using the fact that nFn(t)nF_n(t) has a binomial distribution;

(ii)

using the fact that Fn−1(p)=cnpX(mp)+(1−cnp)X(mp+1)F_n^{-1}(p)=c_{np}X_{(m_p)}+(1-c_{np})X_{(m_p+1)}, where X(j)X_{(j)} is the jjth order statistic, mpm_p is the integer part of npnp, cnp=1c_{np}=1 if npnp is an integer, and cnp=0c_{np}=0 if npnp is not an integer.

Remark.
?

Corresponding exercise in Shao (2003): 5.51.

Problem 5.25

Let FnF_n be the empirical distribution based on a random sample from a distribution FF on R\mathbb {R} with a finite mean. Show that the ppth sample quantile Fn−1(p)F_n^{-1}(p) has a finite jjth moment for sufficiently large nn, j=1,2,…j=1,2,\ldots, where p∈(0,1)p\in (0,1).

?
Remark.
?

Corresponding exercise in Shao (2003): 5.52.

Problem 5.26

Let FnF_n be the empirical distribution based on a random sample from a distribution FF on R\mathbb {R} with Lebesgue density ff that is positive and continuous at θ=F−1(p)\theta =F^{-1}(p), p∈(0,1)p\in (0,1). Using Scheffe's theorem (e.g., Proposition 1.18 in Shao, 2003), prove that

n[Fn−1(p)−θ]→dN(0,p(1−p)[f(θ)]2). \sqrt n[F_n^{-1}(p)-\theta ]\to _d\mathscr {N}\left(0, \frac{p(1-p)}{[f(\theta )]^2}\right).
?
Remark.
?

Corresponding exercise in Shao (2003): 5.54.

Problem 5.27

Let {kn}\left\{ k_n\right\} be a sequence of integers satisfying kn/n=p+o(n−1/2)k_n/n=p+o(n^{-1/2}) with p∈(0,1)p\in (0,1), and let (X1,…,Xn)(X_1,\ldots ,X_n) be a random sample from a distribution FF on R\mathbb {R} with F′(θp)>0F'(\theta_p)>0, where θp=F−1(p)\theta_p=F^{-1}(p). Let X(j)X_{(j)} be the jjth order statistic. Show that

n(X(kn)−θp)→dN(0,p(1−p)/[F′(θp)]2). \sqrt n(X_{(k_n)}-\theta _p)\to _d\mathscr {N}\left(0, p(1-p)/[F'(\theta _p)]^2\right).
?
Remark.
?

Corresponding exercise in Shao (2003): 5.55.

Problem 5.28

Let G,G1,G2,…G,G_1,G_2,\ldots be cumulative distributions on R\mathbb {R}. Suppose that lim⁡nsup⁡x∣Gn(x)−G(x)∣=0\lim_n\sup_x\left|G_n(x)-G(x)\right|=0, GG is continuous, and G−1G^{-1} is continuous at p∈(0,1)p\in (0,1).

?
(i)

Show that lim⁡nGn−1(p)=G−1(p)\lim_nG_n^{-1}(p)=G^{-1}(p).

(ii)

Show that the result in (i) holds for any p∈(0,1)p\in (0,1) if G′(x)G'(x) exists and is positive for any x∈Rx\in \mathbb {R}.

Remark.
?

Corresponding exercise in Shao (2003): 5.112.

Problem 5.29

Calculate the asymptotic relative efficiency of the Hodges--Lehmann estimator with respect to the sample mean based on a random sample from FF when

?
(i)

FF is the cumulative distribution of N(μ,σ2)\mathscr {N}\left(\mu , \sigma^2\right);

(ii)

FF is the cumulative distribution of the logistic distribution with location parameter μ\mu and scale parameter σ\sigma;

(iii)

FF is the cumulative distribution of the double exponential distribution with location parameter μ\mu and scale parameter σ\sigma;

(iv)

F(x)=F0(x−μ)F(x)=F_0(x-\mu ), where F0(x)F_0(x) is the cumulative distribution function of the t-distribution tνt_\nu with ν≥3\nu \geq 3.

Remark.
?

The Hodges--Lehmann estimator is a median of the values (Xi+Xj)/2(X_i+X_j)/2, i,j=1,…,ni,j=1,\ldots ,n; if there are two middle values, any point between them may be used.

Corresponding exercise in Shao (2003): 5.47.

Problem 5.30

Consider a random sample from a distribution FF on R\mathbb {R}. In each of the following cases, obtain the asymptotic relative efficiency of the sample median with respect to the sample mean.

?
(i)

FF is the cumulative distribution function of the uniform distribution on the interval (θ−12,θ+12)(\theta -\frac12,\theta +\frac12), θ∈R\theta \in \mathbb {R}.

(ii)

F(x)=F0(x−θ)F(x)=F_0(x-\theta ) and F0F_0 is the cumulative distribution function with Lebesgue density (1+x2)−1  1(−c,c)(x)/∫−cc(1+x2)−1 dt(1+x^2)^{-1}\; \mathbb {1}_{(-c,c)}(x)/\int_{-c}^c(1+x^2)^{-1}\, dt.

(iii)

F(x)=(1−ϵ)Φ(x−μσ)+ϵD(x−μσ)F(x)=(1-\epsilon )\Phi \left(\frac{x-\mu }\sigma \right)+\epsilon D\left(\frac{x-\mu }\sigma \right), where ϵ∈(0,1)\epsilon \in (0,1) is a known constant, Φ\Phi is the cumulative distribution function of the standard normal distribution, DD is the cumulative distribution function of the double exponential distribution with location parameter 0 and scale parameter 1, and μ∈R\mu \in \mathbb {R} and σ>0\sigma >0 are unknown parameters.

Remark.
?

Corresponding exercises in Shao (2003): 5.61, 5.62, 5.63.

Problem 5.31

Let (X1,…,Xn)(X_1,\ldots ,X_n) be a random sample from a distribution on R\mathbb {R} with Lebesgue density 2−1(1−θ2)eθx−∣x∣2^{-1}(1-\theta^2)e^{\theta x-\left|x\right|}, where θ∈(−1,1)\theta \in (-1,1) is unknown.

?
(i)

Show that the median of the distribution of X1X_1 is given by m(θ)=(1−θ)−1log⁡(1+θ)m(\theta )=(1-\theta )^{-1}\log (1+\theta ) when θ≥0\theta \geq 0 and m(θ)=−m(−θ)m(\theta )=-m(-\theta ) when θ<0\theta <0.

(ii)

Show that the mean of the distribution of X1X_1 is μ(θ)=2θ/(1−θ2)\mu (\theta )=2\theta /(1-\theta^2).

(iii)

Show that the inverse functions of m(θ)m(\theta ) and μ(θ)\mu (\theta ) exist. Obtain the asymptotic relative efficiency of m−1(m^)m^{-1}(\widehat{m}) with respect to μ−1(X‾)\mu^{-1}(\overline{X}), where m^\widehat{m} is the sample median and X‾\overline{X} is the sample mean.

(iv)

Is μ−1(X‾)\mu^{-1}(\overline{X}) asymptotically efficient in estimating θ\theta?

Remark.
?

Corresponding exercise in Shao (2003): 5.64.

Problem 5.32

Obtain the asymptotic relative efficiency of the trimmed sample mean X‾α\overline{X}_\alpha (Exercise 5.21) with respect to

?
(i)

the sample mean, based on a random sample of size nn from the double exponential distribution with location parameter θ∈R\theta \in \mathbb {R} and scale parameter 1;

(ii)

the sample median, based on a random sample of size nn from the Cauchy distribution with location parameter θ∈R\theta \in \mathbb {R} and scale parameter 1.

Remark.
?

For 0<α<120<\alpha <\frac12, the trimmed sample mean is X‾α=[(1−2α)n]−1∑j=mα+1n−mαX(j)\overline{X}_\alpha =[(1-2\alpha )n]^{-1}\sum_{j=m_\alpha +1}^{n-m_\alpha }X_{(j)}, where mαm_\alpha is the integer part of nαn\alpha and X(1)≤⋯≤X(n)X_{(1)}\leq \cdots \leq X_{(n)} are the order statistics.

Corresponding exercises in Shao (2003): 5.70, 5.71.

Problem 5.33

Let (X1,…,Xn)(X_1,\ldots ,X_n) be a random sample from a distribution FF on R\mathbb {R} that is symmetric about θ∈R\theta \in \mathbb {R}. Huber's estimator of θ\theta is defined as a solution of ∑i=1nψ(Xi,t)=0\sum_{i=1}^n\psi (X_i,t)=0, where

ψ(x,t)={Ct−x>C,t−x∣x−t∣≤C,−Ct−x<−C \psi (x,t)=\begin{cases} C& t-x>C,\\ t-x& \left|x-t\right|\leq C,\\ -C& t-x<-C\end{cases}

and C>0C>0 is a constant. Assume that FF is continuous at θ−C\theta -C and θ+C\theta +C.

?
(i)

Show that the function

Ψ(γ)=∫γ−Cγ+C(γ−x) dF(x)+CF(γ−C)−C[1−F(γ+C)] \Psi (\gamma )=\int _{\gamma -C}^{\gamma +C}(\gamma -x)\, dF(x)+CF(\gamma -C)-C[1-F(\gamma +C)]

is differentiable at θ\theta and Ψ(θ)=0\Psi (\theta )=0.

(ii)

Show that the asymptotic relative efficiency of Huber's estimator with respect to the sample mean is Var⁡[X1]/σF2\operatorname {Var}\left[X_1\right]/\sigma_F^2, where

σF2=∫θ−Cθ+C(θ−x)2 dF(x)+C2F(θ−C)+C2[1−F(θ+C)][F(θ+C)−F(θ−C)]2. \sigma _F^2=\frac{\int _{\theta -C}^{\theta +C}(\theta -x)^2\, dF(x)+C^2F(\theta -C)+C^2[1-F(\theta +C)]}{[F(\theta +C)-F(\theta -C)]^2}.
Remark.
?

Corresponding exercise in Shao (2003): 5.85.

Problem 5.34

For Huber's estimator θ^\widehat{\theta } in the previous exercise, obtain a formula e(F)e(F) for the asymptotic relative efficiency of θ^\widehat{\theta } with respect to the sample mean, when

F(x)=(1−ϵ)Φ(x−θσ)+ϵΦ(x−θτσ), F(x)=(1-\epsilon )\Phi \left(\frac{x-\theta }\sigma \right)+\epsilon \Phi \left(\frac{x-\theta }{\tau \sigma }\right),

where Φ\Phi is the cumulative distribution function of N(0,1)\mathscr {N}\left(0, 1\right), σ>0\sigma >0, τ>0\tau >0, and 0<ϵ<10<\epsilon <1. Show that lim⁡τ→∞e(F)=∞\lim_{\tau \to \infty }e(F)=\infty. Find the value of e(F)e(F) when ϵ=0\epsilon =0, σ=1\sigma =1, and C=1.5C=1.5.

?
Remark.
?

Huber's estimator is a root t=θ^t=\widehat{\theta } of ∑i=1nψ(Xi,t)=0\sum_{i=1}^n\psi (X_i,t)=0, where ψ(x,t)=max⁡{−C,min⁡{C,t−x}}\psi (x,t)=\max \left\{ -C,\min \left\{ C,t-x\right\} \right\} and C>0C>0. The observations are a random sample from FF, symmetric about θ\theta, with FF continuous at θ−C\theta -C and θ+C\theta +C.

Corresponding exercise in Shao (2003): 5.86.

Problem 5.35

Consider the L-functional TT defined in Exercise 5.17. Let FnF_n be the empirical distribution based on a random sample of size nn from a distribution FF, σ2=∫[ϕF(x)]2 dF(x)\sigma^2=\int [\phi_F(x)]^2\, dF(x), and σFn2=∫[ϕFn(x)]2 dFn(x)\sigma_{F_n}^2=\int [\phi_{F_n}(x)]^2\, dF_n(x), where ϕG\phi_G denotes the influence function of TT at distribution GG. Show that lim⁡nσFn2=σF2\lim_n\sigma_{F_n}^2=\sigma_F^2 a.s., under one of the following two conditions:

  1. JJ is bounded, J(t)=0J(t)=0 when t∈[0,α]∪[β,1]t\in [0,\alpha ]\cup [\beta ,1] for some constants α<β\alpha <\beta, and the set D={x:J is discontinuous at F(x)}D=\left\{ x:J\text{ is discontinuous at }F(x)\right\} has Lebesgue measure 0.

  2. JJ is continuous on [0,1][0,1] and ∫x2 dF(x)<∞\int x^2\, dF(x)<\infty.

?
Remark.
?

An L-functional is T(G)=∫xJ(G(x)) dG(x)T(G)=\int xJ(G(x))\, dG(x) on the distributions for which the integral is well defined, with JJ a Borel function on [0,1][0,1]. Its influence function is ϕF(x)=LF(δx−F)\phi_F(x)=L_F(\delta_x-F), where LFL_F is the Gateaux differential and δx\delta_x is the point mass at xx.

Corresponding exercise in Shao (2003): 5.99.

Problem 5.36

Let (X1,…,Xn)(X_1,\ldots ,X_n) be a random sample from a distribution FF on R\mathbb {R} and let UnU_n be a U-statistic (see Exercise 3.25) with kernel h(x1,…,xm)h(x_1,\ldots ,x_m) satisfying E[[h(X1,…,Xm)]2]<∞\mathbb {E}\left[[h(X_1,\ldots ,X_m)]^2\right]<\infty, where m<nm<n. Assume that ζ1=Var⁡[h1(X1)]>0\zeta_1=\operatorname {Var}\left[h_1(X_1)\right]>0, where h1(x)=E[h(x,X2,…,Xm)]h_1(x)=\mathbb {E}\left[h(x,X_2,\ldots ,X_m)\right]. Derive a consistent variance estimator for UnU_n.

?
Remark.
?

For a symmetric kernel hh, the U-statistic in Exercise 3.25 is Un=(nm)−1∑1≤i1<⋯<im≤nh(Xi1,…,Xim)U_n=\binom nm^{-1}\sum_{1\leq i_1<\cdots <i_m\leq n}h(X_{i_1},\ldots ,X_{i_m}). That exercise gives Var⁡[Un]=m2ζ1/n+O(n−2)\operatorname {Var}\left[U_n\right]=m^2\zeta_1/n+O(n^{-2}) when ζ1>0\zeta_1>0.

Corresponding exercise in Shao (2003): 5.100.

Problem 5.37

For Huber's estimator defined in Exercise 5.33, derive a consistent estimator of its asymptotic variance σF2\sigma_F^2.

?
Remark.
?

Huber's estimator is a root t=θ^t=\widehat{\theta } of ∑i=1nψ(Xi,t)=0\sum_{i=1}^n\psi (X_i,t)=0, where ψ(x,t)=max⁡{−C,min⁡{C,t−x}}\psi (x,t)=\max \left\{ -C,\min \left\{ C,t-x\right\} \right\} and C>0C>0. The observations are a random sample from FF, symmetric about θ\theta, with FF continuous at θ−C\theta -C and θ+C\theta +C.

Corresponding exercise in Shao (2003): 5.101.

Problem 5.38

Let X1,…,XnX_1,\ldots ,X_n be random variables. For any estimator θ^\widehat{\theta }, its jackknife variance estimator is defined as

vJ=n−1n∑i=1n(θ^−i−1n∑j=1nθ^−j)2, v_J=\frac{n-1}n\sum _{i=1}^n\left(\widehat{\theta }_{-i}-\frac1n\sum _{j=1}^n\widehat{\theta }_{-j}\right)^2,

where θ^−i\widehat{\theta }_{-i} is the same as θ^\widehat{\theta } but is based on n−1n-1 observations X1,…,Xi−1,Xi+1,…,XnX_1,\ldots ,X_{i-1},X_{i+1},\ldots ,X_n, i=1,…,ni=1,\ldots ,n. Let X‾\overline{X} be the sample mean and θ^=X‾2\widehat{\theta }=\overline{X}^2. Show that

vJ=4X‾2c^2n−1−4X‾c^3(n−1)2+c^4−c^22(n−1)3, v_J=\frac{4\overline{X}^2\widehat{c}_2}{n-1}-\frac{4\overline{X}\widehat{c}_3}{(n-1)^2}+\frac{\widehat{c}_4-\widehat{c}_2^2}{(n-1)^3},

where c^k=n−1∑i=1n(Xi−X‾)k\widehat{c}_k=n^{-1}\sum_{i=1}^n(X_i-\overline{X})^k, k=2,3,4k=2,3,4.

?
Remark.
?

Corresponding exercise in Shao (2003): 5.104.

Problem 5.39

Let X1,…,XnX_1,\ldots ,X_n be random variables and X1∗,…,Xn∗X_1^*,\ldots ,X_n^* be a random sample (i.e., a simple random sample with replacement) from X1,…,XnX_1,\ldots ,X_n. For any estimator θ^\widehat{\theta }, its bootstrap variance estimator is vB=Var⁡∗[θ^∗]v_B=\operatorname {Var}_{*}\left[\widehat{\theta }^*\right], where θ^∗\widehat{\theta }^* is the same as θ^\widehat{\theta } but is based on X1∗,…,Xn∗X_1^*,\ldots ,X_n^* and Var⁡∗[ ⋅ ]\operatorname {Var}_{*}\left[\, \cdot \, \right] is the variance with respect to the distribution of X1∗,…,Xn∗X_1^*,\ldots ,X_n^*, given X1,…,XnX_1,\ldots ,X_n. Let X‾\overline{X} be the sample mean and θ^=X‾2\widehat{\theta }=\overline{X}^2. Show that

vB=4X‾2c^2n+4X‾c^3n2+c^4−c^22n3, v_B=\frac{4\overline{X}^2\widehat{c}_2}n+\frac{4\overline{X}\widehat{c}_3}{n^2}+\frac{\widehat{c}_4-\widehat{c}_2^2}{n^3},

where c^k=n−1∑i=1n(Xi−X‾)k\widehat{c}_k=n^{-1}\sum_{i=1}^n(X_i-\overline{X})^k, k=2,3,4k=2,3,4.

?
Remark.
?

Corresponding exercise in Shao (2003): 5.111.

Problem 5.40

Let X1,…,XnX_1,\ldots ,X_n be a random sample from a distribution on Rk\mathbb {R}^{k} with a finite Var⁡[X1]\operatorname {Var}\left[X_1\right]. Let X1∗,…,Xn∗X_1^*,\ldots ,X_n^* be a random sample from X1,…,XnX_1,\ldots ,X_n. Show that for almost all given sequences X1,X2,…X_1,X_2,\ldots, n(X‾∗−X‾)→dNk(0,Var⁡[X1])\sqrt n(\overline{X}^*-\overline{X})\to_dN_k(0,\operatorname {Var}\left[X_1\right]), where X‾\overline{X} is the sample mean based on X1,…,XnX_1,\ldots ,X_n and X‾∗\overline{X}^* is the sample mean based on X1∗,…,Xn∗X_1^*,\ldots ,X_n^*.

?
Remark.
?

Corresponding exercise in Shao (2003): 5.113.