首页 > 解决方案 > 将独立截距/斜率混合模型公式转换为数学公式

问题描述

在 R 会话中,我将有以下公式:

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

我将如何在数学公式中写这个?顺便问一下,它是一个随机斜率随机截距模型吗?

标签: rmixed-modelsmultilevel-analysis

解决方案


tl;博士

方程公式

这只是“困难的部分”;未参与模型随机效应部分的协变量将作为第一个方程 ( 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")

推荐阅读