Несмещенная оценка минимальной дисперсии

редактировать
Несмещенная статистическая оценка, минимизирующая дисперсию

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

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

Хотя сочетание ограничения непредвзятости с метрикой желательности наименьшей дисперсии приводит к хорошим результатам в большинстве практических ситуаций, что делает MVUE естественной отправной точкой для широкого диапазона анализа - целевая спецификация может работать лучше для данной проблемы; таким образом, MVUE не всегда является лучшей точкой остановки.

Содержание
  • 1 Определение
  • 2 Выбор оценщика
  • 3 Пример
  • 4 Другие примеры
  • 5 См. Также
    • 5.1 Байесовские аналоги
  • 6 Ссылки
Определение

Рассмотрим оценку g (θ) {\ displaystyle g (\ theta)}g (\ theta) на основе данных X 1, X 2,…, X n {\ displaystyle X_ {1 }, X_ {2}, \ ldots, X_ {n}}X_ {1}, X_ {2}, \ ldots, X_ {n} iid от некоторого члена семейства плотностей p θ, θ ∈ Ω {\ displaystyle p _ {\ theta}, \ theta \ in \ Omega}p _ {\ theta}, \ theta \ in \ Omega , где Ω {\ displaystyle \ Omega }\ Omega - это пространство параметров. Несмещенная оценка δ (X 1, X 2,…, X n) {\ displaystyle \ delta (X_ {1}, X_ {2}, \ ldots, X_ {n})}\ delta (X_ {1}, X_ {2}, \ ldots, X_ {n}) из g (θ) {\ displaystyle g (\ theta)}g (\ theta) является UMVUE, если ∀ θ ∈ Ω {\ displaystyle \ forall \ theta \ in \ Omega}\ forall \ theta \ in \ Omega ,

var ⁡ (δ (Икс 1, Икс 2,…, Икс n)) ≤ var ⁡ (δ ~ (X 1, X 2,…, X n)) {\ displaystyle \ operatorname {var} (\ delta (X_ {1}, X_ {2}, \ ldots, X_ {n})) \ leq \ operatorname {var} ({\ tilde {\ delta}} (X_ {1}, X_ {2}, \ ldots, X_ {n}))}{\ displaystyle \ operatorname {var} (\ delta (X_ {1}, X_ {2}, \ ldots, X_ {n})) \ leq \ operatorname {var} ({\ tilde {\ delta}} (X_ {1}, X_ {2}, \ ldots, X_ {n}))}

для любой другой несмещенной оценки δ ~. {\ displaystyle {\ tilde {\ delta}}.}{\ tilde {\ delta}}.

Если существует несмещенная оценка g (θ) {\ displaystyle g (\ theta)}g (\ theta) , то там можно доказать по сути уникальный MVUE. Используя теорему Рао – Блэквелла, можно также доказать, что определение MVUE - это просто вопрос поиска полной достаточной статистики для семейства p θ, θ ∈ Ω {\ displaystyle p _ {\ theta}, \ theta \ in \ Omega}p _ {\ theta}, \ theta \ in \ Omega и обусловливает на нем любую несмещенную оценку.

Далее, согласно теореме Лемана – Шеффе, несмещенная оценка, которая является функцией полной, достаточной статистики, является оценкой UMVUE.

Формально предположим, что δ (X 1, X 2,…, X n) {\ displaystyle \ delta (X_ {1}, X_ {2}, \ ldots, X_ {n}) }\ delta (X_ {1}, X_ {2}, \ ldots, X_ {n}) несмещен для g (θ) {\ displaystyle g (\ theta)}g (\ theta) , и что T {\ displaystyle T}T является полной достаточной статистикой для семейства плотностей. Тогда

η (Икс 1, Икс 2,…, Икс n) = E ⁡ (δ (X 1, X 2,…, X n) ∣ T) {\ displaystyle \ eta (X_ {1}, X_ { 2}, \ ldots, X_ {n}) = \ operatorname {E} (\ delta (X_ {1}, X_ {2}, \ ldots, X_ {n}) \ mid T) \,}{\ displaystyle \ eta (X_ {1}, X_ {2}, \ ldots, X_ {n}) = \ имя оператора {E} (\ delta (X_ {1}, X_ {2}, \ ldots, X_ {n}) \ mid T) \,}

MVUE для g (θ). {\ displaystyle g (\ theta).}g (\ theta).

A Байесовский аналог - это байесовская оценка, особенно с минимальной среднеквадратичной ошибкой (MMSE).

Выбор оценщика

эффективный оценщик может не существовать, но если он существует и если он несмещен, то это MVUE. Поскольку среднеквадратичная ошибка (MSE) оценщика δ равна

MSE ⁡ (δ) = var ⁡ (δ) + [bias ⁡ (δ)] 2 {\ displaystyle \ operatorname {MSE} (\ delta) = \ operatorname {var} (\ delta) + [\ operatorname {bias} (\ delta)] ^ {2} \}{\ displaystyle \ operatorname {MSE} (\ delta) = \ operatorname {var} (\ delta) + [\ operatorname {bias} (\ delta)] ^ {2} \}

MVUE минимизирует MSE среди несмещенных оценок. В некоторых случаях смещенные оценки имеют более низкую MSE, потому что они имеют меньшую дисперсию, чем любая несмещенная оценка; см. погрешность оценки.

Пример

Считайте данные единым наблюдением из абсолютно непрерывного распределения на R {\ displaystyle \ mathbb {R}}{\ mathbb {R}} с плотностью

п θ (x) = θ e - x (1 + e - x) θ + 1 {\ displaystyle p _ {\ theta} (x) = {\ frac {\ theta e ^ {-x}} {(1 + e ^ {- x}) ^ {\ theta +1}}}}p _ {\ theta} (x) = {\ frac {\ theta e ^ {{- x}}} {(1 + e ^ {{- x}}) ^ {{\ theta +1}}} }

, и мы хотим найти оценку UMVU для

g (θ) = 1 θ 2 { \ displaystyle g (\ theta) = {\ frac {1} {\ theta ^ {2}}}}{ \ displaystyle g (\ theta) = {\ frac {1} {\ theta ^ {2}}}}

Сначала мы понимаем, что плотность можно записать как

e - x 1 + e - x exp ⁡ (- θ журнал ⁡ (1 + e - x) + журнал ⁡ (θ)) {\ displaystyle {\ frac {e ^ {- x}} {1 + e ^ {- x}}} \ exp (- \ theta \ log (1 + e ^ {- x}) + \ log (\ theta))}{\ frac {e ^ {{- x}}} {1 + e ^ {{- x}}}}} \ exp (- \ theta \ log ( 1 + е ^ {{- x}}) + \ log (\ theta))

Это экспоненциальное семейство с достаточной статистикой T = log ⁡ (1 + e - х) {\ displaystyle T = \ log (1 + e ^ {- x})}{ \ displaystyle T = \ log (1 + e ^ {- x})} . Фактически это экспоненциальное семейство полного ранга, и поэтому T {\ displaystyle T}T вполне достаточно. См. экспоненциальное семейство для вывода, который показывает

E ⁡ (T) = 1 θ, var ⁡ (T) = 1 θ 2 {\ displaystyle \ operatorname {E} (T) = {\ frac {1} {\ theta}}, \ quad \ operatorname {var} (T) = {\ frac {1} {\ theta ^ {2}}}}{\ displaystyle \ operatorname {E } (T) = {\ frac {1} {\ theta}}, \ quad \ operatorname {var} (T) = {\ frac {1} {\ theta ^ {2}}}}

Следовательно,

E ⁡ (T 2) = 2 θ 2 {\ displaystyle \ operatorname {E} (T ^ {2}) = {\ frac {2} {\ theta ^ {2}}}}{\ displaystyle \ operatorname {E} (T ^ {2}) = {\ гидроразрыва {2} {\ theta ^ {2}}}}

Здесь мы используем теорему Лемана – Шеффе, чтобы получить MVUE

Очевидно δ (X) = T 2 2 {\ displaystyle \ delta (X) = {\ frac {T ^ {2}} {2}}}{\ displaystyle \ delta (X) = {\ frac {T ^ {2}} {2}}} беспристрастно и T = log ⁡ (1 + e - x) {\ displaystyle T = \ log (1 + e ^ {- x})}{ \ displaystyle T = \ log (1 + e ^ {- x})} достаточно полно, поэтому оценка UMVU

η (Икс) знак равно Е ⁡ (δ (Икс) ∣ T) знак равно E ⁡ (T 2 2 | T) = T 2 2 = журнал ⁡ (1 + е - X) 2 2 {\ Displaystyle \ eta (X) = \ operatorname {E} (\ delta (X) \ mid T) = \ operatorname {E} \ left (\ left. {\ frac {T ^ {2}} {2}} \, \ right | \, T \ right) = {\ frac {T ^ {2}} {2}} = {\ frac {\ log (1 + e ^ {- X}) ^ {2}} {2}}}{\ displaystyle \ eta (X) = \ operatorname {E} (\ delta (X) \ mid T) = \ operatorname {E} \ left (\ left. {\ frac {T ^ {2}} {2}} \, \ right | \, T \ right) = {\ frac {T ^ {2}} {2}} = {\ frac {\ log (1 + e ^ {- X}) ^ {2}} {2 }}}

Этот пример показывает, что несмещенная функция полной достаточной статистики будет UMVU, поскольку Lehm ann – теорема Шеффе утверждает.

Другие примеры
k + 1 км - 1, { \ displaystyle {\ frac {k + 1} {k}} m-1,}{\ frac {k + 1} {k} } м-1,
, где m - максимум выборки. Это масштабированное и сдвинутое (так несмещенное) преобразование максимума выборки, которое является достаточной и полной статистикой. Подробнее см. Проблема немецкого танка.
См. Также

Байесовские аналоги

Ссылки
  • Кинер, Роберт В. (2006). Статистическая теория: Заметки к курсу теоретической статистики. Springer. С. 47–48, 57–58.
  • Воинов В. Г., Никулин М. С. (1993). Беспристрастные оценки и их приложения, Том 1: Одномерный случай. Kluwer Academic Publishers. стр. 521p.
Последняя правка сделана 2021-05-30 13:16:32
Содержание доступно по лицензии CC BY-SA 3.0 (если не указано иное).
Обратная связь: support@alphapedia.ru
Соглашение
О проекте