VGAM (version 0.9-4)

Coef.qrrvglm-class: Class ``Coef.qrrvglm''

Description

The most pertinent matrices and other quantities pertaining to a QRR-VGLM (CQO model).

Objects from the Class

Objects can be created by calls of the form Coef(object, ...) where object is an object of class "qrrvglm" (created by cqo).

In this document, $R$ is the rank, $M$ is the number of linear predictors and $n$ is the number of observations.

References

Yee, T. W. (2004) A new technique for maximum-likelihood canonical Gaussian ordination. Ecological Monographs, 74, 685--701.

See Also

Coef.qrrvglm, cqo, print.Coef.qrrvglm.

Examples

x2 <- rnorm(n <- 100)
x3 <- rnorm(n)
x4 <- rnorm(n)
latvar1 <- 0 + x3 - 2*x4
lambda1 <- exp(3 - 0.5 * ( latvar1-0)^2)
lambda2 <- exp(2 - 0.5 * ( latvar1-1)^2)
lambda3 <- exp(2 - 0.5 * ((latvar1+4)/2)^2)
y1 <- rpois(n, lambda1)
y2 <- rpois(n, lambda2)
y3 <- rpois(n, lambda3)
yy <- cbind(y1, y2, y3)
# vvv p1 <- cqo(yy ~ x2 + x3 + x4, fam = poissonff, trace = FALSE)
lvplot(p1, y = TRUE, lcol = 1:3, pch = 1:3, pcol = 1:3)
# vvv print(Coef(p1), digits = 3)