大佬教程收集整理的这篇文章主要介绍了如何计算多个 lm() 的 box cox \lambda,大佬教程大佬觉得挺不错的,现在分享给大家,也给大家做个参考。
library(tIDyversE)
获取和清理数据@H_675_7@
data(lme4::sleepstudy)
sleep <- as_tibble(sleepstudy) %>%
rename_all(tolower) %>%
rename(rt = reaction)
我想应用到多个lm()
的函数@H_675_7@
sleepBC <- function(tbl) {
Boxcox(tbl,plotit = falSE) %>%
as_tibble() %>%
filter(y == max(y)) %>%
pull(X) %>%
round_any(.,0.5)
}
round_any
复制自 plyr
:@H_675_7@
round_any = function(x,accuracy,f = round){f(x / accuracy) * accuracy}
测试一个主题:@H_675_7@
d308 <- sleep %>% filter(subject == "308")
m308 <- lm(rt ~ days,data = d308)
sleepBC(m308)
现在介绍多个模型。定义一个函数:@H_675_7@
sleepLM <- function(tbl) {
lm(rt ~ days,data = tbl)
}
我想要这样的东西:@H_675_7@
sleepMod <- sleep %>%
group_by(subject) %>%
nest() %>%
mutate(mods = map(data,sleepLM),tIDy = map(mods,broom::tIDy))
除了 sleepBC()
而不是 broom:tIDy()
@H_675_7@
暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!@H_675_7@
如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。@H_675_7@
小编邮箱:dio#foxmail.com (将#修改为@)@H_675_7@
以上是大佬教程为你收集整理的如何计算多个 lm() 的 box cox \lambda全部内容,希望文章能够帮你解决如何计算多个 lm() 的 box cox \lambda所遇到的程序开发问题。
如果觉得大佬教程网站内容还不错,欢迎将大佬教程推荐给程序员好友。
本图文内容来源于网友网络收集整理提供,作为学习参考使用,版权属于原作者。
如您有任何意见或建议可联系处理。小编QQ:384754419,请注明来意。