Метод моментов и максимальное правдоподобие
Два способа построить оценку параметров из статистической модели.
Перед чтением
Что вы разберёте в этой статье
Два способа построить оценку параметров из статистической модели.
Статистический вывод использует выборку, чтобы осторожно судить о всей совокупности. Любой такой вывод содержит неопределённость, которую нужно явно оценить.
- Метод моментов и максимальное правдоподобие
- estimation
- expectation-integral
Под каждой формулой дана расшифровка величин и условий. Сначала определите, что известно, затем проверьте ограничения и только после этого подставляйте числа.
После теории решите три задачи: простую, среднюю и сложную. Подсказка, шаги решения и ответ открываются отдельно.
Раздел: Статистические выводы · университетский уровень
Метод моментов
Если модель связывает параметры с математическими ожиданиями, замените теоретические моменты соответствующими выборочными средними и решите уравнения. Для одного параметра часто достаточно одного момента; для нескольких нужны достаточно информативные уравнения. Решение может отсутствовать, быть неединственным или выйти за область параметров.
Максимальное правдоподобие
Правдоподобие рассматривает фиксированные наблюдения как функцию параметра. В независимой выборке это произведение масс или плотностей наблюдений. Ищут допустимый параметр, при котором функция максимальна. Логарифм удобен, поскольку превращает произведение в сумму и не меняет расположение максимума.
Правдоподобие не является автоматически вероятностной плотностью параметра. Для байесовского апостериорного закона нужны априорное распределение и нормировка. Нельзя интерпретировать отношение значений плотности данных как уже готовую вероятность гипотезы.
Набросок вывода для Бернулли
При k успехах в n испытаниях логарифм правдоподобия равен k·log p+(n−k)·log(1−p). Внутренняя критическая точка p=k/n максимальна; при k=0 или k=n максимум находится на границе [0,1]. Метод моментов даёт то же решение, потому что E X=p.
Условия и границы
Обещания состоятельности и асимптотической нормальности MLE требуют идентифицируемости и других условий регулярности. Граница параметра или зависящий от параметра носитель может менять выводы. Например, у равномерного [0,θ] нельзя искать максимум только по производной внутри области: ограничение θ≥max xi существенно.
Три примера
- В 20 опытах Бернулли 7 успехов: и метод моментов, и MLE дают p̂=7/20=0,35.
- Для экспоненциальных ожиданий со средним выборки 4 оба метода дают интенсивность λ̂=1/4, а не 4.
- Для равномерного [0,θ] и данных 1,2,4 метод моментов даёт 2·7/3=14/3, а MLE — максимум выборки 4. Методы могут различаться.
Формула, смысл и ограничения
Главные формулы с расшифровкой
Формула словами: Произведение относится к независимой выборке; fθ обозначает плотность либо массу. Для зависимых данных нужна совместная функция.
Как применять: Разделите наблюдаемые данные, неизвестный параметр и выбранный уровень уверенности; проверьте предпосылки метода.
От формулы к наблюдению
Попробуйте самостоятельно
Измените параметры и сопоставьте расчёт с результатами случайного эксперимента.
От простого к сложному
Проверьте себя
Сначала запишите решение самостоятельно. Затем можно открыть намёк, сравнить каждый переход с пошаговым разбором и только после этого посмотреть ответ. Прогресс сохраняется в этом браузере.
В 80 испытаниях Бернулли 28 успехов. Найдите MLE вероятности успеха.
Принимаются равные дроби и десятичная запись. Бесконечные дроби округляйте минимум до 6 знаков. Для процентов следуйте условию.
В независимой экспоненциальной выборке среднее равно 2,5. Найдите MLE интенсивности λ.
Принимаются равные дроби и десятичная запись. Бесконечные дроби округляйте минимум до 6 знаков. Для процентов следуйте условию.
Наблюдения 2,3,7 независимы и равномерны на [0,θ]. Найдите MLE θ.
Принимаются равные дроби и десятичная запись. Бесконечные дроби округляйте минимум до 6 знаков. Для процентов следуйте условию.
Источники и соглашения
Объяснения и задачи — авторские. Источники помогают проверить определения, условия и школьную привязку. Обозначения и параметризация указаны в статье.