VGAM (version 0.7-1)

lgammaff: Log-gamma Distribution Family Function

Description

Estimation of the parameter of the standard and nonstandard log-gamma distribution.

Usage

lgammaff(link = "loge", earg=list(), init.k = NULL)
lgamma3ff(llocation="identity", lscale="loge", lshape="loge",
          elocation=list(), escale=list(), eshape=list(),
          ilocation=NULL, iscale=NULL, ishape=1, zero=NULL)

Value

• An object of class "vglmff" (see vglmff-class). The object is used by modelling functions such as vglm, and vgam.

Details

The probability density function of the standard log-gamma distribution is given by $$f(y)=\exp[ky - \exp(y)] / \Gamma(k),$$ for parameter $k>0$ and all real $y$. The mean of $Y$ is digamma(k) (returned as the fitted values) and its variance is trigamma(k).

For the non-standard log-gamma distribution, one replaces $y$ by $(y-a)/b$, where $a$ is the location parameter and $b$ is the positive scale parameter. Then the density function is $$f(y)=\exp[k(y-a)/b - \exp((y-a)/b)] / (b \Gamma(k)).$$ The mean and variance of $Y$ are a + b*digamma(k) (returned as the fitted values) and b^2 * trigamma(k), respectively.

References

Kotz, S. and Nadarajah, S. (2000) Extreme Value Distributions: Theory and Applications, pages 48--49, London: Imperial College Press.

Johnson, N. L. and Kotz, S. and Balakrishnan, N. (1995) Continuous Univariate Distributions, 2nd edition, Volume 2, p.89, New York: Wiley.

See Also

rlgamma, ggamma, prentice74, gamma1, lgamma.

Examples

ldat = data.frame(y = rlgamma(nn <- 100, k=exp(1)))
fit = vglm(y ~ 1, lgammaff, ldat, trace=TRUE, crit="c")
summary(fit)
coef(fit, matrix=TRUE)
Coef(fit)

ldat = data.frame(x = runif(nn <- 5000))     # Another example
ldat = transform(ldat, loc = -1 + 2*x,
                       Scale = exp(1+x))
ldat = transform(ldat, y = rlgamma(nn, loc=loc, scale=Scale, k=exp(0)))
fit = vglm(y ~ x, lgamma3ff(zero=3), ldat, trace=TRUE, crit="c")
coef(fit, matrix=TRUE)