Results in lme4 and nlme differ from glmfit_multilevel

[liznecka]

Hi there,

I’m trying to verify that the results I get from glmfit_multilevel and that I get from lme4 and nlme in R are the same, but I am finding that they are not. My model has one continuous predictor and one categorical predictor, both on level 1, and I can include a continuous interaction as well. I’m new to using glmfit_multilevel, but have been using R for years, so am wondering how glmfit_multilevel computes betas and SE that is different from the R packages. I’ve attached a reproducible example. Any clarification on this would help.

Thanks in advance,

Liz

Glmfit_multilevelvsR.zip