6.7

Обобщения

[34/0%]
Показать
LaTeX
Задача 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]=α+βti\mathbb {E}\left[X_{i}\right]=\alpha +\beta t_{i} и Var⁡[Xi]=1\operatorname {Var}\left[X_{i}\right]=1, и

(b)

E[Xi]=α+βti\mathbb {E}\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)]=0\mathbb {E}_{\theta_{0}}\left[\psi (X, t)\right]=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]\mathbb {E}_{\theta_{0}}\left[\left.\frac{\partial }{\partial t} \psi (X, t)\right|_{t=t_{0}}\right] конечно и отлично от нуля,

  2. Eθ0[ψ2(X,t0)]<∞\mathbb {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}=\mathbb {E}_{\theta_{0}}\left[\psi^{2}\left(X, t_{0}\right)\right] /\left(\mathbb {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.

?