Dvouparametrická rodina spojitých rozdělení pravděpodobnosti
v teorie pravděpodobnosti a statistika , inverzní rozdělení gama je dvojparametrová rodina spojitých rozdělení pravděpodobnosti na pozitivní skutečná linie , což je distribuce reciproční proměnné distribuované podle gama distribuce . Možná je hlavní využití inverzní gama distribuce Bayesovské statistiky , kde distribuce vzniká jako marginální zadní distribuce pro neznámé rozptyl a normální distribuce , pokud neinformativní předchozí se používá a je analyticky zpracovatelný před konjugátem , pokud je vyžadován informativní předchozí.
U Bayesianů je však běžné uvažovat o alternativě parametrizace z normální distribuce z hlediska přesnost , definovaný jako převrácená hodnota rozptylu, která umožňuje přímé použití gama distribuce jako konjugátu před. Jiní Bayesians preferují parametrizovat inverzní gama distribuci jinak, jak a škálované inverzní rozdělení chí-kvadrát .
Charakterizace Funkce hustoty pravděpodobnosti Inverzní rozdělení gama funkce hustoty pravděpodobnosti je definována přes Podpěra, podpora X > 0 { displaystyle x> 0}
F ( X ; α , β ) = β α Γ ( α ) ( 1 / X ) α + 1 exp ( − β / X ) { displaystyle f (x; alpha, beta) = { frac { beta ^ { alpha}} { Gamma ( alpha)}} (1 / x) ^ { alpha +1} exp vlevo (- beta / x vpravo)} s parametr tvaru α { displaystyle alpha} a parametr měřítka β { displaystyle beta} .[1] Tady Γ ( ⋅ ) { displaystyle gama ( cdot)} označuje funkce gama .
Na rozdíl od Distribuce gama , který obsahuje poněkud podobný exponenciální člen, β { displaystyle beta} je parametr měřítka, protože distribuční funkce splňuje:
F ( X ; α , β ) = F ( X / β ; α , 1 ) β { displaystyle f (x; alfa, beta) = { frac {f (x / beta; alfa, 1)} { beta}}} Funkce kumulativní distribuce The kumulativní distribuční funkce je regularizovaná funkce gama
F ( X ; α , β ) = Γ ( α , β X ) Γ ( α ) = Q ( α , β X ) { displaystyle F (x; alpha, beta) = { frac { gama vlevo ( alpha, { frac { beta} {x}} vpravo)} { gama ( alfa)}} = Q left ( alpha, { frac { beta} {x}} right) !} kde čitatel je horní neúplná funkce gama a jmenovatelem je funkce gama . Mnoho matematických balíčků umožňuje přímý výpočet Q { displaystyle Q} , legalizovaná funkce gama.
Okamžiky The n -tý okamžik inverzní distribuce gama je dán vztahem[2]
E [ X n ] = β n ( α − 1 ) ⋯ ( α − n ) . { displaystyle mathrm {E} [X ^ {n}] = { frac { beta ^ {n}} {( alpha -1) cdots ( alpha -n)}}} Charakteristická funkce K. α ( ⋅ ) { displaystyle K _ { alpha} ( cdot)} ve výrazu charakteristická funkce je upravený Besselova funkce druhého druhu.
Vlastnosti Pro α > 0 { displaystyle alpha> 0} a β > 0 { displaystyle beta> 0} ,
E [ ln ( X ) ] = ln ( β ) − ψ ( α ) { displaystyle mathbb {E} [ ln (X)] = ln ( beta) - psi ( alpha) ,} a
E [ X − 1 ] = α β , { displaystyle mathbb {E} [X ^ {- 1}] = { frac { alpha} { beta}}, ,} The informační entropie je
H ( X ) = E [ − ln ( p ( X ) ) ] = E [ − α ln ( β ) + ln ( Γ ( α ) ) + ( α + 1 ) ln ( X ) + β X ] = − α ln ( β ) + ln ( Γ ( α ) ) + ( α + 1 ) ln ( β ) − ( α + 1 ) ψ ( α ) + α = α + ln ( β Γ ( α ) ) − ( α + 1 ) ψ ( α ) . { displaystyle { begin {aligned} operatorname {H} (X) & = operatorname {E} [- ln (p (X))] & = operatorname {E} left [- alpha ln ( beta) + ln ( Gamma ( alpha)) + ( alpha +1) ln (X) + { frac { beta} {X}} vpravo] & = - alpha ln ( beta) + ln ( Gamma ( alpha)) + ( alpha +1) ln ( beta) - ( alpha +1) psi ( alpha) + alpha & = alpha + ln ( beta gama ( alpha)) - ( alpha +1) psi ( alpha). end {zarovnáno}}} kde ψ ( α ) { displaystyle psi ( alpha)} je funkce digamma .
The Kullback-Leiblerova divergence inverzní gama (αp , βp ) z inverzní gama (αq , βq ) je stejná jako divergence KL gama (αp , βp ) z Gamma (αq , βq ):
D K. L ( α p , β p ; α q , β q ) = E [ log ρ ( X ) π ( X ) ] = E [ log ρ ( 1 / Y ) π ( 1 / Y ) ] = E [ log ρ G ( Y ) π G ( Y ) ] , { displaystyle D _ { mathrm {KL}} ( alpha _ {p}, beta _ {p}; alpha _ {q}, beta _ {q}) = mathbb {E} vlevo [ 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],}
kde ρ , π { displaystyle rho, pi} jsou soubory PDF distribucí inverzní gama a ρ G , π G { displaystyle rho _ {G}, pi _ {G}} jsou soubory PDF distribucí gama, Y { displaystyle Y} je Gamma (αp , βp ) distribuován.
D K. L ( α p , β p ; α q , β q ) = ( α p − α q ) ψ ( α p ) − log Γ ( α p ) + log Γ ( α q ) + α q ( log β p − log β q ) + α p β q − β p β p . { displaystyle { begin {aligned} 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 {zarovnáno}}} Související distribuce Li X ∼ Inv-gama ( α , β ) { displaystyle X sim { mbox {Inv-gama}} ( alfa, beta)} pak k X ∼ Inv-gama ( α , k β ) { displaystyle kX sim { mbox {Inv-gama}} ( alpha, k beta) ,} Li X ∼ Inv-gama ( α , 1 2 ) { displaystyle X sim { mbox {Inv-gama}} ( alpha, { tfrac {1} {2}})} pak X ∼ Inv- χ 2 ( 2 α ) { displaystyle X sim { mbox {Inv -}} chi ^ {2} (2 alpha) ,} (inverzní-chi-kvadrát distribuce ) Li X ∼ Inv-gama ( α 2 , 1 2 ) { displaystyle X sim { mbox {Inv-Gamma}} ({ tfrac { alpha} {2}}, { tfrac {1} {2}})} pak X ∼ Scaled Inv- χ 2 ( α , 1 α ) { displaystyle X sim { mbox {Scaled Inv -}} chi ^ {2} ( alpha, { tfrac {1} { alpha}}) ,} (distribuce v měřítku-inverze-chi-kvadrát ) Li X ∼ Inv-gama ( 1 2 , C 2 ) { displaystyle X sim { textrm {Inv-Gamma}} ({ tfrac {1} {2}}, { tfrac {c} {2}})} pak X ∼ Levy ( 0 , C ) { displaystyle X sim { textrm {Levy}} (0, c) ,} (Lévyho distribuce ) Li X ∼ Inv-gama ( 1 , C ) { displaystyle X sim { textrm {Inv-gama}} (1, c)} pak 1 X ∼ Exp ( C ) { displaystyle { tfrac {1} {X}} sim { textrm {Exp}} (c) ,} (Exponenciální rozdělení ) Li X ∼ Gama ( α , β ) { displaystyle X sim { mbox {Gamma}} ( alfa, beta) ,} (Distribuce gama s hodnotit parametr β { displaystyle beta} ) pak 1 X ∼ Inv-gama ( α , β ) { displaystyle { tfrac {1} {X}} sim { mbox {Inv-gama}} ( alpha, beta) ,} (podrobnosti viz odvození v dalším odstavci) Všimněte si, že pokud X ~ Gamma (k , θ ) (Distribuce gama s parametrem měřítka θ ) pak 1 /X ~ Inv-gama (k , θ −1 ) Inverzní gama distribuce je speciální případ typu 5 Pearsonova distribuce A vícerozměrný zobecnění inverzní gama distribuce je inverzní Wishartova distribuce . Rozdělení součtu nezávislých invertovaných gama proměnných viz Witkovsky (2001) Odvození od gama distribuce Nechat X ∼ Gama ( α , β ) { displaystyle X sim { mbox {Gamma}} ( alfa, beta)} a připomenout, že pdf souboru gama distribuce je
F X ( X ) = β α Γ ( α ) X α − 1 E − β X { displaystyle f_ {X} (x) = { frac { beta ^ { alpha}} { Gamma ( alpha)}} x ^ { alpha -1} e ^ {- beta x}} , X > 0 { displaystyle x> 0} .Všimněte si, že β { displaystyle beta} je parametr rychlosti z hlediska rozdělení gama.
Definujte transformaci Y = G ( X ) = 1 X { displaystyle Y = g (X) = { tfrac {1} {X}}} . Pak pdf Y { displaystyle Y} je
F Y ( y ) = F X ( G − 1 ( y ) ) | d d y G − 1 ( y ) | = β α Γ ( α ) ( 1 y ) α − 1 exp ( − β y ) 1 y 2 = β α Γ ( α ) ( 1 y ) α + 1 exp ( − β y ) = β α Γ ( α ) ( y ) − α − 1 exp ( − β y ) { displaystyle { begin {aligned} 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 {zarovnáno}}} Všimněte si, že β { displaystyle beta} je parametr měřítka z hlediska inverzní distribuce gama.
Výskyt Viz také Reference Hoff, P. (2009). "První kurz bayesovských statistických metod". Springer. Witkovsky, V. (2001). "Výpočet distribuce lineární kombinace invertovaných gama proměnných". Kybernetika . 37 (1): 79–90. PAN 1825758 . Zbl 1263.62022 . Diskrétní univariate s konečnou podporou Diskrétní univariate s nekonečnou podporou Kontinuální univariate podporováno v omezeném intervalu Kontinuální univariate podporováno v poloneomezeném intervalu Kontinuální univariate podporováno na celé reálné linii Kontinuální univariate s podporou, jejíž typ se liší Smíšené spojité diskrétní univariate Vícerozměrný (společný) Směrový Degenerovat a jednotné číslo Rodiny