我使用 dplyr
和 ifelse
根据两个条件和以下数据创建一个新列。
dat <- structure(list(GenIndID = c("BHS_034", "BHS_034", "BHS_068",
"BHS_068", "BHS_068", "BHS_068", "BHS_068", "BHS_068", "BHS_068",
"BHS_068", "BHS_068"), IndID = c("BHS_034_A", "BHS_034_A", "BHS_068_A",
"BHS_068_A", "BHS_068_A", "BHS_068_A", "BHS_068_A", "BHS_068_A",
"BHS_068_A", "BHS_068_A", "BHS_068_A"), Fate = c("Mort", "Mort",
"Alive", "Alive", "Alive", "Alive", "Alive", "Alive", "Alive",
"Alive", "Alive"), Status = c("Alive", "Mort", "Alive", "Alive",
"MIA", "Alive", "MIA", "Alive", "MIA", "Alive", "Alive"), Type = c("Linked",
"Linked", "SOB", "SOB", "SOB", "SOB", "SOB", "SOB", "SOB", "SOB",
"SOB"), SurveyID = c("GYA13-1", "GYA14-1", "GYA13-1", "GYA14-1",
"GYA14-2", "GYA15-1", "GYA16-1", "GYA16-2", "GYA17-1", "GYA17-3",
"GYA15-2"), SurveyDt = structure(c(1379570400, 1407477600, 1379570400,
1407477600, 1409896800, NA, 1462946400, 1474351200, 1495519200,
1507010400, 1441951200), tzone = "", class = c("POSIXct", "POSIXt"
))), row.names = c(NA, 11L), .Names = c("GenIndID", "IndID",
"Fate", "Status", "Type", "SurveyID", "SurveyDt"), class = "data.frame")
> dat
GenIndID IndID Fate Status Type SurveyID SurveyDt
1 BHS_034 BHS_034_A Mort Alive Linked GYA13-1 2013-09-19
2 BHS_034 BHS_034_A Mort Mort Linked GYA14-1 2014-08-08
3 BHS_068 BHS_068_A Alive Alive SOB GYA13-1 2013-09-19
4 BHS_068 BHS_068_A Alive Alive SOB GYA14-1 2014-08-08
5 BHS_068 BHS_068_A Alive MIA SOB GYA14-2 2014-09-05
6 BHS_068 BHS_068_A Alive Alive SOB GYA15-1 <NA>
7 BHS_068 BHS_068_A Alive MIA SOB GYA16-1 2016-05-11
8 BHS_068 BHS_068_A Alive Alive SOB GYA16-2 2016-09-20
9 BHS_068 BHS_068_A Alive MIA SOB GYA17-1 2017-05-23
10 BHS_068 BHS_068_A Alive Alive SOB GYA17-3 2017-10-03
11 BHS_068 BHS_068_A Alive Alive SOB GYA15-2 2015-09-11
更具体地说,按 GenIndID
分组,我想根据 Type
的两个条件创建一个新的日期字段,该字段是最大 SurveyDt
和命运
。此外,我希望最大日期仅在 Status == Alive
时评估 SurveyDt
。下面的代码生成所有 NA
值,而不是满足所有指定条件的 BHS_068
所描述的日期字段。
我最近看到 case_when
这可能适合这里,但我无法正确实现它。
dat %>% group_by(GenIndID) %>%
mutate(NewDat = as.POSIXct(ifelse(Type == "SOB" & Fate == "Alive", max(SurveyDt[Status == "Alive"], na.rm = F), NA),
origin='1970-01-01', na.rm=T)) %>%
as.data.frame()
如有任何建议,我们将不胜感激。
最佳答案
如果您想坚持使用 dplyr
并使用 case_when
,您必须确保每个 case 语句的值都是相同的类型。
在这种情况下,您的 TRUE 值将是日期时间,因此您必须将默认值包装在 as.POSIXct
中,将其设置为日期时间。
dat %>%
group_by(GenIndID) %>%
mutate(NewDat = case_when(Type == "SOB" & Fate == "Alive" ~ max(SurveyDt[Status == "Alive"], na.rm = TRUE),
TRUE ~ as.POSIXct(NA, origin = "1970-01-01")))
使用ifelse
dat %>%
group_by(GenIndID) %>%
mutate(NewDat = ifelse(Type == "SOB" & Fate == "Alive",
max(SurveyDt[Status == "Alive"], na.rm = TRUE),
as.POSIXct(NA, origin = "1970-01-01")))
关于r - dplyr ifelse 语句中的嵌套条件,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/46696348/