6.1

Оценка качества в больших выборках

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

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

?
(a)

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

(b)

Для положительного целого kk покажите, что E[(Xˉ−ξ)2k−1]\mathbb {E}\left[(\bar{X}-\xi )^{2 k-1}\right] и E[(Xˉ−ξ)2k]\mathbb {E}\left[(\bar{X}-\xi )^{2 k}\right], если они существуют, оба являются 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]=ξ\mathbb {E}\left[X_1\right]=\xi, Var⁡[X1]=σ2\operatorname {Var}\left[X_1\right]=\sigma^2 и конечным четвёртым моментом, и пусть hh — функция вещественной переменной, первые четыре производные которой h′,h′′,h′′′,h(iv)h', h'', h''', h^{(iv)} существуют для всех x∈Ix \in I, где II — интервал с P(X1∈I)=1\mathbb {P}\left(X_1 \in I\right)=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, \mathbb {E}\left[h(\bar{X})\right]=h(\xi )+\frac{\sigma ^{2}}{2 n} h''(\xi )+R_{n},

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

Var⁡[h(Xˉ)]=σ2n[h′(ξ)]2+Rn, \operatorname {Var}\left[h(\bar{X})\right]=\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−θ\mathbb {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)\mathbb {E}\left[\left|\bar{X}-\xi \right|^{2 k-1}\right]=O\left(n^{-k+1 / 2}\right).

?
Задача 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]=∞\mathbb {E}\left[\left(\sqrt{n}[h(\bar{X}) -h(\xi )]\right)^{2}\right]=\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ˉ]\mathbb {E}\left[\sqrt{\bar{X}}\right] и Var⁡[Xˉ]\operatorname {Var}\left[\sqrt{\bar{X}}\right].

(b)

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

Var⁡[Xˉ]=[1−1n(Γ(n+1/2)Γ(n))2]θ \operatorname {Var}\left[\sqrt{\bar{X}}\right]=\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}\left[\sqrt{\bar{X}}\right]=\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) \mathbb {E}\left[h(T)\right]=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}\left[h(T)\right]=\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) \mathbb {E}\left[h(\bar{X})\right]=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}\left[h(\bar{X})\right]=\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}=\mathbb {E}\left[(X-\xi )^{3}\right].

?
Задача 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/n\mathbb {E}\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)\mathbb {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, \mathbb {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)→1\mathbb {P}\left(\left|h(T_{n})-h(\xi )\right|<a\right) \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=\mathbb {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.

?