2.5

Информационное неравенство

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

При условиях Задачи 1.3 для каждого p1p_{1} найдите значение LV(p1)L_{V}\left(p_{1}\right) LMVU-оценки pp в точке p1p_{1} и сравните функцию LV(p),0<p<1L_{V}(p), 0<p<1, с дисперсией Vp0(p)V_{p_{0}}(p) оценки, являющейся LMVU в точке

?
(a)

p0=1/3p_{0}=1 / 3 и

(b)

p0=1/2p_{0}=1 / 2.

Задача 2.5.2

Найдите условия, при которых в (5.1) достигается равенство.

?
Задача 2.5.3

Проверьте I(θ)I(\theta ) для распределений из Таблицы 5.1.

?
Задача 2.5.4

Если XX имеет нормальное распределение с нулевым средним и стандартным отклонением σ\sigma, найдите I(σ)I(\sigma ).

?
Задача 2.5.5

Найдите I(p)I(p) для отрицательного биномиального распределения.

?
Задача 2.5.6

Если XX имеет распределение P(λ)P(\lambda ), покажите, что содержащаяся в ней информация о λ\sqrt{\lambda } не зависит от λ\lambda.

?
Задача 2.5.7

Проверьте следующие утверждения, приведённые у Basu (1988, глава 1), которые иллюстрируют связь между информацией, достаточностью и вспомогательностью (ancillarity). Пусть I(θ)=Eθ[−∂2/∂θ2log⁡f(x∣(θ)]I(\theta )=E_{\theta }\left[-\partial^{2} / \partial \theta^{2} \log f(x \mid (\theta )]\right. — информация о θ\theta, содержащаяся в XX, а J(θ)=Eθ[−∂2/∂θ2log⁡g(T∣θ)]J(\theta )=E_{\theta }\left[-\partial^{2} / \partial \theta^{2} \log g(T \mid \theta )\right] — информация о θ\theta, содержащаяся в статистике TT, где g(⋅∣θ)g(\cdot \mid \theta ) — функция плотности TT. Определим λ(θ)=I(θ)−J(θ)\lambda (\theta )=I(\theta )-J(\theta ) — меру информации, теряемой при использовании TT вместо XX. При подходящих условиях регулярности покажите, что

?
(a)

λ(θ)≥0\lambda (\theta ) \geq 0 для всех θ\theta

(b)

λ(θ)=0\lambda (\theta )=0 тогда и только тогда, когда TT достаточна для θ\theta.

(c)

Если YY — вспомогательная статистика, но (T,Y)(T, Y) достаточна, то I(θ)=Eθ[J(θ∣Y)]I(\theta )=E_{\theta }[J(\theta \mid Y)], где

J(θ∣y)=Eθ[−∂2∂θ2log⁡h(T∣y,θ)∣ Y=y] J(\theta \mid y)=E_{\theta }\left[\left.-\frac{\partial ^{2}}{\partial \theta ^{2}} \log h(T \mid y, \theta ) \right\rvert \, Y=y\right]

а h(t∣y,θ)h(t \mid y, \theta ) — условная плотность TT при заданном Y=yY=y.

(«Условия регулярности» Basu в основном касаются перестановки интегрирования и дифференцирования. Считайте, что все подобные перестановки корректны.)

Задача 2.5.8

Найдите функцию от θ\theta, для которой количество информации не зависит от θ\theta:

?
(a)

для гамма-распределения Γ(α,β)\Gamma (\alpha , \beta ) с известным α\alpha и θ=β\theta =\beta;

(b)

для биномиального распределения b(p,n)b(p, n) с θ=p\theta =p.

Задача 2.5.9

Для обратного биномиального выборочного наблюдения (см. Пример 3.2):

?
(a)

Покажите, что наилучшая несмещённая оценка pp задаётся формулой δ∗(Y)=(m−1)/(Y+m−1)\delta^{*}(Y)=(m-1) /(Y+m-1).

(b)

Покажите, что информация о PP, содержащаяся в YY, равна I(p)=m/p2(1−p)I(p)=m / p^{2}(1-p).

(c)

Покажите, что var⁡δ∗>1/I(p)\operatorname {var} \delta^{*}>1 / I(p).

(Оценку δ∗\delta^{*} можно интерпретировать как долю успехов, если не учитывать последнее испытание, которое, как известно, обязательно является успехом.)

Задача 2.5.10

Покажите, что (5.13) можно продифференцировать, дифференцируя под знаком интеграла, когда pθ(x)p_{\theta }(x) задаётся формулой (5.24), для каждого из распределений Таблицы 5.2. [Указание: Составьте разностное отношение и примените теорему о мажорируемой сходимости.]

?
Задача 2.5.11

Проверьте записи Таблицы 5.2.

?
Задача 2.5.12

Вычислите (5.25), когда ff — плотность распределения Стьюдента (tt-распределения) с vv степенями свободы. [Указание: Используйте тот факт, что

∫−∞∞dx(1+x2)k=Γ(1/2)Γ(k−1/2)Γ(k).] \left.\int _{-\infty }^{\infty } \frac{d x}{\left(1+x^{2}\right)^{k}}=\frac{\Gamma (1 / 2) \Gamma (k-1 / 2)}{\Gamma (k)}.\right]
?
Задача 2.5.13

Для распределения с плотностью (5.24) покажите, что I(θ)I(\theta ) не зависит от θ\theta.

?
Задача 2.5.14

Проверьте

?
(a)

формулу (5.25) и

(b)

формулу (5.27).

Задача 2.5.15

Для tt-плотности сдвига вычислите границу информационного неравенства для несмещённых оценок θ\theta.

?
Задача 2.5.16
?
(a)

Для семейства масштаба с плотностью (1/θ)f(x/θ),θ>0(1 / \theta ) f(x / \theta ), \theta >0, количество информации, которое одно наблюдение XX содержит о θ\theta, равно

1θ2∫[yf′(y)f(y)+1]2f(y)dy \frac{1}{\theta ^{2}} \int \left[\frac{y f^{\prime }(y)}{f(y)}+1\right]^{2} f(y) d y
(b)

Покажите, что информация, которую XX содержит о ξ=log⁡θ\xi =\log \theta, не зависит от θ\theta.

(c)

Для распределения Коши C(0,θ)\mathrm{C}(0, \theta ) имеем I(θ)=1/(2θ2)I(\theta )=1 /\left(2 \theta^{2}\right).

Задача 2.5.17

Если pθ(x)p_{\theta }(x) задаётся формулой 1.5.1 с s=1s=1 и T(x)=δ(x)T(x)=\delta (x), покажите, что var⁡[δ(X)]\operatorname {var}[\delta (X)] достигает нижней границы (5.31) и является единственной оценкой, для которой это выполняется. [Указание: Используйте (5.18) и (1.5.15).]

?
Задача 2.5.18

Покажите, что если данная функция g(θ)g(\theta ) имеет несмещённую оценку, то существует несмещённая оценка δ\delta, которая при всех значениях θ\theta достигает нижней границы (5.1) для некоторой ψ(x,θ)\psi (x, \theta ), удовлетворяющей (5.2), тогда и только тогда, когда g(θ)g(\theta ) имеет UMVU-оценку δ0\delta_{0}. [Указание: По Теореме 5.1, ψ(x,θ)=δ0(x)\psi (x, \theta )=\delta_{0}(x) удовлетворяет (5.2). Для любой другой несмещённой оценки δ\delta имеем cov⁡(δ−δ0,δ0)=0\operatorname {cov}\left(\delta -\delta_{0}, \delta_{0}\right)=0, и, следовательно, var⁡(δ0)=[cov⁡(δ,δ0)]2/var⁡(δ0)\operatorname {var}\left(\delta_{0}\right)=\left[\operatorname {cov}\left(\delta , \delta_{0}\right)\right]^{2} / \operatorname {var}\left(\delta_{0}\right), так что ψ=δ0\psi =\delta_{0} даёт достижимую границу.] (Blyth 1974).

?
Задача 2.5.19

Покажите, что если Eθδ=g(θ)E_{\theta } \delta =g(\theta ) и var⁡(δ)\operatorname {var}(\delta ) достигает границы информационного неравенства (5.31), то

δ(x)=g(θ)+g′(θ)I(θ)∂∂θpθ(x). \delta (x)=g(\theta )+\frac{g^{\prime }(\theta )}{I(\theta )} \frac{\partial }{\partial \theta } p_{\theta }(x).
?
Задача 2.5.20

Если Eθδ=g(θ)E_{\theta } \delta =g(\theta ), то нижняя граница информационного неравенства равна IB(θ)=[g′(θ)]2/I(θ)I B(\theta )=\left[g^{\prime }(\theta )\right]^{2} / I(\theta ). Если θ=h(ξ)\theta =h(\xi ), где hh дифференцируема, покажите, что IB(ξ)=IB(θ)I B(\xi )=I B(\theta ).

?
Задача 2.5.21

Пусть XX — наблюдение из смеси нормальных распределений с плотностью

pθ(x)=122π{e−(1/2)(x−θ)2+e−(1/2)(x+θ)2},θ∈Ω, p_{\theta }(x)=\frac{1}{2 \sqrt{2 \pi }}\left\{ e^{-(1 / 2)(x-\theta )^{2}}+e^{-(1 / 2)(x+\theta )^{2}}\right\} , \theta \in \Omega ,

где Ω\Omega — произвольная окрестность нуля. Таким образом, случайная величина XX имеет либо распределение N(θ,1)N(\theta , 1), либо N(−θ,1)N(-\theta , 1), каждое с вероятностью 1/21 / 2. Покажите, что θ=0\theta =0 является особой точкой, то есть, если несмещённая оценка θ\theta существует, она будет иметь бесконечную дисперсию при θ=0\theta =0.

?
Задача 2.5.22

Пусть X1,…,XnX_{1}, \ldots , X_{n} — выборка из распределения Пуассона (λ)(\lambda ), усечённого слева в 0, то есть с выборочным пространством X={1,2,3,…}\mathcal{X}=\left\{ 1,2,3, \ldots \right\} (см. Задачу 3.20). Покажите, что нижняя граница Крамера—Рао для дисперсии несмещённых оценок λ\lambda равна

λ(1−e−λ)2n(1−e−λ−λe−λ) \frac{\lambda \left(1-e^{-\lambda }\right)^{2}}{n\left(1-e^{-\lambda }-\lambda e^{-\lambda }\right)}

и не достигается UMVU-оценкой. (Однако она является асимптотической дисперсией ML-оценки.)

?
Задача 2.5.23

Пусть X1,…,XnX_{1}, \ldots , X_{n} — независимые одинаково распределённые случайные величины с плотностью p(x,θ)p(x, \theta ), положительной при всех xx. Тогда дисперсия любой несмещённой оценки δ\delta величины θ\theta удовлетворяет

var⁡θ0(δ)≥(θ−θ0)2{∫−∞∞[p(x,θ)]2p(x,θ0)}n−1,θ≠θ0 \operatorname {var}_{\theta _{0}}(\delta ) \geq \frac{\left(\theta -\theta _{0}\right)^{2}}{\left\{ \int _{-\infty }^{\infty } \frac{[p(x, \theta )]^{2}}{p\left(x, \theta _{0}\right)}\right\} ^{n}-1}, \quad \theta \neq \theta _{0}
?
Задача 2.5.24

Если X1,…,XnX_{1}, \ldots , X_{n} — независимые одинаково распределённые случайные величины с распределением N(θ,σ2)N\left(\theta , \sigma^{2}\right), где σ\sigma известна, а про θ\theta известно, что она принимает одно из значений 0, ±1, ±2, . . ., то неравенство предыдущей задачи показывает, что любая несмещённая оценка δ\delta ограниченного параметра θ\theta удовлетворяет

var⁡θ0(δ)≥Δ2enΔ2/σ2−1,Δ≠0 \operatorname {var}_{\theta _{0}}(\delta ) \geq \frac{\Delta ^{2}}{e^{n \Delta ^{2} / \sigma ^{2}}-1}, \quad \Delta \neq 0

где Δ=θ−θ0\Delta =\theta -\theta_{0}, и, следовательно, sup⁡Δ≠0var⁡θ0(δ)≥1/[en/σ2−1]\sup_{\Delta \neq 0} \operatorname {var}_{\theta_{0}}(\delta ) \geq 1 /\left[e^{n / \sigma^{2}}-1\right].

?
Задача 2.5.25

При условиях предыдущей задачи пусть Xˉ∗\bar{X}^{*} — ближайшее к Xˉ\bar{X} целое число.

?
(a)

Оценка Xˉ∗\bar{X}^{*} несмещённа для ограниченного параметра θ\theta.

(b)

Существуют такие положительные константы aa и bb, что при всех достаточно больших nvar⁡θ(Xˉ∗)≤ae−bnn \operatorname {var}_{\theta }\left(\bar{X}^{*}\right) \leq a e^{-b n} для всех целых θ\theta.

Задача 2.5.26

Неравенство Кифера.

?
(a)

Пусть XX имеет плотность (относительно μ)p(x,θ)\mu ) p(x, \theta ), которая > 0 при всех xx, и пусть Λ1\Lambda_{1} и Λ2\Lambda_{2} — два распределения на вещественной прямой с конечными первыми моментами. Тогда любая несмещённая оценка δ\delta величины θ\theta удовлетворяет

var⁡(δ)≥[∫ΔdΛ1(Δ)−∫ΔdΛ2(Δ)]2∫ψ2(x,θ)p(x,θ)dμ(x) \operatorname {var}(\delta ) \geq \frac{\left[\int \Delta d \Lambda _{1}(\Delta )-\int \Delta d \Lambda _{2}(\Delta )\right]^{2}}{\int \psi ^{2}(x, \theta ) p(x, \theta ) d \mu (x)}

где

ψ(x,θ)=∫Ωθp(x,θ+Δ)[dΛ1(Δ)−dΛ2(Δ)]p(x,θ) \psi (x, \theta )=\frac{\int _{\Omega _{\theta }} p(x, \theta +\Delta )\left[d \Lambda _{1}(\Delta )-d \Lambda _{2}(\Delta )\right]}{p(x, \theta )}

с Ωθ={Δ:θ+ΔεΩ}\Omega_{\theta }=\left\{ \Delta : \theta +\Delta_{\varepsilon } \Omega \right\}.

(b)

Если Λ1\Lambda_{1} и Λ2\Lambda_{2} приписывают вероятность 1 значениям Δ=0\Delta =0 и Δ\Delta соответственно, неравенство сводится к (5.6) с g(θ)=θg(\theta )=\theta. [Указание: Примените (5.1).] (Kiefer 1952.)

Задача 2.5.27

Проверьте непосредственно, что следующие семейства плотностей удовлетворяют (5.38).

?
(a)

Экспоненциальное семейство из (1.5.1),

pη(x)=h(x)eηT(x)−A(η). p_{\eta }(x)=h(x) e^{\eta T(x)-A(\eta )}.
(b)

tt-семейство сдвига из Примера 5.16.

(c)

Логистическая плотность из Таблицы 1.4.1.

Задача 2.5.28

Распространите условие (5.38) на векторные параметры и покажите, что оно выполняется для экспоненциального семейства (1.5.1) при s>1s>1.

?
Задача 2.5.29

Покажите, что предположение (5.36(b)) влечёт (5.38), так что Теорема 5.15 фактически является следствием Теоремы 5.10.

?
Задача 2.5.30

Покажите, что (5.38) выполняется, если верно одно из следующих условий:

?
(a)

∣∂log⁡pθ/∂θ∣\left|\partial \log p_{\theta } / \partial \theta \right| ограничена.

(b)

[pθ+Δ(x)−pθ(x)]/Δ→∂log⁡pθ/∂θ\left[p_{\theta +\Delta }(x)-p_{\theta }(x)\right] / \Delta \rightarrow \partial \log p_{\theta } / \partial \theta равномерно.

Задача 2.5.31
?
(a)

Покажите, что если (5.38) выполняется, то семейство плотностей сильно дифференцируемо (см. Замечание 8.6).

(b)

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

Задача 2.5.32

Brown and Gajek (1990) приводят два различных достаточных условия для выполнения (8.2), которые указаны ниже. Покажите, что каждое из них влечёт (8.2). (Заметим, что при переходе от (a) к (b) условия становятся слабее, а значит, более широко применимыми и труднее проверяемыми.)

?
(a)

Для некоторого B<∞B<\infty

Eθ0[∂2∂θ2pθ(X)/pθ0(X)]2<B E_{\theta _{0}}\left[\frac{\partial ^{2}}{\partial \theta ^{2}} p_{\theta }(X) / p_{\theta _{0}}(X)\right]^{2}<B

для всех θ\theta из окрестности θ\theta.

(b)

Если pt∗(x)=∂/∂θpθ(x)∣θ=tp_{t}^{*}(x)=\partial /\left.\partial \theta p_{\theta }(x)\right|_{\theta =t}, то

lim⁡Δ→0Eθ0[pθ0+Δ∗(X)−pθ0∗(X)pθ0(X)]2=0. \lim _{\Delta \rightarrow 0} E_{\theta _{0}}\left[\frac{p_{\theta _{0}+\Delta }^{*}(X)-p_{\theta _{0}}^{*}(X)}{p_{\theta _{0}}(X)}\right]^{2}=0.
Задача 2.5.33

Пусть F\mathcal{F} — класс всех унимодальных симметричных плотностей или, более общо, плотностей, симметричных относительно нуля и удовлетворяющих f(x)≤f(0)f(x) \leq f(0) при всех xx. Покажите, что

min⁡f∈F∫x2f(x)dx=112, \min _{f \in \mathcal{F}} \int x^{2} f(x) d x=\frac{1}{12},

и что минимум достигается равномерным распределением (−12,12)\left(-\frac{1}{2}, \frac{1}{2}\right). Таким образом, равномерное распределение имеет минимальную дисперсию среди симметричных унимодальных распределений. (О свойствах масштабного равномерного распределения при больших выборках см. Пример 4.8.6.) [Указание: Побочное условие ∫f(x)dx=1\int f(x) d x=1 вместе с методом неопределённых множителей приводит к эквивалентной задаче — минимизации ∫(x2−a2)f(x)dx\int \left(x^{2}-a^{2}\right) f(x) d x, где aa выбирается так, чтобы удовлетворить ограничению. Здесь сработает рассуждение типа Неймана—Пирсона.]

?