Оценка параметров нормального распределения: медиана вместо среднего?

15

Общий подход для оценки параметров нормального распределения заключается в использовании среднего значения и стандартного отклонения / дисперсии выборки.

Однако, если есть некоторые выбросы, медиана и срединное отклонение от медианы должны быть намного более устойчивыми, верно?

На некоторых наборах данных я пытался, нормальное распределение , оцененное N(median(x),median|xmedian(x)|) , кажется, производит гораздо лучшее прилегание , чем классический N(μ^,σ^) , используя среднее и Среднеквадратичное отклонение

Есть ли причина не использовать медиану, если вы предполагаете, что в наборе данных есть некоторые выбросы? Знаете ли вы некоторые ссылки для этого подхода? Быстрый поиск в Google не нашел мне полезных результатов, которые обсуждают преимущества использования медиан здесь (но очевидно, что «медиана оценки параметров нормального распределения» не очень специфический набор поисковых терминов).

Срединное отклонение, это предвзято? Должен ли я умножить его на чтобы уменьшить смещение?n1n

Знаете ли вы похожие подходы надежной оценки параметров для других распределений, таких как гамма-распределение или экспоненциально модифицированное гауссово распределение (которому требуется асимметрия в оценке параметров, а выбросы действительно портят это значение)?

Эрих Шуберт
источник
2
Если у вас есть выбросы, возможно, ваше распределение не является гауссовским нормальным. Конечно, это не отвечает на ваш вопрос, но, IMO, это возможность, которую всегда нужно развлекать.
SdS
2
У меня нет простого, чистого, математического распределения. У меня есть реальные данные, которые по природе грязны. Никакое распределение не будет идеально подходить, потому что вы не можете больше справляться с ситуацией аналитически. И выбросы на самом деле мой интерес. :-)
Эрих Шуберт

Ответы:

15

Замечание, что в примере с данными, полученными из загрязненного гауссовского распределения, вы получите более точные оценки параметров, описывающих объем данных, используя вместо med | х - мед ( х ) | где безумный ( х ) это:madmed|xmed(x)|mad(x)

mad=1.4826×med|xmed(x)|

- где, - это коэффициент согласованности, разработанный для обеспечения того, чтобы E ( mad ( x ) 2 ) = Var ( x ), когда х не загрязнен, - изначально было сделано Гауссом (Уокер) H. (1931)).(Φ1(0.75))1=1.4826

E(mad(x)2)=Var(x)
x

Я не могу думать о какой - либо причине не использовать вместо выборочного среднего в данном случае. Более низкая эффективность (по Гауссу!) Сумасшедших может быть причиной не использовать сумасшедшие в вашем примере. Однако существуют такие же надежные и высокоэффективные альтернативы безумным . Одним из них является Q nmedmadmadmadQn, Этот оценщик имеет много других преимуществ, кроме. Он также очень нечувствителен к выбросам (фактически почти так же нечувствителен, как и к безумным). Вопреки безумному, он не построен вокруг оценки местоположения и не предполагает, что распределение незагрязненной части данных симметрично. Как и сумасшедший, он основан на статистике заказов, так что он всегда четко определен, даже если в исходном распределении вашей выборки нет моментов. Как и безумный, он имеет простую явную форму. Даже больше, чем для сумасшедших, я не вижу причин использовать стандартное отклонение образца вместо в примере, который вы описываете (см. Rousseeuw и Croux 1993 для получения дополнительной информации о Q n ).QnQn

xΓ(ν,λ)

med(x)λ(ν1/3)

и

mad(x)λν

ν>1.5

ν^=(med(x)mad(x))2

и

λ^=mad(x)2med(x)

См. Чен и Рубин (1986) для полного вывода.

  • Дж. Чен и Х. Рубин, 1986. Границы для различия между медианой и средней величиной гамма-распределения и распределения Пуассона, Statist. Вероятно. Lett., 4, 281–283.
  • PJ Rousseeuw и C. Croux, 1993. Альтернативы журналу «Среднее абсолютное отклонение» Американской статистической ассоциации, Vol. 88, № 424, с. 1273-1283
  • Уокер, Х. (1931). Исследования по истории статистического метода. Балтимор, MD: Williams & Wilkins Co. стр. 24–25.
user603
источник
1
- это значение, которое нужно использовать, или одна из двух дополнительных инверсий? Φ1(0.75)11.4826
Эрих Шуберт
@ErichSchubert: вы правы: я забыл второе обратное .. исправлено.
user603
2
n/(n1)
1
@whuber: спасибо за это, теперь я понимаю, что мое предложение «это похоже по духу » может быть легко понято. Я удалил это.
user603
1
Я сделал часть ExNormal отдельным вопросом: stats.stackexchange.com/questions/48907/… Но у меня есть еще один вопрос: распределение LogNormal - обработать, применив журнал, затем продолжить как обычное распределение?
Эрих Шуберт
7

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

Если бы вы знали некоторые априорные ограничения на долю выбросов, вы могли бы обрезать эту пропорцию до среднего и Winsorize стандартное отклонение. Альтернативой, которая не требует таких знаний, будет использование M-оценок для определения местоположения и связанных величин для дисперсии. Повышение эффективности, если ваши предположения верны (например, данные действительно нормальные, за исключением небольшого процента выбросов), могут в некоторых случаях быть существенными.

Среднее отклонение смещено как оценка стандартного отклонения - но не какnn1

Glen_b - Восстановить Монику
источник