Section 1.6

Sufficient Statistics

[38/0%]
Show
LaTeX
Problem 1.6.1

Extend Example 6.2 to the case that X1,…,XrX_{1}, \ldots , X_{r} are independently distributed with Poisson distributions P(λi)P\left(\lambda_{i}\right) where λi=aiλ(ai>0, known)\lambda_{i}=a_{i} \lambda \left(a_{i}>0 \text{, known}\right).

?
Problem 1.6.2

Let X1,…,XnX_{1}, \ldots , X_{n} be iid according to a distribution FF and probability density ff. Show that the conditional distribution given X(i)=aX_{(i)}=a of the i−1i-1 values to the left of aa and the n−in-i values to the right of aa is that of i−1i-1 variables distributed independently according to the probability density f(x)/F(a)f(x) / F(a) and n−in-i variables distributed independently with density f(x)/[1−F(a)]f(x) /[1-F(a)], respectively, with the two sets being (conditionally) independent of each other.

?
Problem 1.6.3

Let ff be a positive integrable function over (0,∞)(0, \infty ), and let pθ(x)p_{\theta }(x) be the density over (0,θ)(0, \theta ) defined by pθ(x)=c(θ)f(x)p_{\theta }(x)=c(\theta ) f(x) if 0<x<θ0<x<\theta, and 0 otherwise. If X1,…,XnX_{1}, \ldots , X_{n} are iid with density pθp_{\theta }, show that X(n)X_{(n)} is sufficient for θ\theta.

?
Problem 1.6.4

Let ff be a positive integrable function defined over (−∞,∞)(-\infty , \infty ) and let pξ,η(x)p_{\xi , \eta }(x) be the probability density defined by pξ,η(x)=c(ξ,η)f(x)p_{\xi , \eta }(x)=c(\xi , \eta ) f(x) if ξ<x<η\xi <x<\eta, and 0 otherwise. If X1,…,XnX_{1}, \ldots , X_{n} are iid with density pξ,ηp_{\xi , \eta }, show that (X(1),X(n))\left(X_{(1)}, X_{(n)}\right) is sufficient for (ξ,n)(\xi , n).

?
Problem 1.6.5

Show that each of the statistics T1−T4T_{1}-T_{4} of Example 6.11 is sufficient.

?
Problem 1.6.6

Theorem 6.5 (Factorization Criterion) states: a necessary and sufficient condition for a statistic TT to be sufficient for a family P={Pθ,θ∈Ω}\mathcal{P}=\left\{ P_{\theta }, \theta \in \Omega \right\} of distributions of XX dominated by a σ\sigma-finite measure μ\mu is that there exist non-negative functions gθg_{\theta } and hh such that the densities pθp_{\theta } of PθP_{\theta } satisfy pθ(x)=gθ[T(x)]h(x)p_{\theta }(x)=g_{\theta }[T(x)] h(x) (a.e. μ\mu).

Corollary 6.13 states: under the assumptions of Theorem 6.5, a necessary and sufficient condition for a statistic UU to be sufficient is that for any fixed θ\theta and θ0\theta_{0}, the ratio pθ(x)/pθ0(x)p_{\theta }(x) / p_{\theta_{0}}(x) is a function only of U(x)U(x).

Prove Corollary 6.13.

?
Problem 1.6.7

Let X1,…,XmX_{1}, \ldots , X_{m} and Y1,…,YnY_{1}, \ldots , Y_{n} be independently distributed according to N(ξ,σ2)N\left(\xi , \sigma^{2}\right) and N(η,τ2)N\left(\eta , \tau^{2}\right), respectively. Find the minimal sufficient statistics for these cases:

?
(a)

ξ,η,σ,τ\xi , \eta , \sigma , \tau are arbitrary: −∞<ξ,η<∞,0<σ,τ-\infty <\xi , \eta <\infty , 0<\sigma , \tau.

(b)

σ=τ\sigma =\tau and ξ,η,σ\xi , \eta , \sigma are arbitrary.

(c)

ξ=η\xi =\eta and ξ,σ,τ\xi , \sigma , \tau are arbitrary.

Problem 1.6.8

Let X1,…,XnX_{1}, \ldots , X_{n} be iid according to N(σ,σ2),0<σN\left(\sigma , \sigma^{2}\right), 0<\sigma. Find a minimal set of sufficient statistics.

?
Problem 1.6.9
?
(a)

If (x1,…,xn)\left(x_{1}, \ldots , x_{n}\right) and (y1,…,yn)\left(y_{1}, \ldots , y_{n}\right), have the same elementary symmetric functions Σxi=Σyi,Σi≠jxiyj=Σi≠jyiyj,…,x1⋯xn=y1⋯yn\Sigma x_{i}=\Sigma y_{i}, \Sigma_{i \neq j} x_{i} y_{j}=\Sigma_{i \neq j} y_{i} y_{j}, \ldots , x_{1} \cdots x_{n}=y_{1} \cdots y_{n}, then the yy 's are a permutation of the xx 's.

(b)

In the notation of Example 6.10, show that UU is equivalent to VV. [Hint: Compare the coefficients and the roots of the polynomials P(x)=Π(x−ui)P(x)=\Pi \left(x-u_{i}\right) and Q(x)=Π(x−vi)Q(x)= \Pi \left(x-v_{i}\right).]

Problem 1.6.10

Show that the order statistics are minimal sufficient for the location family (6.7) when ff is the density of

?
(a)

the double exponential distribution D(0,1)D(0,1).

(b)

the Cauchy distribution C(0,1)C(0,1).

Problem 1.6.11

Prove the following generalization of Theorem 6.12 to families without common support.

Theorem 9.1 Let P\mathcal{P} be a finite family with densities pi,i=0,…,kp_{i}, i=0, \ldots , k, and for any xx, let S(x)S(x) be the set of pairs of subscripts (i,j)(i, j) for which pi(x)+pj(x)>0p_{i}(x)+p_{j}(x)>0. Then, the statistic

T(X)={pj(X)pi(X),i<j and (i,j)∈S(X)} T(X)=\left\{ \frac{p_{j}(X)}{p_{i}(X)}, \quad i<j \text{ and }(i, j) \in S(X)\right\}

is minimal sufficient. Here, pj(x)/pi(x)=∞p_{j}(x) / p_{i}(x)=\infty if pi(x)=0p_{i}(x)=0 and pj(x)>0p_{j}(x)>0.

?
Problem 1.6.12

In Problem 6.11 it is not enough to replace pi(X)p_{i}(X) by p0(X)p_{0}(X). To see this let k=2k=2 and p0=U(−1,0),p1=U(0,1)p_{0}=U(-1,0), p_{1}=U(0,1), and p2(x)=2x,0<x<1p_{2}(x)=2 x, 0<x<1.

?
Problem 1.6.13

Let k=1k=1 and Pi=U(i,i+1),i=0,1P_{i}=U(i, i+1), i=0,1.

?
(a)

Show that a minimal sufficient statistic for P={P0,P1}\mathcal{P}=\left\{ P_{0}, P_{1}\right\} is T(X)=iT(X)=i if i<X<i+1i<X<i+1, i=0,1i=0,1.

(b)

Let X1X_{1} and X2X_{2} be iid according to a distribution from P\mathcal{P}. Show that each of the two statistics T1=T(X1)T_{1}=T\left(X_{1}\right) and T2=T(X2)T_{2}=T\left(X_{2}\right) is sufficient for (X1,X2)\left(X_{1}, X_{2}\right).

(c)

Show that T(X1)T\left(X_{1}\right) and T(X2)T\left(X_{2}\right) are equivalent.

Problem 1.6.14

In Lemma 6.14, show that the assumption of common support can be replaced by the weaker assumption that every P0\mathcal{P}_{0}-null set is also a P\mathcal{P}-null set so that (a.e. P0\mathcal{P}_{0} ) is equivalent to (a.e. P\mathcal{P} ).

?
Problem 1.6.15

Let X1,…,XnX_{1}, \ldots , X_{n} be iid according to a distribution from P={U(0,θ),θ>0}\mathcal{P}=\left\{ U(0, \theta ), \theta >0\right\}, and let P0\mathcal{P}_{0} be the subfamily of P\mathcal{P} for which θ\theta is rational. Show that every P0\mathcal{P}_{0}-null set in the sample space is also a P\mathcal{P}-null set.

?
Problem 1.6.16

Let X1,…,XnX_{1}, \ldots , X_{n} be iid according to a distribution from a family P\mathcal{P}. Show that TT is minimal sufficient in the following cases:

?
(a)

P={U(0,θ),θ>0};T=X(n)\mathcal{P}=\left\{ U(0, \theta ), \theta >0\right\} ; T=X_{(n)}.

(b)

P={U(θ1,θ2),−∞<θ1<θ2<∞};T=(X(1),X(n))\mathcal{P}=\left\{ U\left(\theta_{1}, \theta_{2}\right),-\infty <\theta_{1}<\theta_{2}<\infty \right\} ; T=\left(X_{(1)}, X_{(n)}\right).

(c)

P={U(θ−1/2,θ+1/2),−∞<θ<∞};T=(X(1),X(n))\mathcal{P}=\left\{ U(\theta -1 / 2, \theta +1 / 2),-\infty <\theta <\infty \right\} ; T=\left(X_{(1)}, X_{(n)}\right).

Problem 1.6.17

Solve the preceding problem for the following cases:

?
(a)

P={E(θ,1),−∞<θ<∞};T=X(1)\mathcal{P}=\left\{ E(\theta , 1),-\infty <\theta <\infty \right\} ; T=X_{(1)}.

(b)

P={E(0,b),0<b};T=ΣXi\mathcal{P}=\left\{ E(0, b), 0<b\right\} ; T=\Sigma X_{i}.

(c)

P={E(a,b),−∞<a<∞,0<b};T=(X(1),Σ[Xi−X(1)])\mathcal{P}=\left\{ E(a, b),-\infty <a<\infty , 0<b\right\} ; T=\left(X_{(1)}, \Sigma \left[X_{i}-X_{(1)}\right]\right).

Problem 1.6.18

Show that the statistics X(1)X_{(1)} and Σ[Xi−X(1)]\Sigma \left[X_{i}-X_{(1)}\right] of Problem 6.17(c) are independently distributed as E(a,b/n)E(a, b / n) and bb Gamma (n−2,1)(n-2,1) respectively.

?
Problem 1.6.19

Show that the sufficient statistics of

  1. Problem 6.3 and

  2. Problem 6.4

are minimal sufficient.

?
Problem 1.6.20
?
(a)

Show that in the N(θ,θ)N(\theta , \theta ) curved exponential family, the sufficient statistic T=(∑xi,∑xi2)T= (\sum x_{i}, \sum x_{i}^{2}) is not minimal.

(b)

For the density of Example 6.19, show that T=(∑xi,∑xi2,∑xi3)T=\left(\sum x_{i}, \sum x_{i}^{2}, \sum x_{i}^{3}\right) is a minimal sufficient statistic.

Problem 1.6.21

For the situation of Example 6.25(ii), find an unbiased estimator of ξ\xi based on ∑Xi\sum X_{i}, and another based on ∑Xi2\sum X_{i}^{2} ); hence, deduce that T=(∑Xi,∑Xi2)T=\left(\sum X_{i}, \sum X_{i}^{2}\right) is not complete.

?
Problem 1.6.22

For the situation of Example 6.26, show that X\mathbf{X} is minimal sufficient and complete.

?
Problem 1.6.23

For the situation of Example 6.27:

?
(a)

Show that X=(X1,X2)\mathbf{X}=\left(X_{1}, X_{2}\right) is minimal sufficient for the family (6.16) with restriction (6.17).

(b)

Establish (6.18), and hence that the minimal sufficient statistic of part (a) is not complete.

Problem 1.6.24

Show that for the general dose-response model

pθ(x)=∏i=1m(nixi)[ηθ(di)]xi[1−ηθ(di)]ni−xi, p_{\theta }(\mathbf{x})=\prod _{i=1}^{m}\binom {n_{i}}{x_{i}}\left[\eta _{\theta }\left(d_{i}\right)\right]^{x_{i}}\left[1-\eta _{\theta }\left(d_{i}\right)\right]^{n_{i}-x_{i}},

the statistic X=(X1,X2,…,Xm)\mathbf{X}=\left(X_{1}, X_{2}, \ldots , X_{m}\right) is minimal sufficient if there exist vectors θ1,θ2,⋯ ,θm)\left.\theta_{1}, \theta_{2}, \cdots , \theta_{m}\right) such that the m×mm \times m matrix

P={log⁡(ηθj(di)[1−ηθ0(di)]ηθ0(di)[1−ηθj(di)])} P=\left\{ \log \left(\frac{\eta _{\theta _{j}}\left(d_{i}\right)\left[1-\eta _{\theta _{0}}\left(d_{i}\right)\right]}{\eta _{\theta _{0}}\left(d_{i}\right)\left[1-\eta _{\theta _{j}}\left(d_{i}\right)\right]}\right)\right\}

is invertible. (Hint: Theorem 6.12.)

?
Problem 1.6.25

Let (Xi,Yi),i=1,…,n\left(X_{i}, Y_{i}\right), i=1, \ldots , n, be iid according to the uniform distribution over a set RR in the (x,y)(x, y) plane and let P\mathcal{P} be the family of distributions obtained by letting RR range over a class R\mathcal{R} of sets RR. Determine a minimal sufficient statistic for the following cases:

?
(a)

R\mathcal{R} is the set of all rectangles a1<x<a2,b1<y<b2,−∞<a1<a2<∞a_{1}<x<a_{2}, b_{1}<y<b_{2},-\infty <a_{1}<a_{2}<\infty, −∞<b1<b2<∞-\infty <b_{1}<b_{2}<\infty.

(b)

R′\mathcal{R}^{\prime } is the subset of R\mathcal{R}, for which a2−a1=b2−b1a_{2}-a_{1}=b_{2}-b_{1}.

(c)

R′′\mathcal{R}^{\prime \prime } is the subset of R′\mathcal{R}^{\prime } for which a2−a1=b2−b1=1a_{2}-a_{1}=b_{2}-b_{1}=1.

Problem 1.6.26

Solve the preceding problem if

?
(a)

R\mathcal{R} is the set of all triangles with sides parallel to the xx axis, the yy axis, and the line y=xy=x, respectively.

(b)

R′\mathcal{R}^{\prime } is the subset of R\mathcal{R} in which the sides parallel to the xx and yy axes are equal.

Problem 1.6.27

Formulate a general result of which Problems 6.25(a) and 6.26(a) are special cases.

?
Problem 1.6.28

If YY is distributed as E(η,1)E(\eta , 1), the distribution of X=e−YX=e^{-Y} is U(0,e−η)U\left(0, e^{-\eta }\right). (This result is useful in the computer generation of random variables; see Problem 4.4.14.)

?
Problem 1.6.29

If a minimal sufficient statistic exists, a necessary condition for a sufficient statistic to be complete is for it to be minimal.

?
Problem 1.6.30

Show that the minimal sufficient statistics T=(X(1),X(n))T=\left(X_{(1)}, X_{(n)}\right) of Problem 6.16(b) are complete.

?
Problem 1.6.31

For each of the following problems, determine whether the minimal sufficient statistic is complete:

?
(a)

Problem 6.7(a)-(c);

(b)

Problem 6.25(a)-(c);

(c)

Problem 6.26(a) and (b).

Problem 1.6.32
?
(a)

Show that if P0,P1\mathcal{P}_{0}, \mathcal{P}_{1} are two families of distributions such that P0∈P1\mathcal{P}_{0} \in \mathcal{P}_{1} and every null set of P0\mathcal{P}_{0} is also a null set of P1\mathcal{P}_{1}, then a sufficient statistic TT that is complete for P0\mathcal{P}_{0} is also complete for P1\mathcal{P}_{1}.

(b)

Let P0\mathcal{P}_{0} be the class of binomial distributions b(p,n),0<p<1,n=b(p, n), 0<p<1, n= fixed, and let P1=P0∪{Q}\mathcal{P}_{1}=\mathcal{P}_{0} \cup \left\{ Q\right\} where QQ is the Poisson distribution with expectation 1. Then P0\mathcal{P}_{0} is complete but P1\mathcal{P}_{1} is not.

Problem 1.6.33

Let X1,…,XnX_{1}, \ldots , X_{n} be iid each with density f(x)f(x) (with respect to Lebesgue measure), which is unknown. Show that the order statistics are complete.

?
Problem 1.6.34

Suppose that X1,…,XnX_{1}, \ldots , X_{n} are an iid sample from a location-scale family with distribution function F((x−a)/b)F((x-a) / b).

?
(a)

If bb is known, show that the differences (X1−Xi)/b,i=2,…,n\left(X_{1}-X_{i}\right) / b, i=2, \ldots , n, are ancillary.

(b)

If aa is known, show that the ratios (X1−a)/(Xi−a),i=2,…,n\left(X_{1}-a\right) /\left(X_{i}-a\right), i=2, \ldots , n, are ancillary.

(c)

If neither aa or bb are known, show that the quantities (X1−Xi)/(X2−Xi),i=3,…,n\left(X_{1}-X_{i}\right) /\left(X_{2}-X_{i}\right), i= 3, \ldots , n, are ancillary.

Problem 1.6.35

Use Basu's theorem to prove independence of the following pairs of statistics:

?
(a)

Xˉ\bar{X} and Σ(Xi−Xˉ)2\Sigma \left(X_{i}-\bar{X}\right)^{2} where the XX 's are iid as N(ξ,σ2)N\left(\xi , \sigma^{2}\right).

(b)

X(1)X_{(1)} and Σ[Xi−X(1)]\Sigma \left[X_{i}-X_{(1)}\right] in Problem 6.18.

Problem 1.6.36
?
(a)

Under the assumptions of Problem 6.18, the ratios Zi=[X(n)−X(i)]/[X(n)−X(n−1)],i=1,…,n−2Z_{i}=\left[X_{(n)}-X_{(i)}\right] / \left[X_{(n)}-X_{(n-1)}\right], i=1, \ldots , n-2, are independent of {X(1),Σ[Xi−X(1)]}\left\{ X_{(1)}, \Sigma \left[X_{i}-X_{(1)}\right]\right\}.

(b)

Under the assumptions of Problems 6.16(b) and 6.30 the ratios Zi=[X(i)−X(1)]/[X(n)−X(1)],i=2,…,n−1Z_{i}=\left[X_{(i)}-X_{(1)}\right] / \left[X_{(n)}-X_{(1)}\right], i=2, \ldots , n-1, are independent of {X(1),X(n)}\left\{ X_{(1)}, X_{(n)}\right\}.

Problem 1.6.37

Under the assumptions of Theorem 6.5, let AA be any fixed set in the sample space, Pθ∗P_{\theta }^{*} the distribution PθP_{\theta } truncated on AA, and P∗={Pθ∗,θ∈Ω}\mathcal{P}^{*}=\left\{ P_{\theta }^{*}, \theta \in \Omega \right\}. Then prove

?
(a)

if TT is sufficient for P\mathcal{P}, it is sufficient for P∗\mathcal{P}^{*}.

(b)

if, in addition, TT is complete for P\mathcal{P}, it is also complete for P∗\mathcal{P}^{*}.

Generalizations of this result were derived by Tukey in the 1940s and also by Smith (1957). The analogous problem for observations that are censored rather than truncated is discussed by Bhattacharyya, Johnson, and Mehrotra (1977).

Problem 1.6.38

If X1,…,XnX_{1}, \ldots , X_{n} are iid as B(a,b)B(a, b),

?
(a)

Show that [ΠXi,Π(1−Xi)]\left[\Pi X_{i}, \Pi \left(1-X_{i}\right)\right] is minimal sufficient for (a,b)(a, b).

(b)

Determine the minimal sufficient statistic when a=ba=b.