[R-sig-ME] Converting lme4 model syntax into corresponding regression equation

Simon Harmel @|m@h@rme| @end|ng |rom gm@||@com
Sun Oct 4 02:43:41 CEST 2020


Hello all,

I was wondering if there is any way (e.g., from an existing R package or
function) to obtain the combined regression equation from an R formula used
to fit an lme4::lmer() model?

For example, if I used `lmer(math ~ ses*sector + (ses | sch.id), data =
data)`, then the expected output in R expression would be:

ex <- expression(math == gamma[00] + gamma[01]*sector + gamma[10]*ses +
gamma[11]*(sector *ses)+U[0]+U1*ses+e)

plot(1, main = ex, cex.main = .7)

	[[alternative HTML version deleted]]



More information about the R-sig-mixed-models mailing list