Задачи: максимальное правдоподобие
Статистический анализ данных
Задачи
Задача 1. Вероятность и правдоподобие
Одно наблюдение \(x\) получено из экспоненциального распределения \[ f(x\mid\lambda)=\lambda e^{-\lambda x}, \qquad x\ge0. \]
- Нарисуйте \(f(x\mid\lambda)\) как функцию \(x\) при нескольких фиксированных \(\lambda\).
- После наблюдения \(x=x_{\mathrm{obs}}\) запишите \(L(\lambda;x_{\mathrm{obs}})\).
- Вычислите \(\int_0^\infty L(\lambda;x_{\mathrm{obs}})\,d\lambda\). Почему этот интеграл не является условием нормировки вероятности параметра?
- Найдите \(\widehat\lambda_{\mathrm{ММП}}\).
- Объясните словами различие между плотностью вероятности и функцией правдоподобия.
Задача 2. Биномиальная оценка
В \(N\) независимых испытаниях наблюдалось \(k\) успехов.
- Запишите биномиальное правдоподобие \(L(p)\).
- Найдите \(\widehat p_{\mathrm{ММП}}\).
- Вычислите \(\mathbb E[\widehat p]\) и \(\operatorname{Var}(\widehat p)\).
- Покажите, что оценка несмещённая и состоятельная.
- Исследуйте форму \(L(p)\) для случаев \(k=0\), \(k=N/2\) и \(k=N\).
Задача 3. Пуассоновский счёт с известным фоном
Пусть \[ N\sim\operatorname{Pois}(s+b), \] где \(s\ge0\) — неизвестный сигнал, а фон \(b\) известен точно.
- Запишите \(L(s)\) и \(\ell(s)\).
- Найдите неограниченную оценку максимального правдоподобия сигнала.
- Найдите оценку максимального правдоподобия при физическом ограничении \(s\ge0\).
- Объясните, что происходит при \(n<b\).
Задача 4. Гауссово среднее
Пусть \[ X_i\sim\mathrm N(\mu,\sigma^2), \] где \(\sigma\) известна.
- Запишите функцию правдоподобия независимой выборки.
- Покажите, что \(\widehat\mu_{\mathrm{ММП}}=\langle X\rangle\).
- Найдите \(\operatorname{Var}(\widehat\mu_{\mathrm{ММП}})\).
Задача 5. Взвешенное среднее
Имеются независимые измерения \[ X_i\sim\mathrm N(\mu,\sigma_i^2), \] где \(\sigma_i\) известны и различаются.
- Запишите совместный логарифм функции правдоподобия.
- Найдите оценку максимального правдоподобия параметра \(\mu\).
- Покажите, что результат является взвешенным средним с весами \(w_i=1/\sigma_i^2\).
- Найдите дисперсию оценки.
- Объясните, почему измерения с меньшей ошибкой получают больший вес.