J u m p t o c o n t e n t
M a i n m e n u
M a i n m e n u
N a v i g a t i o n
● M a i n p a g e
● C o n t e n t s
● C u r r e n t e v e n t s
● R a n d o m a r t i c l e
● A b o u t W i k i p e d i a
● C o n t a c t u s
● D o n a t e
C o n t r i b u t e
● H e l p
● L e a r n t o e d i t
● C o m m u n i t y p o r t a l
● R e c e n t c h a n g e s
● U p l o a d f i l e
S e a r c h
Search
A p p e a r a n c e
● C r e a t e a c c o u n t
● L o g i n
P e r s o n a l t o o l s
● C r e a t e a c c o u n t
● L o g i n
P a g e s f o r l o g g e d o u t e d i t o r s l e a r n m o r e
● C o n t r i b u t i o n s
● T a l k
( T o p )
1
C h a r a c t e r i s t i c s
2
R e l a t e d d i s t r i b u t i o n s
3
M o m e n t s
4
P r o p e r t i e s
5
K u l l b a c k - L e i b l e r d i v e r g e n c e
6
S o f t w a r e i m p l e m e n t a t i o n
7
S e e a l s o
8
R e f e r e n c e s
T o g g l e t h e t a b l e o f c o n t e n t s
G e n e r a l i z e d g a m m a d i s t r i b u t i o n
3 l a n g u a g e s
● ا ل ع ر ب ي ة
● ف ا ر س ی
● F r a n ç a i s
E d i t l i n k s
● A r t i c l e
● T a l k
E n g l i s h
● R e a d
● E d i t
● V i e w h i s t o r y
T o o l s
T o o l s
A c t i o n s
● R e a d
● E d i t
● V i e w h i s t o r y
G e n e r a l
● W h a t l i n k s h e r e
● R e l a t e d c h a n g e s
● U p l o a d f i l e
● S p e c i a l p a g e s
● P e r m a n e n t l i n k
● P a g e i n f o r m a t i o n
● C i t e t h i s p a g e
● G e t s h o r t e n e d U R L
● D o w n l o a d Q R c o d e
● W i k i d a t a i t e m
P r i n t / e x p o r t
● D o w n l o a d a s P D F
● P r i n t a b l e v e r s i o n
A p p e a r a n c e
F r o m W i k i p e d i a , t h e f r e e e n c y c l o p e d i a
Generalized gamma
Probability density function
Parameters
a
>
0
{\displaystyle a>0}
(scale),
d
,
p
>
0
{\displaystyle d,p>0}
Support
x
∈
(
0
,
∞
)
{\displaystyle x\;\in \;(0,\,\infty )}
PDF
p
/
a
d
Γ
(
d
/
p
)
x
d
−
1
e
−
(
x
/
a
)
p
{\displaystyle {\frac {p/a^{d}}{\Gamma (d/p)}}x^{d-1}e^{-(x/a)^{p}}}
CDF
γ
(
d
/
p
,
(
x
/
a
)
p
)
Γ
(
d
/
p
)
{\displaystyle {\frac {\gamma (d/p,(x/a)^{p})}{\Gamma (d/p)}}}
Mean
a
Γ
(
(
d
+
1
)
/
p
)
Γ
(
d
/
p
)
{\displaystyle a{\frac {\Gamma ((d+1)/p)}{\Gamma (d/p)}}}
Mode
a
(
d
−
1
p
)
1
p
f
o
r
d
>
1
,
o
t
h
e
r
w
i
s
e
0
{\displaystyle a\left({\frac {d-1}{p}}\right)^{\frac {1}{p}}\mathrm {for} \;d>1,\mathrm {otherwise} \;0}
Variance
a
2
(
Γ
(
(
d
+
2
)
/
p
)
Γ
(
d
/
p
)
−
(
Γ
(
(
d
+
1
)
/
p
)
Γ
(
d
/
p
)
)
2
)
{\displaystyle a^{2}\left({\frac {\Gamma ((d+2)/p)}{\Gamma (d/p)}}-\left({\frac {\Gamma ((d+1)/p)}{\Gamma (d/p)}}\right)^{2}\right)}
Entropy
ln
a
Γ
(
d
/
p
)
p
+
d
p
+
(
1
p
−
d
p
)
ψ
(
d
p
)
{\displaystyle \ln {\frac {a\Gamma (d/p)}{p}}+{\frac {d}{p}}+\left({\frac {1}{p}}-{\frac {d}{p}}\right)\psi \left({\frac {d}{p}}\right)}
The generalized gamma distribution is a continuous probability distribution with two shape parameters (and a scale parameter ). It is a generalization of the gamma distribution which has one shape parameter (and a scale parameter). Since many distributions commonly used for parametric models in survival analysis (such as the exponential distribution , the Weibull distribution and the gamma distribution ) are special cases of the generalized gamma, it is sometimes used to determine which parametric model is appropriate for a given set of data.[1] Another example is the half-normal distribution .
Characteristics [ edit ]
The generalized gamma distribution has two shape parameters ,
d
>
0
{\displaystyle d>0}
and
p
>
0
{\displaystyle p>0}
, and a scale parameter ,
a
>
0
{\displaystyle a>0}
. For non-negative x from a generalized gamma distribution, the probability density function is [2]
f
(
x
;
a
,
d
,
p
)
=
(
p
/
a
d
)
x
d
−
1
e
−
(
x
/
a
)
p
Γ
(
d
/
p
)
,
{\displaystyle f(x;a,d,p)={\frac {(p/a^{d})x^{d-1}e^{-(x/a)^{p}}}{\Gamma (d/p)}},}
where
Γ
(
⋅
)
{\displaystyle \Gamma (\cdot )}
denotes the gamma function .
The cumulative distribution function is
F
(
x
;
a
,
d
,
p
)
=
γ
(
d
/
p
,
(
x
/
a
)
p
)
Γ
(
d
/
p
)
,
or
P
(
d
p
,
(
x
a
)
p
)
;
{\displaystyle F(x;a,d,p)={\frac {\gamma (d/p,(x/a)^{p})}{\Gamma (d/p)}},{\text{or}}\,P\left({\frac {d}{p}},\left({\frac {x}{a}}\right)^{p}\right);}
where
γ
(
⋅
)
{\displaystyle \gamma (\cdot )}
denotes the lower incomplete gamma function ,
and
P
(
⋅
,
⋅
)
{\displaystyle P(\cdot ,\cdot )}
denotes the regularized lower incomplete gamma function .
The quantile function can be found by noting that
F
(
x
;
a
,
d
,
p
)
=
G
(
(
x
/
a
)
p
)
{\displaystyle F(x;a,d,p)=G((x/a)^{p})}
where
G
{\displaystyle G}
is the cumulative distribution function of the gamma distribution with parameters
α
=
d
/
p
{\displaystyle \alpha =d/p}
and
β
=
1
{\displaystyle \beta =1}
. The quantile function is then given by inverting
F
{\displaystyle F}
using known relations about inverse of composite functions , yielding:
F
−
1
(
q
;
a
,
d
,
p
)
=
a
⋅
[
G
−
1
(
q
)
]
1
/
p
,
{\displaystyle F^{-1}(q;a,d,p)=a\cdot {\big [}G^{-1}(q ){\big ]}^{1/p},}
with
G
−
1
(
q
)
{\displaystyle G^{-1}(q )}
being the quantile function for a gamma distribution with
α
=
d
/
p
,
β
=
1
{\displaystyle \alpha =d/p,\,\beta =1}
.
Related distributions [ edit ]
Alternative parameterisations of this distribution are sometimes used; for example with the substitution α = d/p .[3] In addition, a shift parameter can be added, so the domain of x starts at some value other than zero.[3] If the restrictions on the signs of a , d and p are also lifted (but α = d /p remains positive), this gives a distribution called the Amoroso distribution , after the Italian mathematician and economist Luigi Amoroso who described it in 1925.[4]
Moments [ edit ]
If X has a generalized gamma distribution as above, then[3]
E
(
X
r
)
=
a
r
Γ
(
d
+
r
p
)
Γ
(
d
p
)
.
{\displaystyle \operatorname {E} (X^{r})=a^{r}{\frac {\Gamma ({\frac {d+r}{p}})}{\Gamma ({\frac {d}{p}})}}.}
Properties [ edit ]
Denote GG(a,d,p) as the generalized gamma distribution of parameters a , d , p .
Then, given
c
{\displaystyle c}
and
α
{\displaystyle \alpha }
two positive real numbers, if
f
∼
G
G
(
a
,
d
,
p
)
{\displaystyle f\sim GG(a,d,p)}
, then
c
f
∼
G
G
(
c
a
,
d
,
p
)
{\displaystyle cf\sim GG(ca,d,p)}
and
f
α
∼
G
G
(
a
α
,
d
α
,
p
α
)
{\displaystyle f^{\alpha }\sim GG\left(a^{\alpha },{\frac {d}{\alpha }},{\frac {p}{\alpha }}\right)}
.
Kullback-Leibler divergence [ edit ]
If
f
1
{\displaystyle f_{1}}
and
f
2
{\displaystyle f_{2}}
are the probability density functions of two generalized gamma distributions, then their Kullback-Leibler divergence is given by
D
K
L
(
f
1
∥
f
2
)
=
∫
0
∞
f
1
(
x
;
a
1
,
d
1
,
p
1
)
ln
f
1
(
x
;
a
1
,
d
1
,
p
1
)
f
2
(
x
;
a
2
,
d
2
,
p
2
)
d
x
=
ln
p
1
a
2
d
2
Γ
(
d
2
/
p
2
)
p
2
a
1
d
1
Γ
(
d
1
/
p
1
)
+
[
ψ
(
d
1
/
p
1
)
p
1
+
ln
a
1
]
(
d
1
−
d
2
)
+
Γ
(
(
d
1
+
p
2
)
/
p
1
)
Γ
(
d
1
/
p
1
)
(
a
1
a
2
)
p
2
−
d
1
p
1
{\displaystyle {\begin{aligned}D_{KL}(f_{1}\parallel f_{2})&=\int _{0}^{\infty }f_{1}(x;a_{1},d_{1},p_{1})\,\ln {\frac {f_{1}(x;a_{1},d_{1},p_{1})}{f_{2}(x;a_{2},d_{2},p_{2})}}\,dx\\&=\ln {\frac {p_{1}\,a_{2}^{d_{2}}\,\Gamma \left(d_{2}/p_{2}\right)}{p_{2}\,a_{1}^{d_{1}}\,\Gamma \left(d_{1}/p_{1}\right)}}+\left[{\frac {\psi \left(d_{1}/p_{1}\right)}{p_{1}}}+\ln a_{1}\right](d_{1}-d_{2})+{\frac {\Gamma {\bigl (}(d_{1}+p_{2})/p_{1}{\bigr )}}{\Gamma \left(d_{1}/p_{1}\right)}}\left({\frac {a_{1}}{a_{2}}}\right)^{p_{2}}-{\frac {d_{1}}{p_{1}}}\end{aligned}}}
where
ψ
(
⋅
)
{\displaystyle \psi (\cdot )}
is the digamma function .[5]
Software implementation [ edit ]
In the R programming language, there are a few packages that include functions for fitting and generating generalized gamma distributions. The gamlss package in R allows for fitting and generating many different distribution families including generalized gamma (family=GG). Other options in R, implemented in the package flexsurv , include the function dgengamma , with parameterization:
μ
=
ln
a
+
ln
d
−
ln
p
p
{\displaystyle \mu =\ln a+{\frac {\ln d-\ln p}{p}}}
,
σ
=
1
p
d
{\displaystyle \sigma ={\frac {1}{\sqrt {pd}}}}
,
Q
=
p
d
{\displaystyle Q={\sqrt {\frac {p}{d}}}}
, and in the package ggamma with parametrisation:
a
=
a
{\displaystyle a=a}
,
b
=
p
{\displaystyle b=p}
,
k
=
d
/
p
{\displaystyle k=d/p}
.
In the python programming language, it is implemented in the SciPy package, with parametrisation:
c
=
p
{\displaystyle c=p}
,
a
=
d
/
p
{\displaystyle a=d/p}
, and scale of 1.
See also [ edit ]
References [ edit ]
^ Box-Steffensmeier, Janet M.; Jones, Bradford S. (2004) Event History Modeling: A Guide for Social Scientists . Cambridge University Press. ISBN 0-521-54673-7 (pp. 41-43)
^ Stacy, E.W. (1962). "A Generalization of the Gamma Distribution." Annals of Mathematical Statistics 33(3 ): 1187-1192. JSTOR 2237889
^ a b c Johnson, N.L.; Kotz, S; Balakrishnan, N. (1994) Continuous Univariate Distributions, Volume 1 , 2nd Edition. Wiley. ISBN 0-471-58495-9 (Section 17.8.7)
^ Gavin E. Crooks (2010), The Amoroso Distribution , Technical Note, Lawrence Berkeley National Laboratory.
^ C. Bauckhage (2014), Computing the Kullback-Leibler Divergence between two Generalized Gamma Distributions, arXiv :1401.6853 .
R e t r i e v e d f r o m " https://en.wikipedia.org/w/index.php?title=Generalized_gamma_distribution&oldid=1225701504 "
C a t e g o r i e s :
● C o n t i n u o u s d i s t r i b u t i o n s
● G a m m a a n d r e l a t e d f u n c t i o n s
H i d d e n c a t e g o r i e s :
● A r t i c l e s w i t h s h o r t d e s c r i p t i o n
● S h o r t d e s c r i p t i o n m a t c h e s W i k i d a t a
● T h i s p a g e w a s l a s t e d i t e d o n 2 6 M a y 2 0 2 4 , a t 0 4 : 4 4 ( U T C ) .
● T e x t i s a v a i l a b l e u n d e r t h e C r e a t i v e C o m m o n s A t t r i b u t i o n - S h a r e A l i k e L i c e n s e 4 . 0 ;
a d d i t i o n a l t e r m s m a y a p p l y . B y u s i n g t h i s s i t e , y o u a g r e e t o t h e T e r m s o f U s e a n d P r i v a c y P o l i c y . W i k i p e d i a ® i s a r e g i s t e r e d t r a d e m a r k o f t h e W i k i m e d i a F o u n d a t i o n , I n c . , a n o n - p r o f i t o r g a n i z a t i o n .
● P r i v a c y p o l i c y
● A b o u t W i k i p e d i a
● D i s c l a i m e r s
● C o n t a c t W i k i p e d i a
● C o d e o f C o n d u c t
● D e v e l o p e r s
● S t a t i s t i c s
● C o o k i e s t a t e m e n t
● M o b i l e v i e w