Математическая статистика. Учебное пособие
.pdf
Глава 9. Эффективность и -эффективность |
161 |
Это означает, что ни оценка максимального правдоподобия ^ =
θ
= ( ), ни какая-либо другая оценка параметра θ не могут быть
-эффективными. Это, однако, не является недостатком. Легко проверить, что среднеквадратическое отклонение оценки ( ) от параметра θ убывает с ростом на порядок быстрее, чем допускает неравенство Рао — Крамера для регулярных семейств распределений (см. задачу ). Такие оценки (построение которых возможно лишь при нарушении условий регулярности) иногда называют супер-
эффективными.
Задача 9.11. Пусть 1, . . . , — выборка из нормального распределения со средним и единичной дисперсией. Улучшить оценку* = 1 усреднением при фиксированном значении достаточной статистики . Найти распределение, математическое ожидание, смещение и дисперсию улучшенной оценки. Является ли улучшенная оценка эффективной?
Решение |
|
Из задачи |
следует, что |
E( 1 | ) = .
Так как эта оценка является несмещённой, полной и достаточной (см. задачу ), то она является эффективной. Как известно (задача ), она имеет нормальное распределение с параметрами и σ2/ .
Задача 9.12. Пусть 1, . . . , — выборка из нормального распределения с параметрами и σ2, где R, σ > 0. Проверить, является ли -эффективной несмещённая оценка * = параметра .
Решение
Найдём информацию Фишера относительно параметра (считая, что имеется один неизвестный параметр ). Выпишем плотность распределения и её логарифм:
|
|
2 ( ) = |
1 |
−( − )2/2σ2 |
, ln |
|
2 ( ) = |
|
1 |
ln(2πσ2) |
|
( − )2 |
. |
|
|
|
|
|
|
|
|
|
|||||||
|
, σ |
|
√2πσ2 |
|
|
,σ |
|
− |
2 |
|
− |
2σ2 |
||
162 |
Глава 9. Эффективность и -эффективность |
|
Соответственно, |
ln ,σ2 ( ) = − . |
|
|
∂ |
|
|
∂ |
σ2 |
Найдя второй момент этого выражения при = 1, получим информацию Фишера
( ) = E ( |
∂ |
2 |
E( 1 |
)2 |
|
D 1 |
|
1 |
|
|
|
||||||
∂ ln ,σ2 ( 1)) |
= |
σ4− |
|
= |
σ4 |
= |
σ2 . |
|
Найдём дисперсию оценки : |
|
|
|
|
|
|
||
|
D 1 |
σ2 |
|
|
|
|
|
|
|
D = |
|
= . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Сравнивая левую и правую части в неравенстве Рао — Крамера, получаем равенство
σ2 |
1 |
D = = |
. |
|
( ) |
Итак, оценка * = является -эффективной оценкой.
Задача 9.13. Найти эффективную несмещённую оценку для параметра α показательного распределения по выборке объёма 2.
Решение
Мы будем использовать полную и достаточную статистику (читатель должен был убедиться в этом, решив задачу ). Усредним
по ней оценку максимального правдоподобия α = 1/ (см. зада- |
|
^ |
|
чу ): |
|
E(1/ | ) = 1/ . |
|
Она является эффективной в классе со смещением |
|
α |
α |
(α) = E(1/ ) − α = − 1 − α = |
− 1. |
Здесь мы воспользовались результатом задачи .
|
Глава 9. Эффективность и -эффективность |
163 |
||||||||
Смещение линейно по параметру, а значит, можно применить |
||||||||||
задачу |
и выписать несмещённую эффективную оценку |
|
||||||||
|
α** = |
1 |
|
= |
− 1 |
. |
|
|||
|
|
|
|
|
|
|
|
|||
|
|
(1 + ( − 1)−1) |
|
|
|
|||||
Задача 9.14. Пусть 1, . . . , — выборка из показательного распределения с параметром α, 2. Будет ли -эффективной оценка
α* = − 1?
Решение
Проверим регулярность этого семейства. Плотность этого распределения имеет вид
α |
|
α , |
если > 0, |
|
|
|
|
√ |
|
|
|
α /2, |
если > 0, |
|||||||||
|
|
|
α |
|
||||||||||||||||||
α( ) = {0, |
− |
|
если 0; |
√ α( ) = {0, |
− |
|
если 0. |
|||||||||||||||
Для любого > 0 существует производная функции |
|
|
по |
|||||||||||||||||||
√ α( ) |
||||||||||||||||||||||
α, и эта производная непрерывна во всех точках α > 0: |
|
|||||||||||||||||||||
|
|
∂ |
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
√ α( ) = |
2√ |
|
−α /2 − √α |
|
−α /2. |
|
|
|
|||||||||||
|
|
∂α |
2 |
|
|
|
||||||||||||||||
|
|
α |
|
|
|
|||||||||||||||||
Условие (RR) проверим непосредственным вычислением (α):
α( 1) = α −α 1 , |
ln α( 1) = ln α−α 1, |
|
|
|
∂ |
|
ln α( 1) = |
1 |
|
− 1, |
||||||
|
|
∂α |
α |
|||||||||||||
(α) = E ( |
∂ |
ln α( 1)) |
2 |
( 1 − |
1 |
|
) |
2 |
|
1 |
|
|
|
|
||
|
= E |
|
|
|
|
= D 1 = |
|
. |
|
|||||||
∂α |
α |
|
|
α2 |
|
|||||||||||
Итак, информация Фишера (α) = 1/α2 существует, положительна и непрерывна по α при всех α > 0, т. е. условие (RR) выполнено.
164 |
Глава 9. Эффективность и -эффективность |
В задаче |
мы нашли математическое ожидание оценки макси- |
мального правдоподобия параметра α:
|
|
|
|
|
α |
|
E(1/ ) |
= |
|
|
. |
||
− 1 |
||||||
Следовательно, α* — несмещённая оценка:
Eα* = − 1 E(1/ ) = α.
Итак, оценка α* принадлежит классу 0. Информацию Фишера относительно параметра α мы вычислили: (α) = 1/α2. Найдём второй момент и дисперсию оценки α* . Так как сумма независимых случайных величин, имеющих показательное распределение с параметром α, имеет гамма-распределение с параметрами α и , то
E(α* )2 |
= E(( − 1)2/ |
|
|
|
|
|
|
|
|
2) |
|
= |
|
|
|
|
|
|
||||||||||||
|
|
=1 |
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
∞ |
( |
∑ |
|
|
|
) |
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
1 |
|
|
|
|
|
α |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
= ( − 1)2 |
Z0 |
|
|
|
|
|
|
−1 −α = |
||||||||||||||||||||
|
|
2 |
|
( − 1)! |
||||||||||||||||||||||||||
|
|
|
|
|
|
|
2 |
∞ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
= α2 |
( − |
1) |
Z0 |
|
(α ) |
− |
3 |
− |
α (α ) = |
|
|
|||||||||||||||||
|
|
|
|
|
|
|||||||||||||||||||||||||
|
|
( − 1)! |
|
|
|
|||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
= α2 |
( − 1) |
|
( |
|
|
3)! = |
|
− 1 |
α2 |
|
|
|
||||||||||||||||
|
|
|
( |
− |
2)! |
· |
|
|
|
|
− |
|
|
|
|
|
|
|
|
− |
2 |
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
и |
|
|
|
|
|
|
|
|
|
|
|
|
− 1 |
|
|
|
|
|
|
|
|
|
|
α2 |
|
|||||
Dα* |
= E(α* )2 |
− |
(Eα* )2 |
= |
|
|
α2 |
− |
α2 |
= |
|
. |
||||||||||||||||||
|
|
|
|
|
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
− |
2 |
|
|
|
|
|
|
|
− |
2 |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Подставив дисперсию и информацию Фишера в неравенство Рао — Крамера, получим, что при любом 2 есть строгое неравенство
Dα* |
|
α2 |
> |
α2 |
1 |
|
||
= |
|
|
|
= |
|
. |
||
|
|
|
||||||
|
|
− 2 |
|
|
|
(α) |
||
|
|
|
|
|||||
|
Глава 9. Эффективность и -эффективность |
165 |
Вместе с тем эта оценка является полной и достаточной (см. за- |
||
дачу |
), а значит, эффективной в 0. В силу единственности эф- |
|
фективной оценки в данном классе никакая другая оценка не может быть -эффективной, иначе мы бы имели две различных несмещённых эффективных оценки. А это, в свою очередь, означает, что для α вообще не существует -эффективных несмещённых оценок. Такой же факт имел место, как мы видели, и для равномерного распределения, но здесь ситуация качественно иная, так как параметрическое семейство распределений является регулярным.
Заметим также, что оценка α* |
= − 1 не является -эффектив- |
ной в силу теоремы , так как |
|
|
Ψ (α) = α −α .
Из этого представления функции правдоподобия видно, что -эффек- тивной оценкой для параметра α является оценка , но эта оценка
лежит в очень «плохом» классе оценок со смещением (α) = 1 −α,
α
которое может принимать сколь угодно большие значения; более того, эта оценка не является даже состоятельной оценкой параметра α.
Задача 9.15. Пусть 1, . . . , — выборка из смещённого показательного распределения с плотностью
|
β− , |
если β, |
|
β( ) = {0, |
если < β. |
Найти эффективную несмещённую оценку для параметра β R. |
||
Решение |
|
|
В задаче |
мы показали, что статистика (1) является пол- |
|
ной и достаточной, а значит, и эффективной в классе со смещением,
равным 1/ , как было показано в задаче |
. Очевидно, оценка |
(1) − 1/ будет эффективной в 0. |
|
166 |
Глава 9. Эффективность и -эффективность |
Задача 9.16. Пусть 1, . . . , — выборка из распределения Пуассона с параметром λ. В качестве оценки параметра θ = −λ = = P( 1 = 0) рассматривается θ* = I( 1 = 0). Вычислить смещение(θ) = Eθ* − θ этой оценки и построить эффективную оценку в классе оценок со смещением (θ) усреднением по полной и достаточной для параметра θ статистике.
Решение
Смещение оценки θ* , очевидно, равно нулю:
(θ) = EI( 1 = 0) − P( 1 = 0) = 0.
Статистика является полной и достаточной (см. задачи
и). Заметим, что θ* принимает значения 0 и 1. Поэтому
E(θ* |
= ) = 0 |
· |
P(θ* = 0 |
| |
= )+1 |
· |
P(θ* = 1 |
| |
= ) = |
|
| |
|
|
|
|
|
|
||||
|
|
|
|
|
|
= P( 1 = 0 | = ). |
||||
Вычислив последнюю вероятность по определению условной ве- |
||||||||||
роятности (см. также решение задачи |
), получим |
|
|
|||||||
E(θ* | = ) = (1 − 1/ ) .
Следовательно, оценка θ** = (1 − 1/ ) эффективна в классе несмещённых оценок.
Задача 9.17. Пусть 1, . . . , — выборка из биномиального распределения с параметрами 5 и , где 0 < < 1. Проверить, выполнены ли условия неравенства Рао — Крамера, и доказать, что оценка
* = + 1 неизвестного параметра эффективна в некотором
5
классе. В каком именно?
Решение
Для биномиального распределения плотность относительно считающей меры (она же вероятность попасть в точку ) имеет вид
( ) = P( 1 = ) = 5 (1 − )5− при = 0, 1, 2, 3, 4, 5.
Глава 9. Эффективность и -эффективность |
167 |
Условие (R), очевидно, выполнено: функция
√ √
( ) = /2 (1 − )(5− )/2
5
при всех возможных = 0, 1, . . . , 5 непрерывно дифференцируема по при 0 < < 1.
Чтобы проверить выполнение условия (RR), вычислим информацию Фишера:
|
|
ln |
( ) |
= |
ln |
|
+ |
|
ln + |
|
5 − |
ln(1 |
− |
), |
||||||||
|
|
|
2 |
|
||||||||||||||||||
|
|
|
|
|
|
|
|
5 |
|
|
|
|
|
|
|
2 |
|
|
|
|
||
∂ |
|
ln |
( ) |
= |
|
|
|
|
5 − |
= |
|
− 5 |
. |
|
|
|||||||
∂ |
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
− (1 |
− |
) |
|
(1 |
− |
) |
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Вспомним, что E 1 = 5 и D 1 = 5 (1 − ). Тогда
|
|
|
∂ |
|
|
|
|
2 |
|
|
|
|
|
|
|
1 |
|
5 |
2 |
|
|||||
( ) = |
E ( |
|
|
ln ( 1)) |
|
= E ( |
(1 |
− |
) |
) |
|
|
|||||||||||||
∂ |
|
− |
|
|
|||||||||||||||||||||
|
|
5 (1 − ) |
|
|
|
5 |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
= |
|
= |
|
|
|
|
|
. |
|
|
|
|
|
|
|
|
|||||||||
|
|
2(1 − )2 |
|
(1 − ) |
|
|
|
|
|
|
|
||||||||||||||
Итак, информация Фишера ( ) |
= |
|
5 |
|
|
|
|||||||||||||||||||
|
|
|
|||||||||||||||||||||||
(1 − ) |
|||||||||||||||||||||||||
тельна и непрерывна по при всех 0 < < 1. |
|
|
|||||||||||||||||||||||
Найдём смещение оценки *: |
|
|
|
|
|
|
|
||||||||||||||||||
( ) = E * − |
|
+ 1 |
|
|
|
|
|
+ 1 |
|||||||||||||||||
= |
|
|
|
|
E − = |
|
|
||||||||||||||||||
|
5 |
|
|
|
5 |
||||||||||||||||||||
|
|
|
+ 1 |
5 − |
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
= |
|
|
|
= |
|
|
, |
|
|
|
|
|
|
|
||||||||||
|
|
5 |
|
|
|
|
|
|
|
|
|||||||||||||||
D 1
= 2(1 − )2 =
конечна, положи-
E 1 − =
т. е. оценка * лежит в классе / оценок со смещением / . Вычислим её дисперсию:
D * = |
( + 1)2 |
|
|
|
( + 1)2 |
|
D 1 |
|
( + 1)2 |
|
5 (1 − ) |
|
||
D = |
= |
= |
||||||||||||
|
|
25 2 |
|
|
25 2 |
|
||||||||
|
25 2 |
|
|
|
|
|
||||||||
= |
( + 1)2 (1 − ) |
. |
|
|
|
|
|
|
|
|||||
|
5 3 |
|
|
|
|
|
|
|
|
|
||||
168 |
Глава 9. Эффективность и -эффективность |
Итак, все условия неравенства Рао — Крамера выполнены. Посмотрим, обращается ли для нашей оценки неравенство Рао — Крамера в равенство:
( + 1)2 (1 − ) |
= D * |
(1 + ′( ))2 = |
|
|
5 3 |
|
( ) |
|
|
|
|
= (1 + 1/ )2 |
= |
( + 1)2 (1 − ). |
|
|
5 |
|
5 3 |
|
|
· (1− ) |
|
|
В неравенстве Рао — Крамера достигается равенство, следовательно, оценка * является -эффективной оценкой в классе / .
Задача 9.18. Пусть 1, . . . , — выборка из распределения Коши с параметром сдвига . Исследовать -эффективность выборочной медианы как оценки параметра .
Решение
Без ограничения общности будем считать, что σ2 = 1. В задаче было показано, что Eζ* = . Таким образом, неравенство Рао — Крамера в данном случае принимает следующий вид:
E(ζ* − )2 (1 ).
|
Глава 9. Эффективность и -эффективность |
169 |
|||||||||||||||||
|
Найдём информацию Фишера. Так как ( ) = |
|
|
1 |
, |
||||||||||||||
то |
π(1 + ( − )2) |
||||||||||||||||||
|
( ) = |
E |
∂ ln ( 1) |
) |
2 |
= E |
2( 1 − ) |
|
2 |
= |
|
||||||||
|
|
( |
|
|
∂ |
|
|
|
|
(1 + ( 1 − )2 ) |
|
|
|
||||||
|
= 4E ( |
1 + ( − |
|
)2 2 ) = |
|
|
|
|
|
|
|||||||||
|
|
∞ ( |
( 1 |
|
|
)2 |
) |
|
|
|
|
|
|
|
|
||||
|
|
|
1 −2 |
|
|
|
|
|
|
|
|
|
|||||||
|
= |
4 |
|
|
( − ) |
|
2 · |
|
1 |
|
|
= |
|
||||||
|
|
Z |
|
1 + ( |
|
|
)2 |
π(1 + ( − )2) |
|
|
|
||||||||
|
|
8 ∞ |
( 2 |
− |
|
|
) |
|
|
|
|
|
|
|
|
|
|||
|
|
−∞ |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
= |
π Z0 |
(1 + 2)3 . |
|
|
|
|
|
|
|
|
|
|
||||||
В последнем интеграле сделаем замену = |
1 |
. Тогда |
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
1 + 2 |
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
π |
|
|
|
|
|
|
|
8 |
|
1 |
Z0 1/2(1 − |
)1/2 |
4 |
|
|
|
|
|
|||||||
|
( ) = π |
· |
2 |
= π |
· 8 = 1/2. |
|
|||||||||||||
Теперь найдём E(ζ* − )2.
Распределение Коши симметрично относительно . Значит, является медианой распределения Коши. Тогда, как показано в задаче , выборочная медиана ζ* является асимптотически нормальной оценкой для параметра с коэффициентом асимптотической нормальности (2 ( ))−2 = π2/4, т. е. √ (ζ* − ) η при → ∞, где η имеет нормальное распределение с параметрами 0 и π2/4. Тогда(ζ* − )2 η2 при → ∞ по теореме . Далее воспользуемся следующим вариантом леммы Фату: если последовательность неотрицательных случайных величин {ξ } 1 слабо сходится к случайной величине ξ, то
Eξ lim inf Eξ .
→∞
170 |
Глава 9. Эффективность и -эффективность |
Применяя лемму Фату к последовательности ξ = (ζ* − )2, получаем:
lim inf E(ζ* − )2 Eη2 = π2/4.
→∞
Так как π2/4 > 2 = 1/ ( ), то ζ* не является -эффективной оценкой, по крайней мере, при достаточно больших значениях .
Задача 9.19. Пусть 1, . . . , — выборка из следующего трёхточечного распределения, зависящего от параметра θ (0, 1/3):
P( 1 = 1) = θ, |
P( 1 = 2) = 2θ, |
P( 1 = 3) = 1 − 3θ. |
Проверить -эффективность оценки максимального правдоподобия параметра θ.
