Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Математическая статистика. Учебное пособие

.pdf
Скачиваний:
0
Добавлен:
12.08.2026
Размер:
1 Мб
Скачать

Глава 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

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θ.

Проверить -эффективность оценки максимального правдоподобия параметра θ.

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]