将独立的截距/斜率混合模型公式转换为数学公式

问题描述

在R会话中，我将具有以下公式：

members ~ face + class + gender + number + (day || country)

我该如何用数学公式写出来？顺便问一下，它是一个随机斜率随机拦截模型吗？

解决方法

tl; dr

这只是“困难的部分”；您未纳入模型随机效应部分的协变量将作为第一个方程（beta_2*face + beta_3*class + ...）的一部分被添加。 （顺便说一句，这很不寻常[即除非您有某种特定的原因，否则可能是错误的，如果days是一个连续变量，则将其包含在随机效应模型中而不是固定变量中，效果模型...）

使用equatiomatic的开发版本（请注意，这里的大多数复杂性是将extract_eq()的结果转换为PDF，然后将PNG张贴到此处-需要一些命令行工具（ LaTeX，pdfcrop，ImageMagick）。也许可以使用tinytex包以一种更加独立的方式来实现此目的：

while (!require("equatiomatic")) {
    remotes::install_github("datalorax/equatiomatic")
}
library(lme4)
m1 <- lmer(Reaction ~ Days + (Days||Subject),sleepstudy)
unlink("tmp.tex")

writeLines(c("\\documentclass{article}","\\usepackage{amsmath}","\\begin{document}","\\thispagestyle{empty}",format(extract_eq(m1)),"\\end{document}"),con="tmp.tex")

system("pdflatex tmp.tex")
system("pdfcrop tmp.pdf")
system("convert tmp-crop.pdf tmp.png")