Template:Generalized gamma probability density function: Difference between revisions

From ReliaWiki
Jump to navigation Jump to search
No edit summary
Line 22: Line 22:


While this makes the distribution converge much more easily in computations, it does not facilitate manual manipulation of the equation. By allowing <math>\lambda </math> to become negative, the ''pdf'' of the reparameterized distribution is given by:
While this makes the distribution converge much more easily in computations, it does not facilitate manual manipulation of the equation. By allowing <math>\lambda </math> to become negative, the ''pdf'' of the reparameterized distribution is given by:
::<math>f(t)=\left\{ \begin{matrix}
  \tfrac{|\lambda |}{\sigma \cdot t}\cdot \tfrac{1}{\Gamma \left( \tfrac{1}{{{\lambda }^{2}}} \right)}\cdot {{e}^{\left[ \tfrac{\lambda \cdot \tfrac{\text{ln}(t)-\mu }{\sigma }+\text{ln}\left( \tfrac{1}{{{\lambda }^{2}}} \right)-{{e}^{\lambda \cdot \tfrac{\text{ln}(t)-\mu }{\sigma }}}}{{{\lambda }^{2}}} \right]}}\text{ if }\lambda \ne 0  \\
  \tfrac{1}{t\cdot \sigma \sqrt{2\pi }}{{e}^{-\tfrac{1}{2}{{\left( \tfrac{\text{ln}(t)-\mu }{\sigma } \right)}^{2}}}}\text{                            if }\lambda =0  \\
\end{matrix} \right.</math>

Revision as of 20:58, 14 February 2012

Generalized Gamma Probability Density Function

The generalized gamma function is a three-parameter distribution. One version of the generalized gamma distribution uses the parameters k, [math]\displaystyle{ \beta }[/math], and [math]\displaystyle{ \theta }[/math]. The pdf for this form of the generalized gamma distribution is given by:

[math]\displaystyle{ f(t)=\frac{\beta }{\Gamma (k)\cdot \theta }{{\left( \frac{t}{\theta } \right)}^{k\beta -1}}{{e}^{-{{\left( \frac{t}{\theta } \right)}^{\beta }}}} }[/math]

where [math]\displaystyle{ \theta \gt 0 }[/math] is a scale parameter, [math]\displaystyle{ \beta \gt 0 }[/math] and [math]\displaystyle{ k\gt 0 }[/math] are shape parameters and [math]\displaystyle{ \Gamma (x) }[/math] is the gamma function of x, which is defined by:

[math]\displaystyle{ \Gamma (x)=\int_{0}^{\infty }{{{s}^{x-1}}}\cdot {{e}^{-s}}ds }[/math]

With this version of the distribution, however, convergence problems arise that severely limit its usefulness. Even with data sets containing 200 or more data points, the MLE methods may fail to converge. Further adding to the confusion is the fact that distributions with widely different values of k, [math]\displaystyle{ \beta }[/math], and [math]\displaystyle{ \theta }[/math] may appear almost identical [21]. In order to overcome these difficulties, Weibull++ uses a reparameterization with parameters k, [math]\displaystyle{ \beta }[/math], and [math]\displaystyle{ \theta }[/math] [21] where:

[math]\displaystyle{ \begin{align} & \mu =\ln (\theta )+\frac{1}{\beta }\cdot \ln \left( \frac{1}{{{\lambda }^{2}}} \right) \\ & \\ & \sigma =\frac{1}{\beta \sqrt{k}} \\ & \\ & \lambda =\frac{1}{\sqrt{k}} \\ \end{align} }[/math]

where [math]\displaystyle{ -\infty \lt \mu \lt \infty ,\text{ }\sigma \gt 0\text{, }0\lt \lambda . }[/math]

While this makes the distribution converge much more easily in computations, it does not facilitate manual manipulation of the equation. By allowing [math]\displaystyle{ \lambda }[/math] to become negative, the pdf of the reparameterized distribution is given by:


[math]\displaystyle{ f(t)=\left\{ \begin{matrix} \tfrac{|\lambda |}{\sigma \cdot t}\cdot \tfrac{1}{\Gamma \left( \tfrac{1}{{{\lambda }^{2}}} \right)}\cdot {{e}^{\left[ \tfrac{\lambda \cdot \tfrac{\text{ln}(t)-\mu }{\sigma }+\text{ln}\left( \tfrac{1}{{{\lambda }^{2}}} \right)-{{e}^{\lambda \cdot \tfrac{\text{ln}(t)-\mu }{\sigma }}}}{{{\lambda }^{2}}} \right]}}\text{ if }\lambda \ne 0 \\ \tfrac{1}{t\cdot \sigma \sqrt{2\pi }}{{e}^{-\tfrac{1}{2}{{\left( \tfrac{\text{ln}(t)-\mu }{\sigma } \right)}^{2}}}}\text{ if }\lambda =0 \\ \end{matrix} \right. }[/math]