R 将独立截距/斜率混合模型公式转换为数学公式
在R会话中,我将使用以下公式:R 将独立截距/斜率混合模型公式转换为数学公式,r,mixed-models,multilevel-analysis,R,Mixed Models,Multilevel Analysis,在R会话中,我将使用以下公式: members ~ face + class + gender + number + (day || country) 我该如何在数学公式中写出这个?顺便问一下,这是一个随机斜率随机截距模型吗;博士 这只是“难的部分”;模型随机效应部分中未涉及的协变量将作为第一个方程的一部分添加(beta_2*face+beta_3*class+…)。(顺便说一句,如果days是一个连续变量,将其置于随机效应模型中而非固定效应模型中是不寻常的[即,除非您有特定的理由这样做,否
members ~ face + class + gender + number + (day || country)
我该如何在数学公式中写出这个?顺便问一下,这是一个随机斜率随机截距模型吗;博士
这只是“难的部分”;模型随机效应部分中未涉及的协变量将作为第一个方程的一部分添加(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")
while(!require(“equatiomatic”)){
remotes::install_github(“datalorax/equatiomatic”)
}
图书馆(lme4)
m1查看equatiomatic软件包,看看这是否有帮助:它通常被称为“具有不相关斜率和截距的随机斜率模型”