fam.link.mer | R Documentation |
This function extracts the distribution family and link function of a
generalized linear mixed model fit with glmer
or lmer
.
fam.link.mer(mod)
mod |
an object of |
This utility function extracts the information from an mer
or
merMod
object resulting from glmer
or lmer
. The
function is called by modavg
, modavgEffect
,
modavgPred
, and predictSE
.
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 |
Marc J. Mazerolle
modavg
, modavgPred
,
predictSE
, glmer
, lmer
##modified example from ?glmer
## Not run:
if(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)
}
## End(Not run)
##example with linear mixed model with Orthodont data from
##Pinheiro and Bates (2000)
## Not run:
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)
## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.