首页
学习
活动
专区
圈层
工具
发布
首页
学习
活动
专区
圈层
工具
MCP广场
社区首页 >问答首页 >等式: lmer模型“超越边界”

等式: lmer模型“超越边界”
EN

Stack Overflow用户
提问于 2021-06-14 07:24:01
回答 1查看 217关注 0票数 1

我正试着用等价物来画我的模型。

代码语言:javascript
复制
c1<-lmer(sqrt(co22)~t.air+t.soil+tdr+ph+TC+TN+wtd+Fe+0+(1|treatment/site), data=data)

然而,当我试图提取方程时,我总是得到一个错误:

代码语言:javascript
复制
extract_eq(c1,wrap = TRUE, ital_vars = FALSE)

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

但是,如果我将模型转换为一个没有随机效应的简单线性模型,效果会很好:

代码语言:javascript
复制
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,有什么建议吗?

EN

回答 1

Stack Overflow用户

回答已采纳

发布于 2021-06-14 20:36:29

我是那个包的开发者。它应该适用于lme4::lmer()方程。这里的问题是放弃拦截,但在更高的层次上,它的变化是随机的。有关详细信息,请参阅问题。

如果您有一些建议,您可以期望方程式呈现,我愿意做一个修正。但就目前而言,equatiomatic::extract_eq()假设无论有什么随机效应,也有相应的固定效应。

票数 1
EN
页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接:

https://stackoverflow.com/questions/67966270

复制
相关文章

相似问题

领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档