5.2

Допустимость и минимаксность в экспоненциальных семействах

[25/0%]
Показать
LaTeX
Задача 5.2.1

Лемма 2.1 была расширена Berger (1990a) для случая, когда оцениваемая функция не обязана быть ограничена конечным интервалом, а вместо этого достигает максимума или минимума в конечной точке параметра.

Лемма 8.1. Пусть оцениваемая функция g(θ)g(\theta ) непостоянна и достигает глобального максимума или минимума в точке θ∗∈Ω\theta^{*} \in \Omega, для которой f(x∣θ∗)>0f\left(x \mid \theta^{*}\right)>0 п.в. (относительно доминирующей меры μ\mu), и пусть функция потерь L(θ,d)L(\theta , d) удовлетворяет предположениям Леммы 2.1. Тогда любая оценка δ\delta, принимающая значения выше максимума g(θ)g(\theta ) или ниже минимума, недопустима.

?
(a)

Покажите, что если θ∗\theta^{*} минимизирует g(θ)g(\theta ) и g^(x)\hat{g}(x) является несмещённой оценкой g(θ)g(\theta ), то существует ϵ>0\epsilon >0 такое, что множество Aϵ={x∈X:g^(x)<g(θ)−ϵ}A_{\epsilon }=\left\{ x \in \mathcal{X}: \hat{g}(x)<g(\theta )-\epsilon \right\} удовлетворяет P(Aϵ)>0P\left(A_{\epsilon }\right)>0. Аналогичный вывод верен, если g(θ∗)g\left(\theta^{*}\right) является максимумом.

(b)

Предположим, что g(θ∗)g\left(\theta^{*}\right) является минимумом. (Случай максимума рассматривается аналогично.) Покажите, что оценка

δ(x)={g^(x) если g^(x)≥g(θ∗)g(θ∗) если g^(x)<g(θ∗) \delta (x)= \begin{cases} \hat{g}(x) & \text{ если } \hat{g}(x) \geq g\left(\theta ^{*}\right) \\ g\left(\theta ^{*}\right) & \text{ если } \hat{g}(x)<g\left(\theta ^{*}\right)\end{cases}

удовлетворяет R(δ,θ)−R(g^(x),θ)<0R(\delta , \theta )-R(\hat{g}(x), \theta )<0.

(c)

Для ситуации из Примера 2.3 примените Лемму 8.1, чтобы установить недопустимость UMVU-оценки σA2\sigma_{A}^{2}. Также объясните, почему предположения Леммы 8.1 не выполняются для оценивания σ2\sigma^{2}.

Задача 5.2.2

Определите байесовский риск оценки (2.4), когда θ\theta имеет априорное распределение N(μ,τ2)N\left(\mu , \tau^{2}\right).

?
Задача 5.2.3

Докажите утверждение (d) во втором доказательстве Примера 2.8: что существует последовательность значений θi→−∞\theta_{i} \rightarrow -\infty, для которой b′(θi)→0b^{\prime }\left(\theta_{i}\right) \rightarrow 0.

?
Задача 5.2.4

Покажите, что оценка aX+b(0≤a≤1)a X+b(0 \leq a \leq 1) величины Eθ(X)E_{\theta }(X) недопустима (при квадратичной функции потерь) при каждом из следующих условий:

?
(a)

если Eθ(X)≥0E_{\theta }(X) \geq 0 для всех θ\theta и b<0b<0;

(b)

если Eθ(X)≤kE_{\theta }(X) \leq k для всех θ\theta и ak+b>ka k+b>k.

Задача 5.2.5

Покажите, что оценка [1/(1+λ)+ε]X[1 /(1+\lambda )+\varepsilon ] X величины Eθ(X)E_{\theta }(X) недопустима (при квадратичной функции потерь) при каждом из следующих условий:

?
(a)

если var⁡θ(X)/Eθ2(X)>λ>0\operatorname {var}_{\theta }(X) / E_{\theta }^{2}(X)>\lambda >0 и ε>0\varepsilon >0,

(b)

если var⁡θ(X)/Eθ2(X)<λ\operatorname {var}_{\theta }(X) / E_{\theta }^{2}(X)<\lambda и ε<0\varepsilon <0.

Задача 5.2.6

Покажите, что если var⁡θ(X)/Eθ2(X)>λ>0\operatorname {var}_{\theta }(X) / E_{\theta }^{2}(X)>\lambda >0, то оценка [1/(1+λ)+ε]X+b[1 /(1+\lambda )+\varepsilon ] X+b недопустима (при квадратичной функции потерь) при каждом из следующих условий:

?
(a)

если Eθ(X)>0E_{\theta }(X)>0 для всех θ,b>0\theta , b>0 и ε>0\varepsilon >0;

(b)

если Eθ(X)<0E_{\theta }(X)<0 для всех θ,b<0\theta , b<0 и ε>0\varepsilon >0 (Gupta 1966).

Задача 5.2.7

Brown (1986a) указывает на связь между информационным неравенством и несмещённой оценкой риска оценок стейновского типа.

?
(a)

Покажите, что из (2.7) следует

R(θ,δ)≥[1+b′(θ)]2n+b2(θ)≥1n+2b′(θ)n+b2(θ) R(\theta , \delta ) \geq \frac{\left[1+b^{\prime }(\theta )\right]^{2}}{n}+b^{2}(\theta ) \geq \frac{1}{n}+\frac{2 b^{\prime }(\theta )}{n}+b^{2}(\theta )

и, следовательно, если R(θ,δ)≤R(θ,Xˉ)R(\theta , \delta ) \leq R(\theta , \bar{X}), то 2b′(θ)n+b2≤0\frac{2 b^{\prime }(\theta )}{n}+b^{2} \leq 0.

(b)

Покажите, что нетривиальное решение b(θ)b(\theta ) привело бы к улучшенной оценке x−g(x)x-g(x) при p=1p=1 в Следствии 4.7.2.

Задача 5.2.8

Функция плотности f(x∣θ)f(x \mid \theta ) называется понижающей вариацию порядка n+1(VRn+1)n+1\left(V R_{n+1}\right), если для любой функции g(x)g(x) с k(k≤n)k(k \leq n) переменами знака (без учёта нулей) математическое ожидание Eθg(X)=∫g(x)f(x∣θ)dxE_{\theta } g(X)= \int g(x) f(x \mid \theta ) d x имеет не более kk перемен знака. Если Eθg(X)E_{\theta } g(X) имеет ровно kk перемен знака, они происходят в том же порядке. Покажите, что f(x∣θ)f(x \mid \theta ) обладает свойством VR2V R_{2} тогда и только тогда, когда она имеет монотонное отношение правдоподобия. (См. TSH2, Лемма 2, Раздел 3.3, для доказательства импликации «если».) Brown и др. (1981) дают подробное введение в эту тему, включая характеризации VRV R для многих семейств распределений (экспоненциальное семейство обладает свойством VR∞V R_{\infty }, как и χv2\chi_{v}^{2} с параметром ν\nu, и нецентральное χν2(λ)\chi_{\nu }^{2}(\lambda ) по параметру λ\lambda). Существует эквивалентность между VRnV R_{n} и TPnT P_{n} — полной положительностью порядка nn по Karlin (1968), а именно VRn=TPnV R_{n}=T P_{n}.

?
Задача 5.2.9

Для ситуации из Примера 2.9 покажите, что:

?
(a)

без потери общности ограничение θ∈[a,b]\theta \in [a, b] можно свести к θ∈[−m,m]\theta \in [-m, m], m>0m>0.

(b)

Если Λ\Lambda — априорное распределение, помещающее массу 1/2 в каждую из точек ±m\pm m, то байесовская оценка при квадратичной функции потерь есть

δΛ(xˉ)=memnxˉ−e−mnxˉemnxˉ+e−mnxˉ=mtanh⁡(mnxˉ). \delta ^{\Lambda }(\bar{x})=m \frac{e^{m n \bar{x}}-e^{-m n \bar{x}}}{e^{m n \bar{x}}+e^{-m n \bar{x}}}=m \tanh (m n \bar{x}).
(c)

Для m<1/nm<1 / \sqrt{n}

max⁡θ∈[−m,m]R(θ,δ(Xˉ))=max⁡{R(−m,δΛ(Xˉ)),R(m,δΛ(Xˉ))} \max _{\theta \in [-m, m]} R(\theta , \delta (\bar{X}))=\max \left\{ R\left(-m, \delta ^{\Lambda }(\bar{X})\right), R\left(m, \delta ^{\Lambda }(\bar{X})\right)\right\}

и, следовательно, по Следствию 1.6, δΛ\delta^{\Lambda } является минимаксной. [Указание: Задача 2.8 может быть использована, чтобы показать, что производная функции риска может иметь не более одной перемены знака, с отрицательного на положительный, и, следовательно, любой внутренний экстремум может быть только минимумом.]

(d)

Для m>1.05/nm>1.05 / \sqrt{n} оценка δΛ\delta^{\Lambda } из пункта (b) больше не является минимаксной. Объясните, почему это так, и предложите альтернативную оценку в этом случае. [Указание: Рассмотрите R(0,δΛ)R\left(0, \delta^{\Lambda }\right).]

Задача 5.2.10

Для ситуации из Примера 2.10 покажите, что:

?
(a)

max⁡θ∈[−m,m]R(θ,aXˉ+b)=max⁡{R(−m,aXˉ+b),R(m,aXˉ+b)}\max_{\theta \in [-m, m]} R(\theta , a \bar{X}+b)=\max \left\{ R(-m, a \bar{X}+b), R(m, a \bar{X}+b)\right\}.

(b)

Оценка a∗Xˉa^{*} \bar{X} с a∗=m2/(1n+m2)a^{*}=m^{2} /\left(\frac{1}{n}+m^{2}\right) является линейной минимаксной оценкой для всех mm с минимаксным риском a∗/na^{*} / n.

(c)

Xˉ\bar{X} является линейной минимаксной оценкой для m=∞m=\infty.

Задача 5.2.11

Предположим, что XX имеет распределение FξF_{\xi }, а YY имеет распределение GηG_{\eta }, где ξ\xi и η\eta изменяются независимо. Если известно, что η=η0\eta =\eta_{0}, то любая оценка δ(X,Y)\delta (X, Y) может быть улучшена с помощью

δ∗(x)=EYδ(x,Y)=∫δ(x,y)dGη0(y) \delta ^{*}(x)=E_{Y} \delta (x, Y)=\int \delta (x, y) d G_{\eta _{0}}(y)
?
Задача 5.2.12

В Примере 2.13 докажите, что оценка aY+ba Y+b недопустима при a>1/(r+1)a>1 /(r+1). [Указание: Задачи 2.4--2.6]

?
Задача 5.2.13

Пусть X1,…,XnX_{1}, \ldots , X_{n} независимые одинаково распределённые случайные величины с плотностью N(0,σ2)N\left(0, \sigma^{2}\right), и пусть S2=∑Xi2S^{2}=\sum X_{i}^{2}. Нас интересует оценивание σ2\sigma^{2} при квадратичной функции потерь с помощью линейных оценок cS2+dc S^{2}+d, где cc и dd — константы. Покажите, что:

?
(a)

допустимость оценки aY+ba Y+b в Примере 2.13 эквивалентна допустимости cS2+dc S^{2}+d при подходящем выборе cc и dd.

(b)

риск cS2+dc S^{2}+d задаётся формулой R(cS2+d,σ2)=2nc2σ2+[(nc−1)σ2+d]2R\left(c S^{2}+d, \sigma^{2}\right)=2 n c^{2} \sigma^{2}+\left[(n c-1) \sigma^{2}+d\right]^{2}

(c)

при d=0,R(cS2,σ2)<R(0,σ2)d=0, R\left(c S^{2}, \sigma^{2}\right)<R\left(0, \sigma^{2}\right), когда c<2/(n+2)c<2 /(n+2), и, следовательно, оценка aY+ba Y+b в Примере 2.13 недопустима при a=b=0a=b=0.

Примечание.
?

Это упражнение иллюстрирует тот факт, что константы не обязательно являются допустимыми оценками.

Задача 5.2.14

Для ситуации из Примера 2.15 пусть Z=Xˉ/SZ=\bar{X} / S.

?
(a)

Покажите, что риск оценки δ=φ(z)s2\delta =\varphi (z) s^{2} при квадратичной функции потерь минимизируется при выборе

φ(z)=φμ,σ∗(z)=E(S2/σ2∣z)/E((S2/σ2)2∣z). \varphi (z)=\varphi _{\mu , \sigma }^{*}(z)=E\left(S^{2} / \sigma ^{2} \mid z\right) / E\left(\left(S^{2} / \sigma ^{2}\right)^{2} \mid z\right).
(b)

Stein (1964) показал, что φμ,σ∗(z)≤φ0,1∗(z)\varphi_{\mu , \sigma }^{*}(z) \leq \varphi_{0,1}^{*}(z) для всех μ,σ\mu , \sigma. Предполагая это, выведите, что φs(Z)S2\varphi_{s}(Z) S^{2} доминирует над [1/(n+1)]S2[1 /(n+1)] S^{2} при квадратичной функции потерь, где

φs(z)=min⁡{φ0,1∗(z),1n+1}. \varphi _{s}(z)=\min \left\{ \varphi _{0,1}^{*}(z), \frac{1}{n+1}\right\} .
(c)

Покажите, что φ0,1∗(z)=(1+z2)/(n+2)\varphi_{0,1}^{*}(z)=\left(1+z^{2}\right) /(n+2), и, следовательно, φs(Z)S2\varphi_{s}(Z) S^{2} задаётся формулой (2.31).

(d)

Наилучшая эквивариантная оценка вида φ(Z)S2\varphi (Z) S^{2} была получена Brewster и Zidek (1974) и задаётся формулой

φBZ(z)=E(S2∣Z≤z)E(S4∣Z≤z), \varphi _{B Z}(z)=\frac{E\left(S^{2} \mid Z \leq z\right)}{E\left(S^{4} \mid Z \leq z\right)},

где математическое ожидание вычисляется в предположении μ=0\mu =0 и σ=1\sigma =1. Покажите, что φBZ(Z)S2\varphi_{B Z}(Z) S^{2} является обобщённой байесовской оценкой относительно априорного распределения

π(μ,σ)=1σ∫0∞u−1/2(1+u)−1e−unμ2/σ2dudμdσ \pi (\mu , \sigma )=\frac{1}{\sigma } \int _{0}^{\infty } u^{-1 / 2}(1+u)^{-1} e^{-u n \mu ^{2} / \sigma ^{2}} d u d \mu d \sigma

[Brewster и Zidek первоначально получили свою оценку не как байесовскую, а, напротив, сначала нашли оценку, а затем нашли априорное распределение. Brown (1968) рассмотрел семейство оценок, аналогичных оценкам Stein (1964), которые принимали различные значения в зависимости от порогового значения для z2z^{2}. Brewster и Zidek (1974) показали, что число пороговых значений может быть сколь угодно большим. Они построили последовательность оценок с убывающими рисками и всё более плотными пороговыми значениями, предел которой являлся наилучшей эквивалентной оценкой.]

Задача 5.2.15

Покажите эквивалентность следующих соотношений:

?
(a)

(2.26) и (2.27),

(b)

(2.34) и (2.35) при c=(n−1)/(n+1)c=\sqrt{(n-1) /(n+1)}, и

(c)

(2.38) и (2.39).

Задача 5.2.16

В Примере 2.17 покажите, что оценка aX/n+ba X / n+b недопустима для всех (a,b)(a, b) вне треугольника (2.39).

?
Задача 5.2.17

Докажите допустимость оценок, соответствующих внутренности треугольника (2.39), применив Теорему 2.4 и используя результаты Примера 4.1.5.

?
Задача 5.2.18

Используя Теорему 2.14, приведите альтернативное доказательство допустимости оценки aXˉ+ba \bar{X}+b, удовлетворяющей (2.6), в Примере 2.5.

?
Задача 5.2.19

Определите, какие оценки aX+ba X+b допустимы для оценивания E(X)E(X) при квадратичной функции потерь в следующих ситуациях:

?
(a)

XX имеет распределение Пуассона.

(b)

XX имеет отрицательное биномиальное распределение (Gupta 1966).

Задача 5.2.20

Пусть XX имеет распределение Poisson⁡(λ)\operatorname {Poisson}(\lambda ), и рассмотрим оценивание λ\lambda при функции потерь (d−λ)2/λ(d-\lambda )^{2} / \lambda с ограничением 0≤λ≤m0 \leq \lambda \leq m, где mm известно.

?
(a)

Используя рассуждение, аналогичное рассуждению из Примера 2.9, покажите, что XX не является минимаксной оценкой и что наименее благоприятное априорное распределение должно иметь множество w∧w_{\wedge } [из (1.5)], состоящее из конечного числа точек.

(b)

Пусть Λa\Lambda_{a} — априорное распределение, помещающее массу ai,i=1,…,ka_{i}, i=1, \ldots , k, в точки параметра bi,i=1,…,kb_{i}, i=1, \ldots , k. Покажите, что байесовская оценка, соответствующая этому априорному распределению, есть

δΛa(x)=1E(λ−1∣x)=∑i=1kaibixe−bi∑i=1kaibix−1e−bi. \delta ^{\Lambda a}(x)=\frac{1}{E\left(\lambda ^{-1} \mid x\right)}=\frac{\sum _{i=1}^{k} a_{i} b_{i}^{x} e^{-b_{i}}}{\sum _{i=1}^{k} a_{i} b_{i}^{x-1} e^{-b_{i}}}.
(c)

Пусть m0m_{0} — решение уравнения m=e−m(m0≈.57)m=e^{-m}\left(m_{0} \approx .57\right). Покажите, что при 0≤λ≤m0 \leq \lambda \leq m, m≤m0m \leq m_{0}, одноточечное априорное распределение (ai=1,b1=m)(a_{i}=1, b_{1}=m) даёт минимаксную оценку. Вычислите минимаксный риск и сравните его с риском XX.

(d)

Пусть m1m_{1} — первый положительный нуль уравнения (1+δΛ(m))2=2+m2/2\left(1+\delta^{\Lambda }(m)\right)^{2}=2+m^{2} / 2, где Λ\Lambda — двухточечное априорное распределение (a1=a,b1=0;a2=1−a,b2=m)(a_{1}=a, b_{1}=0 ; a_{2}=1-a, b_{2}=m). Покажите, что при 0≤λ≤m0 \leq \lambda \leq m, m0<m≤m1m_{0}<m \leq m_{1}, двухточечное априорное распределение даёт минимаксную оценку (используйте Следствие 1.6). Вычислите минимаксный риск и сравните его с риском XX. [С ростом mm ситуация усложняется, и точные минимаксные решения становятся труднодостижимыми. В этих случаях линейные приближения могут быть вполне удовлетворительными. См. Johnstone и MacGibbon 1992, 1993.]

Задача 5.2.21

Покажите, что условия (2.41) и (2.42) из Примера 2.22 являются не только достаточными, но и необходимыми для допустимости (2.40).

?
Задача 5.2.22

Пусть XX и YY независимо распределены по закону Пуассона с E(X)=ξE(X)=\xi и E(Y)=ηE(Y)=\eta соответственно. Покажите, что aX+bY+ca X+b Y+c допустима для оценивания ξ\xi при квадратичной функции потерь тогда и только тогда, когда либо 0≤a<1,b≥0,c≥00 \leq a<1, b \geq 0, c \geq 0, либо a=1,b=c=0a=1, b=c=0 (Makani 1972).

?
Задача 5.2.23

Пусть XX имеет плотность 12β(θ)eθxe−∣x∣,∣θ∣<1\frac{1}{2} \beta (\theta ) e^{\theta x} e^{-\left|x\right|},\left|\theta \right|<1.

?
(a)

Покажите, что β(θ)=1−θ2\beta (\theta )=1-\theta^{2}.

(b)

Покажите, что aX+ba X+b допустима для оценивания Eθ(X)E_{\theta }(X) при квадратичной функции потерь тогда и только тогда, когда 0≤a≤1/20 \leq a \leq 1 / 2.

Задача 5.2.24

Пусть XX распределена как N(θ,1)N(\theta , 1), и пусть θ\theta имеет несобственную априорную плотность π(θ)=eθ\pi (\theta )=e^{\theta } (−∞<θ<∞)(-\infty <\theta <\infty ). При квадратичной функции потерь формальная байесовская оценка θ\theta равна X+1X+1, и она не является ни минимаксной, ни допустимой. (См. также Задачу 2.15.) Условия, при которых формальная байесовская оценка, соответствующая несобственному априорному распределению для θ\theta в Примере 3.4, является допустимой, даны в Zidek (1970).

?
Задача 5.2.25

Покажите, что естественное параметрическое пространство семейства (2.16) есть (−∞,∞)(-\infty , \infty ) для нормального (с известной дисперсией), биномиального распределений и распределения Пуассона, но не в случае гамма-распределения или отрицательного биномиального распределения.

?