Логарифм функции правдоподобия

 

Система уравнений для нахождения оценок параметров

 

Из первого уравнения следует:

или окончательно

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

Из второго уравнения можно найти .

Эмпирическая дисперсия является смещенной. После устранения смещения

Фактические значения оценок параметров: m =27,51, s² = 0,91.

Для проверки того, что полученные оценки максимизируют значение функции правдоподобия, возьмем вторые производные

 

Вторые производные от функции ln L(m,S) независимо от значений параметров меньше нуля, следовательно, найденные значения параметров являются оценками максимального правдоподобия.

Метод максимального правдоподобия позволяет получить состоятельные, эффективные (если таковые существуют, то полученное решение даст эффективные оценки), достаточные, асимптотически нормально распределенные оценки. Этот метод может давать как смещенные, так и несмещенные оценки. Смещение удается устранить введением поправок. Метод особенно полезен при малых выборках.

 

Метод моментов

 

Метод предложен К. Пирсоном в 1894 г. Сущность метода:

· выбирается столько эмпирических моментов, сколько требуется оценить неизвестных параметров распределения. Желательно применять моменты младших порядков, так как погрешности вычисления оценок резко возрастают с увеличением порядка момента;

· вычисленные по ЭД оценки моментов приравниваются к теоретическим моментам;

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

Пример 1.4. Предположим, что случайная величина Х имеет гамма-распределение. Необходимо найти оценки параметров этого распределения (можно отметить, что нормальное распределение является частным случаем гамма-распределения).

Решение. Функция плотности гамма-распределения имеет вид

Распределение характеризуется двумя параметрами ν и λ, поэтому следует выразить один параметр через оценку математического ожидания, а другой – через оценку дисперсии. Математическое ожидание и дисперсия этого распределения равны ν/ λ и ν/ λ ² соответственно. Пусть их оценки определены равны:

α₁= 27,51, μ₂ = 0,91.

Составим систему уравнений для оцениваемых параметров

Разделив оценку математического ожидания на оценку дисперсии, получим

λ =α₁/μ₂ =30,12,

Метод моментов позволяет получить состоятельные, достаточные оценки, они при довольно общих условиях распределены асимптотически нормально. Смещение удается устранить введением поправок. Эффективность оценок невысокая, т.е. даже при больших объемах выборок дисперсия оценок относительно велика (за исключением нормального распределения, для которого метод моментов дает эффективные оценки). В реализации метод моментов проще метода максимального правдоподобия. Напомним, что метод целесообразно применять для оценки не более чем четырех параметров, так как точность выборочных моментов резко падает с увеличением их порядка.

 

Метод квантилей

 

Сущность метода квантилей схожа с методом моментов: выбирается столько квантилей, сколько требуется оценить параметров; неизвестные теоретические квантили, выраженные через параметры распределения, приравниваются к эмпирическим квантилям. Решение полученной системы уравнений дает искомые оценки параметров.

Дисперсия D(x_g) выборочной квантили обратно пропорциональна квадрату плотности распределения

D(x_g)=[g(1–g)]/[nf ²(x_g)]

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

Пример 1.5. Оценить методом квантилей параметры нормального распределения случайной величины.

Решение. Так как требуется определить два параметра распределения m и S, то выберем из вариационного ряда две эмпирические квантили. Например, можно взять

g₁ =5/44 =0,114; х_g1 = 26,13;

a_g =31/44=0,705; х_g2 = 28,01.

Используя стандартные функции математических пакетов, для выбранных значений g₁ и g₂ определим значения аргументов теоретической функции распределения для стандартизованной переменной

U_g1 = – 1, 207; u_g2 = 0,538.