r - group_modify() 下的 case_when() 中的 %within% 不起作用
问题描述
我有以下类型的数据:
library(tidyverse)
library(lubridate)
data <- tibble(a = c(1, 1, 2, 3, 3, 3, 3),
b = c('x', 'y', 'z', 'z', 'z', 'z', 'z'),
c = c('ps', 'ps', 'qs', 'rs', 'rs', 'rs', 'rs'),
d = c(100, 200, 300, 400, 500, 600, 700),
strt = ymd(c('2019-03-20', '2020-01-01', '2018-01-02', '2020-05-01', '2016-01-01', '2020-03-01', '2020-01-01')),
fnsh = ymd(c('3019-03-20', '3020-01-01', '3018-01-02', '2020-06-01', '2016-05-01', '2020-04-01', '2020-06-10')))
我正在使用基于变量 a、b 和 c(即data %>% group_by(a, b, c)
)的分组操作group_modify()
。对于每个组,我需要找到在过去一年内具有真正开始日期的行。如果 strt 不在组中任何其他行的 strt 和 fnsh 之间,则它是真实的。我目前的做法是:
test <- data %>%
group_by(a, b, c) %>%
group_modify(function(.x, .y) {
.x %>%
mutate(startLatestYear = case_when(strt > today(tzone = 'CET') - years(1) &
strt <= today(tzone = 'CET') &
!strt %within% (.x %>%
mutate(pushInterval = interval(strt + days(1), fnsh)) %>%
select(pushInterval)) ~ 1,
TRUE ~ 0))}) %>%
ungroup()
这种方法给出:
data <- tibble(a = c(1, 1, 2, 3, 3, 3, 3),
b = c('x', 'y', 'z', 'z', 'z', 'z', 'z'),
c = c('ps', 'ps', 'qs', 'rs', 'rs', 'rs', 'rs'),
d = c(100, 200, 300, 400, 500, 600, 700),
strt = ymd(c('2019-03-20', '2020-01-01', '2018-01-02', '2020-05-01', '2016-01-01', '2020-03-01', '2020-01-01')),
fnsh = ymd(c('3019-03-20', '3020-01-01', '3018-01-02', '2020-06-01', '2016-05-01', '2020-04-01', '2020-06-10')),
startLatestYear = c(0, 1, 0, 1, 0, 1, 1))
需要的是:
data <- tibble(a = c(1, 1, 2, 3, 3, 3, 3),
b = c('x', 'y', 'z', 'z', 'z', 'z', 'z'),
c = c('ps', 'ps', 'qs', 'rs', 'rs', 'rs', 'rs'),
d = c(100, 200, 300, 400, 500, 600, 700),
strt = ymd(c('2019-03-20', '2020-01-01', '2018-01-02', '2020-05-01', '2016-01-01', '2020-03-01', '2020-01-01')),
fnsh = ymd(c('3019-03-20', '3020-01-01', '3018-01-02', '2020-06-01', '2016-05-01', '2020-04-01', '2020-06-10')),
startLatestYear = c(0, 1, 0, 0, 0, 0, 1))
该组基于a == 3
,b == 'z'
并且c == 'rs'
有一行(最后一行)应该是该组中 startLatestYear 中具有 1 的唯一行。最后一行是组中唯一在最近一年内具有 strt 并且 strt 在组中其他行的间隔之外的行。
目前使用的前两个条件case_when()
似乎有效。使用的第三个条件%within%
似乎不起作用。使用条件如何%within%
发挥作用?或者如何实施替代解决方案?
PS:我曾尝试在分组 tibble 之前制作 pushInterval。这样做会为 startLatestYear 生成相同的列,但该操作会导致bind_rows_()
剥离间隔属性的“问题”。因此,当前产生 pushInterval 的解决方案是动态的。
解决方案
我认为您不需要使用group_modify
,这在一个简单的组中有效mutate
:
data %>%
group_by(a, b, c) %>%
mutate(x = +(purrr::map_lgl(strt, ~ sum(strt <= .x & .x <= fnsh) < 2) &
difftime(Sys.time(), strt, "days") < 365)) %>%
ungroup()
# # A tibble: 7 x 7
# a b c d strt fnsh x
# <dbl> <chr> <chr> <dbl> <date> <date> <int>
# 1 1 x ps 100 2019-03-20 3019-03-20 0
# 2 1 y ps 200 2020-01-01 3020-01-01 1
# 3 2 z qs 300 2018-01-02 3018-01-02 0
# 4 3 z rs 400 2020-05-01 2020-06-01 0
# 5 3 z rs 500 2016-01-01 2016-05-01 0
# 6 3 z rs 600 2020-03-01 2020-04-01 0
# 7 3 z rs 700 2020-01-01 2020-06-10 1
.x
是作为第一个参数传递给 的参数的占位符map_lgl
。在这种情况下,它也是strt
,但让我们暂时忘记它。
在波浪号函数内部,strt
指的是整个向量,并且.x
指的是每个单独的strt
值(长度始终为 1)。strt <= .x
第一次是有效的strt <= strt[1]
。只是计算有sum
多少事件是真实的。(应该总是有一个,因为一个数字总是在它自己的范围内。)
推荐阅读
- javascript - 针对父属性的猫鼬子模式验证
- python - sqlalchemy:从表中选择查询中的列
- c# - 以编程方式从 forms.office.com 检索数据
- android - 错误:不兼容的类型:PesoR_Male 无法转换为上下文(AlertDialog 自定义)
- html - 不会在 AngularJs 中访问特定的 URL this.router.navigate()
- mongodb - 无法使用 monstache 使用 URL 连接到 mongodb
- r - SublimeREPL: OSError: [WinError 6] 句柄无效
- ruby - Ruby 2.3.1 如何在config.rb 中设置项目根目录
- polymer - 为什么我的 this._setPropertyName 函数不起作用?
- javascript - 用于替换 Kendo Grid 中数据的 Javascript 函数