Displaying 1 result from an estimated 1 matches for "fixep".
Did you mean:
fixed
2023 Dec 02
1
Try reproduce glmm by hand
...x=30), runif(40,
min=40, max=60))
x <- (x-min(x))/(max(x)-min(x))
# In g0, I have the results of the glmm
library(lme4)
g0 <- glmer(formula = df ~ x + (1 | ID), family =
binomial(link="logit"), nAGQ=1)
-logLik(g0) # 'log Lik.' 268.0188 (df=3)
# I get the fitted parameters
fixep <- fixef(g0)
par <- getME(g0, c("theta","beta"))
# _______________________________________________________________________
# Question 1: how theta is converted into the specific effect on
(intercept) for the random effect ?
# Then how a theta parameter is converted into i...