[R-sig-ME] family() in probit-glmer
Alessandro Moscatelli
moskante at gmail.com
Fri Sep 5 12:22:01 CEST 2014
Dear all,
I have a problem using the function family(). I am using lme4 version 1.1-7
running on R 3.1.0
I fit a glmer model with a probit link function:
formula.mod = cbind(Longer, Total - Longer) ~ X * condition + (1 + X|
Subject)
mod1 <- glmer(formula = formula.mod, family = binomial(link = "probit"),
data = datafr)
> mod1
Generalized linear mixed model fit by maximum likelihood (Laplace
Approximation) ['glmerMod']
Family: binomial ( probit )
Formula: cbind(Longer, Total - Longer) ~ X * condition + (1 + X | Subject)
Data: datafr
AIC BIC logLik deviance df.resid
330.1026 348.8775 -158.0513 316.1026 101
Random effects:
Groups Name Std.Dev. Corr
Subject (Intercept) 1.452885
X 0.001128 -0.99
Number of obs: 108, groups: Subject, 6
Fixed Effects:
(Intercept) X conditionB X:conditionB
-7.1760859 0.0090841 -0.2190119 0.0001522
However, using the family() function:
> family(mod1)
Family: binomial
Link function: logit
Does anyone has the same problem? Is it a bug? For Now I solved the problem
using:
> summary(mod1)$link
[1] "probit"
Best
Alessandro
[[alternative HTML version deleted]]
More information about the R-sig-mixed-models
mailing list