从 R 中的线性模型列表映射 emmeans

Question

我有一个超过 100 个线性模型的列表，我想为每个模型获取估计均值和标准误差。

我们以mtcars为例。

library(tidyverse); library(magrittr); library(emmeans)

mtcars %<>% 
  mutate(
    cyl = as.factor(cyl)
  )

df <- mtcars %>% select(cyl, hp, mpg)

我可以使用 emmeans:

轻松获得每个模型的估计均值和标准误差

mod <- lm(hp ~ cyl, data = df)
emmeans(mod, "cyl")

但是如果我有一个模型列表呢？

list_lm <- df %>% 
   select(-c(cyl)) %>%
   map(function(dv) lm(dv ~ df$cyl, data = .))

我不能使用：

emmeans(list_lm$hp, "cyl")
Error in ref_grid(object, ...) : Perhaps a 'data' or 'params' argument is needed

理想情况下，我想要一些东西可以为我提供所有模型的这些统计数据。模型系数类似于 broom::tidy，但对于 emmeans:

list_lm %>% 
   map(broom::tidy)

Answer 1

你的直觉是对的。该解决方案需要将中间结果保存在列表列中，然后将它们解包，但考虑到 emmeans 输出的结构，broom::tidy() 不是必需的。只需将 emmeans 输出转换为 data.frames 和 unnest() 的列表列。

library(dplyr)
library(purrr)
library(tidyr)
library(emmeans)

ds_mtcars <- 
  mtcars %>% 
  mutate(cyl = as.factor(cyl)) 

ds_nest <- 
  ds_mtcars %>% 
  group_by(am) %>% 
  nest() 

foo_model <- function(data){
  lm(hp ~ cyl, data = data)
}

ds_nest <- ds_nest %>% mutate(model = map(.x = data, .f = foo_model))

ds_temp <- 
  ds_nest %>% 
  mutate(
    emmeans = pmap(
      .l = list(
        object = model, 
        specs = "cyl"
      ),
      .f = emmeans
    )
  ) 

ds_temp %>% 
  mutate(emm2 = map(emmeans, data.frame)) %>% 
  unnest(emm2)

另外，purrr 函数 map() 和 pmap() 可能会让人费解，但我尽我所能在我的 [=16= 上引导我自己和新用户使用这些函数].

从 R 中的线性模型列表映射 emmeans

map emmeans from a list of linear models in R

dictionary

r

linear-regression

tidyverse

emmeans