vcdExtra (version 0.5-2)

glmlist: Create a Model List Object

Description

glmlist creates a glmlist object containing a list of fitted glm objects with their names. loglmlist does the same for loglm objects. The intention is to provide object classes to facilitate model comparison, extraction, summary and plotting of model components, etc., perhaps using lapply or similar.

Usage

glmlist(...)
loglmlist(...)

Arguments

...

One or more model objects, as appropriate to the function, optionally assigned names as in list.

Value

• An object of class glmlist loglmlist, just like a list, except that each is given a name attribute.

Details

The arguments to glmlist or loglmlist are of the form value or name=value. Any objects which do not inherit the appropriate class glm or loglm are excluded, with a warning.

See Also

The function llist in package Hmisc is similar, but perplexingly more general. The function anova.glm also handles glmlist objects

Examples

data(Mental)
indep <- glm(Freq ~ mental+ses,
                family = poisson, data = Mental)
Cscore <- as.numeric(Mental$ses)
Rscore <- as.numeric(Mental$mental)

coleff <- glm(Freq ~ mental + ses + Rscore:ses,
                family = poisson, data = Mental)
roweff <- glm(Freq ~ mental + ses + mental:Cscore,
                family = poisson, data = Mental)
linlin <- glm(Freq ~ mental + ses + Rscore:Cscore,
                family = poisson, data = Mental)
                
# use object names
mods <- glmlist(indep, coleff, roweff, linlin)
names(mods)

# assign new names
mods <- glmlist(Indep=indep, Col=coleff, Row=roweff, LinxLin=linlin)
names(mods)

summarise(mods)

#extract model components
unlist(lapply(mods, deviance))

res <- lapply(mods, residuals)
boxplot(as.data.frame(res), main="Residuals from various models")