方程式：lmer 模型 "subsript out of bounds"

Question

我正在尝试使用 equatiomatic 包来绘制我的 lmer 模型。

c1<-lmer(sqrt(co22)~t.air+t.soil+tdr+ph+TC+TN+wtd+Fe+0+(1|treatment/site), data=data)

但是，在尝试提取方程式时，我不断收到错误消息：

extract_eq(c1,wrap = TRUE, ital_vars = FALSE)

Error in `[<-`(`*tmp*`, v[1], v[2], value = greek_vcov) : 
  subscript out of bounds

但是，如果我将模型转换为没有随机效应的简单线性模型，效果很好：

c1<-lm(sqrt(co22)~t.air+t.soil+tdr+ph+TC+TN+wtd+Fe+0, data=data)

extract_eq(c1,wrap = TRUE, ital_vars = FALSE)

$$
\begin{aligned}
\operatorname{sqrt(co22)} &= \beta_{0}(\operatorname{t.air}) + \beta_{1}(\operatorname{t.soil}) + \beta_{2}(\operatorname{tdr}) + \beta_{3}(\operatorname{ph})\ + \
&\quad \beta_{4}(\operatorname{TC}) + \beta_{5}(\operatorname{TN}) + \beta_{6}(\operatorname{wtd}) + \beta_{7}(\operatorname{Fe})\ + \
&\quad \epsilon
\end{aligned}
$$

关于如何解决“下标越界”错误的任何建议？

Answer 1

我是那个包的开发者。它应该适用于 lme4::lmer() 方程。这里的问题是丢弃截距但让它在更高级别随机变化。有关详细信息，请参阅 this 问题。

如果您对方程式的呈现方式有任何建议，我愿意制定解决方案。但是现在，equatiomatic::extract_eq() 假定您拥有的任何随机效应也具有相应的固定效应。

方程式：lmer 模型 "subsript out of bounds"

Equatiomatic: lmer model "subsript out of bounds"

r

equation

lme4