extractfit_lmer <-
function(x)
{
mu <- fixef(x)
sfitx <- summary(x)
ni <- c(sfitx$ngrps, "total"=nrow(x@frame))
vari <- c(unlist(VarCorr(x)), sfitx$sigma^2)
return(list(mu=mu, ni=ni, vari=vari))
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.