2.1

Равномерно минимальные по дисперсии несмещённые оценки (UMVU)

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

Проверьте

?
(a)

что (1.4) определяет вероятностное распределение, и

(b)

условие (1.5).

Задача 2.1.2

В Примере 1.5 покажите, что ai∗a_{i}^{*} минимизирует (1.6) при i=0,1i=0,1, и упростите выражение для a0∗a_{0}^{*}.

?
Задача 2.1.3

Пусть XX принимает значения -1, 0, 1, 2, 3 с вероятностями P(X=−1)=2pq\mathbb {P}\left(X=-1\right)=2 p q и P(X=k)=pkq3−k\mathbb {P}\left(X=k\right)=p^{k} q^{3-k} при k=0,1,2,3k=0,1,2,3.

?
(a)

Проверьте, что это вероятностное распределение.

(b)

Найдите LMVU-оценку в точке p0p_{0} для (i) pp и (ii) pqp q, и для каждой определите, является ли она UMVU.

Задача 2.1.4

Для выборки объёма nn предположим, что оценка T(x)T(\mathbf{x}) параметра τ(θ)\tau (\theta ) имеет математическое ожидание

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

где aka_{k} может зависеть от θ\theta, но не от nn.

?
(a)

Покажите, что математическое ожидание джекнайф-оценки TJT_{J} из (1.3) равно

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)

Покажите, что если var T∼c/nT \sim c / n для некоторой константы cc, то var TJ∼c′/nT_{J} \sim c^{\prime } / n для некоторой константы c′c^{\prime }. Таким образом, джекнайф уменьшает смещение, не увеличивая дисперсию.

Джекнайф второго порядка можно определить, применив джекнайф к TJT_{J}, что приведёт к дальнейшему уменьшению смещения, но может не сохранить дисперсию того же порядка (Robson and Whitlock 1964; см. также Thorburn 1976 и Замечание 8.3).

Задача 2.1.5
?
(a)

Любые две случайные величины XX и YY с конечными вторыми моментами удовлетворяют неравенству для ковариации [cov⁡(X,Y)]2≤var⁡(X)⋅var⁡(Y)[\operatorname {cov}(X, Y)]^{2} \leq \operatorname {var}(X) \cdot \operatorname {var}(Y).

(b)

Неравенство из пункта (a) обращается в равенство тогда и только тогда, когда существуют константы aa и bb, для которых P(X=aY+b)=1\mathbb {P}\left(X=a Y+b\right)=1.

Задача 2.1.6

Альтернативное доказательство неравенства Шварца получается из наблюдения, что

∫(f+λg)2dP=∫f2dP+2λ∫fgdP+λ2∫g2dP≥0 для всех λ, \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{ для всех } \lambda ,

так что этот квадратный трёхчлен относительно λ\lambda имеет не более одного корня.

?
Задача 2.1.7

Предположим, что XX распределена на ( 0,1 ) с плотностью вероятности pθ(x)=(1−θ)+θ/2xp_{\theta }(x)=(1-\theta )+\theta / 2 \sqrt{x} для всех 0<x<1,0≤θ≤10<x<1,0 \leq \theta \leq 1. Покажите, что не существует LMVU-оценки θ\theta.

?
Задача 2.1.8

Если δ\delta и δ′\delta^{\prime } имеют конечную дисперсию, то δ′−δ\delta^{\prime }-\delta также имеет конечную дисперсию.

?
Задача 2.1.9

В Примере 1.9,

?
(a)

найдите все несмещённые оценки нуля;

(b)

покажите, что ни одна непостоянная оценка не является UMVU.

Задача 2.1.10

Если оценки ограничены классом линейных оценок, характеризация наилучших несмещённых оценок несколько упрощается. Хотя приведённый ниже результат является следствием Теоремы 1.7, его следует установить, не используя эту теорему.

Пусть Xp×1\mathbf{X}_{p \times 1} удовлетворяет E[X]=Bψ\mathbb {E}\left[\mathbf{X}\right]=B \psi и var⁡(X)=I\operatorname {var}(\mathbf{X})=I, где Bp×rB_{p \times r} известна, а ψr×1\psi_{r \times 1} неизвестен. Линейная оценка — это оценка вида a′Xa^{\prime } \mathbf{X}, где ap×1a_{p \times 1} — известный вектор. Мы рассматриваем класс оценок

D={δ(x):δ(x)=a′x, для некоторого известного вектора a}. \mathcal{D}=\left\{ \delta (\mathbf{x}): \delta (\mathbf{x})=a^{\prime } \mathbf{x}, \text{ для некоторого известного вектора } a\right\} .
?
(a)

Для известного вектора cc покажите, что оценки из D\mathcal{D}, являющиеся несмещёнными оценками c′ψc^{\prime } \psi, удовлетворяют a′B=c′a^{\prime } B=c^{\prime }.

(b)

Пусть 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\} — класс линейных несмещённых оценок c′ψc^{\prime } \psi. Покажите, что наилучшая линейная несмещённая оценка (BLUE) величины c′ψc^{\prime } \psi, то есть линейная несмещённая оценка с минимальной дисперсией, есть δ∗(x)=a∗′x\delta^{*}(\mathbf{x})=a^{* \prime } \mathbf{x}, где a∗′=a′B(B′B)−1B′a^{* \prime }=a^{\prime } B\left(B^{\prime } B\right)^{-1} B^{\prime } и a∗′B=ca^{* \prime } B=c, с дисперсией var⁡(δ∗)=c′c\operatorname {var}\left(\delta^{*}\right)=c^{\prime } c.

(c)

Пусть 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\} — класс линейных несмещённых оценок нуля. Покажите, что если δ∈D0\delta \in \mathcal{D}_{0}, то cov⁡(δ,δ∗)=0\operatorname {cov}\left(\delta , \delta^{*}\right)=0.

(d)

Следовательно, установите аналог Теоремы 1.7 для линейных оценок: Теорема. Оценка δ∗∈Dc\delta^{*} \in \mathcal{D}_{c} удовлетворяет var⁡(δ∗)=min⁡δ∈Dcvar⁡(δ)\operatorname {var}\left(\delta^{*}\right)=\min_{\delta \in \mathcal{D}_{c}} \operatorname {var}(\delta ) тогда и только тогда, когда cov⁡(δ∗,U)=0\operatorname {cov}\left(\delta^{*}, U\right)=0, где UU — произвольная оценка из D0\mathcal{D}_{0}.

(e)

Покажите, что полученные здесь результаты можно непосредственно распространить на случай var⁡(X)=Σ\operatorname {var}(\mathbf{X})=\Sigma, где Σp×p\Sigma_{p \times p} — известная матрица, рассмотрев преобразованную задачу с X∗=Σ1/2X\mathbf{X}^{*}=\Sigma^{1 / 2} \mathbf{X} и B∗=Σ1/2BB^{*}=\Sigma^{1 / 2} B.

Задача 2.1.11

Используя Теорему 1.7, найдите UMVU-оценки некоторых ηθ(di)\eta_{\theta }\left(d_{i}\right) в модели «доза-эффект» (1.6.16) с ограничением (1.6.17) (Messig and Strawderman 1993). Пусть классы Δ\Delta и U\mathcal{U} определены так же, как в Теореме 1.7.

?
(a)

Покажите, что оценка U∈UU \in \mathcal{U} тогда и только тогда, когда 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] для произвольной константы a<∞a<\infty.

(b)

Используя пункт (a) и (1.7), покажите, что оценка δ\delta является UMVU для своего математического ожидания, только если она имеет вид δ(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), где a,b,ca, b, c и dd — произвольные константы.

(c)

Покажите, что UMVU-оценки ηθ(d1)=1−e−θ\eta_{\theta }\left(d_{1}\right)=1-e^{-\theta } не существует, но UMVU-оценкой ηθ(d2)=1−e−2θ\eta_{\theta }\left(d_{2}\right)=1-e^{-2 \theta } является δ(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)

Покажите, что LMVU-оценкой 1−e−θ1-e^{-\theta } является δ(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].

Задача 2.1.12

Покажите, что если δ(X)\delta (X) — UMVU-оценка g(θ)g(\theta ), то она является единственной UMVU-оценкой g(θ)g(\theta ). (Не предполагайте полноту, а вместо этого используйте неравенство для ковариации и условия, при которых оно обращается в равенство.)

?
Задача 2.1.13

Если δ1\delta_{1} и δ2\delta_{2} принадлежат Δ\Delta и являются UMVU-оценками g1(θ)g_{1}(\theta ) и g2(θ)g_{2}(\theta ) соответственно, то a1δ1+a2δ2a_{1} \delta_{1}+a_{2} \delta_{2} также принадлежит Δ\Delta и является UMVU-оценкой для a1g1(θ)+a2g2(θ)a_{1} g_{1}(\theta )+a_{2} g_{2}(\theta ) при любых вещественных a1a_{1} и a2a_{2}.

?
Задача 2.1.14

Полнота TT не только достаточна, но и необходима для того, чтобы каждая функция g(θ)g(\theta ), допускающая несмещённую оценку, имела единственную несмещённую оценку, являющуюся функцией от TT.

?
Задача 2.1.15

Пусть X1,…,XnX_{1}, \ldots , X_{n} — независимые одинаково распределённые случайные величины, имеющие распределение Пуассона (λ)(\lambda ).

?
(a)

Покажите, что Xˉ\bar{X} — UMVU-оценка для λ\lambda.

(b)

Для S2=∑i=1n(Xi−Xˉ)2/(n−1)S^{2}=\sum_{i=1}^{n}\left(X_{i}-\bar{X}\right)^{2} /(n-1) имеем E[S2]=E[Xˉ]=λ\mathbb {E}\left[S^{2}\right]=\mathbb {E}\left[\bar{X}\right]=\lambda. Чтобы напрямую установить, что var S2>S^{2}> var Xˉ\bar{X}, докажите, что E[S2∣Xˉ]=Xˉ\mathbb {E}\left[S^{2} \mid \bar{X}\right]=\bar{X}.

Замечание: Тождество E[S2∣Xˉ]=Xˉ\mathbb {E}\left[S^{2} \mid \bar{X}\right]=\bar{X} показывает, как полноту можно использовать при вычислении условных математических ожиданий.

Задача 2.1.16
?
(a)

Если X1,…,XnX_{1}, \ldots , X_{n} — независимые одинаково распределённые случайные величины (не обязательно нормальные) с var⁡(Xi)=σ2<∞\operatorname {var}\left(X_{i}\right)=\sigma^{2}<\infty, покажите, что δ=Σ(Xi−Xˉ)2/(n−1)\delta =\Sigma \left(X_{i}-\bar{X}\right)^{2} /(n-1) является несмещённой оценкой σ2\sigma^{2}.

(b)

Если XiX_{i} принимают значения 1 и 0 с вероятностями pp и q=1−pq=1-p, то оценка δ\delta из пункта (a) зависит только от T=ΣXiT=\Sigma X_{i} и, следовательно, является UMVU-оценкой для σ2=pq\sigma^{2}=p q. Сравните этот результат с результатом Примера 1.13.

Задача 2.1.17

Если TT имеет биномиальное распределение b(p,n)b(p, n) с n>3n>3, используйте Метод 1, чтобы найти UMVU-оценку p3p^{3}.

?
Задача 2.1.18

Пусть X1,…,XnX_{1}, \ldots , X_{n} — независимые одинаково распределённые случайные величины, имеющие распределение Пуассона P(λ)P(\lambda ). Используйте Метод 1, чтобы найти UMVU-оценку

?
(a)

λk\lambda^{k} для произвольного положительного целого kk и

(b)

e−λe^{-\lambda }.

Задача 2.1.19

Пусть X1,…,XnX_{1}, \ldots , X_{n} распределены так же, как в Примере 1.14. Используйте Метод 1, чтобы найти UMVU-оценку θk\theta^{k} для произвольного целого k>−nk>-n.

?
Задача 2.1.20

Решите Задачу 1.18(b) с помощью Метода 2, используя тот факт, что несмещённой оценкой e−λe^{-\lambda } является δ=1\delta =1, если X1=0X_{1}=0, и δ=0\delta =0 в противном случае.

?
Задача 2.1.21

В nn испытаниях Бернулли пусть Xi=1X_{i}=1 или 0 в зависимости от того, является ли ii-е испытание успехом или неудачей, и пусть T=ΣXiT=\Sigma X_{i}. Решите Задачу 1.17 с помощью Метода 2, используя тот факт, что несмещённой оценкой p3p^{3} является δ=1\delta =1, если X1=X2=X3=1X_{1}=X_{2}=X_{3}=1, и δ=0\delta =0 в противном случае.

?
Задача 2.1.22

Пусть XX принимает значения 1 и 0 с вероятностями pp и qq соответственно, и предположим, что 1/4<p<3/41 / 4<p<3 / 4. Рассмотрим задачу оценивания pp с функцией потерь L(p,d)=L(p, d)= 1, если ∣d−p∣≥1/4\left|d-p\right| \geq 1 / 4, и 0 в противном случае. Пусть δ∗\delta^{*} — рандомизированная оценка, равная Y0Y_{0} или Y1Y_{1} при X=0X=0 или 1, где Y0Y_{0} и Y1Y_{1} распределены как U(−1/2,1/2)U(-1 / 2,1 / 2) и U(1/2,3/2)U(1 / 2,3 / 2) соответственно.

?
(a)

Покажите, что δ∗\delta^{*} несмещённа.

(b)

Сравните функцию риска δ∗\delta^{*} с функцией риска XX.