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

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

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

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

71

будет несмещённой оценкой? Является ли эта оценка состоятельной оценкой для того же параметра θ?

Решение

В силу независимости случайных величин 1, 2, . . . , имеем равенство

Eθ* = E( 1 · 2 · . . . · 5) = E 1 · E 2 · . . . · E 5 = λ5.

Итак, оценка θ* является несмещённой оценкой для параметра θ = = θ(λ) = λ5. Состоятельной она не является, поскольку не зависит от. Соответственно, при → ∞ оценка θ* стремится по вероятности разве что к себе самой, не равняясь при этом λ5.

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

оценок θ*, = ( + 1) параметра θ, = 1, 2, . . .

Решение

Для доказательства состоятельности, как и в предыдущих задачах, воспользуемся усиленным законом больших чисел. Так как при каждом = 1, 2, . . . выборочный -й момент почти наверное схо-

дится к E 1 = θ /( + 1) при → ∞, а функция ( ) = ( + 1) непрерывна, то и θ*, = ( ) сходится почти наверное к θ =

= (E 1 ) при → ∞, т. е. для каждого = 1, 2, . . . оценка θ*, является сильно состоятельной.

Вообще, справедлива следующая теорема, которую читатель может легко доказать, используя усиленный закон больших чисел и теорему непрерывности (см. теорему ): если функция (θ) являет-

ся непрерывной и строго монотонной для некоторой функции , то (одномерная) оценка метода моментов θ* = −1( ( )) является

сильно состоятельной.

Проверить несмещённость θ*, , непосредственно вычислив математические ожидания этих оценок при произвольном объёме выборки, как мы это сделали для оценки максимального правдоподобия параметра показательного распределения в задаче , очень за-

72

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

труднительно, так как выписать плотность случайной величины

вявном виде для произвольных и технически довольно сложно.

Вситуациях, когда оценка является строго выпуклой (или строго вогнутой) функцией от того или иного выборочного момента, полезным оказывается неравенство Йенсена.

Пусть случайная величина ξ имеет невырожденное распределение, а функция ( ) строго выпукла. Тогда E (ξ) > (Eξ). Для строго

вогнутых функций справедливо противоположное неравенство.

Так как при каждом = 2, 3, . . . функция ( ) = ( + 1) является строго вогнутой в области > 0, то неравенство Йенсена гарантирует, что

Eθ*, < ( + 1)E = θ,

и, значит, оценки θ*, , 2, являются смещёнными при любом

1.

При = 1 имеем несмещённую оценку:

Eθ*,1 = E(2 ) = 2E 1 = θ.

Проверим, являются ли оценки θ*, , 2, хотя бы асимптотически несмещёнными. Это возможно и без знания точного вида смещения. Так как все эти оценки (сильно) состоятельны, то для сходимости последовательности величин Eθ*, к параметру θ при → ∞ достаточно, чтобы при любом 2 существовала равномерная поинтегрируемая мажоранта для θ*, (см. теорему ). Так как распределение элементов выборки имеет ограниченный носитель, то

такая мажоранта существует. Например, в качестве такой мажоран-

ты можно выбрать вырожденную случайную величину + 1θ: так как P(0 θ) = 1 для всех = 1, 2, . . ., то для всех 1

0 θ*, =

( + 1)

 

( + 1)θ =

( + 1)θ п. н.

 

 

 

 

 

 

 

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

73

Задача 4.8. Пусть элементы выборки 1, . . . , имеют равномерное распределение на отрезке [ , ], где < . Проверить, явля-

ется ли оценка θ* = ( ) (1):

а) несмещённой; б) состоятельной оценкой параметра θ = − .

Решение

а) Функции распределения случайных величин ( ) и (1) равны

(см. задачу ), соответственно,

 

 

 

 

 

 

 

 

0,

 

 

 

 

 

< ,

( )

( ) =

 

 

 

 

 

 

,

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

)

 

 

 

 

 

 

(

 

 

 

 

 

> ;

 

 

1,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,

 

 

 

 

 

 

 

< ,

(1)

 

 

 

 

 

 

 

 

 

( ) =

1

 

 

 

 

, ,

 

 

 

 

 

 

 

 

 

 

 

 

 

(

 

 

 

 

)

 

 

 

 

 

 

 

 

> .

 

 

1,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Обозначим

= ( − )/( − ), = 1, 2, . . . , . Тогда, очевидно,

случайные величины 1, . . . , независимы и имеют равномерное распределение на отрезке [0, 1].

Так как (1) = +( − ) (1) и ( ) = +( − ) ( ), то, используя

результат задачи

, получаем

 

 

 

 

 

 

 

 

 

 

 

E

(1)

= + (

)E

(1)

= +

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ 1

 

 

 

 

 

 

E

 

= + (

)E

 

= +

( − )

=

,

 

( )

 

 

 

 

 

 

 

( )

 

 

+ 1

 

 

+ 1

 

 

поэтому

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Eθ* = E

( )

E

 

=

 

2( − )

=

− 1

(

) =

,

 

 

 

(1)

 

 

 

 

+ 1

 

 

+ 1

 

̸

 

т. е. θ* не является несмещённой оценкой параметра θ.

74

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

Этот вывод можно получить и не вычисляя математических ожиданий: заметим, что θ* < θ почти наверное, так что Eθ* < θ (в силу соответствующего свойства математического ожидания).

б) Проверим состоятельность. Докажем, что ( ) сходится по вероятности к и (1) сходится по вероятности к , откуда сразу же будет следовать, что θ* сходится по вероятности к θ. Пусть ε > 0. Согласно (например) неравенству Маркова,

P(| ( ) − | > ε)

= P( ( ) > ε)

E( −ε ( )) =

 

=

0 при

→ ∞

;

 

 

 

( + 1)ε

 

 

 

 

P(| (1) − | > ε)

 

P( (1) − > ε)

E( (1)

)

=

 

ε

=

 

=

0 при

→ ∞

.

 

 

 

( + 1)ε

 

 

 

 

Исчезновение модулей объясняется тем, что ( ) < и (1) > с вероятностью 1.

Итак, θ* — состоятельная оценка параметра θ = − .

Задача 4.9. Пусть 1, . . . , — выборка из распределения Пуассона с параметром λ. Для какого параметра θ = θ(λ) статистика

θ* = будет состоятельной оценкой? Является ли θ* несмещённой оценкой того же параметра?

Решение

Так как является состоятельной оценкой параметра λ (в силу

закона больших чисел), то θ*

=

p

θ(λ) = λ λ при

→ ∞

в

 

 

 

 

силу теоремы .

 

 

 

 

 

Вычислим математическое ожидание оценки θ* . В силу устойчивости распределения Пуассона по суммированию случайная величина ξ = = 1 + . . . + имеет распределение Пуассона

 

 

( λ)

с параметром λ, т. е. P(ξ

= ) =

 

λ для = 0, 1, 2, . . .

 

 

 

!

 

 

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

75

Поэтому

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ξ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Eθ*

= E

 

 

 

 

 

 

 

= E(

ξ/ ) = =0

− / P(ξ = ) =

 

 

(

 

 

)

 

 

 

 

 

 

 

 

 

 

 

 

( λ)

 

λ

λ −1/

 

=

− /

λ =

=

 

=0

!

 

=1

( (

 

1)!)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= λ −1/

 

 

λ

 

 

 

( λ −1/ ) −1

= λ λ( −1/ −1)−1/ .

 

 

 

 

 

·

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1=0

 

(

1)!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Итак, Eθ*

= λ λ( −1/ −1)−1/ = θ(λ) = λ λ, т. е. оценка θ*

не

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

̸

 

 

 

 

 

 

 

 

является несмещённой оценкой параметра θ = θ(λ) = λ λ.

 

Так как −1/ = 1 − 1/ + ( −2) при → ∞, то

 

 

 

 

Eθ*

 

 

λ( −1/

1)

1/

 

 

λ

1/ + (1/ )

λ

λ

= θ

 

 

= λ

 

 

 

 

 

 

 

= λ

 

 

 

 

 

при

→ ∞

, т. е. θ*

 

является асимптотически несмещённой оценкой

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

параметра θ = λ λ.

Асимптотическую несмещённость θ* можно доказать и без непосредственного вычисления Eθ* . В силу закона больших чисел сходится по вероятности (а значит, и слабо) к λ. Но в силу эквивалентного определения слабой сходимости это означает, что для любой ограниченной и непрерывной функции ( ) имеет место сходимость E ( ) → (λ) при → ∞. Взяв ( ) = , получим асимптотическую несмещённость θ* .

Задача 4.10. Будет ли статистика состоятельной оценкой параметра сдвига θ распределения Коши?

Решение

Найдём распределение выборочного среднего . Для этого воспользуемся методом характеристических функций. Характеристическая функция распределения Коши с параметрами θ и σ2 имеет вид

ϕ( ) = θ σ| |.

76

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

По свойствам характеристических функций

ϕ ( ) = ϕ ( / ) = ( θ/ −|σ / |) = ϕ( ).

Таким образом, случайная величина при любом 1 имеет распределение Коши с теми же параметрами, что и элементы выборки, а вовсе не вырожденное в точке θ распределение. Оценка не является состоятельной.

Задача 4.11. Проверить несмещённость выборочной медианы ζ* как оценки параметра сдвига нормального распределения.

Решение

Пусть сначала выборка имеет нечётный объём = 2 − 1. Тогда выборочная медиана равна ( ). Обозначим = ( − )/σ,= 1, 2, . . . , . Тогда 1, . . . , — выборка из стандартного нормального распределения.

Очевидно, E ( ) = (E ( ) − )/σ. Пользуясь результатом задачи , мы можем выписать математическое ожидание ( ) по опреде-

лению:

ZZ

E

=

( ) =

−1

Φ −1( )(1

Φ( ))

ϕ( ) ,

( )

( )

 

−1

 

 

 

RR

где через Φ( ) обозначена функция распределения стандартного нормального закона, а через ϕ( ) — соответствующая плотность. Подынтегральная функция является нечётной, так как стандартное нормальное распределение симметрично, а − = − 1:

Φ −1( )(1 − Φ( )) = Φ −1( )Φ(− ) = (Φ( )Φ(− )) −1

и ϕ(− ) = ϕ( ). Кроме того, интеграл сходится абсолютно:

Z | | 1

Φ −1( )(1 − Φ( ))

ϕ( )

1

Z

| |ϕ( ) =

1

 

 

 

1

 

 

 

 

 

 

R

 

 

 

 

R

 

 

 

 

 

 

 

=

 

−1E

 

1|

<

.

 

 

 

 

−1

|

 

 

 

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

77

Сходящийся интеграл от нечётной функции по симметричной области равен нулю. Но тогда

E ( ) = + σE ( ) = .

Пусть теперь = 2 . Тогда ζ* = ( ( ) + ( +1))/2 и

(Eζ* − )/σ = (E ( ) + E ( +1))/2 =

1

ZR ( ) ϕ( ) = 0,

2

где через ( ) мы обозначили ограниченную чётную функцию

( ) =

−1

Φ −1

( )(1

Φ( ))

+

 

 

Φ ( )(1

Φ( )) −1 =

 

−1

 

 

 

 

 

 

−1

 

 

 

 

 

=

−1

Φ −1

( )(1

Φ( )) −1(1

Φ( ) + Φ( )) =

 

 

 

−1

 

 

 

 

 

 

 

 

(Φ( )Φ(− ))

 

= 1

(Φ( )(1 − Φ( ))

 

= 1

 

 

1

 

 

 

 

 

1

 

 

 

 

1

 

 

 

 

1

(мы использовали равенство −1

=

 

, справедливое при =

= 2 ). Окончательно,

 

 

−1

 

−1

 

 

 

 

 

 

 

Eζ* = . Оценка является несмещённой.

Задача 4.12. Проверить несмещённость выборочной медианы ζ* как оценки параметра сдвига распределения Коши.

Решение

Пусть и — функция распределения и плотность распределения Коши с параметрами 0 и 1 соответственно:

( ) =

1

+ 1 arctan( )

и ( ) =

1 .

 

2

π

 

π(1 + 2)

Несложно заметить, что это распределение симметрично, т. е.

(− ) = 1 − ( ) и (− ) = ( ) для всех R.

Обозначим = ( − )/σ, = 1, 2, . . . , . Тогда, используя, например, формулу плотности линейного преобразования из реше-

ния задачи , получаем, что 1, . . . , — выборка из распределения Коши с параметрами 0 и 1. Далее заметим, что при решении задачи по существу использовалась только симметричность

78

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

стандартного нормального распределения. Повторяя решение задачи с заменой Φ на и ϕ на соответственно, получим равенство Eζ* = и для выборки из распределения Коши. Единственный момент, который требует проверки, это сходимость интеграла (при

= 2 − 1 и = 2 )

E ( ) = Z 1

−1

1

 

R

 

Пусть = 2 − 1. Тогда надо

Z

| | −1( )(1− ( )) ( )

( )(1 − ( )) ( ) .

проверить, что

Z

= | |( ( )(1− ( ))) −1 ( ) =

R R

Z1

=| −1( )|( (1 − )) −1 < ∞.

0

Читателю предлагается убедиться, что последнее условие выполнено при 2 (см. также задачу ).

Случай = 2 сводится к проверке сходимости того же интегра-

ла.

Итак, выборочная медиана ζ* является несмещённой оценкой параметра сдвига распределения Коши.

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

Задача 4.13. Пусть дана выборка из распределения Вейбулла с параметрами α и θ, причём значение параметра α известно. Проверить на несмещённость и состоятельность оценку 1/ α параметра θ.

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

79

Решение

Так как значение параметра α известно, то можно рассмотреть преобразованную выборку 1, . . . , , где = α. Заметим, что при

> 0

 

1/α

θ

P( 1 < ) = P( 1 < 1/α) =

Z0

θ α θα α−1 = Z0 = 1− θ .

Здесь мы использовали замену переменной = θ α (см. также решение задачи ). Таким образом, преобразованная выборка имеет показательное распределение с параметром θ, а оценка 1/ α = = 1/ , как показано в задаче , является состоятельной и смещённой (но асимптотически несмещённой) оценкой параметра θ.

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

β,

если β,

β( ) = {0,

если < β.

Выяснить, являются ли несмещёнными и состоятельными следующие оценки параметра сдвига β:

а) (1); б) − 1.

Решение

а) Обозначим = − β, = 1, 2, . . . , . Из решения задачи следует, что 1, . . . , — выборка из показательного распределения с параметром 1. Тогда (1) = (1) + β и для любого ε > 0

P(| (1) β| ε) =

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

= P ( 1 ε) = ε → 0 при → ∞.

p

при → ∞, т. е. оценка максимального

Это означает, что (1) β

правдоподобия (см. задачу

) β = (1) является состоятельной.

 

^

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

Отметим также, что равенство P( (1) ε) = ε для всех ε > 0 означает, что (1) = (1) β имеет показательное распределение

с параметром . Отсюда, в частности, следует, что E (1) = E (1) +

+ = + 1/ , т. е. оценка ^ = (1) является смещённой, но асимп-

β β β

тотически несмещённой.

б) Так как E 1 = E 1 + β = β + 1, то из свойств математического ожидания и закона больших чисел очевидным образом следует,

что оценка метода моментов (см. задачу

) β*

=

1 является

несмещённой и состоятельной:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Eβ* = E

1 = E

1

1 = β,

β*

=

1

 

p

E

1

1 = β.

 

 

 

 

 

 

 

 

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

α−1 −( −β)/α,

если β,

α,β( ) = {0,

если < β.

Проверить состоятельность оценок параметров масштаба α > 0 и сдвига β R, построенных по методу моментов и по методу максимального правдоподобия.

Решение

Оценки метода моментов α* , β* и максимального правдоподо-

бия ^ , ^ мы нашли в задачах и соответственно:

α β

α* = = 2 − ( )2, β* = − ,

 

 

 

^

 

 

^

 

 

 

 

 

 

 

 

 

 

 

 

β = (1).

 

 

 

 

 

 

 

α = − (1),

 

 

 

 

 

В задаче

мы также нашли первый и второй моменты нашего

распределения:

,

2

2

+(α+β)

2. Так как

 

p

и

 

 

 

 

p

E 1 = α+β

E 1 = α

 

→ E 1

 

 

 

 

 

 

 

 

 

 

2

→ E 12 при → ∞ в силу закона больших чисел, то при → ∞

α* = →

 

= α,

β* = − → E 1 α = β,

 

E 12 − (E 1)2

 

 

p

 

 

 

 

 

 

 

p

 

 

 

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