6.3

Эффективное оценивание методом максимального правдоподобия

[29/0%]
Показать
LaTeX
Задача 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^{*}, — измерим.