Математическая статистика. Воскобойников Ю.Е - 18 стр.

UptoLike

35 36
Эта вероятность есть функция от
n
xxx ...,,,
21
,
которая называется
функцией максимального правдоподобия и обозначается
KK )xX(P),x,,x,x(L
n 1121
=
=
θ
)xX(P
nn
= .
Учитывая, что значение
i
y встречается в выборке
j
n
раз, получаем
)(...)(),...,,(
1
1
1
θθθ
m
n
m
n
n
ppxxL = .
Как уже было сказано, суть метода максимального правдоподобия
состоит в том, что в качестве параметра
θ
берется такое значение, ко-
торое максимизирует функцию
),...,,(
1
θ
n
xxL . Полученное значение,
если оно существует, является функцией от
n
xxx ...,,,
21
, т.е.
)...,,,(
21
*
nМП
xxx
θθ
=
. Заменяя элементы
n
x...,,x,x
21
случайными
величинами
n
XX ...,,
1
, получаем оценку максимального правдоподо-
бия
)X...,,X,X(
n
*
МП 21
θ
.
Точка максимума функции
),...,,(
1
θ
n
xxL удовлетворяет нелиней-
ному (в общем случае) уравнению
0
),,...,(
1
=
θ
θ
n
xxL
(3.20)
и поэтому конкретное значение оценки
)...,,,(
21
*
nМП
xxx
θ
опреде-
ляют как корень уравнения (3.20).
Функции
),...,,(
1
θ
n
xxL ) и ),...,,(ln
1
θ
n
xxL достигают макси-
мума при одном и том же значении
θ
. Поэтому вместо отыскания мак-
симума функции
),...,,(
1
θ
n
xxL находят максимум функции
),...,,(ln
1
θ
n
xxL . Эта функция получила название логарифмической
функции правдоподобия.
Построение оценки максимального правдоподобия можно разбить на
следующие этапы:
1. Определяем производную логарифмической функции правдопо-
добия по параметру
θ
.
2. Приравнивая производную к нулю, находим критическую точку
кр
θ
корень уравнения правдоподобия 0
),,...,(
1
=
θ
θ
n
xxL
.
3. Находим вторую производную
2
2
ln
θ
L
и ее значение в точке
кр
θ
. Если вторая производная в точке
кр
θ
меньше нуля, то в точке
кр
θ
функция ),x,...,x(L
n
θ
1
достигает максимума.
Найденная таким образом
*
МП
θ
является функцией случайных ве-
личин
n
XXX ...,,
21
и, следовательно, сама является случайной вели-
чиной. Конкретное значение оценки
*
МП
θ
получается при подстановке
в
)...,,(
1
*
nМП
XX
θ
вместо
n
XXX ...,,
21
значений выборки
n
xxx ...,,,
21
.
Непрерывная генеральная совокупность. Рассмотрим случай, ко-
гда генеральная совокупность имеет непрерывный ряд распределения.
Функцию максимального правдоподобия определим по правилу
),x(p),x(p),x...,,x(L
nn
θ
θ
θ
L
11
=
,
где
),(
θ
xp плотность распределения генеральной совокупности. Все
остальное, изложенное для дискретного случая, переносится на непре-
рывный.
Пример 3.1. Проводится
п независимых опытов, в каждом из кото-
рых событие
А повторяется с неизвестной вероятностью р. Рассмотрим
генеральную совокупность
Х количество появлений события А в од-
ном опыте. По выборке
n
xx ...,,
1
из генеральной совокупности Х необ-
ходимо оценить параметр
р.
Решение. Выборка
n
xx ...,,
1
состоит из нулей и единиц, причем
1=
i
x , если в i -м опыте событие А произошло, и 0
=
i
x , если собы-
тие не произошло. Предположим, что
т частота появления события А
в п опытах. Тогда выборка
n
xx ...,,
1
содержит m единиц и )( mn
нулей. Так как
pXPpXP
=
=
=
=
1)0(,)1( , то
mnm
n
)p(p),x...,,x(L
= 1
1
θ
.
Найдем точку максимума логарифмической функции максимального
правдоподоподобия
Lln
)1ln()(ln),...,,(ln
1
pmnpmxxL
n
+
=
θ
.
Определим из уравнения
Эта вероятность есть функция от x1 , x 2 , ..., x n , которая называется                                                                                ∂ 2 ln L
                                                                                                      3. Находим вторую производную                                 и ее значение в точке
функцией         максимального        правдоподобия    и    обозначается
                                                                                                                                                         ∂θ 2
L( x1 , x2 ,K, xn ,θ ) = P( X1 = x1 )K P( X n = xn ) .
                                                                                                  θ кр .    Если вторая производная в точке              θ кр   меньше нуля, то в точке
   Учитывая, что значение y i встречается в выборке n                        раз, получаем
                                                                         j
                                                                                                  θ кр     функция L( x1 ,..., x n ,θ ) достигает максимума.

                    L( x1 , ..., x n ,θ ) = p1n1 (θ ) ... p mnm (θ ) .                                Найденная таким образом
                                                                                                                                              *
                                                                                                                                            θ МП   является функцией случайных ве-
   Как уже было сказано, суть метода максимального правдоподобия                                  личин X 1 , X 2 ..., X n и, следовательно, сама является случайной вели-
состоит в том, что в качестве параметра θ берется такое значение, ко-                                                                                *
                                                                                                  чиной. Конкретное значение оценки                θ МП     получается при подстановке
торое максимизирует функцию L( x1 , ..., x n ,θ ) . Полученное значение,                                *
                                                                                                  в   θ МП ( X 1 , ..., X n )        вместо     X 1 , X 2 ..., X n значений выборки
если оно         существует, является           функцией        от x1 , x 2 , ..., x n ,   т.е.
      *
                                                                                                  x1 , x 2 , ..., x n .
θ = θ МП ( x1 , x 2 , ..., x n ) . Заменяя элементы x1 , x2 ,..., xn случайными                      Непрерывная генеральная совокупность. Рассмотрим случай, ко-
величинами X 1 , ..., X n , получаем оценку максимального правдоподо-                             гда генеральная совокупность имеет непрерывный ряд распределения.
                                                                                                  Функцию максимального правдоподобия определим по правилу
         *
бия    θ МП ( X 1 , X 2 , ..., X n ) .                                                                                    L( x1 , ..., xn ,θ ) = p( x1 ,θ )L p( xn ,θ ) ,
   Точка максимума функции L( x1 , ..., x n ,θ ) удовлетворяет нелиней-                           где p ( x,θ ) – плотность распределения генеральной совокупности. Все
ному (в общем случае) уравнению                                                                   остальное, изложенное для дискретного случая, переносится на непре-
                                 ∂L( x1 ,..., x n ,θ )                                            рывный.
                                                       =0                              (3.20)        Пример 3.1. Проводится п независимых опытов, в каждом из кото-
                                        ∂θ
                                                                                                  рых событие А повторяется с неизвестной вероятностью р. Рассмотрим
                                                        *
и поэтому конкретное значение оценки                  θ МП ( x1 , x 2 , ..., x n )   опреде-      генеральную совокупность Х – количество появлений события А в од-
ляют как корень уравнения (3.20).                                                                 ном опыте. По выборке x1 , ..., x n из генеральной совокупности Х необ-
   Функции L( x1 , ..., x n ,θ ) ) и ln L( x1 , ..., x n ,θ ) достигают макси-                    ходимо оценить параметр р.
мума при одном и том же значении θ . Поэтому вместо отыскания мак-                                   Решение. Выборка x1 , ..., x n состоит из нулей и единиц, причем
симума функции L( x1 , ..., x n ,θ ) находят максимум функции                                     xi = 1 , если в i -м опыте событие А произошло, и xi = 0 , если собы-
ln L( x1 , ..., x n ,θ ) . Эта функция получила название логарифмической                          тие не произошло. Предположим, что т – частота появления события А
функции правдоподобия.                                                                            в п опытах. Тогда выборка x1 , ..., x n содержит m единиц и ( n − m)
   Построение оценки максимального правдоподобия можно разбить на
следующие этапы:                                                                                  нулей. Так как P ( X = 1) = p, P ( X = 0) = 1 − p , то
   1. Определяем производную логарифмической функции правдопо-                                                                                                      n−m
                                                                                                                             L( x1 , ..., xn ,θ ) = p m ( 1 − p )         .
добия по параметру θ .
                                                                                                     Найдем точку максимума логарифмической функции максимального
   2. Приравнивая производную к нулю, находим критическую точку
                                                                                                  правдоподоподобия ln L
                                                       ∂L( x1 ,..., x n ,θ )
θ кр   – корень уравнения правдоподобия                                      = 0.                                ln L( x1 , ..., x n ,θ ) = m ln p + (n − m) ln(1 − p)        .
                                                              ∂θ                                  Определим из уравнения
                                               35                                                                                                  36