lindley: 1-parameter Gamma Distribution

Description Usage Arguments Details Value Note Author(s) References See Also Examples

View source: R/family.aunivariate.R

Description

Estimates the (1-parameter) Lindley distribution by maximum likelihood estimation.

Usage

1
lindley(link = "loge", itheta = NULL, zero = NULL)

Arguments

link

Link function applied to the (positive) parameter. See Links for more choices.

itheta, zero

See CommonVGAMffArguments for information.

Details

The density function is given by

f(y; theta) = theta^2 * (1 + y) * exp(-theta * y) / (1 + theta)

for theta > 0 and y > 0. The mean of Y (returned as the fitted values) is mu = (theta + 2) / (theta * (theta + 1)). The variance is (theta^2 + 4 * theta + 2) / (theta * (theta + 1))^2.

Value

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

Note

This VGAM family function can handle multiple responses (inputted as a matrix). Fisher scoring is implemented.

Author(s)

T. W. Yee

References

Lindley, D. V. (1958) Fiducial distributions and Bayes' theorem. Journal of the Royal Statistical Society, Series B, Methodological, 20, 102–107.

Ghitany, M. E. and Atieh, B. and Nadarajah, S. (2008) Lindley distribution and its application. Math. Comput. Simul., 78, 493–506.

See Also

dlind, gammaR, simulate.vlm.

Examples

1
2
3
4
5
ldata <- data.frame(y = rlind(n = 1000, theta = exp(3)))
fit <- vglm(y ~ 1, lindley, data = ldata, trace = TRUE, crit = "coef")
coef(fit, matrix = TRUE)
Coef(fit)
summary(fit)

VGAM documentation built on May 19, 2017, 10:21 a.m.

Search within the VGAM package
Search all R packages, documentation and source code

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

Please suggest features or report bugs in the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.