Глава 6

Асимптотическая оптимальность

[169/0%]
Показать
LaTeX
§
Задача 6.1.1

Пусть X1,…,XnX_{1}, \ldots , X_{n} — независимые одинаково распределённые случайные величины с E(Xi)=ξE\left(X_{i}\right)=\xi.

?
(a)

Если Xi sX_{i} \mathrm{~ s} имеют конечный четвёртый момент, установите (1.3)

(b)

Для положительного целого kk покажите, что E(Xˉ−ξ)2k−1E(\bar{X}-\xi )^{2 k-1} и E(Xˉ−ξ)2kE(\bar{X}-\xi )^{2 k}, если они существуют, оба являются O(1/nk)O\left(1 / n^{k}\right).

Задача 6.1.2

Для фиксированного nn опишите относительную ошибку в Примере 1.3 как функцию от pp.

?
Задача 6.1.3

Теорема 1.1 утверждает: пусть X1,…,XnX_1, \ldots , X_n — независимые одинаково распределённые случайные величины с E(X1)=ξE(X_1)=\xi, var⁡(X1)=σ2\operatorname {var}(X_1)=\sigma^2 и конечным четвёртым моментом, и пусть hh — функция вещественной переменной, первые четыре производные которой h′,h′′,h′′′,h(iv)h', h'', h''', h^{(iv)} существуют для всех x∈Ix \in I, где II — интервал с P(X1∈I)=1P(X_1 \in I)=1. Предположим далее, что ∣h(iv)(x)∣≤M\left|h^{(iv)}(x)\right| \leq M для всех x∈Ix \in I при некотором M<∞M<\infty. Тогда,

E[h(Xˉ)]=h(ξ)+σ22nh′′(ξ)+Rn, E[h(\bar{X})]=h(\xi )+\frac{\sigma ^{2}}{2 n} h''(\xi )+R_{n},

и если, кроме того, четвёртая производная от h2h^2 также ограничена,

var⁡[h(Xˉ)]=σ2n[h′(ξ)]2+Rn, \operatorname {var}[h(\bar{X})]=\frac{\sigma ^{2}}{n}\left[h'(\xi )\right]^{2}+R_{n},

где остаточный член RnR_n в обоих случаях есть O(1/n2)O(1/n^2).

Теорема 1.5 утверждает: предположим, что предположения Теоремы 1.1 выполнены, и что cnc_n — последовательность констант, удовлетворяющая cn=1+an+O ⁣(1n2)c_{n}=1+\dfrac {a}{n}+O\! \left(\dfrac {1}{n^{2}}\right). Тогда дисперсия δn(Xˉ)=h(cnXˉ)\delta_{n}(\bar{X})=h\left(c_{n} \bar{X}\right) удовлетворяет

var⁡[δn(Xˉ)]=σ2n[h′(ξ)]2+O ⁣(1n2). \operatorname {var}\left[\delta _{n}(\bar{X})\right]=\frac{\sigma ^{2}}{n}\left[h'(\xi )\right]^{2}+O\! \left(\frac{1}{n^{2}}\right).

Докажите Теорему 1.5.

?
Задача 6.1.4

Пусть X1,…,XnX_{1}, \ldots , X_{n} — независимые одинаково распределённые случайные величины с распределением N(ξ,σ2)N\left(\xi , \sigma^{2}\right), σ2\sigma^{2} известно, и пусть g(ξ)=ξr,r=2,3,4g(\xi )=\xi^{r}, r=2,3,4. Определите, с точностью до членов порядка 1/n1 / n,

?
(a)

дисперсию UMVU-оценки g(ξ)g(\xi );

(b)

смещение MLE для g(ξ)g(\xi ).

Задача 6.1.5

Пусть X1,…,XnX_{1}, \ldots , X_{n} — независимые одинаково распределённые случайные величины с распределением N(ξ,σ2)N\left(\xi , \sigma^{2}\right), ξ\xi известно. Для чётного rr определите дисперсию UMVU-оценки (2.2.4) величины σr\sigma^{r} с точностью до членов порядка rr.

?
Задача 6.1.6

Решите предыдущую задачу для случая, когда ξ\xi неизвестно.

?
Задача 6.1.7

Для оценивания pmp^{m} в Примере 3.3.1 определите, с точностью до порядка 1/n1 / n,

?
(a)

дисперсию UMVU-оценки (2.3.2);

(b)

смещение MLE.

Задача 6.1.8

Решите предыдущую задачу, если pmp^{m} заменить оцениваемым параметром из Задачи 2.3.3.

?
Задача 6.1.9

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

?
(a)

Определите UMVU-оценку P(Xi=0)=e−θP\left(X_{i}=0\right)=e^{-\theta }.

(b)

Вычислите дисперсию оценки из пункта (a) с точностью до членов порядка 1/n1 / n.

Задача 6.1.10

Решите пункт (b) предыдущей задачи для оценки (2.3.22).

?
Задача 6.1.11

В условиях Задачи 1.1 покажите, что E∣Xˉ−ξ∣2k−1=O(n−k+1/2)E|\bar{X}-\xi |^{2 k-1}=O\left(n^{-k+1 / 2}\right). [Указание: Используйте тот факт, что E∣Xˉ−ξ∣2k−1≤[E(Xˉ−ξ)4k−2]1/2E|\bar{X}-\xi |^{2 k-1} \leq \left[E(\bar{X}-\xi )^{4 k-2}\right]^{1 / 2}, вместе с результатом Задачи 1.1.]

?
Задача 6.1.12

Получите вариант Теоремы 1.1, который требует существования и ограниченности только h′′′h^{\prime \prime \prime } вместо h(iv)h^{(i v)}, но где RnR_{n} есть лишь O(n−3/2)O\left(n^{-3 / 2}\right).

?
Задача 6.1.13

Чтобы убедиться, что Теорема 1.1 не обязательно верна без ограниченности четвёртой (или некоторой более высокой) производной, предположим, что XX распределены как N(ξ,σ2)N\left(\xi , \sigma^{2}\right), и пусть h(X)=ex4h(X)=e^{x^{4}}. Тогда все моменты XX и все производные hh существуют.

?
(a)

Покажите, что математическое ожидание h(Xˉ)h(\bar{X}) не существует ни при каком nn, и, следовательно, что E{n[h(Xˉ)−h(ξ)]}2=∞E\left\{ \sqrt{n}[h(\bar{X}) -h(\xi )]\right\}^{2}=\infty при всех значениях nn.

(b)

С другой стороны, покажите, что n[h(Xˉ)−h(ξ)]\sqrt{n}[h(\bar{X})-h(\xi )] имеет предельное нормальное распределение с конечной дисперсией, и определите эту дисперсию.

Задача 6.1.14

Пусть X1,…,XnX_{1}, \ldots , X_{n} — независимые одинаково распределённые случайные величины из показательного распределения с плотностью (1/θ)e−x/θ,x>0(1 / \theta ) e^{-x / \theta }, x>0, и θ>0\theta >0.

?
(a)

Используя Теорему 1.1, найдите приближения для E(Xˉ)E(\sqrt{\bar{X}}) и var⁡(Xˉ)\operatorname {var}(\sqrt{\bar{X}}).

(b)

Проверьте точное вычисление

var⁡(Xˉ)=[1−1n(Γ(n+1/2)Γ(n))2]θ \operatorname {var}(\sqrt{\bar{X}})=\left[1-\frac{1}{n}\left(\frac{\Gamma (n+1 / 2)}{\Gamma (n)}\right)^{2}\right] \theta

и покажите, что lim⁡n→∞nvar⁡(Xˉ)=θ/4\lim_{n \rightarrow \infty } n \operatorname {var}(\sqrt{\bar{X}})=\theta / 4.

(c)

Согласуйте результаты пунктов (a) и (b). Объясните, почему, хотя Теорема 1.1 неприменима, она дала верный ответ.

(d)

Покажите, что аналогичный вывод верен для h(x)=1/xh(x)=1 / x.

Задача 6.1.15

Пусть X1,…,XnX_{1}, \ldots , X_{n} — независимые одинаково распределённые случайные величины с распределением U(0,θ)U(0, \theta ). Определите дисперсию UMVU-оценки θk\theta^{k}, где kk — целое число, k>−nk>-n.

?
Задача 6.1.16

В условиях Задачи 1.15 найдите MLE для θk\theta^{k} и сравните её ожидаемую квадратичную ошибку с дисперсией UMVU-оценки.

?
Задача 6.1.17

Пусть X1,…,XnX_{1}, \ldots , X_{n} — независимые одинаково распределённые случайные величины с распределением U(0,θ)U(0, \theta ), пусть T=max⁡(X1,…,Xn)T=\max \left(X_{1}, \ldots , X_{n}\right), и пусть hh — функция, удовлетворяющая условиям Теоремы 1.1. Покажите, что

E[h(T)]=h(θ)−θnh′(θ)+1n2[θh′(θ)+θ2h′′(θ)]+O(1n3) E[h(T)]=h(\theta )-\frac{\theta }{n} h^{\prime }(\theta )+\frac{1}{n^{2}}\left[\theta h^{\prime }(\theta )+\theta ^{2} h^{\prime \prime }(\theta )\right]+O\left(\frac{1}{n^{3}}\right)

и

var⁡[h(T)]=θ2n2[h′(θ)]2+O(1n3). \operatorname {var}[h(T)]=\frac{\theta ^{2}}{n^{2}}\left[h^{\prime }(\theta )\right]^{2}+O\left(\frac{1}{n^{3}}\right).
?
Задача 6.1.18

Примените результаты Задачи 1.17, чтобы получить приближённые ответы к Задачам 1.15 и 1.16, и сравните ответы с точными решениями.

?
Задача 6.1.19

Если XX таковы, как в Теореме 1.1, и если первые пять производных hh существуют, а пятая производная ограничена, покажите, что

E[h(Xˉ)]=h(ξ)+12h′′σ2n+124n2[4h′′′μ3+3h(iv)σ4]+O(n−5/2) E[h(\bar{X})]=h(\xi )+\frac{1}{2} h^{\prime \prime } \frac{\sigma ^{2}}{n}+\frac{1}{24 n^{2}}\left[4 h^{\prime \prime \prime } \mu _{3}+3 h^{(i v)} \sigma ^{4}\right]+O\left(n^{-5 / 2}\right)

и, если пятая производная h2h^{2} также ограничена,

var⁡[h(Xˉ)]=(h′2)σ2n+1n2[h′h′′μ3+(h′h′′′+12h′′2)σ4]+O(n−5/2) \operatorname {var}[h(\bar{X})]=\left(h^{\prime 2}\right) \frac{\sigma ^{2}}{n}+\frac{1}{n^{2}}\left[h^{\prime } h^{\prime \prime } \mu _{3}+\left(h^{\prime } h^{\prime \prime \prime }+\frac{1}{2} h^{\prime \prime 2}\right) \sigma ^{4}\right]+O\left(n^{-5 / 2}\right)

где μ3=E(X−ξ)3\mu_{3}=E(X-\xi )^{3}.

?
Задача 6.1.20

В условиях предыдущей задачи проведите вычисление дисперсии (1.16) до членов порядка 1/n21 / n^{2} и сравните результат с результатом предыдущей задачи.

?
Задача 6.1.21

Проведите вычисление из Задачи 1.4 до членов порядка 1/n21 / n^{2}.

?
Задача 6.1.22

Для оцениваемых параметров из Задачи 1.4 вычислите ожидаемую квадратичную ошибку MLE до членов порядка 1/n21 / n^{2} и сравните её с дисперсией, вычисленной в Задаче 1.21.

?
Задача 6.1.23

Вычислите дисперсию (1.18) до членов порядка 1/n21 / n^{2} и сравните её с ожидаемой квадратичной ошибкой MLE, вычисленной до того же порядка.

?
Задача 6.1.24

Найдите дисперсию оценки (2.3.17) с точностью до членов порядка 1/n31 / n^{3}.

?
Задача 6.1.25

Для ситуации Примера 1.12 покажите, что UMVU-оценка δ1n\delta_{1 n} является MLE с поправкой на смещение, где MLE есть δ3n\delta_{3 n}.

?
Задача 6.1.26

Для оценок из Примера 1.13:

?
(a)

Вычислите их точные дисперсии.

(b)

Используя результат пункта (a), проверьте (1.27).

Задача 6.1.27
?
(a)

В условиях Теоремы 1.5, если все четвёртые моменты XiνX_{i \nu } конечны, покажите, что E(Xˉi−ξi)(Xˉj−ξj)=σij/nE\left(\bar{X}_{i}-\xi_{i}\right)\left(\bar{X}_{j}-\xi_{j}\right)=\sigma_{i j} / n и что все третьи и четвёртые моменты E(Xˉi−ξi)(Xˉj−ξj)(Xˉk−ξk)E\left(\bar{X}_{i}-\xi_{i}\right)\left(\bar{X}_{j}-\xi_{j}\right)\left(\bar{X}_{k}-\xi_{k}\right) и так далее имеют порядок 1/n21 / n^{2}.

(b)

Если, кроме того, все производные hh суммарного порядка ≤4\leq 4 существуют, а производные порядка 4 равномерно ограничены, то

E[h(Xˉ1,…,Xˉs)]=h(ξ1,…,ξs)+12n∑i=1s∑j=1sσij∂2h(ξ1,…,ξs)∂ξi∂ξj+Rn, E\left[h\left(\bar{X}_{1}, \ldots , \bar{X}_{s}\right)\right]=h\left(\xi _{1}, \ldots , \xi _{s}\right)+\frac{1}{2 n} \sum _{i=1}^{s} \sum _{j=1}^{s} \sigma _{i j} \frac{\partial ^{2} h\left(\xi _{1}, \ldots , \xi _{s}\right)}{\partial \xi _{i} \partial \xi _{j}}+R_{n},

и если производные h2h^{2} порядка 4 также ограничены,

var⁡[h(Xˉ1,…,Xˉs)]=1nΣΣσij∂h∂ξi∂h∂ξj+Rn \operatorname {var}\left[h\left(\bar{X}_{1}, \ldots , \bar{X}_{s}\right)\right]=\frac{1}{n} \Sigma \Sigma \sigma _{i j} \frac{\partial h}{\partial \xi _{i}} \frac{\partial h}{\partial \xi _{j}}+R_{n}

где остаточный член RnR_{n} в обоих случаях есть O(1/n2)O\left(1 / n^{2}\right).

Задача 6.1.28

На основе выборки из N(ξ,σ2)N\left(\xi , \sigma^{2}\right) пусть Pn(ξ,σ)P_{n}(\xi , \sigma ) — вероятность того, что UMVU-оценка Xˉ2−σ2/n\bar{X}^{2}-\sigma^{2} / n величины ξ2\xi^{2} (σ\sigma известно) отрицательна.

?
(a)

Покажите, что Pn(ξ,σ)P_{n}(\xi , \sigma ) — убывающая функция от n∣ξ∣/σ\sqrt{n}\left|\xi \right| / \sigma.

(b)

Покажите, что Pn(ξ,σ)→0P_{n}(\xi , \sigma ) \rightarrow 0 при n→∞n \rightarrow \infty для любых фиксированных ξ≠0\xi \neq 0 и σ\sigma.

(c)

Определите значение Pn(0,σ)P_{n}(0, \sigma ).

Задача 6.1.29

Используя tt-распределение, найдите значение Pn(0,σ)P_{n}(0, \sigma ) из предыдущей задачи для UMVU-оценки ξ2\xi^{2}, когда σ\sigma неизвестно, для характерных значений nn.

?
Задача 6.1.30

Теорема 1.9 утверждает: если последовательность случайных величин TnT_{n} стремится к ξ\xi по вероятности, и если hh непрерывна в точке ξ\xi, то h(Tn)→Ph(ξ)h(T_{n}) \xrightarrow {P} h(\xi ). Чтобы показать, что P(∣h(Tn)−h(ξ)∣<a)→1P(\left|h(T_{n})-h(\xi )\right|<a) \rightarrow 1, достаточно заметить, что в силу непрерывности hh разность ∣h(Tn)−h(ξ)∣\left|h(T_{n})-h(\xi )\right| будет сколь угодно мала, если TnT_{n} близко к ξ\xi, и что для любого aa вероятность того, что ∣Tn−ξ∣<a\left|T_{n}-\xi \right|<a, стремится к 11 при n→∞n \rightarrow \infty.

Заполните детали доказательства Теоремы 1.9. (См. также Задачу 1.8.8.)

?
Задача 6.1.31

В Примере 8.13 при θ=0\theta =0 покажите, что δ2n\delta_{2 n} не имеет в точности распределение σ2(χ12−1)/n\sigma^{2}\left(\chi_{1}^{2}-1\right) / n.

?
Задача 6.1.32

В Примере 8.13 пусть δ4n=max⁡(0,Xˉ2−σ2/n)\delta_{4 n}=\max \left(0, \bar{X}^{2}-\sigma^{2} / n\right), что является улучшением по сравнению с δ1n\delta_{1 n}.

?
(a)

Покажите, что n(δ4n−θ2)\sqrt{n}\left(\delta_{4 n}-\theta^{2}\right) имеет то же предельное распределение, что и n(δ1n−θ2)\sqrt{n}\left(\delta_{1 n}-\theta^{2}\right), при θ≠0\theta \neq 0.

(b)

Опишите предельное распределение nδ4nn \delta_{4 n} при θ=0\theta =0.

Задача 6.1.33

Пусть XX имеет биномиальное распределение b(p,n)b(p, n), и пусть g(p)=pqg(p)=p q. UMVU-оценкой g(p)g(p) является δ=X(n−X)/n(n−1)\delta =X(n-X) / n(n-1). Определите предельное распределение n(δ−pq)\sqrt{n}(\delta - p q) и n(δ−pq)n(\delta -p q), когда g′(p)≠0g^{\prime }(p) \neq 0 и g′(p)=0g^{\prime }(p)=0 соответственно.

?
Задача 6.1.34

Пусть X1,…,XnX_{1}, \ldots , X_{n} — независимые одинаково распределённые случайные величины с распределением N(ξ,1)N(\xi , 1). Определите предельное поведение распределения UMVU-оценки p=P[∣Xi∣≤u]p=P\left[\left|X_{i}\right| \leq u\right].

?
Задача 6.1.35

Определите предельное поведение оценки (2.3.22) при n→∞n \rightarrow \infty.

?
Задача 6.1.36

Пусть X1,…,XnX_{1}, \ldots , X_{n} — независимые одинаково распределённые случайные величины с распределением PθP_{\theta }, и предположим, что δn\delta_{n} является UMVU-оценкой g(θ)g(\theta ) на основе X1,…,XnX_{1}, \ldots , X_{n}. Если существует n0n_{0} и несмещённая оценка δ0(X1,…,Xn0)\delta_{0}\left(X_{1}, \ldots , X_{n_{0}}\right), имеющая конечную дисперсию при всех θ\theta, то δn\delta_{n} состоятельна для g(θ)g(\theta ).

?
Задача 6.1.37

Пусть YnY_{n} имеет распределение N(0,1)N(0,1) с вероятностью πn\pi_{n} и распределение N(0,τn2)N\left(0, \tau_{n}^{2}\right) с вероятностью 1−πn1-\pi_{n}. Если τn→∞\tau_{n} \rightarrow \infty и πn→π\pi_{n} \rightarrow \pi, определите, при каких значениях π\pi последовательность {Yn}\left\{ Y_{n}\right\} имеет предельное распределение, а при каких — нет.

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

В Задаче 1.37 определите, к каким значениям может стремиться var⁡(Yn)\operatorname {var}\left(Y_{n}\right) при n→∞n \rightarrow \infty, если πn→1\pi_{n} \rightarrow 1 и τn→∞\tau_{n} \rightarrow \infty, но в остальном обе последовательности произвольны.

(b)

Используя (a), покажите, что предел дисперсии не обязан совпадать с дисперсией предельного распределения.

Задача 6.1.39

Пусть bm,n,m,n=1,2,…b_{m, n}, m, n=1,2, \ldots — двойная последовательность вещественных чисел, которая при каждом фиксированном mm не убывает по nn. Покажите, что lim⁡n→∞lim⁡m→∞bm,n=lim⁡m,n→∞inf⁡bm,n\lim_{n \rightarrow \infty } \lim_{m \rightarrow \infty } b_{m, n}=\lim_{m, n \rightarrow \infty } \inf b_{m, n} и lim⁡m→∞lim⁡n→∞bm,n=lim⁡m,n→∞sup⁡bm,n\lim_{m \rightarrow \infty } \lim_{n \rightarrow \infty } b_{m, n}=\lim_{m, n \rightarrow \infty } \sup b_{m, n}, при условии что указанные пределы существуют (они могут быть бесконечны), и где lim⁡inf⁡bm,n\lim \inf b_{m, n} и lim⁡sup⁡bm,n\lim \sup b_{m, n} обозначают, соответственно, наименьшую и наибольшую предельные точки, достижимые последовательностью bmk,nk,k=1,2,…b_{m_{k}, n_{k}}, k=1,2, \ldots, с mk→∞m_{k} \rightarrow \infty и nk→∞n_{k} \rightarrow \infty.

?
§
Задача 6.2.1

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

Tn={Xˉn если Sn≤ann если Sn>an, T_{n}= \begin{cases} \bar{X}_{n} & \text{ если } S_{n} \leq a_{n} \\ n & \text{ если } S_{n}>a_{n},\end{cases}

где Sn=Σ(Xi−Xˉ)2,P(Sn>an)=1/nS_{n}=\Sigma \left(X_{i}-\bar{X}\right)^{2}, P\left(S_{n}>a_{n}\right)=1 / n, и Tn′=(X1+⋯+Xkn)/knT_{n}^{\prime }=\left(X_{1}+\cdots +X_{k_{n}}\right) / k_{n} с knk_{n} — наибольшим целым числом ≤n\leq \sqrt{n}.

?
(a)

Покажите, что асимптотическая эффективность Tn′T_{n}^{\prime } относительно TnT_{n} равна нулю.

(b)

Покажите, что для любого фиксированного ε>0,P[∣Tn−θ∣>ε]=1n+o(1n)\varepsilon >0, P\left[\left|T_{n}-\theta \right|>\varepsilon \right]=\frac{1}{n}+o\left(\frac{1}{n}\right), но P[∣Tn′−θ∣>ε]=o(1n)P\left[\left|T_{n}^{\prime }-\theta \right|>\varepsilon \right]= o\left(\frac{1}{n}\right).

(c)

Для больших значений nn что можно сказать о двух вероятностях из пункта (b), когда ε\varepsilon заменяется на a/na / \sqrt{n}? (Basu 1956).

Задача 6.2.2

Если kn[δn−g(θ)]→LHk_{n}\left[\delta_{n}-g(\theta )\right] \xrightarrow {\mathcal{L}} H для некоторой последовательности knk_{n}, покажите, что тот же результат верен, если knk_{n} заменить на kn′k_{n}^{\prime }, где kn/kn′→1k_{n} / k_{n}^{\prime } \rightarrow 1.

?
Задача 6.2.3

Предположим, что распределение Yn=n(δn−g(θ))Y_{n}=\sqrt{n}\left(\delta_{n}-g(\theta )\right) сходится к распределению со средним 0 и дисперсией v(θ)v(\theta ). Используя лемму Фату (Лемма 1.2.6), установите, что var⁡θ(δn)→0\operatorname {var}_{\theta }\left(\delta_{n}\right) \rightarrow 0 для всех θ\theta.

?
Задача 6.2.4

Если X1,…,XnX_{1}, \ldots , X_{n} — выборка из однопараметрического экспоненциального семейства (1.5.2), то ΣT(Xi)\Sigma T\left(X_{i}\right) является минимальной достаточной статистикой, и E[(1/n)ΣT(Xi)]=(∂/∂η)A(η)=τE\left[(1 / n) \Sigma T\left(X_{i}\right)\right]=(\partial / \partial \eta ) A(\eta )=\tau. Покажите, что для любой функции g(⋅)g(\cdot ), для которой верна Теорема 1.8.12, g((1/n)ΣT(Xi))g\left((1 / n) \Sigma T\left(X_{i}\right)\right) асимптотически несмещена для g(τ)g(\tau ).

?
Задача 6.2.5

Если X1,…,XnX_{1}, \ldots , X_{n} — независимые одинаково распределённые случайные величины с распределением n(μ,σ2)n\left(\mu , \sigma^{2}\right), покажите, что Sr=[1/(n−1)Σ(xi−xˉ)2]r/2S^{r}=\left[1 /(n-1) \Sigma \left(x_{i}-\bar{x}\right)^{2}\right]^{r / 2} является асимптотически несмещённой оценкой σr\sigma^{r}.

?
Задача 6.2.6

Пусть X1,…,XnX_{1}, \ldots , X_{n} — независимые одинаково распределённые случайные величины с распределением U(0,θ)U(0, \theta ). Согласно Примеру 2.1.14, δn=(n+1)X(n)/n\delta_{n}=(n+1) X_{(n)} / n является UMVU-оценкой θ\theta, тогда как MLE есть X(n)X_{(n)}. Определите предельное распределение

?
(a)

n[θ−δn]n\left[\theta -\delta_{n}\right] и

(b)

n[θ−X(n)]n\left[\theta -X_{(n)}\right].

Прокомментируйте асимптотическое смещение этих оценок.

Задача 6.2.7

Для ситуации Задачи 2.6:

?
(a)

Вычислите средние квадратичные ошибки δn\delta_{n} и X(n)X_{(n)} как оценок θ\theta.

(b)

Покажите, что

lim⁡n→∞E(X(n)−θ)2E(δn−θ)2=2. \lim _{n \rightarrow \infty } \frac{E\left(X_{(n)}-\theta \right)^{2}}{E\left(\delta _{n}-\theta \right)^{2}}=2.
Задача 6.2.8

Проверьте асимптотическое распределение, заявленное для δn\delta_{n} в Примере 2.5.

?
Задача 6.2.9

Пусть δn\delta_{n} — произвольная оценка, удовлетворяющая (2.2) с g(θ)=θg(\theta )=\theta. Постройте последовательность δn′\delta_{n}^{\prime }, такую что n(δn′−θ)→LN[0,w2(θ)]\sqrt{n}\left(\delta_{n}^{\prime }-\theta \right) \xrightarrow {\mathcal{L}} N\left[0, w^{2}(\theta )\right] с w(θ)=v(θ)w(\theta )=v(\theta ) при θ≠θ0\theta \neq \theta_{0} и w(θ0)=0w\left(\theta_{0}\right)=0.

?
Задача 6.2.10

В предыдущей задаче постройте δn′\delta_{n}^{\prime }, такую что w(θ)=v(θ)w(\theta )=v(\theta ) для всех θ≠θ0\theta \neq \theta_{0} и θ1\theta_{1}, и <v(θ)<v(\theta ) при θ=θ0\theta =\theta_{0} и θ1\theta_{1}.

?
Задача 6.2.11

Постройте последовательность {δn}\left\{ \delta_{n}\right\}, удовлетворяющую (2.2), но для которой смещение bn(θ)b_{n}(\theta ) не стремится к нулю.

?
Задача 6.2.12

В Примере 2.7 с Rn(θ)R_{n}(\theta ), заданной формулой (2.11), покажите, что Rn(θ)→1R_{n}(\theta ) \rightarrow 1 при θ≠0\theta \neq 0 и что Rn(0)→a2R_{n}(0) \rightarrow a^{2}.

?
Задача 6.2.13

Пусть bn(θ)=Eθ(δn)−θb_{n}(\theta )=E_{\theta }\left(\delta_{n}\right)-\theta — смещение оценки δn\delta_{n} из Примера 2.5.

?
(a)

Покажите, что

bn(θ)=−(1−a)n∫−n4n4xϕ(x−nθ)dx b_{n}(\theta )=\frac{-(1-a)}{\sqrt{n}} \int _{-\sqrt[4]{n}}^{\sqrt[4]{n}} x \phi (x-\sqrt{n} \theta ) d x
(b)

Покажите, что bn′(θ)→0b_{n}^{\prime }(\theta ) \rightarrow 0 при любом θ≠0\theta \neq 0 и bn′(0)→(1−a)b_{n}^{\prime }(0) \rightarrow (1-a).

(c)

Используя (b), объясните, как оценка Ходжеса δn\delta_{n} может нарушать (2.7), не нарушая при этом информационное неравенство.

Задача 6.2.14

В Примере 2.7 покажите, что если θn=c/n\theta_{n}=c / \sqrt{n}, то Rn(θn)→a2+c2(1−a)2R_{n}\left(\theta_{n}\right) \rightarrow a^{2}+c^{2}(1-a)^{2}.

?
§
Задача 6.3.1

Пусть XX имеет биномиальное распределение b(p,n),0≤p≤1b(p, n), 0 \leq p \leq 1. Определите MLE для pp

?
(a)

обычным методом математического анализа, определяя максимум функции;

(b)

показав, что pxqn−x≤(x/n)x[(n−x)/n]n−xp^{x} q^{n-x} \leq (x / n)^{x}[(n-x) / n]^{n-x}.

Задача 6.3.2

В предыдущей задаче покажите, что MLE не существует, когда pp ограничено условием 0<p<10<p<1, а x=0x=0 или =n=n.

?
Задача 6.3.3

Пусть X1,…,XnX_{1}, \ldots , X_{n} — независимые одинаково распределённые случайные величины с распределением N(ξ,σ2)N\left(\xi , \sigma^{2}\right). Определите MLE для

?
(a)

ξ\xi, когда σ\sigma известно,

(b)

σ\sigma, когда ξ\xi известно, и

(c)

(ξ,σ)(\xi , \sigma ), когда оба неизвестны.

Задача 6.3.4

Предположим, что X1,…,XnX_{1}, \ldots , X_{n} — независимые одинаково распределённые случайные величины с распределением N(ξ,1)N(\xi , 1), где ξ>0\xi >0. Покажите, что MLE есть Xˉ\bar{X}, когда Xˉ>0\bar{X}>0, и не существует, когда Xˉ≤0\bar{X} \leq 0.

?
Задача 6.3.5

Пусть XX принимает значения 0 и 1 с вероятностями pp и qq соответственно. Если известно, что 1/3≤p≤2/31 / 3 \leq p \leq 2 / 3,

?
(a)

найдите MLE и

(b)

покажите, что ожидаемая квадратичная ошибка MLE равномерно больше, чем у δ(x)=1/2\delta (x)=1 / 2.

[Аналогичная задача оценивания возникает в опросах с рандомизированным ответом. См. Пример 5.2.2.]

Задача 6.3.6

Когда Ω\Omega конечно, покажите, что MLE состоятельна тогда и только тогда, когда она удовлетворяет (3.2).

?
Задача 6.3.7

Покажите, что Теорема 3.2 остаётся верной, если предположение A1 ослабить до A1': Существует непустое множество Ω0∈Ω\Omega_{0} \in \Omega, такое что θ0∈Ω0\theta_{0} \in \Omega_{0}, и Ω0\Omega_{0} содержится в носителе каждого PθP_{\theta }.

?
Задача 6.3.8

Докажите существование и единственность 0<ak<ak−1,k=1,2,…0<a_{k}<a_{k-1}, k=1,2, \ldots, удовлетворяющих (3.4).

?
Задача 6.3.9

Докажите (3.9).

?
Задача 6.3.10

В Примере 3.6 с 0<c<1/20<c<1 / 2 определите состоятельную оценку kk.

?
Задача 6.3.11

Проверьте характер корней в Примере 3.9.

?
Задача 6.3.12

Пусть XX имеет распределение N(θ,1)N(\theta , 1). Покажите, что при условии a<X<ba<X<b переменная XX стремится по вероятности к bb при θ→∞\theta \rightarrow \infty.

?
Задача 6.3.13

Рассмотрим выборку X1,…,XnX_{1}, \ldots , X_{n} из распределения Пуассона, обусловленного положительностью, так что P(Xi=x)=θxe−θ/x!(1−e−θ)P\left(X_{i}=x\right)=\theta^{x} e^{-\theta } / x!\left(1-e^{-\theta }\right) при x=1,2,…x=1,2, \ldots. Покажите, что уравнение правдоподобия имеет единственный корень при всех значениях xx.

?
Задача 6.3.14

Пусть XX имеет отрицательное биномиальное распределение (2.3.3). Найдите ELE для pp.

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

Плотность называется строго унимодальной, или, что то же самое, логарифмически вогнутой, если log⁡f(x)\log f(x) является вогнутой функцией. Покажите, что такая плотность имеет единственную моду.

(b)

Пусть X1,…,XnX_{1}, \ldots , X_{n} — независимые одинаково распределённые случайные величины с плотностью f(x−θ)f(x-\theta ). Покажите, что функция правдоподобия имеет единственный корень, если f′(x)/f(x)f^{\prime }(x) / f(x) монотонна, и что этот корень является максимумом, если f′(x)/f(x)f^{\prime }(x) / f(x) убывает. Следовательно, логарифмически вогнутые плотности дают единственные MLE.

(c)

Пусть X1,…,XnX_{1}, \ldots , X_{n} — положительные случайные величины (или симметрично распределённые относительно нуля) с совместной плотностью anΠf(axi),a>0a^{n} \Pi f\left(a x_{i}\right), a>0. Покажите, что уравнение правдоподобия имеет единственный максимум, если xf′(x)/f(x)x f^{\prime }(x) / f(x) строго убывает при x>0x>0.

(d)

Если X1,…,XnX_{1}, \ldots , X_{n} — независимые одинаково распределённые случайные величины с плотностью f(xi−θ)f\left(x_{i}-\theta \right), где ff унимодальна, и если уравнение правдоподобия имеет единственный корень, покажите, что уравнение правдоподобия также имеет единственный корень, когда плотность каждой XiX_{i} равна af[a(xi−θ)]a f\left[a\left(x_{i}-\theta \right)\right], при известном aa.

Задача 6.3.16

Для каждой из следующих плотностей f(⋅)f(\cdot ) определите,

(a) является ли она строго унимодальной, и

(b) является ли xf′(x)/f(x)x f^{\prime }(x) / f(x) строго убывающей при x>0x>0.

Соответственно прокомментируйте, имеют ли соответствующие параметры сдвига и масштаба единственные MLE:

?
(a)

f(x)=12πe−12x2,−∞<x<∞f(x)=\frac{1}{\sqrt{2 \pi }} e^{-\frac{1}{2} x^{2}},-\infty <x<\infty

(b)

f(x)=12π1xe−12(log⁡x)2,0≤x<∞f(x)=\frac{1}{\sqrt{2 \pi }} \frac{1}{x} e^{-\frac{1}{2}(\log x)^{2}}, 0 \leq x<\infty

(c)

f(x)=e−x/(1+e−x)2,−∞<x<∞f(x)=e^{-x} /\left(1+e^{-x}\right)^{2},-\infty <x<\infty

(d)

f(x)=Γ(v+1/2)Γ(v/2)1vπ1[1+(x/v)2]v+12,−∞<x<∞f(x)=\frac{\Gamma (v+1 / 2)}{\Gamma (v / 2)} \frac{1}{\sqrt{v \pi }} \frac{1}{\left[1+(x / v)^{2}\right]^{\frac{v+1}{2}}},-\infty <x<\infty

Задача 6.3.17

Если X1,…,XnX_{1}, \ldots , X_{n} — независимые одинаково распределённые случайные величины с плотностью f(xi−θ)f\left(x_{i}-\theta \right) или af(axi)a f\left(a x_{i}\right), и ff — логистическая плотность L(0,1)L(0,1), уравнение правдоподобия имеет единственные решения θ^\hat{\theta } и a^\hat{a} как в случае сдвига, так и в случае масштаба. Определите предельное распределение n(θ^−θ)\sqrt{n}(\hat{\theta }-\theta ) и n(a^−a)\sqrt{n}(\hat{a}-a).

?
Задача 6.3.18

В Задаче 3.15(b), при ff — плотности Коши C(0,a)C(0, a), уравнение правдоподобия имеет единственный корень a^\hat{a}, и n(a^−a)→LN(0,2a2)\sqrt{n}(\hat{a}-a) \xrightarrow {\mathcal{L}} N\left(0,2 a^{2}\right).

?
Задача 6.3.19

Если X1,…,XnX_{1}, \ldots , X_{n} — независимые одинаково распределённые случайные величины с распределением C(θ,1)C(\theta , 1), то при любом фиксированном nn существует положительная вероятность того,

?
(a)

что уравнение правдоподобия имеет 2n−12 n-1 корней, и

(b)

что уравнение правдоподобия имеет единственный корень.

Задача 6.3.20

Если X1,…,XnX_{1}, \ldots , X_{n} — независимые одинаково распределённые случайные величины с гамма-распределением Γ(θ,1)\Gamma (\theta , 1), уравнение правдоподобия имеет единственный корень.

?
Задача 6.3.21

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

fθ(x)=θxθ−1e−xθ,x>0,θ>0, f_{\theta }(x)=\theta x^{\theta -1} e^{-x^{\theta }}, \quad x>0, \theta >0,

которое не принадлежит ни экспоненциальному семейству, ни семейству сдвига, ни семейству масштаба. Тем не менее покажите, что существует единственный внутренний максимум функции правдоподобия.

?
Задача 6.3.22

В условиях Теоремы 3.2 покажите, что

[L(θ0+1n)−L(θ0)+12I(θ0)]/I(θ0) \left[L\left(\theta _{0}+\frac{1}{\sqrt{n}}\right)-L\left(\theta _{0}\right)+\frac{1}{2} I\left(\theta _{0}\right)\right] / \sqrt{I\left(\theta _{0}\right)}

стремится по распределению к N(0,1)N(0,1).

?
Задача 6.3.23

Пусть X1,…,XnX_{1}, \ldots , X_{n} — независимые одинаково распределённые случайные величины с распределением N(θ,aθ2),θ>0N\left(\theta , a \theta^{2}\right), \theta >0, где aa — известная положительная константа.

?
(a)

Найдите явное выражение для ELE параметра θ\theta.

(b)

Определите, существует ли MRE-оценка при подходящей группе преобразований.

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

Этот случай рассматривался в Berk (1972).

Задача 6.3.24

Проверьте, что предположения Теоремы 3.10 выполнены в Примере 3.12.

?
Задача 6.3.25

Для X1,…,XnX_{1}, \ldots , X_{n} — независимых одинаково распределённых случайных величин с распределением DE(θ,1)D E(\theta , 1) покажите, что

?
(a)

выборочная медиана является MLE для θ\theta, и

(b)

выборочная медиана асимптотически нормальна с дисперсией 1/n1 / n — границей информационного неравенства.

Задача 6.3.26

В Примере 3.12 покажите непосредственно, что (1/n)ΣT(Xi)(1 / n) \Sigma T\left(X_{i}\right) является асимптотически эффективной оценкой θ=Eη[T(X)]\theta =E_{\eta }[T(X)], рассмотрев её предельное распределение.

?
Задача 6.3.27

Пусть X1,…,XnX_{1}, \ldots , X_{n} — независимые одинаково распределённые случайные величины с распределением θg(x)+(1−θ)h(x)\theta g(x)+(1-\theta ) h(x), где (g,h)(g, h) — пара заданных плотностей вероятности относительно μ\mu, и 0<θ<10<\theta <1.

?
(a)

Приведите один пример (g,h)(g, h), для которого предположения Теоремы 3.10 выполнены, и один, для которого они не выполнены.

(b)

Обсудите существование и характер корней уравнения правдоподобия при n=1n=1, 2, 3.

Задача 6.3.28

В условиях Теоремы 3.7 предположим, что θ^1n\hat{\theta }_{1 n} и θ^2n\hat{\theta }_{2 n} — две состоятельные последовательности корней уравнения правдоподобия. Докажите, что Pθ0(θ^1n=θ^2n)→1P_{\theta_{0}}\left(\hat{\theta }_{1 n}=\hat{\theta }_{2 n}\right) \rightarrow 1 при n→∞n \rightarrow \infty. [Указание:

?
(a)

Пусть Sn={x:x=(x1,…,xn) такие, что θ^1n(x)≠θ^2n(x)}S_{n}=\left\{ \mathbf{x}: \mathbf{x}=\left(x_{1}, \ldots , x_{n}\right) \text{ такие, что } \hat{\theta }_{1 n}(\mathbf{x}) \neq \hat{\theta }_{2 n}(\mathbf{x})\right\}. Для всех x∈Sn\mathbf{x} \in S_{n} существует θn∗\theta_{n}^{*} между θ^1n\hat{\theta }_{1 n} и θ^2n\hat{\theta }_{2 n}, такое что L′′(θn∗)=0L^{\prime \prime }\left(\theta_{n}^{*}\right)=0. Для всех x∉Sn\mathbf{x} \notin S_{n} пусть θn∗\theta_{n}^{*} — общее значение θ^1n\hat{\theta }_{1 n} и θ^2n\hat{\theta }_{2 n}. Тогда θn∗\theta_{n}^{*} — состоятельная последовательность корней уравнения правдоподобия.

(b)

(1/n)L′′(θn∗)−(1/n)L′′(θ0)→0(1 / n) L^{\prime \prime }\left(\theta_{n}^{*}\right)-(1 / n) L^{\prime \prime }\left(\theta_{0}\right) \rightarrow 0 по вероятности, и, следовательно, (1/n)L′′(θn∗)→−I(θ0)(1 / n) L^{\prime \prime }\left(\theta_{n}^{*}\right) \rightarrow -I\left(\theta_{0}\right) по вероятности.

(c)

Пусть 0<ε<I(θ0)0<\varepsilon <I\left(\theta_{0}\right), и пусть

Sn′={x:1nL′′(θn∗)<−I(θ0)+ε}. S_{n}^{\prime }=\left\{ \mathbf{x}: \frac{1}{n} L^{\prime \prime }\left(\theta _{n}^{*}\right)<-I\left(\theta _{0}\right)+\varepsilon \right\} .

Тогда Pθ0(Sn′)→1P_{\theta_{0}}\left(S_{n}^{\prime }\right) \rightarrow 1. С другой стороны, L′′(θn∗)=0L^{\prime \prime }\left(\theta_{n}^{*}\right)=0 на SnS_{n}, так что SnS_{n} содержится в дополнении к Sn′S_{n}^{\prime } (Huzurbazar 1948).]

Задача 6.3.29

Чтобы установить измеримость последовательности корней θ^n∗\hat{\theta }_{n}^{*} из Теоремы 3.7, можно следовать доказательству Serfling (1980, Раздел 4.2.2), где доказана измеримость аналогичной последовательности.

?
(a)

Для определённости определим θ^n(a)\hat{\theta }_{n}(a) как значение, минимизирующее ∣θ^−θ0∣\left|\hat{\theta }-\theta_{0}\right| при условии

θ0−a≤θ^≤θ0+a и ∂∂θl(θ∣x)∣θ=θ^=0. \theta _{0}-a \leq \hat{\theta } \leq \theta _{0}+a \quad \text{ и }\left.\quad \frac{\partial }{\partial \theta } l(\theta \mid \mathbf{x})\right|_{\theta =\hat{\theta }}=0.

Покажите, что θ~n(a)\tilde{\theta }_{n}(a) измеримо.

(b)

Покажите, что θn∗\theta_{n}^{*} — корень, ближайший к θ∗\theta^{*}, — измерим.

§
Задача 6.4.1

Пусть

u(t)={c∫0te−1/x(1−x)dx при 0<t<10 при t≤01 при t≥1 u(t)= \begin{cases} c \int _{0}^{t} e^{-1 / x(1-x)} d x & \text{ при } 0<t<1 \\ 0 & \text{ при } t \leq 0 \\ 1 & \text{ при } t \geq 1\end{cases}

Покажите, что при подходящем cc функция uu непрерывна и бесконечно дифференцируема при −∞<t<∞-\infty <t<\infty.

?
Задача 6.4.2

Покажите, что плотность (4.1) с Ω=(0,∞)\Omega =(0, \infty ) удовлетворяет всем условиям Теоремы 3.10, за исключением условия (d) Теоремы 2.6.

?
Задача 6.4.3

Покажите, что плотность (4.4) с Ω=(0,∞)\Omega =(0, \infty ) удовлетворяет всем условиям Теоремы 3.10.

?
Задача 6.4.4

В Примере 4.5 вычислите оценки (4.8) и (4.14) для случая Коши, используя в качестве θ~n\tilde{\theta }_{n} выборочную медиану.

?
Задача 6.4.5

В Примере 4.7 покажите, что l(θ)l(\theta ) вогнута.

?
Задача 6.4.6

В Примере 4.7, если η=ξ\eta =\xi, покажите, как получить n\sqrt{n}-состоятельную оценку, приравнивая выборочный и теоретический вторые моменты.

?
Задача 6.4.7

В Теореме 4.8 покажите, что σ11=σ12\sigma_{11}=\sigma_{12}.

?
Задача 6.4.8

Не используя Теорему 4.8, покажите в Примере 4.13, что EM-последовательность сходится к MLE.

?
Задача 6.4.9

Рассмотрим следующие 12 наблюдений из двумерного нормального распределения с параметрами μ1=μ2=0,σ12,σ22,ρ\mu_{1}=\mu_{2}=0, \sigma_{1}^{2}, \sigma_{2}^{2}, \rho:

x1x_{1}11-1-122-2-2****
x2x_{2}1-11-1****22-2-2

где «∗*» обозначает пропущенное значение.

?
(a)

Покажите, что функция правдоподобия имеет глобальные максимумы при ρ=±1/2,σ12=σ22=8/3\rho = \pm 1 / 2, \sigma_{1}^{2}=\sigma_{2}^{2}=8 / 3, и седловую точку при ρ=0,σ12=σ22=5/2\rho =0, \sigma_{1}^{2}=\sigma_{2}^{2}=5 / 2.

(b)

Покажите, что если EM-последовательность начинается с ρ=0\rho =0, то она остаётся при ρ=0\rho =0 на всех последующих итерациях.

(c)

Покажите, что если EM-последовательность начинается с ρ\rho, отделённого от нуля, то она сходится к максимуму.

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

Эта задача принадлежит Murray (1977) и обсуждается в Wu (1983).

Задача 6.4.10

Покажите, что если плотность полных данных EM f(y,z∣θ)f(\mathbf{y}, \mathbf{z} \mid \theta ) из (4.21) принадлежит искривлённому экспоненциальному семейству, то предположения Теоремы 4.12 выполнены.

?
Задача 6.4.11

В EM-алгоритме вычисление E-шага, то есть вычисление математического ожидания, может быть сложным. В таких случаях E-шаг иногда можно заменить методом Монте-Карло, что приводит к MCEM-алгоритму (Wei and Tanner 1990). Рассмотрим следующее MCEM-вычисление Q(θ∣θ^(j),y)Q\left(\theta \mid \hat{\theta }_{(j)}, \mathbf{y}\right): При заданном θ^(j)(k)\hat{\theta }_{(j)}^{(k)}

  1. Сгенерируйте Z1,⋯ ,ZkZ_{1}, \cdots , Z_{k} — независимые одинаково распределённые случайные величины из k(z∣θ^(j)(k),y)k\left(\mathbf{z} \mid \hat{\theta }_{(j)}^{(k)}, \mathbf{y}\right),

  2. Положите Q^(θ∣θ^(j)(k),y)=1k∑i=1klog⁡L(θ∣y,z)\hat{Q}\left(\theta \mid \hat{\theta }_{(j)}^{(k)}, \mathbf{y}\right)=\frac{1}{k} \sum_{i=1}^{k} \log L(\theta \mid \mathbf{y}, \mathbf{z})

а затем вычислите θ^(j+1)(k)\hat{\theta }_{(j+1)}^{(k)} как значение, максимизирующее Q^(θ∣θ^(j)(k),y)\hat{Q}\left(\theta \mid \hat{\theta }_{(j)}^{(k)}, \mathbf{y}\right).

?
(a)

Покажите, что Q^(θ∣θ^(j)(k),y)→Q^(θ∣θ^(j),y)\hat{Q}\left(\theta \mid \hat{\theta }_{(j)}^{(k)}, \mathbf{y}\right) \rightarrow \hat{Q}\left(\theta \mid \hat{\theta }_{(j)}, \mathbf{y}\right) при k→∞k \rightarrow \infty.

(b)

Какие условия обеспечат, что L(θ^(j+1)(k)∣y)≥L(θ^(j)(k)∣y)L\left(\hat{\theta }_{(j+1)}^{(k)} \mid \mathbf{y}\right) \geq L\left(\hat{\theta }_{(j)}^{(k)} \mid \mathbf{y}\right) при достаточно больших kk? Достаточны ли для этого предположения Теоремы 4.12?

Задача 6.4.12

Для смеси распределений из Примера 4.7, то есть

Xi∼θg(x)+(1−θ)h(x),i=1,…,n, независимые  X_{i} \sim \theta g(x)+(1-\theta ) h(x), \quad i=1, \ldots , n, \text{ независимые }

где g(⋅)g(\cdot ) и h(⋅)h(\cdot ) известны, для нахождения ML-оценки θ\theta можно использовать EM-алгоритм. Пусть Z1,⋯ ,ZnZ_{1}, \cdots , Z_{n}, где ZiZ_{i} указывает, из какого распределения получена XiX_{i}, так что

Xi∣Zi=1∼g(x)Xi∣Zi=0∼h(x). \begin{align} & X_{i} \mid Z_{i}=1 \sim g(x) \\ & X_{i} \mid Z_{i}=0 \sim h(x). \end{align}
?
(a)

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

L(θ∣x,z)=∏i=1n[zig(xi)+(1−zi)h(xi)]θzi(1−θ)1−zi. L(\theta \mid \mathbf{x}, \mathbf{z})=\prod _{i=1}^{n}\left[z_{i} g\left(x_{i}\right)+\left(1-z_{i}\right) h\left(x_{i}\right)\right] \theta ^{z_{i}}(1-\theta )^{1-z_{i}}.
(b)

Покажите, что E(Zi∣θ,xi)=θg(xi)/[θg(xi)+(1−θ)h(xi)]E\left(Z_{i} \mid \theta , x_{i}\right)=\theta g\left(x_{i}\right) /\left[\theta g\left(x_{i}\right)+(1-\theta ) h\left(x_{i}\right)\right], и, следовательно, EM-последовательность задаётся формулой

θ^(j+1)=1n∑i=1nθ^(j)g(xi)θ^(j)g(xi)+(1−θ^(j))h(xi). \hat{\theta }_{(j+1)}=\frac{1}{n} \sum _{i=1}^{n} \frac{\hat{\theta }_{(j)} g\left(x_{i}\right)}{\hat{\theta }_{(j)} g\left(x_{i}\right)+\left(1-\hat{\theta }_{(j)}\right) h\left(x_{i}\right)}.
(c)

Покажите, что θ^(j)→θ^\hat{\theta }_{(j)} \rightarrow \hat{\theta } — ML-оценка θ\theta.

Задача 6.4.13

Для ситуации Примера 4.10:

?
(a)

Покажите, что M-шаг EM-алгоритма задаётся формулами

μ^=(∑i=14∑j=1niyij+z1+z2)/12, \hat{\mu }=\left(\sum _{i=1}^{4} \sum _{j=1}^{n_{i}} y_{i j}+z_{1}+z_{2}\right) / 12, α^i=(∑j=12yij+zi)/3−μ^,i=1,3=(∑j=13yij)/3−μ^,i=2,4. \begin{align} \hat{\alpha }_{i} & =\left(\sum _{j=1}^{2} y_{i j}+z_{i}\right) / 3-\hat{\mu }, \quad i=1,3 \\ & =\left(\sum _{j=1}^{3} y_{i j}\right) / 3-\hat{\mu }, \quad i=2,4. \end{align}
(b)

Покажите, что E-шаг EM-алгоритма задаётся формулой

zi=E[Yi3∣μ=μ^,αi=α^i]=μ^+α^ii=1,3. z_{i}=E\left[Y_{i 3} \mid \mu =\hat{\mu }, \alpha _{i}=\hat{\alpha }_{i}\right]=\hat{\mu }+\hat{\alpha }_{i} \quad i=1,3.
(c)

При ограничении ∑iαi=0\sum_{i} \alpha_{i}=0 покажите, что EM-последовательность сходится к α^i=yˉi.−μ^\hat{\alpha }_{i}= \bar{y}_{i.}-\hat{\mu }, где μ^=∑iyˉi./4\hat{\mu }=\sum_{i} \bar{y}_{i.} / 4.

(d)

При ограничении ∑iniαi=0\sum_{i} n_{i} \alpha_{i}=0 покажите, что EM-последовательность сходится к α^i=yˉi.−μ^\hat{\alpha }_{i}= \bar{y}_{i.}-\hat{\mu }, где μ^=∑ijyij/10\hat{\mu }=\sum_{i j} y_{i j} / 10.

(e)

Для общей однофакторной схемы с aa обработками и nijn_{i j} наблюдениями на обработку покажите, как использовать EM-алгоритм для дополнения данных так, чтобы каждая обработка имела nn наблюдений. Запишите EM-последовательность и покажите, к чему она сходится при ограничениях пунктов (c) и (d).

[Ограничения пунктов (c) и (d) встречались в Примере 3.4.9, где они привели, соответственно, к анализу невзвешенных средних и анализу взвешенных средних.]

Задача 6.4.14

В двухфакторной схеме (см. Пример 3.4.11) EM-алгоритм может быть очень полезен при вычислении ML-оценок в несбалансированном случае. Предположим, что мы наблюдаем

Yijk:N(ξij,σ2),i=1,…,I,j=1,…,J,k=1,…,nij, Y_{i j k}: N\left(\xi _{i j}, \sigma ^{2}\right), \quad i=1, \ldots , I, \quad j=1, \ldots , J, \quad k=1, \ldots , n_{i j},

где ξij=μ+αi+βj+γij\xi_{i j}=\mu +\alpha_{i}+\beta_{j}+\gamma_{i j}. Данные будут дополнены так, чтобы полные данные содержали nn наблюдений на ячейку.

?
(a)

Покажите, как вычислить и E-шаг, и M-шаг EM-алгоритма.

(b)

При ограничении ∑iαi=∑jβj=∑iγij=∑jγij=0\sum_{i} \alpha_{i}=\sum_{j} \beta_{j}=\sum_{i} \gamma_{i j}=\sum_{j} \gamma_{i j}=0 покажите, что EM-последовательность сходится к ML-оценкам, соответствующим анализу невзвешенных средних.

(c)

При ограничении ∑ini.αi=∑jn⋅jβj=∑ini⋅γij=∑j⋅jγij=0\sum_{i} n_{i}. \alpha_{i}=\sum_{j} n_{\cdot j} \beta_{j}=\sum_{i} n_{i \cdot } \gamma_{i j}=\sum_{j} \cdot j \gamma_{i j}=0 покажите, что EM-последовательность сходится к ML-оценкам, соответствующим анализу взвешенных средних.

Задача 6.4.15

Для однофакторной схемы со случайными эффектами (Пример 3.5.1) EM-алгоритм полезен для вычисления ML-оценок. (На самом деле, он весьма полезен во многих смешанных моделях; см. Searle и др. 1992, Глава 8.) Предположим, что имеет место модель

Xij=μ+Ai+Uij(j=1,…,ni,i=1,…,s) X_{i j}=\mu +A_{i}+U_{i j} \quad \left(j=1, \ldots , n_{i}, i=1, \ldots , s\right)

где AiA_{i} и UijU_{i j} — независимые нормальные случайные величины с нулевым средним и известной дисперсией. Для вычисления ML-оценок μ,σU2\mu , \sigma_{U}^{2} и σU2\sigma_{U}^{2} обычно применяют EM-алгоритм, используя ненаблюдаемые AiA_{i} в качестве дополненных данных. Выпишите E-шаг и M-шаг и покажите, что EM-последовательность сходится к ML-оценкам.

?
Задача 6.4.16

Оценивание максимального правдоподобия в пробит-модели из Раздела 3.6 можно реализовать с помощью EM-алгоритма. Мы наблюдаем независимые бернуллиевские величины X1,…,XnX_{1}, \ldots , X_{n}, которые зависят от ненаблюдаемых величин ZiZ_{i}, независимо распределённых как N(ζi,σ2)N\left(\zeta_{i}, \sigma^{2}\right), где

Xi={0 если Zi≤u1 если Zi>u. X_{i}= \begin{cases} 0 & \text{ если } Z_{i} \leq u \\ 1 & \text{ если } Z_{i}>u.\end{cases}

Предполагая, что uu известно, мы хотим получить ML-оценки ζ\zeta и σ2\sigma^{2}.

?
(a)

Покажите, что функция правдоподобия равна pΣxi(1−p)n−Σxip^{\Sigma x_{i}}(1-p)^{n-\Sigma x_{i}}, где

p=P(Zi>u)=Φ(ζ−uσ). p=P\left(Z_{i}>u\right)=\Phi \left(\frac{\zeta -u}{\sigma }\right).
(b)

Если считать Z1,…,ZnZ_{1}, \ldots , Z_{n} полными данными, то правдоподобие полных данных равно

∏i=1n12πσe−12σ2(zi−ζ)2 \prod _{i=1}^{n} \frac{1}{\sqrt{2 \pi } \sigma } e^{-\frac{1}{2 \sigma ^{2}}\left(z_{i}-\zeta \right)^{2}}

а ожидаемый логарифм правдоподобия полных данных равен

−n2log⁡(2πσ2)−12σ2∑i=1n[E(Zi2∣Xi)−2ζE(Zi∣Xi)+ζ2]. -\frac{n}{2} \log \left(2 \pi \sigma ^{2}\right)-\frac{1}{2 \sigma ^{2}} \sum _{i=1}^{n}\left[E\left(Z_{i}^{2} \mid X_{i}\right)-2 \zeta E\left(Z_{i} \mid X_{i}\right)+\zeta ^{2}\right].
(c)

Покажите, что EM-последовательность задаётся формулами

ζ^(j+1)=1n∑i=1nti(ζ^(j),σ^(j)2)σ^(j+1)2=1n[∑i=1nvi(ζ^(j),σ^(j)2)−1n(∑i=1nti(ζ^(j),σ^(j)2))2] \begin{align} & \hat{\zeta }_{(j+1)}=\frac{1}{n} \sum _{i=1}^{n} t_{i}\left(\hat{\zeta }_{(j)}, \hat{\sigma }_{(j)}^{2}\right) \\ & \hat{\sigma }_{(j+1)}^{2}=\frac{1}{n}\left[\sum _{i=1}^{n} v_{i}\left(\hat{\zeta }_{(j)}, \hat{\sigma }_{(j)}^{2}\right)-\frac{1}{n}\left(\sum _{i=1}^{n} t_{i}\left(\hat{\zeta }_{(j)}, \hat{\sigma }_{(j)}^{2}\right)\right)^{2}\right] \end{align}

где

ti(ζ,σ2)=E(Zi∣Xi,ζ,σ2) и vi(ζ,σ2)=E(Zi2∣Xi,ζ,σ2). t_{i}\left(\zeta , \sigma ^{2}\right)=E\left(Z_{i} \mid X_{i}, \zeta , \sigma ^{2}\right) \quad \text{ и } \quad v_{i}\left(\zeta , \sigma ^{2}\right)=E\left(Z_{i}^{2} \mid X_{i}, \zeta , \sigma ^{2}\right).
(d)

Покажите, что

E(Zi∣Xi,ζ,σ2)=ζ+σHi(u−ζσ)E(Zi2∣Xi,ζ,σ2)=ζ2+σ2+σ(u+ζ)Hi(u−ζσ) \begin{align} & E\left(Z_{i} \mid X_{i}, \zeta , \sigma ^{2}\right)=\zeta +\sigma H_{i}\left(\frac{u-\zeta }{\sigma }\right) \\ & E\left(Z_{i}^{2} \mid X_{i}, \zeta , \sigma ^{2}\right)=\zeta ^{2}+\sigma ^{2}+\sigma (u+\zeta ) H_{i}\left(\frac{u-\zeta }{\sigma }\right) \end{align}

где

Hi(t)={φ(t)1−Φ(t) если Xi=1−φ(t)Φ(t) если Xi=0 H_{i}(t)= \begin{cases} \frac{\varphi (t)}{1-\Phi (t)} & \text{ если } X_{i}=1 \\ -\frac{\varphi (t)}{\Phi (t)} & \text{ если } X_{i}=0\end{cases}
(e)

Покажите, что ζ^(j)→ζ^\hat{\zeta }_{(j)} \rightarrow \hat{\zeta } и σ^(j)2→σ^2\hat{\sigma }_{(j)}^{2} \rightarrow \hat{\sigma }^{2} — ML-оценки ζ\zeta и σ2\sigma^{2}.

Задача 6.4.17

Проверьте (4.30).

?
Задача 6.4.18

EM-алгоритм также можно применить в байесовской иерархической модели для нахождения моды апостериорного распределения. Вспомним модель (4.5.5.1),

X∣θ∼f(x∣θ),Θ∣λ∼π(θ∣λ),Λ∼γ(λ), \begin{align} & X \mid \theta \sim f(x \mid \theta ), \\ & \Theta \mid \lambda \sim \pi (\theta \mid \lambda ), \\ & \Lambda \sim \gamma (\lambda ), \end{align}

где интерес представляет оценивание величин, связанных с π(θ∣x)\pi (\theta \mid x). Поскольку

π(θ∣x)=∫π(θ,λ∣x)dλ, \pi (\theta \mid x)=\int \pi (\theta , \lambda \mid x) d \lambda ,

где π(θ,λ∣x)=π(θ∣λ,x)π(λ∣x)\pi (\theta , \lambda \mid x)=\pi (\theta \mid \lambda , x) \pi (\lambda \mid x), EM-алгоритм является одним из возможных методов нахождения моды π(θ∣x)\pi (\theta \mid x), где λ\lambda используется в качестве дополненных данных.

?
(a)

Определите k(λ∣θ,x)=π(θ,λ∣x)/π(θ∣x)k(\lambda \mid \theta , x)=\pi (\theta , \lambda \mid x) / \pi (\theta \mid x) и покажите, что

log⁡π(θ∣x)=∫log⁡π(θ,λ∣x)k(λ∣θ∗,x)dλ−∫log⁡k(λ∣θ,x)k(λ∣θ∗,x)dλ. \log \pi (\theta \mid x)=\int \log \pi (\theta , \lambda \mid x) k\left(\lambda \mid \theta ^{*}, x\right) d \lambda -\int \log k(\lambda \mid \theta , x) k\left(\lambda \mid \theta ^{*}, x\right) d \lambda .
(b)

Если последовательность {θ^(j)}\left\{ \hat{\theta }_{(j)}\right\} удовлетворяет

max⁡θ∫log⁡π(θ,λ∣x)k(λ∣θ(j),x)dλ=∫log⁡π(θ(j+1),λ∣x)k(λ∣θ(j),x)dλ, \max _{\theta } \int \log \pi (\theta , \lambda \mid x) k\left(\lambda \mid \theta _{(j)}, x\right) d \lambda =\int \log \pi \left(\theta _{(j+1)}, \lambda \mid x\right) k\left(\lambda \mid \theta _{(j)}, x\right) d \lambda ,

покажите, что log⁡π(θ(j+1)∣x)≥log⁡π(θ(j)∣x)\log \pi \left(\theta_{(j+1)} \mid x\right) \geq \log \pi \left(\theta_{(j)} \mid x\right). При каких условиях последовательность {θ^(j)}\left\{ \hat{\theta }_{(j)}\right\} будет сходиться к моде π(θ∣x)\pi (\theta \mid x)?

(c)

Для иерархии

X∣θ∼N(θ),1),Θ∣λ∼N(λ,1)),Λ∼Uniform⁡(−∞,∞), \begin{align} & X \mid \theta \sim N(\theta ), 1), \\ & \Theta \mid \lambda \sim N(\lambda , 1)), \\ & \Lambda \sim \operatorname {Uniform}(-\infty , \infty ), \end{align}

покажите, как использовать EM-алгоритм для вычисления моды апостериорного распределения π(θ∣x)\pi (\theta \mid x).

Задача 6.4.19

Между EM-алгоритмом и сэмплированием Гиббса существует связь, поскольку оба метода основаны на теории марковских цепей. Один из способов убедиться в этом — показать, что правдоподобие неполных данных является решением интегрального уравнения выборки методом последовательных подстановок (см. Задачи 4.5.9–4.5.11), и что затем сэмплирование Гиббса можно использовать для вычисления функции правдоподобия. Если L(θ∣y)L(\theta \mid \mathbf{y}) — правдоподобие неполных данных, а L(θ∣y,z)L(\theta \mid \mathbf{y}, \mathbf{z}) — правдоподобие полных данных, определим

L∗(θ∣y)=L(θ∣y)∫L(θ∣y)dθL∗(θ∣y,z)=L(θ∣y,z)∫L(θ∣y,z)dθ \begin{align} L^{*}(\theta \mid \mathbf{y}) & =\frac{L(\theta \mid \mathbf{y})}{\int L(\theta \mid \mathbf{y}) d \theta } \\ L^{*}(\theta \mid \mathbf{y}, \mathbf{z}) & =\frac{L(\theta \mid \mathbf{y}, \mathbf{z})}{\int L(\theta \mid \mathbf{y}, \mathbf{z}) d \theta } \end{align}
?
(a)

Покажите, что L∗(θ∣y)L^{*}(\theta \mid \mathbf{y}) является решением уравнения

L∗(θ∣y)=∫[∫L∗(θ∣y,z)k(z∣θ′,y)dz]L∗(θ′∣y)dθ′ L^{*}(\theta \mid \mathbf{y})=\int \left[\int L^{*}(\theta \mid \mathbf{y}, \mathbf{z}) k\left(\mathbf{z} \mid \theta ^{\prime }, \mathbf{y}\right) d \mathbf{z}\right] L^{*}\left(\theta ^{\prime } \mid \mathbf{y}\right) d \theta ^{\prime }

где, как обычно, k(z∣θ,y)=L(θ∣y,z)/L(θ∣y)k(\mathbf{z} \mid \theta , \mathbf{y})=L(\theta \mid \mathbf{y}, \mathbf{z}) / L(\theta \mid \mathbf{y}).

(b)

Покажите, как последовательность θ(j)\theta_{(j)}, полученная из итерации Гиббса,

θ(j)∼L∗(θ∣y,z(j−1))z(j)∼k(z∣θ(j),y) \begin{align} & \theta _{(j)} \sim L^{*}\left(\theta \mid \mathbf{y}, \mathbf{z}_{(j-1)}\right) \\ & \mathbf{z}_{(j)} \sim k\left(\mathbf{z} \mid \theta _{(j)}, \mathbf{y}\right) \end{align}

будет сходиться к случайной величине с плотностью L∗(θ∣y)L^{*}(\theta \mid \mathbf{y}) при j→∞j \rightarrow \infty. Как это можно использовать для вычисления функции правдоподобия L(θ∣y)L(\theta \mid \mathbf{y})?

[Используя функции L(θ∣y,z)L(\theta \mid \mathbf{y}, \mathbf{z}) и k(z∣θ,y)k(\mathbf{z} \mid \theta , \mathbf{y}), EM-алгоритм даёт нам ML-оценку из L(θ∣y)L(\theta \mid \mathbf{y}), тогда как сэмплер Гиббса даёт нам всю функцию целиком. Эта реализация сэмплера Гиббса для вычисления правдоподобия использовалась Casella и Berger (1994), а также описана у Smith и Roberts (1993). Версия EM-алгоритма, в которой связь с марковскими цепями вполне очевидна, была предложена Baum и Petrie (1966) и Baum и др. (1970).]

§
Задача 6.5.1
?
(a)

Если вектор Yn\mathbf{Y}_{n} в EsE_{s} сходится по вероятности к постоянному вектору a\mathbf{a}, а hh — непрерывная функция, определённая на EsE_{s}, покажите, что h(Yn)→h(a)h\left(\mathbf{Y}_{n}\right) \rightarrow h(\mathbf{a}) по вероятности.

(b)

Используя (a), покажите, что элементы ∥Ajkn∥−1\left\| A_{j k n}\right\|^{-1} сходятся по вероятности к элементам BB, как утверждается в доказательстве Леммы 5.2.

Задача 6.5.2
?
(a)

Покажите, что (5.26) без учёта остаточного члена имеет тот же вид, что и (5.15), и определите AjknA_{j k n}.

(b)

Покажите, что получаемые ajka_{j k} из Леммы 5.2 совпадают с теми, что в (5.23).

(c)

Покажите, что остаточным членом в (5.26) можно пренебречь в доказательстве Теоремы 5.3.

Задача 6.5.3

Пусть X1,…,XnX_{1}, \ldots , X_{n} — независимые одинаково распределённые случайные величины с распределением N(ξ,σ2)N\left(\xi , \sigma^{2}\right).

?
(a)

Покажите, что уравнения правдоподобия имеют единственный корень.

(b)

Покажите непосредственно (т.е. не прибегая к Теореме 5.1), что MLE ξ^\hat{\xi } и σ^\hat{\sigma } асимптотически эффективны.

Задача 6.5.4

Пусть (X0,…,Xs)\left(X_{0}, \ldots , X_{s}\right) имеет полиномиальное распределение M(p0,…,ps;n)M\left(p_{0}, \ldots , p_{s} ; n\right).

?
(a)

Покажите, что уравнения правдоподобия имеют единственный корень.

(b)

Покажите непосредственно, что MLE p^i\hat{p}_{i} асимптотически эффективны.

Задача 6.5.5

Теорема 5.3 утверждает: предположим, что выполнены предположения Теоремы 5.1 (стандартные условия многопараметрической регулярности правдоподобия и состоятельности для независимой одинаково распределённой выборки из регулярного параметрического семейства с параметром θ=(θ1,…,θs)\theta =(\theta_1,\ldots ,\theta_s)), и что θ~jn\tilde\theta_{jn} — n\sqrt{n}-состоятельная оценка θj\theta_j при j=1,…,sj=1,\ldots ,s. Пусть {δkn,k=1,…,s}\left\{ \delta_{kn}, k=1,\ldots ,s\right\} — решение уравнений

∑k=1s(δkn−θ~kn)ljk′′(θ~n)=−lj′(θ~n),j=1,…,s, \sum _{k=1}^{s}\left(\delta _{k n}-\tilde{\theta }_{k n}\right) l_{j k}''\left(\tilde{\theta }_{n}\right)=-l_{j}'\left(\tilde{\theta }_{n}\right), \qquad j=1,\ldots ,s,

где l(θ)l(\theta ) — логарифм правдоподобия, lj′=∂l/∂θjl_j' = \partial l/\partial \theta_j, а ljk′′=∂2l/∂θj∂θkl_{jk}'' = \partial^2 l /\partial \theta_j\partial \theta_k. Тогда δn=(δ1n,…,δsn)\delta_n=(\delta_{1n},\ldots ,\delta_{sn}) асимптотически эффективна.

Следствие 5.4 утверждает: предположим, что выполнены предположения Теоремы 5.3 и что элементы Ijk(θ)I_{jk}(\theta ) информационной матрицы величин XX непрерывны. Тогда решения δkn′\delta_{kn}' уравнений

n(δkn′−θ~kn)Ijk(θ~n)=lj′(θ~n) n\left(\delta _{k n}'-\tilde{\theta }_{k n}\right) I_{j k}\left(\tilde{\theta }_{n}\right)=l_{j}'\left(\tilde{\theta }_{n}\right)

асимптотически эффективны.

Докажите Следствие 5.4.

?
Задача 6.5.6

Покажите, что существует функция ff двух переменных, для которой уравнения ∂f(x,y)/∂x=\partial f(x, y) / \partial x= 0 и ∂f(x,y)/∂y=0\partial f(x, y) / \partial y=0 имеют единственное решение, причём это решение является локальным, но не глобальным максимумом ff.

?
§
Задача 6.6.1

В Примере 6.1 покажите, что уравнения правдоподобия задаются формулами (6.2) и (6.3).

?
Задача 6.6.2

В Примере 6.1 проверьте уравнение (6.4).

?
Задача 6.6.3

Проверьте (6.5).

?
Задача 6.6.4

Если θ=(θ1,…,θr,θr+1,…,θs)\boldsymbol {\theta }=\left(\theta_{1}, \ldots , \theta_{r}, \theta_{r+1}, \ldots , \theta_{s}\right) и если

cov⁡[∂∂θiL(θ),∂∂θjL(θ)]=0 для любых i≤r<j, \operatorname {cov}\left[\frac{\partial }{\partial \theta _{i}} L(\boldsymbol {\theta }), \frac{\partial }{\partial \theta _{j}} L(\boldsymbol {\theta })\right]=0 \quad \text{ для любых } i \leq r<j,

то асимптотическое распределение (θ^1,…,θ^r)\left(\hat{\theta }_{1}, \ldots , \hat{\theta }_{r}\right) при предположениях Теоремы 5.1 не зависит от того, известны ли θr+1,…,θs\theta_{r+1}, \ldots , \theta_{s}.

?
Задача 6.6.5

Пусть X1,…,XnX_{1}, \ldots , X_{n} — независимые одинаково распределённые случайные величины с гамма-распределением Γ(α,β)\Gamma (\alpha , \beta ) и плотностью 1/(Γ(α)βα)×xα−1e−x/β1 /\left(\Gamma (\alpha ) \beta^{\alpha }\right) \times x^{\alpha -1} e^{-x / \beta }.

?
(a)

Вычислите информационную матрицу для обычной параметризации (α,β)(\alpha , \beta ).

(b)

Запишите плотность через параметры (α,μ)=(α,α/β)(\alpha , \mu )=(\alpha , \alpha / \beta ). Вычислите информационную матрицу для параметризации (α,μ)(\alpha , \mu ) и покажите, что она диагональна, а значит, параметры ортогональны.

(c)

Если MLE в пункте (a) равны (α^,β^)(\hat{\alpha }, \hat{\beta }), покажите, что (α^,μ^)=(α^,α^/β^)(\hat{\alpha }, \hat{\mu })=(\hat{\alpha }, \hat{\alpha } / \hat{\beta }). Таким образом, обе модели одинаково хорошо оценивают среднее.

(Теорию и другие примеры ортогональности параметров см. в Cox и Reid 1987.)

Задача 6.6.6

В Примере 6.4 проверьте MLE ξ^i\hat{\xi }_{i} и σ^jk\hat{\sigma }_{j k}, когда ξ\xi неизвестны.

?
Задача 6.6.7

В Примере 6.4 покажите, что SjkS_{j k}, заданные формулой (6.15), не зависят от (X1,…,Xp)(X_{1}, \ldots , X_{p}) и имеют то же совместное распределение, что и статистики (6.13), в которых nn заменено на n−1n-1.

?
Задача 6.6.8

Проверьте матрицы

?
(a)

(6.17) и

(b)

(6.18).

Задача 6.6.9

Рассмотрим ситуацию, приводящую к (6.20), где (Xi,Yi),i=1,…,n\left(X_{i}, Y_{i}\right), i=1, \ldots , n, независимы и одинаково распределены по двумерному нормальному закону с E(Xi)=E(Yi)=0,var⁡(Xi)=var⁡(Yi)=1E\left(X_{i}\right)=E\left(Y_{i}\right)=0, \operatorname {var}\left(X_{i}\right)=\operatorname {var}\left(Y_{i}\right)=1 и неизвестным коэффициентом корреляции ρ\rho.

?
(a)

Покажите, что уравнение правдоподобия является кубическим, для которого вероятность единственного корня стремится к 1 при n→∞n \rightarrow \infty. [Указание: Для кубического уравнения ax3+3bx2+3cx+d=0a x^{3}+3 b x^{2}+3 c x+d=0 положим G=a2d−3abc+2b3G=a^{2} d-3 a b c+2 b^{3} и H=ac−b2H=a c-b^{2}. Тогда условие единственности вещественного корня есть G2+4H3>0G^{2}+4 H^{3}>0.]

(b)

Покажите, что если ρ^n\hat{\rho }_{n} — состоятельное решение уравнения правдоподобия, то оно удовлетворяет (6.20).

(c)

Покажите, что δ=ΣXiYi/n\delta =\Sigma X_{i} Y_{i} / n — состоятельная оценка ρ\rho и что n(δ−ρ)→LN(0,1+ρ2)\sqrt{n}(\delta -\rho ) \xrightarrow {\mathcal{L}} N\left(0,1+\rho^{2}\right), а значит, δ\delta менее эффективна, чем MLE для ρ\rho.

Задача 6.6.10

Проверьте предельное распределение, утверждаемое в (6.21).

?
Задача 6.6.11

Пусть X,…,XnX, \ldots , X_{n} — независимые одинаково распределённые по распределению Пуассона P(λ)P(\lambda ). Найдите ARE δ2n=[\delta_{2 n}=\left[\right. Число Xi=0]/n\left.X_{i}=0\right] / n по отношению к δ1n=e−Xˉn\delta_{1 n}=e^{-\bar{X}_{n}} как оценок e−λe^{-\lambda }.

?
Задача 6.6.12

Покажите, что эффективность (6.27) стремится к 0 при ∣a−θ∣→∞\left|a-\theta \right| \rightarrow \infty.

?
Задача 6.6.13

Для ситуации Примера 6.9 рассмотрим в качестве другого семейства распределений семейство загрязнённых нормальных смесей, предложенное Tukey (1960) в качестве модели для наблюдений, которые обычно подчиняются нормальному распределению, но иногда в эксперименте или при регистрации его результатов что-то идёт не так, так что полученное наблюдение оказывается грубой ошибкой. В модели Tukey функция распределения имеет вид

Fτ,ϵ(t)=(1−ϵ)Φ(t)+ϵΦ(tτ). F_{\tau , \epsilon }(t)=(1-\epsilon ) \Phi (t)+\epsilon \Phi \left(\frac{t}{\tau }\right).

То есть в случаях грубой ошибки предполагается, что наблюдения нормально распределены с тем же средним θ\theta, но с другой (большей) дисперсией τ2.9\tau^{2}.{ }^{9}

?
(a)

Покажите, что если XiX_{i} имеют распределение Fτ,ϵ(x−θ)F_{\tau , \epsilon }(x-\theta ), предельное распределение δ2n\delta_{2 n} не меняется.

(b)

Покажите, что предельное распределение δ1n\delta_{1 n} нормально со средним нуль и дисперсией nn−1{ϕ[nn−1(a−θ)]}2(1−ϵ+ϵτ2)\frac{n}{n-1}\left\{ \phi \left[\sqrt{\frac{n}{n-1}}(a-\theta )\right]\right\}^{2}\left(1-\epsilon +\epsilon \tau^{2}\right).

(c)

Сравните асимптотическую относительную эффективность δ1n\delta_{1 n} и δ2n\delta_{2 n}.

Задача 6.6.14

Пусть X1,…,XnX_{1}, \ldots , X_{n} независимы и одинаково распределены как N(0,σ2)N\left(0, \sigma^{2}\right).

?
(a)

Покажите, что δn=kΣ∣Xi∣/n\delta_{n}=k \Sigma \left|X_{i}\right| / n является состоятельной оценкой σ\sigma тогда и только тогда, когда k=π/2k=\sqrt{\pi / 2}.

(b)

Определите ARE оценки δ\delta при k=π/2k=\sqrt{\pi / 2} по отношению к MLE ΣXi2/n\sqrt{\Sigma X_{i}^{2} / n}.

Задача 6.6.15

Пусть X1,…,XnX_{1}, \ldots , X_{n} независимы и одинаково распределены с E(Xi)=θ,var⁡(Xi)=1E\left(X_{i}\right)=\theta , \operatorname {var}\left(X_{i}\right)=1 и E(Xi−θ)4=μ4E\left(X_{i}-\theta \right)^{4}=\mu_{4}, и рассмотрим несмещённые оценки δ1n=(1/n)ΣXi2−1\delta_{1 n}=(1 / n) \Sigma X_{i}^{2}-1 и δ2n=Xˉn2−1/n\delta_{2 n}=\bar{X}_{n}^{2}-1 / n параметра θ2\theta^{2}.

?
(a)

Определите ARE e2,1e_{2,1} оценки δ2n\delta_{2 n} по отношению к δ1n\delta_{1 n}.

(b)

Покажите, что e2,1≥1e_{2,1} \geq 1, если XiX_{i} симметричны относительно θ\theta.

(c)

Найдите распределение XiX_{i}, для которого e2,1<1e_{2,1}<1.

Задача 6.6.16

Свойство асимптотической относительной эффективности было определено (Определение 6.6) для оценок, сходящихся к нормальному распределению со скоростью n\sqrt{n}. Это определение, а также Теорема 6.7, можно обобщить, включив другие распределения и скорости сходимости.

9{ }^{9} Как отметил Stigler (1973), подобные модели для распределений с тяжёлыми хвостами уже предлагались значительно раньше в забытой работе Newcomb (1882, 1886). Теорема 9.1 Пусть {δin}\left\{ \delta_{i n}\right\} — две последовательности оценок g(θ)g(\theta ), такие что

nα[δin−g(θ)]→LτiT,α>0,τi>0,i=1,2, n^{\alpha }\left[\delta _{i n}-g(\theta )\right] \xrightarrow {\mathcal{L}} \tau _{i} T, \alpha >0, \tau _{i}>0, i=1,2,

где распределение HH величины TT имеет носитель на интервале −∞≤A<B≤∞-\infty \leq A<B \leq \infty со строго возрастающей функцией распределения на (A,B)(A, B). Тогда ARE {δ2n}\left\{ \delta_{2 n}\right\} по отношению к {δ1n}\left\{ \delta_{1 n}\right\} существует и равна

e21=lim⁡n2→∞n1(n2)n2=[τ1τ2]1/α. e_{21}=\lim _{n_{2} \rightarrow \infty } \frac{n_{1}\left(n_{2}\right)}{n_{2}}=\left[\frac{\tau _{1}}{\tau _{2}}\right]^{1 / \alpha }.
?
Задача 6.6.17

В Примере 6.10 покажите, что условия Теоремы 5.1 выполнены.

?
§
Задача 6.7.1

Рассмотрим rr независимых выборок с индексом α=1,…,r\alpha =1,\ldots ,r объёмов nαn_\alpha, таких что при стремлении общего объёма выборки N=∑αnα→∞N=\sum_\alpha n_\alpha \to \infty имеем nα/N→λαn_\alpha /N \to \lambda_\alpha, где ∑αλα=1\sum_\alpha \lambda_\alpha =1 и каждое λα>0\lambda_\alpha >0.

Теорема 7.1 утверждает: предположим, что предположения Теоремы 5.1 выполнены для каждой из плотностей fα,θf_{\alpha ,\theta }. Пусть I(α)(θ)I^{(\alpha )}(\theta ) обозначает информационную матрицу, соответствующую fα,θf_{\alpha ,\theta }, и пусть I(θ)=∑αλαI(α)(θ)I(\theta )=\sum_\alpha \lambda_\alpha I^{(\alpha )}(\theta ). Логарифм правдоподобия l(θ)l(\theta ) задаётся формулой l(θ)=∑α=1r∑j=1nαlog⁡fα,θ(xαj)l(\theta )=\sum_{\alpha =1}^{r}\sum_{j=1}^{n_\alpha } \log f_{\alpha ,\theta }(x_{\alpha j}), а уравнения правдоподобия — ∂∂θjl(θ)=0\dfrac {\partial }{\partial \theta_j} l(\theta )=0 (j=1,…,s)(j=1,\ldots ,s). При этих обозначениях выводы Теоремы 5.1 остаются в силе.

Докажите Теорему 7.1.

?
Задача 6.7.2

Для ситуации Примера 7.3 при m=nm=n:

?
(a)

Покажите, что необходимым условием сходимости (7.5) к N(0,1)N(0,1) является n(λ^−λ)→\sqrt{n}(\hat{\lambda }-\lambda ) \rightarrow 0, где λ^=σ^2/τ^2\hat{\lambda }=\hat{\sigma }^{2} / \hat{\tau }^{2} и λ=σ2/τ2\lambda =\sigma^{2} / \tau^{2}, для σ^2\hat{\sigma }^{2} и τ^2\hat{\tau }^{2} из (7.4).

(b)

Используя тот факт, что λ^/λ\hat{\lambda } / \lambda имеет FF-распределение, покажите, что n(λ^−λ)↛0\sqrt{n}(\hat{\lambda }-\lambda ) \nrightarrow 0.

(c)

Покажите, что полная MLE задаётся решением

ξ=(m/σ2)Xˉ+(n/τ2)Yˉm/σ2+n/τ2,σ2=1mΣ(Xi−ξ)2,τ2=1nΣ(Yj−ξ)2, \xi =\frac{\left(m / \sigma ^{2}\right) \bar{X}+\left(n / \tau ^{2}\right) \bar{Y}}{m / \sigma ^{2}+n / \tau ^{2}}, \quad \sigma ^{2}=\frac{1}{m} \Sigma \left(X_{i}-\xi \right)^{2}, \quad \tau ^{2}=\frac{1}{n} \Sigma \left(Y_{j}-\xi \right)^{2},

и выведите её асимптотическую эффективность из Теоремы 5.1.

Задача 6.7.3

В Примере 7.4 определите совместное распределение

?
(a)

(σ^2,τ^2)(\hat{\sigma }^{2}, \hat{\tau }^{2}) и

(b)

(σ^2,σ^A2)(\hat{\sigma }^{2}, \hat{\sigma }_{A}^{2}).

Задача 6.7.4

Рассмотрим выборки (X1,Y1),…,(Xm,Ym)\left(X_{1}, Y_{1}\right), \ldots ,\left(X_{m}, Y_{m}\right) и (X1′,Y1′),…,(Xn′,Yn′)\left(X_{1}^{\prime }, Y_{1}^{\prime }\right), \ldots ,\left(X_{n}^{\prime }, Y_{n}^{\prime }\right) из двух двумерных нормальных распределений с нулевыми средними и дисперсиями-ковариациями (σ2,τ2,ρστ)\left(\sigma^{2}, \tau^{2}, \rho \sigma \tau \right) и (σ′2,τ′2,ρ′σ′τ′)\left(\sigma^{\prime 2}, \tau^{\prime 2}, \rho^{\prime } \sigma^{\prime } \tau^{\prime }\right) соответственно. Используя Теорему 7.1 и Примеры 6.5 и 6.8, найдите предельное распределение

?
(a)

σ^2\hat{\sigma }^{2} и τ^2\hat{\tau }^{2}, когда известно, что ρ′=ρ\rho^{\prime }=\rho

(b)

ρ^\hat{\rho }, когда известно, что σ′=σ\sigma^{\prime }=\sigma и τ′=τ\tau^{\prime }=\tau.

Задача 6.7.5

В предыдущей задаче найдите выигрыш в эффективности (если он есть)

?
(a)

в пункте (a), обусловленный знанием того, что ρ′=ρ\rho^{\prime }=\rho

(b)

в пункте (b), обусловленный знанием того, что σ′=σ\sigma^{\prime }=\sigma и τ′=τ\tau^{\prime }=\tau.

Задача 6.7.6

Покажите, что уравнения правдоподобия (7.11) имеют не более одного решения.

?
Задача 6.7.7

В Примере 7.6 предположим, что pi=1−F(α+βti)p_{i}=1-F\left(\alpha +\beta t_{i}\right) и что log⁡F(x)\log F(x) и log⁡[1−F(x)]\log [1- F(x)] строго вогнуты. Тогда уравнения правдоподобия имеют не более одного решения.

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

Если функция распределения FF симметрична и log⁡F(x)\log F(x) строго вогнута, то и log⁡[1−F(x)]\log [1-F(x)] строго вогнута.

(b)

Покажите, что log⁡F(x)\log F(x) строго вогнута, когда FF сильно унимодальна, но не когда FF — распределение Коши.

Задача 6.7.9

В Примере 7.7 покажите, что YnY_{n} менее информативна, чем YY.

?
Задача 6.7.10

Покажите, что оценка δn\delta_{n} из Примера 7.7 удовлетворяет (7.14).

?
Задача 6.7.11

Найдите подходящие нормирующие константы для δn\delta_{n} из Примера 7.7, когда

?
(a)

γi=i\gamma_{i}=i,

(b)

γi=i2\gamma_{i}=i^{2}, и

(c)

γi=1/i\gamma_{i}=1 / i.

Задача 6.7.12

Пусть Xi(i=1,…,n)X_{i}(i=1, \ldots , n) независимы, нормальны, с дисперсией 1 и средним βti\beta t_{i} (где tit_{i} известны). Обсудите оценивание β\beta по аналогии с Примером 7.7.

?
Задача 6.7.13

Обобщите предыдущую задачу на ситуацию, в которой

?
(a)

E(Xi)=α+βtiE\left(X_{i}\right)=\alpha +\beta t_{i} и var⁡(Xi)=1\operatorname {var}\left(X_{i}\right)=1, и

(b)

E(Xi)=α+βtiE\left(X_{i}\right)=\alpha +\beta t_{i} и var⁡(Xi)=σ2\operatorname {var}\left(X_{i}\right)=\sigma^{2}, где α,β\alpha , \beta и σ2\sigma^{2} — неизвестные параметры, подлежащие оцениванию.

Задача 6.7.14

Пусть Xj(j=1,…,n)X_{j}(j=1, \ldots , n) независимо распределены с плотностями fj(xj∣θ)f_{j}\left(x_{j} \mid \theta \right) (θ\theta вещественнозначный), пусть Ij(θ)I_{j}(\theta ) — информация о θ\theta, содержащаяся в XjX_{j}, и пусть Tn(θ)=∑j=1nIj(θ)T_{n}(\theta )=\sum_{j=1}^{n} I_{j}(\theta ) — полная информация о θ\theta в выборке. Предположим, что θ^n\hat{\theta }_{n} — состоятельный корень уравнения правдоподобия L′(θ)=0L^{\prime }(\theta )=0 и что, обобщая (3.18)–(3.20),

1Tn(θ0)L′(θ0)→LN(0,1) \frac{1}{\sqrt{T_{n}\left(\theta _{0}\right)}} L^{\prime }\left(\theta _{0}\right) \xrightarrow {\mathcal{L}} N(0,1)

и

−L′′(θ0)Tn(θ0)→P1 и L′′′(θn∗)Tn(θ0) ограничено по вероятности.  -\frac{L^{\prime \prime }\left(\theta _{0}\right)}{T_{n}\left(\theta _{0}\right)} \xrightarrow {P} 1 \quad \text{ и } \quad \frac{L^{\prime \prime \prime }\left(\theta _{n}^{*}\right)}{T_{n}\left(\theta _{0}\right)} \quad \text{ ограничено по вероятности. }

Покажите, что

Tn(θ0)(θ^n−θ0)→LN(0,1). \sqrt{T_{n}\left(\theta _{0}\right)}\left(\hat{\theta }_{n}-\theta _{0}\right) \xrightarrow {\mathcal{L}} N(0,1).
?
Задача 6.7.15

Докажите, что последовательность X1,X2,…X_{1}, X_{2}, \ldots из Примера 7.8 стационарна при условии, что она удовлетворяет (7.17).

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

В Примере 7.8 покажите, что уравнение правдоподобия имеет единственное решение, что оно является MLE и что оно имеет то же асимптотическое распределение, что и δn′=∑i=1nXiXi+1/∑i=1nXi2\delta_{n}^{\prime }=\sum_{i=1}^{n} X_{i} X_{i+1} / \sum_{i=1}^{n} X_{i}^{2}.

(b)

Покажите непосредственно, что δn′\delta_{n}^{\prime } — состоятельная оценка β\beta.

Задача 6.7.17

В Примере 7.8:

?
(a)

Покажите, что при j>1j>1 математическое ожидание условной информации (при заданном Xj−1X_{j-1}), которую XjX_{j} содержит о β\beta, равно 1/(1−β2)1 /\left(1-\beta^{2}\right).

(b)

Определите информацию о β\beta, содержащуюся в X1X_{1}.

Задача 6.7.18

Покажите, что при τ=σ\tau =\sigma в (7.21) MLE существует и состоятельна.

?
Задача 6.7.19

Предположим, что в (7.21) сами ξ\xi являются случайными величинами, независимо одинаково распределёнными как N(μ,γ2)N\left(\mu , \gamma^{2}\right).

?
(a)

Покажите, что совместная плотность (Xi,Yi)\left(X_{i}, Y_{i}\right) является плотностью выборки из двумерного нормального распределения, и определите параметры этого распределения.

(b)

В модели пункта (a) найдите асимптотически эффективные оценки параметров μ,γ,β,σ\mu , \gamma , \beta , \sigma и τ\tau.

Задача 6.7.20

Проверьте корни (7.22).

?
Задача 6.7.21

Покажите, что правдоподобие (7.21) неограничено.

?
Задача 6.7.22

Покажите, что если ρ\rho определяется формулой (7.24), то ρ\rho и ρ′\rho^{\prime } всюду непрерывны.

?
Задача 6.7.23

Пусть FF имеет дифференцируемую плотность ff, и пусть ∫ψ2f<∞\int \psi^{2} f<\infty.

?
(a)

Используя интегрирование по частям, запишите знаменатель (7.27) в виде [∫ψ(x)f′(x)dx]2\left[\int \psi (x) f^{\prime }(x) d x\right]^{2}.

(b)

Покажите, что σ2(F,ψ)≥[∫(f′/f)2f]−1=If−1\sigma^{2}(F, \psi ) \geq \left[\int \left(f^{\prime } / f\right)^{2} f\right]^{-1}=I_{f}^{-1}, применив неравенство Шварца к пункту (a).

Следующие три задачи посвящены исследованию технических условий, необходимых для состоятельности и асимптотической нормальности M-оценок, как отмечено в (7.26).

Задача 6.7.24

Для состоятельности MM-оценок достаточным условием является единственность и изолированность корня оценивающей функции. Установите следующую теорему.

Теорема 9.2 Предположим, что выполнены условия (A0)–(A3). Пусть t0t_{0} — изолированный корень уравнения Eθ0[ψ(X,t)]=0E_{\theta_{0}}[\psi (X, t)]=0, где ψ(⋅,t)\psi (\cdot , t) монотонна по tt и непрерывна в окрестности t0t_{0}. Если T0(x)T_{0}(\mathbf{x}) — решение уравнения ∑i=1nψ(xi,t)=0\sum_{i=1}^{n} \psi \left(x_{i}, t\right)=0, то T0T_{0} сходится к t0t_{0} по вероятности.

?
Задача 6.7.25

Теорема 9.3 При условиях Теоремы 9.2, если, кроме того

  1. Eθ0[∂∂tψ(X,t)∣t=t0]E_{\theta_{0}}\left[\left.\frac{\partial }{\partial t} \psi (X, t)\right|_{t=t_{0}}\right] конечно и отлично от нуля,

  2. Eθ0[ψ2(X,t0)]<∞E_{\theta_{0}}\left[\psi^{2}\left(X, t_{0}\right)\right]<\infty,

то

n(T0−t0)→LN(0,σT02), \sqrt{n}\left(T_{0}-t_{0}\right) \xrightarrow {\mathcal{L}} N\left(0, \sigma _{T_{0}}^{2}\right),

где σT02=Eθ0[ψ2(X,t0)]/(Eθ0[∂∂tψ(X,t)∣t=t0])2\sigma_{T_{0}}^{2}=E_{\theta_{0}}\left[\psi^{2}\left(X, t_{0}\right)\right] /\left(E_{\theta_{0}}\left[\left.\frac{\partial }{\partial t} \psi (X, t)\right|_{t=t_{0}}\right]\right)^{2}. [Заметим, что это небольшое обобщение (7.27).]

?
Задача 6.7.26

Для каждой из следующих оценок выпишите определяющую её функцию ψ\psi и покажите, что оценка состоятельна и асимптотически нормальна при условиях Теорем 9.2 и 9.3.

?
(a)

Оценка наименьших квадратов, минимизирующая ∑(xi−t)2\sum \left(x_{i}-t\right)^{2}.

(b)

Оценка наименьших абсолютных отклонений, минимизирующая ∑∣xi−t∣\sum \left|x_{i}-t\right|.

(c)

Усечённое среднее Хьюбера, минимизирующее (7.24).

Задача 6.7.27

В Примере 7.12 сравните

?
(a)

асимптотические распределения ξ^\hat{\xi } и δn\delta_{n};

(b)

нормированную ожидаемую квадратичную ошибку ξ^\hat{\xi } и δn\delta_{n}.

Задача 6.7.28

В Примере 7.12 покажите, что

?
(a)

n(b^^−b)→LN(0,b2)\sqrt{n}(\hat{\hat{b}}-b) \xrightarrow {\mathcal{L}} N\left(0, b^{2}\right) и

(b)

n(b^−b)→LN(0,b2)\sqrt{n}(\hat{b}-b) \xrightarrow {\mathcal{L}} N\left(0, b^{2}\right).

Задача 6.7.29

В Примере 7.13 покажите, что

?
(a)

c^\hat{c} и a^\hat{a} независимы и имеют указанные распределения;

(b)

X(1)X_{(1)} и Σlog⁡[Xi/X(1)]\Sigma \log \left[X_{i} / X_{(1)}\right] — полные достаточные статистики на основе выборки из (7.33).

Задача 6.7.30

В Примере 7.13 определите UMVU-оценки aa и cc и асимптотические распределения этих оценок.

?
Задача 6.7.31

В предыдущей задаче сравните

?
(a)

асимптотическое распределение MLE и UMVU-оценки cc;

(b)

нормированную ожидаемую квадратичную ошибку этих двух оценок.

Задача 6.7.32

В Примере 7.15,

?
(a)

проверьте уравнение (7.39),

(b)

покажите, что выбор a=−2a=-2 даёт оценку с наилучшей эффективностью второго порядка,

(c)

покажите, что предельное отношение риска MLE⁡(a=0)\operatorname {MLE}(a=0) к δn(a=−2)\delta_{n}(a=-2) равно 2, и

(d)

обсудите поведение этой оценки при малых объёмах выборки.

Задача 6.7.33

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

?
(a)
p∗(x∣ξ)=sup⁡γ,σ2p(x∣ξ,γ,σ2)=c/[σ^(ξ)]nΠ[1/(xi−ξ)] p^{*}(\mathbf{x} \mid \xi )=\sup _{\gamma , \sigma ^{2}} p\left(\mathbf{x} \mid \xi , \gamma , \sigma ^{2}\right)=c /[\hat{\sigma }(\xi )]^{n} \Pi \left[1 /\left(x_{i}-\xi \right)\right]

где

p(x∣ξ,γ,σ2)=∏i=1nf(xi∣ξ,σ2),σ^2(ξ)=1nΣ[log⁡(xi−ξ)−γ^(ξ)]2 и γ^(ξ)=1nΣlog⁡(xi−ξ). \begin{align} p\left(\mathbf{x} \mid \xi , \gamma , \sigma ^{2}\right) & =\prod _{i=1}^{n} f\left(x_{i} \mid \xi , \sigma ^{2}\right), \\ \hat{\sigma }^{2}(\xi ) & =\frac{1}{n} \Sigma \left[\log \left(x_{i}-\xi \right)-\hat{\gamma }(\xi )\right]^{2} \text{ и } \hat{\gamma }(\xi )=\frac{1}{n} \Sigma \log \left(x_{i}-\xi \right). \end{align}
(b)

p∗(x∣ξ)→∞p^{*}(\mathbf{x} \mid \xi ) \rightarrow \infty при ξ→x(1)\xi \rightarrow x_{(1)}. [Указание: (b) Для ξ\xi, достаточно близкого к x(1)x_{(1)},

σ^2(ξ)≤1nΣ[log⁡(xi−ξ)]2≤[log⁡(x(1)−ξ)]2 \hat{\sigma }^{2}(\xi ) \leq \frac{1}{n} \Sigma \left[\log \left(x_{i}-\xi \right)\right]^{2} \leq \left[\log \left(x_{(1)}-\xi \right)\right]^{2}

и, следовательно,

p∗(x∣ξ)≥∣log⁡(x(1)−ξ)∣−nΠ(x(i)−ξ)−1. p^{*}(\mathbf{x} \mid \xi ) \geq \left|\log \left(x_{(1)}-\xi \right)\right|^{-n} \Pi \left(x_{(i)}-\xi \right)^{-1}.

Правая часть стремится к бесконечности при ξ→x(1)\xi \rightarrow x_{(1)} (Hill 1963.]

Задача 6.7.34

Производные всех порядков плотности (7.37) стремятся к нулю при x→ξx \rightarrow \xi.

?
§
Задача 6.8.1

Определите предельное распределение байесовской оценки, соответствующей квадратичной функции потерь, и проверьте, что она асимптотически эффективна, в каждом из следующих случаев:

?
(a)

Наблюдения X1,…,XnX_{1}, \ldots , X_{n} независимы и одинаково распределены как N(θ,σ2)N\left(\theta , \sigma^{2}\right), где σ\sigma известно, а оцениваемый параметр — θ\theta. Априорное распределение для Θ\Theta — сопряжённое нормальное распределение, скажем, N(μ,b2)N\left(\mu , b^{2}\right). (См. Пример 4.2.2.)

(b)

Наблюдения YiY_{i} имеют гамма-распределение Γ(γ,1/τ)\Gamma (\gamma , 1 / \tau ), оцениваемый параметр — 1/τ1 / \tau, а τ\tau имеет сопряжённую априорную плотность Γ(g,α)\Gamma (g, \alpha ).

(c)

Наблюдения и априорное распределение такие же, как в Задаче 4.1.9, а оцениваемый параметр — λ\lambda.

(d)

Наблюдения YiY_{i} имеют отрицательное биномиальное распределение (4.3), pp имеет априорную плотность B(a,b)B(a, b), а оцениваемый параметр — (a) pp и (b) 1/b1 / b.

Задача 6.8.2

Обращаясь к Примеру 8.1, рассмотрим вместо этого минимаксную оценку δn\delta_{n} параметра pp, заданную формулой (1.11), которая соответствует последовательности бета-априорных распределений с a=b=n/2a=b=\sqrt{n} / 2. Тогда,

n[δn−p]=n(Xn−p)+n1+n(12−Xn). \sqrt{n}\left[\delta _{n}-p\right]=\sqrt{n}\left(\frac{X}{n}-p\right)+\frac{\sqrt{n}}{1+\sqrt{n}}\left(\frac{1}{2}-\frac{X}{n}\right).
?
(a)

Покажите, что предельное распределение n[δn−p]\sqrt{n}\left[\delta_{n}-p\right] есть N[12−p,p(1−p)]N\left[\frac{1}{2}-p, p(1-p)\right], так что δn\delta_{n} имеет ту же асимптотическую дисперсию, что и X/nX / n, но при p≠12p \neq \frac{1}{2} является асимптотически смещённой.

(b)

Покажите, что ARE оценки δn\delta_{n} относительно X/nX / n не существует, кроме случая p=12p=\frac{1}{2}, когда она равна 1.

Задача 6.8.3

Условия регулярности Теоремы 2.6 таковы: пусть X1,…,XnX_1, \ldots , X_n независимы и одинаково распределены, каждая с плотностью f(x∣θ)f(x \mid \theta ) относительно σ\sigma-конечной меры μ\mu, где θ\theta вещественнозначный, и предположим:

?
(a)

параметрическое пространство Ω\Omega является открытым интервалом (не обязательно конечным);

(b)

распределения PθP_\theta величин XiX_i имеют общий носитель, так что A={x:f(x∣θ)>0}A=\left\{ x: f(x\mid \theta )>0\right\} не зависит от θ\theta;

(c)

для каждого x∈Ax \in A плотность f(x∣θ)f(x\mid \theta ) дважды дифференцируема по θ\theta, и вторая производная непрерывна по θ\theta;

(d)

интеграл ∫f(x∣θ) dμ(x)\int f(x\mid \theta )\, d\mu (x) можно дважды дифференцировать под знаком интеграла;

(e)

информация Фишера I(θ)I(\theta ) удовлетворяет условию 0<I(θ)<∞0<I(\theta )<\infty;

(f)

для любого заданного θ0∈Ω\theta_0 \in \Omega существуют c>0c>0 и функция M(x)M(x) (возможно, зависящие от θ0\theta_0), такие что ∣∂2log⁡f(x∣θ)/∂θ2∣≤M(x)\left|\partial^{2} \log f(x\mid \theta )/\partial \theta^{2}\right| \leq M(x) для всех x∈Ax \in A, θ0−c<θ<θ0+c\theta_0-c<\theta <\theta_0+c, и Eθ0[M(X)]<∞E_{\theta_0}[M(X)]<\infty.

Если l(θ)=∑i=1nlog⁡f(Xi∣θ)l(\theta )=\sum_{i=1}^{n} \log f(X_i \mid \theta ) — логарифм правдоподобия, а θ0\theta_0 — истинное значение, то (8.1) и (8.2) — это утверждения

l(θ)=l(θ0)+(θ−θ0)l′(θ0)−12(θ−θ0)2[nI(θ0)+Rn(θ)](8.1) l(\theta )=l(\theta _{0})+(\theta -\theta _{0}) l'(\theta _{0})-\frac{1}{2}(\theta -\theta _{0})^{2}\left[n I(\theta _{0})+R_{n}(\theta )\right] \tag {8.1}

где

1nRn(θ)→P0при n→∞.(8.2) \frac{1}{n} R_{n}(\theta ) \xrightarrow {P} 0 \quad \text{при } n \rightarrow \infty . \tag {8.2}

Предположения Теоремы 2.6 влекут (8.1) и (8.2).

Задача 6.8.4

В Примере 8.5 апостериорная плотность θ\theta после одного наблюдения равна f(x1−θ)f\left(x_{1}-\theta \right); это собственная плотность, и она удовлетворяет (B5) при условии Eθ∣X1∣<∞E_{\theta }\left|X_{1}\right|<\infty.

?
Задача 6.8.5

Пусть X1,…,XnX_{1}, \ldots , X_{n} — независимые положительные величины, каждая с плотностью (1/τ)f(xi/τ)(1 / \tau ) f\left(x_{i} / \tau \right), а τ\tau имеет несобственную плотность π(τ)=1/τ(τ>0)\pi (\tau )=1 / \tau (\tau >0). Апостериорная плотность после одного наблюдения является собственной плотностью и удовлетворяет (B5) при условии Eτ(1/X1)<∞E_{\tau }\left(1 / X_{1}\right)<\infty.

?
Задача 6.8.6

Приведите пример, в котором апостериорная плотность является собственной (с вероятностью 1) после двух наблюдений, но не после одного.

?
Задача 6.8.7

Докажите результат, сформулированный перед Примером 8.6.

?
Задача 6.8.8

Пусть X1,…,XnX_{1}, \ldots , X_{n} независимы и одинаково распределены как N(θ,1)N(\theta , 1), и рассмотрим несобственную плотность π(θ)=eθ4\pi (\theta )=e^{\theta^{4}}. Тогда апостериорное распределение будет несобственным при всех nn.

?
Задача 6.8.9

Пусть X1,…,XnX_1, \ldots , X_n независимы и одинаково распределены с плотностью f(xi∣θ)f(x_i \mid \theta ) (относительно μ\mu), где θ\theta вещественнозначный, а параметрическое пространство Ω\Omega — открытый интервал; пусть θ0\theta_0 обозначает истинное значение, l(θ)l(\theta ) — логарифм правдоподобия, а l′(θ0)l'(\theta_0) — score в точке θ0\theta_0. Предположим, что для любой последовательности θ~=θ~n→Pθ0\tilde\theta =\tilde\theta_n \xrightarrow {P} \theta_0,

l(θ)=l(θ0)+(θ−θ0)l′(θ0)−12(θ−θ0)2[nI(θ0)+Rn(θ)], l(\theta )=l(\theta _{0})+(\theta -\theta _{0}) l'(\theta _{0})-\frac{1}{2}(\theta -\theta _{0})^{2}\left[n I(\theta _{0})+R_{n}(\theta )\right],

где 1nRn(θ)→P0\frac{1}{n} R_{n}(\theta ) \xrightarrow {P} 0 при n→∞n \rightarrow \infty. Пусть

Tn=θ0+1nI(θ0)l′(θ0), T_{n}=\theta _{0}+\frac{1}{n I(\theta _{0})} l'(\theta _{0}),

и определим

ω(t)=l(Tn+tn)−l(θ0)−12nI(θ0)[l′(θ0)]2. \omega (t)=l\left(T_{n}+\frac{t}{\sqrt{n}}\right)-l(\theta _{0})-\frac{1}{2 n I(\theta _{0})}\left[l'(\theta _{0})\right]^{2}.

Лемма 8.7 утверждает: величина ω(t)\omega (t) равна

ω(t)=−I(θ0)t22n−12nRn(Tn+tn)[t+1I(θ0)nl′(θ0)]2. \omega (t)=-I(\theta _{0}) \frac{t^{2}}{2 n}-\frac{1}{2 n} R_{n}\left(T_{n}+\frac{t}{\sqrt{n}}\right)\left[t+\frac{1}{I(\theta _{0}) \sqrt{n}} l'(\theta _{0})\right]^{2}.

Докажите Лемму 8.7.

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

Если sup⁡∣Yn(t)∣→P0\sup \left|Y_{n}(t)\right| \xrightarrow {P} 0 и sup⁡∣Xn(t)−c∣→P0\sup \left|X_{n}(t)-c\right| \xrightarrow {P} 0 при n→∞n \rightarrow \infty, то sup⁡∣Xn(t)−ceYn(t)∣→P0\sup \mid X_{n}(t)- c e^{Y_{n}(t)} \mid \xrightarrow {P} 0, где супремум берётся по общему множеству t∈Tt \in T.

(b)

Используя (a), покажите, что (8.22) и (8.23) влекут (8.21).

Задача 6.8.11

Условия регулярности Теоремы 2.6 таковы: пусть X1,…,XnX_1, \ldots , X_n независимы и одинаково распределены, каждая с плотностью f(x∣θ)f(x \mid \theta ) относительно σ\sigma-конечной меры μ\mu, где θ\theta вещественнозначный, и предположим:

?
(a)

параметрическое пространство Ω\Omega является открытым интервалом (не обязательно конечным);

(b)

распределения PθP_\theta величин XiX_i имеют общий носитель, так что A={x:f(x∣θ)>0}A=\left\{ x: f(x\mid \theta )>0\right\} не зависит от θ\theta;

(c)

для каждого x∈Ax \in A плотность f(x∣θ)f(x\mid \theta ) дважды дифференцируема по θ\theta, и вторая производная непрерывна по θ\theta;

(d)

интеграл ∫f(x∣θ) dμ(x)\int f(x\mid \theta )\, d\mu (x) можно дважды дифференцировать под знаком интеграла;

(e)

информация Фишера I(θ)I(\theta ) удовлетворяет условию 0<I(θ)<∞0<I(\theta )<\infty;

(f)

для любого заданного θ0∈Ω\theta_0 \in \Omega существуют c>0c>0 и функция M(x)M(x) (возможно, зависящие от θ0\theta_0), такие что ∣∂2log⁡f(x∣θ)/∂θ2∣≤M(x)\left|\partial^{2} \log f(x\mid \theta )/\partial \theta^{2}\right| \leq M(x) для всех x∈Ax \in A, θ0−c<θ<θ0+c\theta_0-c<\theta <\theta_0+c, и Eθ0[M(X)]<∞E_{\theta_0}[M(X)]<\infty.

Предположение (B1) таково: логарифм правдоподобия l(θ)=∑i=1nlog⁡f(Xi∣θ)l(\theta )=\sum_{i=1}^{n}\log f(X_i\mid \theta ) удовлетворяет этим условиям.

Пусть θ0\theta_0 обозначает истинное значение, l′(θ0)l'(\theta_0) — score в точке θ0\theta_0, I(θ0)I(\theta_0) — информацию Фишера, и

Tn=θ0+1nI(θ0)l′(θ0). T_{n}=\theta _{0}+\frac{1}{n I(\theta _{0})} l'(\theta _{0}).

Тогда (8.24) — это утверждение Tn→Pθ0T_{n} \xrightarrow {P} \theta_{0}, а (8.26) — утверждение, что 1I(θ0)nl′(θ0)\dfrac {1}{I(\theta_{0}) \sqrt{n}} l'(\theta_{0}) ограничена по вероятности.

Покажите, что (B1) влечёт

  1. (8.24) и

  2. (8.26).