Двухпараметрическое семейство непрерывных распределений вероятностей
В теория вероятности и статистика , то обратное гамма-распределение является двухпараметрическим семейством непрерывных распределения вероятностей на позитиве реальная линия , которое является распределением взаимный переменной, распределенной согласно гамма-распределение . Возможно, основное использование обратного гамма-распределения заключается в Байесовская статистика , где распределение возникает как маргинальное апостериорное распределение для неизвестного отклонение из нормальное распределение , если малоинформативный приор используется, и как аналитически поддающийся обработке сопряженный предшествующий , если требуется информативная предварительная информация.
Однако среди байесовцев принято рассматривать альтернативу. параметризация из нормальное распределение с точки зрения точность , определяемый как величина, обратная дисперсии, что позволяет использовать гамма-распределение непосредственно в качестве сопряженного априорного значения. Другие байесовцы предпочитают иначе параметризовать обратное гамма-распределение, так как масштабированное обратное распределение хи-квадрат .
Характеристика
Функция плотности вероятности Обратное гамма-распределение функция плотности вероятности определяется над поддержка Икс > 0 { displaystyle x> 0}
ж ( Икс ; α , β ) = β α Γ ( α ) ( 1 / Икс ) α + 1 exp ( − β / Икс ) { Displaystyle е (х; альфа, бета) = { гидроразрыва { бета ^ { альфа}} { гамма ( альфа)}} (1 / х) ^ { альфа +1} ехр слева (- beta / x right)} с участием параметр формы α { displaystyle alpha} и масштабный параметр β { displaystyle beta} .[1] Вот Γ ( ⋅ ) { Displaystyle Gamma ( cdot)} обозначает гамма-функция .
в отличие от Гамма-распределение , который содержит несколько похожий экспоненциальный член, β { displaystyle beta} - масштабный параметр, так как функция распределения удовлетворяет:
ж ( Икс ; α , β ) = ж ( Икс / β ; α , 1 ) β { Displaystyle е (х; альфа, бета) = { гидроразрыва {е (х / бета; альфа, 1)} { бета}}} Кумулятивная функция распределения В кумулятивная функция распределения это регуляризованная гамма-функция
F ( Икс ; α , β ) = Γ ( α , β Икс ) Γ ( α ) = Q ( α , β Икс ) { Displaystyle F (х; альфа, бета) = { гидроразрыва { Gamma left ( alpha, { frac { beta} {x}} right)} { Gamma ( alpha)}} = Q left ( alpha, { frac { beta} {x}} right) !} где числитель - верхний неполная гамма-функция а знаменатель - это гамма-функция . Многие математические пакеты позволяют напрямую вычислять Q { displaystyle Q} , регуляризованная гамма-функция.
Моменты В п -й момент обратного гамма-распределения определяется выражением[2]
E [ Икс п ] = β п ( α − 1 ) ⋯ ( α − п ) . { displaystyle mathrm {E} [X ^ {n}] = { frac { beta ^ {n}} {( alpha -1) cdots ( alpha -n)}}.} Характеристическая функция K α ( ⋅ ) { Displaystyle К _ { альфа} ( cdot)} в выражении характеристическая функция это модифицированный Функция Бесселя 2-го рода.
Свойства
Для α > 0 { displaystyle alpha> 0} и β > 0 { displaystyle beta> 0} ,
E [ пер ( Икс ) ] = пер ( β ) − ψ ( α ) { Displaystyle mathbb {E} [ ln (X)] = ln ( бета) - psi ( альфа) ,} и
E [ Икс − 1 ] = α β , { displaystyle mathbb {E} [X ^ {- 1}] = { frac { alpha} { beta}}, ,} В информационная энтропия является
ЧАС ( Икс ) = E [ − пер ( п ( Икс ) ) ] = E [ − α пер ( β ) + пер ( Γ ( α ) ) + ( α + 1 ) пер ( Икс ) + β Икс ] = − α пер ( β ) + пер ( Γ ( α ) ) + ( α + 1 ) пер ( β ) − ( α + 1 ) ψ ( α ) + α = α + пер ( β Γ ( α ) ) − ( α + 1 ) ψ ( α ) . { displaystyle { begin {align} operatorname {H} (X) & = operatorname {E} [- ln (p (X))] & = operatorname {E} left [- alpha ln ( beta) + ln ( Gamma ( alpha)) + ( alpha +1) ln (X) + { frac { beta} {X}} right] & = - альфа ln ( beta) + ln ( Gamma ( alpha)) + ( alpha +1) ln ( beta) - ( alpha +1) psi ( alpha) + alpha & = alpha + ln ( beta Gamma ( alpha)) - ( alpha +1) psi ( alpha). end {выравнивается}}} где ψ ( α ) { Displaystyle psi ( альфа)} это функция дигаммы .
В Расхождение Кульбака-Лейблера обратной гаммы (αп , βп ) из обратной гаммы (αq , βq ) совпадает с KL-расходимостью Gamma (αп , βп ) из гаммы (αq , βq ):
D K L ( α п , β п ; α q , β q ) = E [ журнал ρ ( Икс ) π ( Икс ) ] = E [ журнал ρ ( 1 / Y ) π ( 1 / Y ) ] = E [ журнал ρ г ( Y ) π г ( Y ) ] , { displaystyle D _ { mathrm {KL}} ( alpha _ {p}, beta _ {p}; alpha _ {q}, beta _ {q}) = mathbb {E} left [ log { frac { rho (X)} { pi (X)}} right] = mathbb {E} left [ log { frac { rho (1 / Y)} { pi (1 / Y)}} right] = mathbb {E} left [ log { frac { rho _ {G} (Y)} { pi _ {G} (Y)}} right],}
где ρ , π { displaystyle rho, pi} являются pdf-распределениями обратного гамма-распределения и ρ г , π г { displaystyle rho _ {G}, pi _ {G}} - pdf-файлы гамма-распределений, Y { displaystyle Y} это Гамма (αп , βп ) распределены.
D K L ( α п , β п ; α q , β q ) = ( α п − α q ) ψ ( α п ) − журнал Γ ( α п ) + журнал Γ ( α q ) + α q ( журнал β п − журнал β q ) + α п β q − β п β п . { displaystyle { begin {align} D _ { mathrm {KL}} ( alpha _ {p}, beta _ {p}; alpha _ {q}, beta _ {q}) = {} & ( alpha _ {p} - alpha _ {q}) psi ( alpha _ {p}) - log Gamma ( alpha _ {p}) + log Gamma ( alpha _ {q} ) + alpha _ {q} ( log beta _ {p} - log beta _ {q}) + alpha _ {p} { frac { beta _ {q} - beta _ {p }} { beta _ {p}}}. end {выравнивается}}} Связанные дистрибутивы
Если Икс ∼ Инв-Гамма ( α , β ) { Displaystyle X sim { mbox {Inv-Gamma}} ( alpha, beta)} тогда k Икс ∼ Инв-Гамма ( α , k β ) { Displaystyle кХ сим { mbox {Inv-Gamma}} ( альфа, к бета) ,} Если Икс ∼ Инв-Гамма ( α , 1 2 ) { Displaystyle X sim { mbox {Inv-Gamma}} ( alpha, { tfrac {1} {2}})} тогда Икс ∼ Инв- χ 2 ( 2 α ) { displaystyle X sim { mbox {Inv -}} chi ^ {2} (2 alpha) ,} (обратное распределение хи-квадрат ) Если Икс ∼ Инв-Гамма ( α 2 , 1 2 ) { displaystyle X sim { mbox {Inv-Gamma}} ({ tfrac { alpha} {2}}, { tfrac {1} {2}})} тогда Икс ∼ Масштабированный инв- χ 2 ( α , 1 α ) { displaystyle X sim { mbox {Scaled Inv -}} chi ^ {2} ( alpha, { tfrac {1} { alpha}}) ,} (масштабированное обратное распределение хи-квадрат ) Если Икс ∼ Инв-Гамма ( 1 2 , c 2 ) { displaystyle X sim { textrm {Inv-Gamma}} ({ tfrac {1} {2}}, { tfrac {c} {2}})} тогда Икс ∼ Леви ( 0 , c ) { Displaystyle X sim { textrm {Леви}} (0, c) ,} (Распределение Леви ) Если Икс ∼ Инв-Гамма ( 1 , c ) { displaystyle X sim { textrm {Inv-Gamma}} (1, c)} тогда 1 Икс ∼ Опыт ( c ) { Displaystyle { tfrac {1} {X}} sim { textrm {Exp}} (с) ,} (Экспоненциальное распределение ) Если Икс ∼ Гамма ( α , β ) { Displaystyle X sim { mbox {Gamma}} ( alpha, beta) ,} (Гамма-распределение с участием показатель параметр β { displaystyle beta} ) тогда 1 Икс ∼ Инв-Гамма ( α , β ) { Displaystyle { tfrac {1} {X}} sim { mbox {Inv-Gamma}} ( alpha, beta) ,} (подробности см. в выводе в следующем абзаце) Обратите внимание, что если Икс ~ Гамма (k , θ ) (Гамма-распределение с масштабным параметром θ ), то 1 /Икс ~ Инв-Гамма (k , θ −1 ) Обратное гамма-распределение - частный случай типа 5. Распределение Пирсона А многомерный обобщением обратного гамма-распределения является обратное распределение Вишарта . О распределении суммы независимых инвертированных гамма-переменных см. Витковский (2001). Вывод из гамма-распределения
Позволять Икс ∼ Гамма ( α , β ) { Displaystyle X sim { mbox {Gamma}} ( alpha, beta)} , и напомним, что PDF-файл гамма-распределение является
ж Икс ( Икс ) = β α Γ ( α ) Икс α − 1 е − β Икс { displaystyle f_ {X} (x) = { frac { beta ^ { alpha}} { Gamma ( alpha)}} x ^ { alpha -1} e ^ {- beta x}} , Икс > 0 { displaystyle x> 0} .Обратите внимание, что β { displaystyle beta} - параметр скорости с точки зрения гамма-распределения.
Определите преобразование Y = г ( Икс ) = 1 Икс { Displaystyle Y = г (X) = { tfrac {1} {X}}} . Затем PDF-файл Y { displaystyle Y} является
ж Y ( y ) = ж Икс ( г − 1 ( y ) ) | d d y г − 1 ( y ) | = β α Γ ( α ) ( 1 y ) α − 1 exp ( − β y ) 1 y 2 = β α Γ ( α ) ( 1 y ) α + 1 exp ( − β y ) = β α Γ ( α ) ( y ) − α − 1 exp ( − β y ) { displaystyle { begin {align} f_ {Y} (y) & = f_ {X} left (g ^ {- 1} (y) right) left | { frac {d} {dy}} g ^ {- 1} (y) right | [6pt] & = { frac { beta ^ { alpha}} { Gamma ( alpha)}} left ({ frac {1} { y}} right) ^ { alpha -1} exp left ({ frac {- beta} {y}} right) { frac {1} {y ^ {2}}} [ 6pt] & = { frac { beta ^ { alpha}} { Gamma ( alpha)}} left ({ frac {1} {y}} right) ^ { alpha +1} exp left ({ frac {- beta} {y}} right) [6pt] & = { frac { beta ^ { alpha}} { Gamma ( alpha)}} left (y right) ^ {- alpha -1} exp left ({ frac {- beta} {y}} right) [6pt] end {align}}} Обратите внимание, что β { displaystyle beta} - масштабный параметр с точки зрения обратного гамма-распределения.
Вхождение
Смотрите также
использованная литература
Хофф, П. (2009). «Первый курс байесовских статистических методов». Springer. Витковский, В. (2001). «Вычисление распределения линейной комбинации инвертированных гамма-переменных». Кибернетика . 37 (1): 79–90. Г-Н 1825758 . Zbl 1263.62022 . Дискретный одномерный с конечной опорой Дискретный одномерный с бесконечной поддержкой Непрерывный одномерный поддерживается на ограниченном интервале Непрерывный одномерный поддерживается на полубесконечном интервале Непрерывный одномерный поддерживается на всей реальной линии Непрерывный одномерный с поддержкой, тип которой варьируется Смешанная непрерывно-дискретная одномерная Многовариантный (совместный) Направленный Вырожденный и единственное число Семьи