我正试着用等价物来画我的模型。
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}
$$对于如何解决“下标出界”-error,有什么建议吗?
发布于 2021-06-14 20:36:29
我是那个包的开发者。它应该适用于lme4::lmer()方程。这里的问题是放弃拦截,但在更高的层次上,它的变化是随机的。有关详细信息,请参阅这问题。
如果您有一些建议,您可以期望方程式呈现,我愿意做一个修正。但就目前而言,equatiomatic::extract_eq()假设无论有什么随机效应,也有相应的固定效应。
https://stackoverflow.com/questions/67966270
复制相似问题