错误-lmerMod类的对象没有tidy方法
我正在使用在线实践教程中的数据集,代码可以在第4页的底部找到() 在本教程中,他们使用下面列出的代码使函数正常工作,但在我的R会话中,我收到一个错误,上面说: 我尝试过使用“欧洲防风草”软件包并重新启动我的R会话,我也尝试过像在其他类似问题的答案中所建议的那样使用扫帚 haggis practice csv文件可从以下位置下载:错误-lmerMod类的对象没有tidy方法,r,tidy,broom,R,Tidy,Broom,我正在使用在线实践教程中的数据集,代码可以在第4页的底部找到() 在本教程中,他们使用下面列出的代码使函数正常工作,但在我的R会话中,我收到一个错误,上面说: 我尝试过使用“欧洲防风草”软件包并重新启动我的R会话,我也尝试过像在其他类似问题的答案中所建议的那样使用扫帚 haggis practice csv文件可从以下位置下载: 库(asreml) 图书馆(nadiv) 图书馆(tidytext) 图书馆(tidyverse) 图书馆(扫帚) 需要(扫帚) 图书馆(lme4) 库(数据表) 图书
库(asreml)
图书馆(nadiv)
图书馆(tidytext)
图书馆(tidyverse)
图书馆(扫帚)
需要(扫帚)
图书馆(lme4)
库(数据表)
图书馆(防风草)
HData%
变异(重复性=ID/(ID+残余))
提供答案(来自评论)
多级/混合型模型的
tidy
方法(例如来自lme4
,brms
,MCMCglmm
,…)被移动到broom.mixed
。您可以安装/加载broom.mixed
软件包,或者使用broomExtra
软件包,这是一个“元软件包”,用于在broom
和broom.mixed
中查找方法。这可能只适用于较旧版本的broom
。似乎此功能在某个时候被移动到了一个单独的包“broom.mixed”,请参见此处:。这为lmer-fits提供了一个tidy()
函数。真管用!非常感谢你!
No tidy method for objects of class lmerMod.
library(asreml)
library(nadiv)
library(tidytext)
library(tidyverse)
library(broom)
require(broom)
library(lme4)
library(data.table)
library(parsnip)
HData<- read_csv("haggis practice.csv")
lmer_b <- lmer(boldness ~ scale(assay_rep, scale=FALSE) +
scale(body_size) +
(1|ID),
data = HData)
plot(lmer_b)
qqnorm(residuals(lmer_b))
hist(residuals(lmer_b))
summary(lmer_b)
rep_bold <- tidy(lmer_b, effects = "ran_pars", scales = "vcov") %>%
select(group, estimate) %>%
spread(group, estimate) %>%
mutate(repeatability = ID/(ID + Residual))