r-赤道:lmer模型;subscript越界"

Question

我正在尝试使用equiomatic包来绘制我的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()方程。这里的问题是放弃拦截，但让它在更高的级别上随机变化。有关更多详细信息，请参阅本期。

如果你对如何渲染方程有建议，我愿意解决。但就目前而言，equatiomatic::extract_eq()假设无论你有什么随机效应，也有相应的固定效应。