Section 3.2

The Principle of Equivariance

[28/0%]
Show
LaTeX
Problem 3.2.1

Show that the class G(C)G(\mathcal{C}) is a group.

?
Problem 3.2.2

In Example 2.2(ii), show that the transformations x′=−x\mathbf{x}^{\prime }=-\mathbf{x} together with the identity transformation form a group.

?
Problem 3.2.3

Let {gX,g∈G}\left\{ g X, g \in G\right\} be a group of transformations that leave the model (2.1) invariant. If the distributions Pθ,θ∈ΩP_{\theta }, \theta \in \Omega are distinct, show that the induced transformations gˉ\bar{g} are 1:11: 1 transformations of Ω\Omega.

?
Problem 3.2.4

Under the assumptions of Problem 2.3, show that

?
(a)

the transformations gˉ\bar{g} satisfy g2g1‾=gˉ2⋅gˉ1\overline{g_{2} g_{1}}=\bar{g}_{2} \cdot \bar{g}_{1} and (gˉ)−1=(g−1‾)(\bar{g})^{-1}=\left(\overline{g^{-1}}\right);

(b)

the transformations gˉ\bar{g} corresponding to g∈Gg \in G form a group.

(c)

establish (2.3) and (2.4).

Problem 3.2.5

Show that a loss function satisfies (2.9) if and only if it is of the form (2.10).

?
Problem 3.2.6
?
(a)

The transformations g∗g^{*} defined by (2.12) satisfy (g2g1)∗=g2∗⋅g1∗\left(g_{2} g_{1}\right)^{*}=g_{2}^{*} \cdot g_{1}^{*} and (g∗)−1=(g−1)∗\left(g^{*}\right)^{-1}= \left(g^{-1}\right)^{*}.

(b)

If GG is a group leaving (2.1) invariant and G∗={g∗,g∈G}G^{*}=\left\{ g^{*}, g \in G\right\}, then G∗G^{*} is a group.

Problem 3.2.7

Let XX be distributed as N(ξ,σ2),−∞<ξ<∞,0<σN\left(\xi , \sigma^{2}\right),-\infty <\xi <\infty , 0<\sigma, and let h(ξ,σ)=σ2h(\xi , \sigma )=\sigma^{2}. The problem is invariant under the transformations x′=ax+c;0<a,−∞<c<∞x^{\prime }=a x+c ; 0<a,-\infty <c<\infty. Show that the only equivariant estimator is δ(X)≡0\delta (X) \equiv 0.

?
Problem 3.2.8

Show that:

?
(a)

If (2.11) holds, the transformations g∗g^{*} defined by (2.12) are 1 : 1 from H\mathcal{H} onto itself.

(b)

If L(θ,d)=L(θ,d′)L(\theta , d)=L\left(\theta , d^{\prime }\right) for all θ\theta implies d=d′d=d^{\prime }, then g∗g^{*} defined by (2.14) is unique, and is a 1: 1 transformation from D\mathcal{D} onto itself.

Problem 3.2.9

If θ\theta is the true temperature in degrees Celsius, then θ′=gˉθ=θ+273\theta^{\prime }=\bar{g} \theta =\theta +273 is the true temperature in degrees Kelvin. Given an observation XX, in degrees Celsius:

?
(a)

Show that an estimator δ(X)\delta (X) is functionally equivariant if it satisfies δ(x)+a=δ(x+a)\delta (x)+a= \delta (x+a) for all aa.

(b)

Suppose our estimator is δ(x)=(ax+bθ0)/(a+b)\delta (x)=\left(a x+b \theta_{0}\right) /(a+b), where xx is the observed temperature in degrees Celsius, θ0\theta_{0} is a prior guess at the temperature, and aa and bb are constants. Show that for a constant K,δ(x+K)≠δ(x)+KK, \delta (x+K) \neq \delta (x)+K, so δ\delta does not satisfy the principle of functional equivariance.

(c)

Show that the estimators of part (b) will not satisfy the principle of formal invariance.

Problem 3.2.10

To illustrate the difference between functional equivariance and formal invariance, consider the following. To estimate the amount of electric power obtainable from a stream, one could use the estimate

δ(x)=cmin⁡{100,x−20} \delta (x)=c \min \left\{ 100, x-20\right\}

where x=x= stream flow in m3/sec,100 m3/sec\mathrm{m}^{3} / \mathrm{sec}, 100 \mathrm{~ m}^{3} / \mathrm{sec} is the capacity of the pipe leading to the turbine, and 20 m3/sec20 \mathrm{~ m}^{3} / \mathrm{sec} is the flow reduction necessary to avoid harming the trout. The constant cc, in kilowatts /m3/sec/ \mathrm{m}^{3} / \mathrm{sec} converts the flow to a kilowatt estimate.

?
(a)

If measurements were, instead, made in liters and watts, so g(x)=1000xg(x)=1000 x and gˉ(θ)=1000θ\bar{g}(\theta )=1000 \theta, show that functional equivariance leads to the estimate

gˉ(δ(x))=cmin⁡{105,g(x)−20,000}. \bar{g}(\delta (x))=c \min \left\{ 10^{5}, g(x)-20,000\right\} .
(b)

The principle of formal invariance leads to the estimate δ(g(x))\delta (g(x)). Show that this estimator is not a reasonable estimate of wattage.

(Communicated by L. LeCam.)

Problem 3.2.11

In an invariant probability model, write X=(T,W)X=(T, W), where TT is sufficient for θ\theta, and WW is ancillary.

?
(a)

If the group operation is transitive, show that any invariant statistic must be ancillary.

(b)

What can you say about the invariance of an ancillary statistic?

Problem 3.2.12

In an invariant estimation problem, write X=(T,W)X=(T, W) where TT is sufficient for θ\theta, and WW is ancillary. If the group of transformations is transitive, show:

?
(a)

The best equivariant estimator δ∗\delta^{*} is the solution to min⁡dEθ[L(θ,d(x))∣W=w]\min_{d} \mathbb {E}_{\theta }\left[L(\theta , d(x)) \mid W=w\right].

(b)

If ee is the identity element of the group (g−1g=e)\left(g^{-1} g=e\right), then δ∗=δ∗(t,w)\delta^{*}=\delta^{*}(t, w) can be found by solving, for each w,min⁡dEe[L[e,d(T,w)]∣W=w]w, \min_{d} \mathbb {E}_{e}\left[L[e, d(T, w)] \mid W=w\right].

Problem 3.2.13

For the situation of Example 2.11:

?
(a)

Show that the class of transformations is a group.

(b)

Show that equivariant estimators must satisfy δ(n−x)=1−δ(x)\delta (n-x)=1-\delta (x).

(c)

Show that, using an invariant loss, the risk of an equivariant estimator is symmetric about p=1/2p=1 / 2.

Problem 3.2.14

For the situation of Example 2.12:

?
(a)

Show that the class of transformations is a group.

(b)

Show that estimators of the form φ(xˉ/s2)s2\varphi \left(\bar{x} / s^{2}\right) s^{2}, where xˉ=1/nΣxi\bar{x}=1 / n \Sigma x_{i} and s2=Σ(xi−xˉ)2s^{2}=\Sigma \left(x_{i}-\bar{x}\right)^{2} are equivariant, where φ\varphi is an arbitrary function.

(c)

Show that, using an invariant loss function, the risk of an equivariant estimator is a function only of τ=μ/σ\tau =\mu / \sigma.

Problem 3.2.15

Theorem 2.7 states: if δ\delta is an equivariant estimator in a problem which is invariant under a transformation group GG, then the risk function of δ\delta satisfies R(gˉθ,δ)=R(θ,δ)R(\bar{g}\theta ,\delta )=R(\theta ,\delta ) for all θ\theta and all g∈Gg\in G.

Corollary 2.13 states: under the assumptions of Theorem 2.7, the risk function of any equivariant estimator is constant on the orbits of GG.

Prove Corollary 2.13.

?
Problem 3.2.16
?
(a)

If gg is the transformation (2.20), determine gˉ\bar{g}.

(b)

In Example 2.12, show that (2.22) is not only sufficient for (2.14) but also necessary.

Problem 3.2.17
?
(a)

In Example 2.12, determine the smallest group GG containing both G1G_{1} and G2G_{2}.

(b)

Show that the only estimator that is invariant under GG is δ(X,Y)≡0\delta (\mathbf{X}, \mathbf{Y}) \equiv 0.

Problem 3.2.18

If δ(X)\delta (X) is an equivariant estimator of h(θ)h(\theta ) under a group GG, then so is g∗δ(X)g^{*} \delta (X) with g∗g^{*} defined by (2.12) and (2.13), provided G∗G^{*} is commutative.

?
Problem 3.2.19

Show that:

?
(a)

In Example 2.14(i), XX is not risk-unbiased.

(b)

The group of transformations ax+ca x+c of the real line (0<a,−∞<c<∞)(0<a,-\infty <c<\infty ) is not commutative.

Problem 3.2.20

In Example 2.14, determine the totality of equivariant estimators of Δ\Delta under the smallest group GG containing G1G_{1} and G2G_{2}.

?
Problem 3.2.21

Let θ\theta be real-valued and hh strictly increasing, so that (2.11) is vacuously satisfied. If L(θ,d)L(\theta , d) is the loss resulting from estimating θ\theta by dd, suppose that the loss resulting from estimating θ′=h(θ)\theta^{\prime }=h(\theta ) by d′=h(d)d^{\prime }=h(d) is M(θ′,d′)=L[θ,h−1(d′)]M\left(\theta^{\prime }, d^{\prime }\right)=L\left[\theta , h^{-1}\left(d^{\prime }\right)\right]. Show that:

?
(a)

If the problem of estimating θ\theta with loss function LL is invariant under GG, then so is the problem of estimating h(θ)h(\theta ) with loss function MM.

(b)

If δ\delta is equivariant under GG for estimating θ\theta with loss function LL, show that h[δ(X)]h[\delta (X)] is equivariant for estimating h(θ)h(\theta ) with loss function MM.

(c)

If δ\delta is MRE for θ\theta with LL, then h[δ(X)]h[\delta (X)] is MRE for h(θ)h(\theta ) with MM.

Problem 3.2.22

If δ(X)\delta (\mathbf{X}) is MRE for estimating ξ\xi in Example 2.2(i) with loss function ρ(d−ξ)\rho (d-\xi ), state an optimum property of eδ(X)e^{\delta (\mathbf{X})} as an estimator of eξe^{\xi }.

?
Problem 3.2.23

Let Xij,j=1,…,ni,i=1,…,sX_{i j}, j=1, \ldots , n_{i}, i=1, \ldots , s, and WW be distributed according to a density of the form

[∏i=1sfi(xi−ξi)]h(w) \left[\prod _{i=1}^{s} f_{i}\left(\mathbf{x}_{i}-\xi _{i}\right)\right] h(w)

where xi−ξi=(xi1−ξi,…,xini−ξi)\mathbf{x}_{i}-\xi_{i}=\left(x_{i 1}-\xi_{i}, \ldots , x_{i n_{i}}-\xi_{i}\right), and consider the problem of estimating θ=Σciξi\theta =\Sigma c_{i} \xi_{i} with loss function L(ξi,…,ξs;d)=ρ(d−θ)L\left(\xi_{i}, \ldots , \xi_{s} ; d\right)=\rho (d-\theta ). Show that:

?
(a)

This problem remains invariant under the transformations

Xij′=Xij+ai,ξi′=ξi+ai,θ′=θ+Σaici,d′=d+Σaici. \begin{align} X_{i j}^{\prime }=X_{i j}+a_{i}, \quad & \xi _{i}^{\prime }=\xi _{i}+a_{i}, \quad \theta ^{\prime }=\theta +\Sigma a_{i} c_{i}, \\ & d^{\prime }=d+\Sigma a_{i} c_{i}. \end{align}
(b)

An estimator δ\delta of θ\theta is equivariant under these transformations if

δ(x1+a1,…,xs+as,w)=δ(x1,…,xs,w)+Σaici. \delta \left(\mathbf{x}_{1}+a_{1}, \ldots , \mathbf{x}_{s}+a_{s}, w\right)=\delta \left(\mathbf{x}_{1}, \ldots , \mathbf{x}_{s}, w\right)+\Sigma a_{i} c_{i}.
Problem 3.2.24

Theorem 1.4 states: let XX be distributed with density (1.9), and let δ\delta be equivariant for estimating ξ\xi with loss function (1.15). Then, the bias, risk, and variance of δ\delta are all constant (i.e., do not depend on ξ\xi).

Generalize Theorem 1.4 to the situation of Problem 2.23.

?
Problem 3.2.25

If δ0\delta_{0} is any equivariant estimator of θ\theta in Problem 2.23, and if yi=(xi1−xini,xi2−xini,…,xini−1−xini)\mathbf{y}_{i}=\left(x_{i 1}-x_{i n_{i}}, x_{i 2}-x_{i n_{i}}, \ldots , x_{i n_{i}-1}-x_{i n_{i}}\right), show that the most general equivariant estimator of θ\theta is of the form

δ(x1,…,xs,w)=δ0(x1,…,xs,w)−v(y1,…,ys,w). \delta \left(\mathbf{x}_{1}, \ldots , \mathbf{x}_{s}, w\right)=\delta _{0}\left(\mathbf{x}_{1}, \ldots , \mathbf{x}_{s}, w\right)-v\left(\mathbf{y}_{1}, \ldots , \mathbf{y}_{s}, w\right).
?
Problem 3.2.26

Theorem 1.10 states: let X=(X1,…,Xn)X=(X_{1}, \ldots , X_{n}) be distributed according to (1.9), let Yi=Xi−XnY_{i}=X_{i}-X_{n} (i=1,…,n−1)(i=1, \ldots , n-1) and Y=(Y1,…,Yn−1)Y=(Y_{1}, \ldots , Y_{n-1}). Suppose that the loss function is given by (1.15) and that there exists an equivariant estimator δ0\delta_{0} of ξ\xi with finite risk. Assume that for each yy there exists a number v(y)=v∗(y)v(y)=v^{*}(y) which minimizes E0[ρ[δ0(X)−v(y)]∣y]\mathbb {E}_{0}\left[\rho [\delta_{0}(X)-v(y)] \mid y\right]. Then, a location equivariant estimator δ\delta of ξ\xi with minimum risk exists and is given by δ∗(X)=δ0(X)−v∗(Y)\delta^{*}(X)=\delta_{0}(X)-v^{*}(Y).

Corollary 1.12 states: under the assumptions of Theorem 1.10: (i) if ρ(d−ξ)=(d−ξ)2\rho (d-\xi )=(d-\xi )^{2}, then v∗(y)=E0[δ0(X)∣y]v^{*}(y)=\mathbb {E}_{0}\left[\delta_{0}(X) \mid y\right]; (ii) if ρ(d−ξ)=∣d−ξ∣\rho (d-\xi )=\left|d-\xi \right|, then v∗(y)v^{*}(y) is any median of δ0(X)\delta_{0}(X) under the conditional distribution of XX given yy.

?
(a)

Generalize Theorem 1.10 and Corollary 1.12 to the situation of Problems 2.23 and 2.25.

(b)

Show that the MRE estimators of (a) can be chosen to be independent of ww.

Problem 3.2.27

Suppose that the variables XijX_{i j} in Problem 2.23 are independently distributed as N(ξi,σ2),σN\left(\xi_{i}, \sigma^{2}\right), \sigma is known. Show that:

?
(a)

The MRE estimator of θ\theta is then ΣciXˉi−v∗\Sigma c_{i} \bar{X}_{i}-v^{*}, where Xˉi=(Xi1+⋯+Xini)/ni\bar{X}_{i}=\left(X_{i 1}+\cdots +X_{i n_{i}}\right) / n_{i}, and where v∗v^{*} minimizes (1.24) with X=ΣciXˉiX=\Sigma c_{i} \bar{X}_{i}.

(b)

If ρ\rho is convex and even, the MRE estimator of θ\theta is ΣciXˉi\Sigma c_{i} \bar{X}_{i}.

(c)

The results of (a) and (b) remain valid when σ\sigma is unknown and the distribution of WW depends on σ\sigma (but not the ξ\xi 's).

Problem 3.2.28

Show that the transformation of Example 2.11 and the identity transformation are the only transformations leaving the family of binomial distributions invariant.

?