我正在使用在线练习教程中的数据集,代码可以在第 4 页的底部找到(https://tomhouslay.files.wordpress.com/2017/02/indivvar_mv_tutorial_asreml.pdf)
在本教程中,他们使用下面列出的代码使该函数正常工作,但在我的 R 会话中,我收到一条错误消息:
Run Code Online (Sandbox Code Playgroud)No tidy method for objects of class lmerMod.
我尝试过使用“parsnip”包并重新启动我的 R 会话,并且我尝试过按照类似问题的其他答案中的建议要求扫帚。
哈吉斯练习 csv 文件可以从这里下载:https://figshare.com/articles/Haggis_data_behavioural_syndromes/4702540
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))
Run Code Online (Sandbox Code Playgroud)
Ben*_*ker 17
提供答案(来自评论)。
tidy多级/混合类型模型的方法(例如,从、lme4、brms、MCMCglmm...)移至broom.mixed。您可以安装/加载该broom.mixed包,也可以使用该包,它是一个“元包”,可以在...broomExtra中查找方法broombroom.mixed