AICcmodavg (version 1.35)

fam.link.mer: Extract Distribution Family and Link Function

Description

This function extracts the distribution family and link function of a generalized linear mixed model fit with 'glmer' or 'lmer'.

Usage

fam.link.mer(mod)

Arguments

mod

an object of 'mer' or 'merMod' class resulting from the fit of 'glmer' or 'lmer'.

Value

• 'fam.link.mer' returns a list with the following components:
• family
  the family of the distribution of the model.
• link
  the link function of the model.
• supp.link
  a character value indicating whether the link function used is supported by 'predictSE.mer', 'predictSE.merMod', 'modavgpred.mer', and 'modavgpred.merMod'.

Details

This utility function extracts the information from an 'mer' or 'merMod' object resulting from 'glmer' or 'lmer'. The function is called by 'modavg.mer', 'modavg.merMod', 'modavg.effect.mer', 'modavg.effect.merMod', 'modavgpred.mer', 'modavgpred.merMod', 'predictSE.mer', and 'predictSE.merMod'.

See Also

modavg, modavgpred, predictSE.mer, glmer, lmer

Examples

##modified example from ?glmer
require(lme4)
##create proportion of incidence
cbpp$prop <- cbpp$incidence/cbpp$size
gm1 <- glmer(prop ~ period + (1 | herd), family = binomial,
             weights = size, data = cbpp)
fam.link.mer(gm1)
gm2 <- glmer(prop ~ period + (1 | herd),
             family = binomial(link = "cloglog"), weights = size,
             data = cbpp)
fam.link.mer(gm2)

##example with linear mixed model with Orthodont data from
##Pinheiro and Bates (2000)
data(Orthodont, package = "nlme")
m1 <- lmer(distance ~ Sex + (1 | Subject), data = Orthodont,
           REML = FALSE)
fam.link.mer(m1)
m2 <- glmer(distance ~ Sex + (1 | Subject),
            family = gaussian(link = "log"), data = Orthodont,
            REML = FALSE)
fam.link.mer(m2)

detach(package:lme4)