R语言如何更方便地检索lm()系数

vngu2lb8 于 2023-02-14 发布在其他

关注(0)|答案(1)|浏览(143)

在dplyr中，使用group_map可以生成lm（）模型，并将系数合并到一个 Dataframe 中。我想将组信息和Adjusted_R_squared添加到最终的data.fame中。
目前我键入值作为向量并添加到data.frame。有什么方便的方法吗？

library(dplyr)
library(data.table)

models <- iris %>% 
  group_by(Species) %>% 
  group_map(~ broom::tidy(lm(Petal.Length ~ Sepal.Length,data=.x))) %>% 
  rbindlist()
    
models$Species <- c('setosa','setosa','versicolor','versicolor','virginica','virginica')
models$Adjusted_R_squared <- c('0.05204','0.05204','0.5596','0.5596','0.7416','0.7416')

来源：https://stackoverflow.com/questions/75432774/how-to-retrieve-the-lm-coefficient-more-easily

1条答案

按热度按时间

l7wslrjt1#

如果你在你的group_map()中打开一个函数，你可以对你的输出进行更多的定制，比如：

models <- iris %>% 
  group_by(Species) %>% 
  group_map(
    function(spe, spe_name) {
      results <- lm(Petal.Length ~ Sepal.Length, data = spe)
      broom::tidy(results) %>%
        mutate(Species = spe_name$Species, Adjusted_R_squared = summary(results)$adj.r.squared)
    }
  ) %>% 
  rbindlist()

赞(0）回复(0）举报 2023-02-14

我来回答

R语言如何更方便地检索lm()系数

1条答案

相关问题

热门标签

最新问答

R语言 如何更方便地检索lm()系数

1条答案

相关问题

热门标签

最新问答

R语言如何更方便地检索lm()系数