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

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

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

Глава 3. Метод максимального правдоподобия

61

ЗАДАЧИ ДЛЯ САМОСТОЯТЕЛЬНОГО РЕШЕНИЯ

Задача 3.17. Найти оценку максимального правдоподобия дисперсии σ2 нормального распределения, если математическое ожидание известно.

Задача 3.18. Найти оценку максимального правдоподобия параметра θ > 0, если распределение выборки имеет нормальную плотность с математическим ожиданием θ и дисперсией:

а) 2θ; б) θ2.

Задача 3.19. Найти оценку максимального правдоподобия двумерного параметра ( , ) для равномерного распределения на отрезке [ , ].

Задача 3.20. Найти оценку максимального правдоподобия параметра сдвига µ R распределения Лапласа с плотностью

µ( ) = 12 −| −µ|.

Задача 3.21. Пусть дана выборка из двухпараметрического распределения Лапласа с плотностью

µ,σ( ) = 1 −| −µ|/σ,

2σ

где µ R и σ > 0. Найти оценку максимального правдоподобия двумерного параметра (µ, σ).

Задача 3.22. Распределение Кэптейна определяется плотностью

( ) =

( )

−(θ− ( ))2/2,

 

 

 

θ

2π

где — неубывающая дифференцируемая функция. Найти оценку максимального правдоподобия параметра θ.

62

Глава 3. Метод максимального правдоподобия

Задача 3.23. Найти оценки максимального правдоподобия параметра θ > 0, если распределение выборки имеет плотность:

а) θ θ−1 при [0, 1];

г) θ (lnθ−1 )/ при [1, ];

б) 2 /θ2 при [0, θ];

д) −| |/2(1 − θ) при | | θ.

в) θ θ2/2 /

 

при 0;

 

2π 3

 

Задача 3.24. Пусть дана выборка из распределения Коши с параметром сдвига θ. Построить оценку максимального правдоподобия параметра θ по выборке:

а) объёма 1;

б) объёма 2.

Задача 3.25. На плоской фольге в неизвестной точке находится источник радиоактивного излучения, посылающий лучи равномерно по всем направлениям пространства. Параллельно фольге на расстоянии 1 имеется экран, на котором наблюдаются вспышки, вызываемые радиоактивным излучением. Требуется по местам вспышек на экране определить положение источника излучения на фольге.

Указание: плоскость экрана принять за координатную плоскость0 , ось 0 направить к фольге и рассмотреть оценку какой-нибудь одной координаты точки; используя метод максимального правдоподобия, выписать уравнение правдоподобия и рассмотреть случай

= 2.

Задача 3.26. Пусть дана выборка из биномиального распределения с параметрами и (0, 1). Найти оценку максимального правдоподобия параметра:

Глава 3. Метод максимального правдоподобия

63

а) , если значение параметра известно; б) по выборке объёма = 1, если значение известно.

Задача 3.27. Найти оценку максимального правдоподобия параметра (0, 1) геометрического распределения.

Задача 3.28. Пусть дана выборка из усечённого на заданном уровне 1 геометрического распределения с параметром

(0, 1):

P( 1 = )

=

(1 − ) , = 0, 1, . . . , − 1,

P( 1 = )

=

1 − P( 1 − 1) = (1 − ) .

Найти оценку максимального правдоподобия для .

Задача 3.29. Найти оценку максимального правдоподобия параметра θ равномерного распределения на конечном множестве {1, . . . , θ}, θ — целый положительный параметр.

Задача 3.30. Пусть дана выборка из показательного распределения с параметром α, где α может принимать лишь значения 1, 2 и 3. Построить оценку максимального правдоподобия параметра α.

Задача 3.31. Пусть распределение выборки имеет плотностьθ( ) = ( − θ), где функция имеет единственный максимум в точке = 0. Построить оценку максимального правдоподобия θ*1 параметра сдвига θ по одному наблюдению 1.

Задача 3.32. Пусть в условиях предыдущей задачи функция убывает с ростом | |. Доказать, что оценка максимального правдоподобия θ* , построенная по выборке объёма , лежит в интервале

[ (1), ( )].

Г Л А В А 4

Несмещённость и состоятельность оценок

В этой главе мы познакомимся с базовыми свойствами оценок неизвестного параметра, их взаимосвязями и основными алгоритмами их проверки.

Пусть { θ, θ Θ} — некоторое параметрическое семейство распределений. Пусть 1, . . . , — выборка из распределения θ и θ* = θ* ( 1, . . . , ) — некоторая оценка параметра θ, построенная по данной выборке.

Оценка θ* называется состоятельной оценкой параметра θ, если при любом θ Θ величина θ* сходится по вероятности к θ при

→ ∞.

Оценка θ* называется сильно состоятельной оценкой параметра

θ, если при любом θ Θ величина θ* почти наверное (п. н.) сходится к θ при → ∞.

Смещением оценки θ* называется величина (θ) = Eθ* θ.

Оценка θ* называется несмещённой оценкой параметра θ, если

(θ) = 0 при любом θ Θ.

Оценка θ* называется асимптотически несмещённой оценкой параметра θ, если при любом θ Θ величина (θ) → 0 при → ∞.

Очевидно, из сильной состоятельности следует состоятельность.

Глава 4. Несмещённость и состоятельность

65

Два важнейших асимптотических свойства — состоятельность и асимптотическая несмещённость — не вытекают одно из другого. Из курса теории вероятностей известно, какое условие нужно добавить к сходимости по вероятности, чтобы получить и сходимость математических ожиданий. Это условие равномерной интегрируемости:

→∞ 1

(|ξ |

I(

|ξ |

)

lim sup E

 

 

) = 0.

Имеет место следующее утверждение.

 

p

Теорема 3. Пусть ξ ξ при → ∞ и E|ξ | < ∞ для всех 1.

Тогда

при → ∞ и E|ξ| < ∞

E|ξ ξ| → 0

тогда и только тогда, когда последовательность случайных величин {ξ } 1 равномерно интегрируема.

Так как |Eξ − Eξ| E|ξ ξ|, то из теоремы следует, что последовательность состоятельных оценок, удовлетворяющая условию равномерной интегрируемости, будет также и асимптотически несмещённой.

Существование интегрируемой мажоранты или равномерная ограниченность моментов порядка больше единицы являются простыми достаточными условиями равномерной интегрируемости.

Теорема 4. Пусть последовательность случайных величин {ξ } 1 удовлетворяет одному из следующих условий:

1.|ξ | ξ почти наверное для всех 1 и Eξ < ∞.

2.sup E|ξ |1+ε < ∞ для некоторого ε > 0.

1

Тогда последовательность {ξ } 1 равномерно интегрируема.

Что касается обратной взаимосвязи состоятельных и асимптотическая несмещённых оценок, то с помощью неравенства Чебышёва легко получить следующее утверждение.

Теорема 5. Пусть θ* — асимптотически несмещённая оценка параметра θ и пусть Dθ* → 0 при → ∞ для любого θ Θ. Тогда оценка θ* является состоятельной.

66

Глава 4. Несмещённость и состоятельность

Для доказательства состоятельности и сильной состоятельности часто оказывается полезной так называемая теорема непрерывности.

p

Теорема 6. Пусть — непрерывная функция. Тогда если ξ ξ

p

при → ∞, то (ξ ) → (ξ) при → ∞.

Последняя теорема останется верной, если в ней сходимость по вероятности заменить на сходимость почти наверное или слабую сходимость.

ПРИМЕРЫ РЕШЕНИЯ ЗАДАЧ

Задача 4.1. Для выборки из равномерного распределения на отрезке [0, θ] проверить состоятельность и несмещённость оценки( ) параметра θ.

Решение

 

 

 

Математическое ожидание E ( ) мы уже нашли в задаче

:

 

 

 

E ( ) =

 

θ.

 

 

 

 

+ 1

 

Следовательно, смещение оценки ( ) равно E ( ) θ = −θ/( +1) и она является смещённой, но асимптотически несмещённой.

Проверим состоятельность. Для любого фиксированного ε (0, θ)

P( θ

 

ε) = P(

( )

θ

ε) =

θ ε

)

 

0 при

→ ∞

.

 

|

( )|

 

 

(

θ

 

 

 

При ε > θ

P(|θ ( )| ε) = 0.

Следовательно, оценка ( ) состоятельная.

Более того, она сильно состоятельна, так как из вышеприведённых оценок следует, что для любого ε > 0

P(|θ ( )| ε) < ∞.

=1

Глава 4. Несмещённость и состоятельность

67

Отсюда следует (в силу леммы Бореля — Кантелли), что ( ) θ п. н. при → ∞.

Задача 4.2. Пусть 1, . . . , — выборка из равномерного распределения на отрезке [0, θ]. Проверить состоятельность и несмещённость следующих оценок параметра θ:

а)

(1)

+

( )

;

 

 

 

б) +1

 

 

 

 

;

 

 

 

 

в) ( +1)

(1)

.

 

 

 

 

 

 

 

 

 

 

 

 

 

( )

 

 

 

 

 

 

 

 

 

Решение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а) Воспользуемся результатами задачи

 

 

:

 

 

 

 

 

 

 

E(

(1)

+

( )

) = E + E

( )

=

θ

+ θ = θ.

 

 

 

 

 

 

 

 

 

 

 

(1)

 

 

 

 

+ 1

+ 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Следовательно, оценка (1) + ( ) является несмещённой.

 

 

В задаче

 

мы доказали, что ( )

p

 

при → ∞. Аналогично

 

θ

 

 

 

 

 

 

 

 

 

 

 

p

при

→ ∞. Действительно, так как

 

 

проверяется, что (1) → 0

 

 

P(| (1)| ε) = P( (1) ε) =

 

 

1

 

ε/θ ,

если 0 < ε < θ,

 

 

{0,

 

 

)

 

если ε

 

θ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

 

то P(| (1)| ε) → 0 при → ∞ для любого ε > 0. Значит, (1) → 0

при → ∞.

 

 

 

 

 

 

 

 

p

 

при → ∞, т. е. оценка (1) +

Таким образом, (1) + ( ) θ

+ ( ) является состоятельной.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

б) Так как E

 

= θ/( +1), то оценка

+1

 

( )

является несме-

щённой.

 

 

 

 

 

( )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Так как

( + 1)/ → 1

при

 

→ ∞p

и, как показано в задаче ,

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( ) θ

при → ∞, то +1 ( )

 

θ при → ∞, т. е. оценка

+1

( )

состоятельна.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

в) Так как E (1)

= θ/( + 1), то оценка ( + 1) (1)

является

несмещённой.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

68

Глава 4. Несмещённость и состоятельность

Найдём слабый предел последовательности ( + 1) (1). Так как

в силу задачи

P(( + 1) (1) < ) = 1 − (1 − ( /( + 1)) и

 

 

( ) =

/ ,

если 0 <

,

 

 

 

 

 

0,

θ

если

 

0,

θ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

то

 

 

 

1,

 

если > θ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,

 

 

 

 

 

 

 

 

 

 

 

P(( +1) (1)

 

 

 

1

 

/θ

 

 

 

,

если 0,

< ) = 1

 

( + 1)

 

 

 

если 0 < ( + 1)θ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(

 

 

 

 

)

 

 

 

 

 

1,

 

 

 

 

 

 

если > ( + 1)θ.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

для всех

 

 

при

 

 

 

Тогда, очевидно,

 

 

 

 

R

 

 

→ ∞

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,

 

 

 

 

 

если 0,

P(( + 1) (1)

< ) → {1 − − /θ,

если > 0.

Предельная функция является функцией показательного распределения с параметром 1/θ. Значит, слабый предел последовательности ( + 1) (1) невырожден, а оценка ( + 1) (1) не является состоятельной.

Задача 4.3. Проверить несмещённость и состоятельность оценки максимального правдоподобия векторного параметра ( , σ2) нормального распределения.

Решение

 

 

 

 

 

 

 

 

 

 

 

 

В задаче

 

мы нашли оценку максимального правдоподобия па-

раметра ( , σ

 

): ^ = , σ =

 

=

 

−( )

 

. Проверим состоятель-

 

2

 

 

^2

2

 

2

 

 

 

2

 

ность. В силу усиленного закона больших чисел → E 1 = п. н. и2 → E 12 = 2 +σ2 п. н. при → ∞, откуда следует (в силу элементарных свойств сходимости почти наверное), что 2 σ2 п. н. при→ ∞. Таким образом, обе координаты нашей векторной оценки являются сильно состоятельными.

Глава 4. Несмещённость и состоятельность

69

Теперь проверим несмещённость. Как мы знаем (см. задачу

),

E = , E 2 = ( − 1)σ2/ .

Следовательно, ^ = — несмещённая оценка, а оценка σ^2 = 2 имеет ненулевое смещение −σ2/ , являясь лишь асимптотически несмещённой. Несмещённой оценкой для дисперсии нормального распределения σ2 будет, очевидно, 02 = 2/( − 1).

Задача 4.4. Проверить состоятельность и несмещённость оценки α* = 1/ параметра α показательного распределения.

Решение

В силу усиленного закона больших чисел → E 1 = 1/α п. н. при → ∞. Так как функция ( ) = 1/ является непрерывной в области > 0, то α* α п. н. при → ∞, что означает сильную состоятельность оценки.

Для проверки несмещённости вычислим

E 1/ = E(

 

)−1.

=1

(

 

)

Воспользуемся устойчивостью гамма-распределения по суммированию (показательное распределение с параметром α есть гамма-

распределение с параметрами α и 1), из которой следует, что

=1

имеет гамма-распределение с параметрами α и с плотностью

γα, ( ) = α −1 −α /Γ( ),

0,

70

 

Глава 4. Несмещённость и состоятельность

а значит,

 

 

 

 

 

 

 

 

 

 

 

 

 

E

 

 

 

1

Z

1

 

 

 

α

=1

 

=

· γα, ( ) = Z

Γ( ) −2 −α =

(

 

)

 

 

R

 

 

 

0

 

 

 

 

 

 

 

 

 

 

α

α −1

 

 

 

 

 

 

 

=

 

 

Z0

 

−2 −α =

 

 

 

 

 

 

− 1

Γ( − 1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

α

Z0 γα, −1( ) =

α

,

− 1

− 1

откуда

 

 

 

 

 

 

α

 

 

 

 

 

 

 

 

 

 

 

 

E(1/ )

=

 

 

.

 

 

 

− 1

 

 

Таким образом, оценка α* является смещённой, но асимптотически несмещённой оценкой параметра α.

Задача 4.5. Привести пример оценки, которая является несмещённой, но не является состоятельной.

Решение

Рассмотрим распределение Пуассона с параметром λ. Пусть λ* = = 1. Эта оценка является несмещённой для параметра λ, но не является состоятельной (так как её слабый предел при → ∞ совпадает с 1 и, следовательно, невырожден).

Вместо распределения Пуассона можно рассмотреть любое невырожденное распределение с конечным математическим ожиданием, в качестве параметра взять θ = E 1 и θ* = 1 в качестве оценки. Эта оценка является несмещённой для θ, но не является состоятельной.

Задача 4.6. Пусть элементы выборки 1, . . . , объёма 5 имеют распределение Пуассона с параметром λ > 0. Для какого параметра θ = θ(λ) оценка

θ* = 1 · 2 · . . . · 5

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