将结果从 map_dfr 传递到后续的 map_dfr 以将自定义函数应用于数据组

Pipe result from map_dfr to a subsequent map_dfr to apply custom fuction to groups of data

我想将 certain function(即下面的 AddLags)应用于数据帧的组。为此,我尝试使用两个连续的 map_dfr(一个连接到另一个),以便应用相应的过滤器。对于最后一步,我正在应用自定义函数(前面提到)- 使用 map_dfr(在新对象中捕获新计算的输出数据)。

我目前的代码如下:

# dummy dataset
df <- data.frame(
  date = seq(today(),length.out=12,by='month'),
  dim1 = c('a','a','a','b','b','b','c','c','c','d','d','d'),
  dim2 = c(1,1,1,1,1,1,2,2,2,2,2,2),
  value = 1:12
  )

# function to apply
AddLags <- function(df,lags_vector,target_col,date_col){
  temp_lags <- map_dfc(lags_vector, 
                       ~ df %>% 
                         arrange({{date_col}}) %>% 
                         transmute(
                           across(contains(target_col), lag, .x, .names = '{col}_lag_{ifelse(.x<10,paste0("0",.x),.x)}')
                         )
  )
  return(temp_lags)
}


# prepare for map_dfr approach
lags_features <- c(1,2)
dims1 <- df %>% pull(dim1) %>% unique %>% sort
dims2 <- df %>% pull(dim2) %>% unique %>% sort

# what I am struggling with
map_dfr(dims1, 
        ~ df %>%
          filter(dim1==.x) %>%
          map_dfr(dims2,
                 ~ . %>% 
                   filter(dim2==.x) %>% 
                   AddLags(lags_features,variable,date)
          )
)

# how the loop version would look like
gather_results <- data.frame()
for(d1 in dims1){
  for(d2 in dims2){
    tempdata <- df %>% filter(dim1==d1,dim2==dim2) %>% arrange(date)
    temp <- AddLags(tempdata)
    gather_results %<>% bind_rows(temp)   
  }
}

本质上,我正在遍历不同的组(通过筛选)并分别应用自定义函数,同时尝试使用 map_dfr 合并新计算的结果。

我想知道如何实现上述目标(假设这是可行的)以及我缺少什么,因为目前我得到的只是一个空数据框。

奖金问题: 在我写这篇文章时,我意识到必须有更好的方法来代替循环——例如使用 group_by——但考虑到问题的性质和函数输出新数据的事实,我不确定这会是什么样子(假设一开始是可行的)。因此,任何类型的 suggestion/alternative/best 实践都将不胜感激。

免责声明: 当涉及到 purrr 功能时,我是一个大菜鸟,而且也不是一个有经验的 dplyr 用户,所以请原谅我的无知。

这是预期的输出吗?

library(tidyverse)
library(lubridate)

group_split(df, dim1, dim2) %>%
  map_dfr(~ .x %>% AddLags(1:2, "value", date))
#> # A tibble: 12 × 2
#>    value_lag_01 value_lag_02
#>           <int>        <int>
#>  1           NA           NA
#>  2            1           NA
#>  3            2            1
#>  4           NA           NA
#>  5            4           NA
#>  6            5            4
#>  7           NA           NA
#>  8            7           NA
#>  9            8            7
#> 10           NA           NA
#> 11           10           NA
#> 12           11           10

数据:

# dummy dataset
df <- data.frame(
  date = seq(today(), length.out = 12, by = "month"),
  dim1 = c("a", "a", "a", "b", "b", "b", "c", "c", "c", "d", "d", "d"),
  dim2 = c(1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2),
  value = 1:12
)

# function to apply
AddLags <- function(df, lags_vector, target_col, date_col) {
  temp_lags <- map_dfc(
    lags_vector,
    ~ df %>%
      arrange({{ date_col }}) %>%
      transmute(
        across(contains(target_col), lag, .x, .names = '{col}_lag_{ifelse(.x<10,paste0("0",.x),.x)}')
      )
  )
  return(temp_lags)
}

reprex package (v2.0.1)

创建于 2022-01-13

正如@Limey 所建议的,一种可能的方法是使用 group_map 函数:

results_df <- data.frame()
results_df <- 
  bind_rows(
    df %>% 
      group_by(dim1,dim2) %>% 
      group_map(~AddLags(.,c(1,2),'value',date))
  )

预期结果将是:

   value_lag_01 value_lag_02
          <int>        <int>
 1           NA           NA
 2            1           NA
 3            2            1
 4           NA           NA
 5            4           NA
 6            5            4
 7           NA           NA
 8            7           NA
 9            8            7
10           NA           NA
11           10           NA
12           11           10

但是,我个人会选择@jpdugo17 方法