WikiDer > Серия Эджворта

Edgeworth series

В Грам – Шарлье серия А (назван в честь Йорген Педерсен Грам и Карл Шарлье), а Серия Эджворта (назван в честь Фрэнсис Исидро Эджворт) находятся серии что приблизительно распределение вероятностей с точки зрения его кумулянты.[1] Серии такие же; но порядок терминов (и, следовательно, точность усечения ряда) различается.[2] Ключевая идея этих расширений - написать характеристическая функция распределения, чьи функция плотности вероятности ж должен быть аппроксимирован характеристической функцией распределения с известными и подходящими свойствами, и восстановить ж через обратный преобразование Фурье.

Грам – Шарлье серия А

Мы исследуем непрерывную случайную величину. Позволять - характеристическая функция его распределения, функция плотности которой ж, и его кумулянты. Разложим по известному распределению с функцией плотности вероятности ψ, характеристическая функция , и кумулянты . Плотность ψ обычно выбирается нормальное распределение, но возможны и другие варианты. По определению кумулянтов мы имеем (см. Wallace, 1958)[3]

и

что дает следующую формальную идентичность:

По свойствам преобразования Фурье - преобразование Фурье , где D это дифференциальный оператор относительно Икс. Таким образом, после изменения с участием по обе стороны уравнения, находим для ж формальное расширение

Если ψ выбрана нормальная плотность

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

поскольку для всех р > 2, поскольку старшие кумулянты нормального распределения равны 0. Разложив экспоненту и собрав члены в соответствии с порядком производных, мы приходим к ряду Грама – Шарлье. Такое разложение можно компактно записать в терминах Полиномы Белла так как

Поскольку n-я производная гауссовой функции дается с точки зрения Многочлен Эрмита так как

это дает нам окончательное выражение ряда Грама – Шарлье А как

Интеграция серии дает нам кумулятивная функция распределения

где - функция распределения нормального распределения.

Если мы включим только первые два поправочных члена в нормальное распределение, мы получим

с участием и .

Обратите внимание, что положительное значение этого выражения не гарантируется, и поэтому оно не является допустимым распределением вероятностей. Ряд Грама – Шарлье расходится во многих интересных случаях - он сходится, только если падает быстрее, чем на бесконечности (Cramér 1957). Когда он не сходится, ряд тоже не соответствует действительности. асимптотическое разложение, потому что оценить погрешность разложения невозможно. По этой причине серия Эджворта (см. Следующий раздел) обычно предпочтительнее серии Грама – Шарлье.

Серия Эджворта

Эджворт разработал подобное расширение как усовершенствование Центральная предельная теорема.[4] Преимущество серии Эджворта в том, что ошибка контролируется, так что это истинное асимптотическое разложение.

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

Позволять обозначить кумулятивные функции распределения переменных . Тогда по центральной предельной теореме

для каждого , пока среднее значение и дисперсия конечны.

Теперь предположим, что помимо среднего и дисперсия , i.i.d. случайные переменные иметь более высокие кумулянты . Из свойств аддитивности и однородности кумулянтов кумулянты с точки зрения кумулянтов для ,

Если мы расширим с точки зрения стандартного нормального распределения, то есть если мы установим

то кумулянтные разности формального выражения характеристической функции из находятся

Ряд Грама – Шарлье для функции плотности сейчас

Серия Эджворта разработана аналогично серии А Грама – Шарлье, только теперь термины собираются в соответствии со степенями . Коэффициенты при п-м / 2 член может быть получен путем сбора мономов многочленов Белла, соответствующих целочисленным разбиениям м. Таким образом, мы имеем характеристическую функцию как

где это многочлен степени . Опять же, после обратного преобразования Фурье функция плотности следует как

Аналогично, интегрируя ряд, получаем функцию распределения

Мы можем явно написать многочлен так как

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

Например, если м = 3, то есть три способа разделить это число: 1 + 1 + 1 = 2 + 1 = 3. Таким образом, нам нужно изучить три случая:

  • 1 + 1 + 1 = 1 · k1, так что у нас есть k1 = 3, л1 = 3 и s = 9.
  • 1 + 2 = 1 · k1 + 2 · k2, так что у нас есть k1 = 1, k2 = 1, л1 = 3, л2 = 4 и s = 7.
  • 3 = 3 · k3, так что у нас есть k3 = 1, л3 = 5 и s = 5.

Таким образом, искомый полином равен

Первые пять условий расширения:[5]

Вот, φ(j)(Икс) это j-я производная от φ (·) в точке Икс. Помня, что производные плотности нормального распределения связаны с нормальной плотностью соотношением , (где это Многочлен Эрмита порядка п), это объясняет альтернативные представления в терминах функции плотности. Блинников и Месснер (1998) предложили простой алгоритм для вычисления членов разложения более высокого порядка.

Обратите внимание, что в случае решетчатых распределений (которые имеют дискретные значения), расширение Эджворта должно быть скорректировано для учета разрывных скачков между точками решетки.[6]

Иллюстрация: плотность выборочного среднего из трех

Плотность выборочного среднего трех переменных chi2. На диаграмме сравниваются истинная плотность, нормальное приближение и два разложения Эджворта.

Взять и выборочное среднее .

Мы можем использовать несколько дистрибутивов для :

  • Точное распределение, которое следует гамма-распределение: .
  • Асимптотическое нормальное распределение: .
  • Два разложения Эджворта степени 2 и 3.

Обсуждение результатов

  • Для конечных выборок не гарантируется, что разложение Эджворта будет правильным. распределение вероятностей поскольку значения CDF в некоторых точках могут выходить за рамки .
  • Они гарантируют (асимптотически) абсолютные ошибки, но относительные ошибки можно легко оценить, сравнив ведущий член Эджворта в остатке с общим ведущим членом. [7]

Смотрите также

использованная литература

  1. ^ Стюарт А. и Кендалл М. Г. (1968). Продвинутая теория статистики. Издательство Hafner.
  2. ^ Коласса, Дж. Э. (2006). Методы приближения рядов в статистике (Том 88). Springer Science & Business Media.
  3. ^ Уоллес, Д. Л. (1958). «Асимптотические приближения распределений». Анналы математической статистики. 29 (3): 635–654. Дои:10.1214 / aoms / 1177706528. JSTOR 2237255.
  4. ^ Холл, П. (2013). Бутстрап и расширение Эджворта. Springer Science & Business Media.
  5. ^ Вайсштейн, Эрик В. "Серия Эджворта". MathWorld.
  6. ^ Коласса, Джон Э .; Маккаллах, Питер (1990). «Ряд Эджворта для решетчатых распределений». Анналы статистики. 18 (2): 981–985. Дои:10.1214 / aos / 1176347637. JSTOR 2242145.
  7. ^ Коласса, Джон Э. (2006). Методы приближения рядов в статистике (3-е изд.). Springer. ISBN 0387322272.

дальнейшее чтение