Returns the rank of reduced-rank regression-type models in the VGAM package.
Rank(object, ...)
Rank.rrvglm(object, ...)
Rank.qrrvglm(object, ...)
Rank.rrvgam(object, ...)Returns an integer value, provided the rank of the model makes sense.
Some VGAM object, for example, having
class rrvglm-class.
The class vglm-class is not included since this
is not based on reduced-rank regression.
Other possible arguments fed into the function later (used for added flexibility for the future).
T. W. Yee.
Regression models based on reduced-rank regression have a quantity
called the rank, which is 1 or 2 or 3 etc.
The smaller the value the more dimension reduction, so that there
are fewer parameters.
This function was not called rank() to avoid conflict
with rank.
RR-VGLMs are described in rrvglm-class;
QRR-VGLMs are described in qrrvglm-class.
pneumo <- transform(pneumo, let = log(exposure.time),
x3 = runif(nrow(pneumo)))
(fit1 <- rrvglm(cbind(normal, mild, severe) ~ let + x3,
acat, data = pneumo))
coef(fit1, matrix = TRUE)
constraints(fit1)
Rank(fit1)
Run the code above in your browser using DataLab