# Fit the model
fit_RML <- lme(fixed = {user_fixed},
random = {user_random}, {corr_struct} {var_struct}
data = dat,
control = lmeControl(msMaxIter = 50, apVar = FALSE, returnObject = TRUE))
summary(fit_RML)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.