如果日期小于 R 中特定行的日期，则返回 NA

发布于08月04日

我的数据是这样的:

dput(head(VI_v2_KRR05,28))
structure(list(ID = c("AUR", "AUR", "AUR", "AUR", "AUR", "AUR", 
"AUR", "AUR", "AUR", "AUR", "AUR", "LAM", "LAM", "LAM", "LAM", 
"LAM", "LAM", "LAM", "LAM", "LAM", "LAM", "P0", "P0", "P0", "P0", 
"P0", "P01", "P01"), EVI_SOS = structure(c(16440, 16805, 17124, 
17421, 17599, 17851, 18216, 18403, NA, NA, NA, 16272, 16406, 
16637, 16771, 17148, 17516, 17725, 18022, 18210, NA, 16692, 16845, 
17058, 17212, NA, 16717, 17077), class = "Date"), NDVI_SOS = structure(c(16436, 
16801, 17110, 17420, 17607, 17841, 18196, 18402, NA, NA, NA, 
16270, 16380, 16635, 16745, 17139, 17274, 17522, 17731, 18027, 
18198, 16480, 16683, 17055, NA, NA, 16712, 17076), class = "Date"), 
    NIRv_SOS = structure(c(16424, 16557, 16789, 16922, 17221, 
    17379, 17607, 17821, 17931, 18214, 18400, 16274, 16404, 16639, 
    16769, 17145, 17519, 17727, 18028, 18208, NA, 16695, 16848, 
    17061, 17218, NA, 16720, 17084), class = "Date"), kNDVI_SOS = structure(c(16542, 
    16799, 17212, 17431, 17607, 17854, 18003, 18219, 18401, NA, 
    NA, 16282, 16647, 17139, 17516, 17733, 18026, 18205, NA, 
    NA, NA, 16474, 16706, 16847, 17073, 17209, 16712, 17075), class = "Date"), 
    EVI_EOS = structure(c(16766, 17084, 17356, 17577, 17812, 
    18098, 18385, 18585, NA, NA, NA, 16375, 16618, 16740, 17099, 
    17452, 17707, 17841, 18183, 18549, NA, 16825, 17012, 17193, 
    17396, NA, 17022, 17398), class = "Date"), NDVI_EOS = structure(c(16768, 
    17075, 17375, 17590, 17812, 18095, 18385, 18568, NA, NA, 
    NA, 16359, 16616, 16724, 16998, 17248, 17490, 17714, 17961, 
    18177, 18460, 16648, 17020, 17393, NA, NA, 17028, 17392), class = "Date"), 
    NIRv_EOS = structure(c(16533, 16768, 16899, 17191, 17358, 
    17590, 17798, 17920, 18104, 18379, 18547, 16373, 16619, 16738, 
    17086, 17441, 17709, 17839, 18182, 18471, NA, 16828, 17008, 
    17198, 17387, NA, 17028, 17398), class = "Date"), kNDVI_EOS = structure(c(16764, 
    17184, 17356, 17592, 17805, 17993, 18088, 18381, 18562, NA, 
    NA, 16620, 17007, 17439, 17715, 17845, 18180, 18457, NA, 
    NA, NA, 16625, 16833, 16991, 17192, 17366, 17026, 17391), class = "Date")), row.names = c(NA, 
-28L), class = c("tbl_df", "tbl", "data.frame"))

我想返回小于2015-09-04的ID AUR、2015-09-04的ID LAM、2016-01-02的ID P0和2016-01-09的ID P01的日期.

任何帮助都将不胜感激.

library(dplyr) library(collapse) out <- df1 %>% mutate(across(where(is_date), ~ case_when(.x > as.Date("2015-09-04") & ID == 'AUR' ~ .x, ID == "LAM" & .x > as.Date("2015-09-04 ") ~ .x, ID == "P0" & .x > as.Date("2016-01-02 ") ~ .x, ID == "P01" & .x > as.Date("2016-01-09") ~ .x)))

> df1 %>% filter(ID == "LAM") # A tibble: 10 × 9 ID EVI_SOS NDVI_SOS NIRv_SOS kNDVI_SOS EVI_EOS NDVI_EOS NIRv_EOS kNDVI_EOS <chr> <date> <date> <date> <date> <date> <date> <date> <date> 1 LAM 2014-07-21 2014-07-19 2014-07-23 2014-07-31 2014-11-01 2014-10-16 2014-10-30 2015-07-04 2 LAM 2014-12-02 2014-11-06 2014-11-30 2015-07-31 2015-07-02 2015-06-30 2015-07-03 2016-07-25 3 LAM 2015-07-21 2015-07-19 2015-07-23 2016-12-04 2015-11-01 2015-10-16 2015-10-30 2017-09-30 4 LAM 2015-12-02 2015-11-06 2015-11-30 2017-12-16 2016-10-25 2016-07-16 2016-10-12 2018-07-03 5 LAM 2016-12-13 2016-12-04 2016-12-10 2018-07-21 2017-10-13 2017-03-23 2017-10-02 2018-11-10 6 LAM 2017-12-16 2017-04-18 2017-12-19 2019-05-10 2018-06-25 2017-11-20 2018-06-27 2019-10-11 7 LAM 2018-07-13 2017-12-22 2018-07-15 2019-11-05 2018-11-06 2018-07-02 2018-11-04 2020-07-14 8 LAM 2019-05-06 2018-07-19 2019-05-12 NA 2019-10-14 2019-03-06 2019-10-13 NA 9 LAM 2019-11-10 2019-05-11 2019-11-08 NA 2020-10-14 2019-10-08 2020-07-28 NA 10 LAM NA 2019-10-29 NA NA NA 2020-07-17 NA NA > out %>% filter(ID == "LAM") # A tibble: 10 × 9 ID EVI_SOS NDVI_SOS NIRv_SOS kNDVI_SOS EVI_EOS NDVI_EOS NIRv_EOS kNDVI_EOS <chr> <date> <date> <date> <date> <date> <date> <date> <date> 1 LAM NA NA NA NA NA NA NA NA 2 LAM NA NA NA NA NA NA NA 2016-07-25 3 LAM NA NA NA 2016-12-04 2015-11-01 2015-10-16 2015-10-30 2017-09-30 4 LAM 2015-12-02 2015-11-06 2015-11-30 2017-12-16 2016-10-25 2016-07-16 2016-10-12 2018-07-03 5 LAM 2016-12-13 2016-12-04 2016-12-10 2018-07-21 2017-10-13 2017-03-23 2017-10-02 2018-11-10 6 LAM 2017-12-16 2017-04-18 2017-12-19 2019-05-10 2018-06-25 2017-11-20 2018-06-27 2019-10-11 7 LAM 2018-07-13 2017-12-22 2018-07-15 2019-11-05 2018-11-06 2018-07-02 2018-11-04 2020-07-14 8 LAM 2019-05-06 2018-07-19 2019-05-12 NA 2019-10-14 2019-03-06 2019-10-13 NA 9 LAM 2019-11-10 2019-05-11 2019-11-08 NA 2020-10-14 2019-10-08 2020-07-28 NA 10 LAM NA 2019-10-29 NA NA NA 2020-07-17 NA NA

library(data.table) df2 <- data.table(ID = c("AUR", "LAM", "PO", "P01"), dates = as.Date(c("2015-09-04", "2015-09-04", "2016-01-02", "2016-01-09"))) out1 <- copy(df1) nm1 <- grep("_SOS$", names(out1), value = TRUE) setDT(out1)[df2, (nm1) := lapply(.SD, \(x) fcase(x > dates, x)), on = .(ID), .SDcols = nm1]

df1 %>% left_join(df2) %>% mutate(across(c(where(is_date), -dates), ~ case_when(.x > dates ~ .x)), dates = NULL) # A tibble: 28 × 9 ID EVI_SOS NDVI_SOS NIRv_SOS kNDVI_SOS EVI_EOS NDVI_EOS NIRv_EOS kNDVI_EOS <chr> <date> <date> <date> <date> <date> <date> <date> <date> 1 AUR NA NA NA NA 2015-11-27 2015-11-29 NA 2015-11-25 2 AUR 2016-01-05 2016-01-01 NA 2015-12-30 2016-10-10 2016-10-01 2015-11-29 2017-01-18 3 AUR 2016-11-19 2016-11-05 2015-12-20 2017-02-15 2017-07-09 2017-07-28 2016-04-08 2017-07-09 4 AUR 2017-09-12 2017-09-11 2016-05-01 2017-09-22 2018-02-15 2018-02-28 2017-01-25 2018-03-02 5 AUR 2018-03-09 2018-03-17 2017-02-24 2018-03-17 2018-10-08 2018-10-08 2017-07-11 2018-10-01 6 AUR 2018-11-16 2018-11-06 2017-08-01 2018-11-19 2019-07-21 2019-07-18 2018-02-28 2019-04-07 7 AUR 2019-11-16 2019-10-27 2018-03-17 2019-04-17 2020-05-03 2020-05-03 2018-09-24 2019-07-11 8 AUR 2020-05-21 2020-05-20 2018-10-17 2019-11-19 2020-11-19 2020-11-02 2019-01-24 2020-04-29 9 AUR NA NA 2019-02-04 2020-05-19 NA NA 2019-07-27 2020-10-27 10 AUR NA NA 2019-11-14 NA NA NA 2020-04-27 NA

如果日期小于 R 中特定行的日期，则返回 NA

推荐答案

R相关问答推荐

如果R上的不同时期之间的值发生了变化，则创建假人

NA仅省略具有NA的 Select 行

R箱形图gplot 2 4组但6个参数

在垂直轴中包含多个ggplot2图中的平均值

格点中指数、双曲和反双曲模型曲线的正确绘制

基于多列将值链接到NA

RStudio中相关数据的分组箱形图

如何将旋转后的NetCDF转换回正常的纬度/经度网格，并使用R？

如何在Chart_Series()中更改轴值的 colored颜色？

如何在ggplot图中找到第二轴的比例

一小时满足条件的日期的 Select

在R中使用Scale_y_Break后更改y轴标签

使用带有OR条件的grepl过滤字符串

是否有新方法来更改Facet_WRAP(Ggplot2)中条文本的文本 colored颜色？

如果COLSUM为>；0，则COLNAME为向量

我如何使用tidyselect来传递一个符号数组，比如Pivot_Long？

使用gt_summary是否有一种方法来限制每个变量集进行配对比较？

如何筛选截止年份之前最后一个测量年度的所有观测值以及截止年份之后所有年份的所有观测值

如何在使用因子时获得Sankey图的Scale_Fill_Viridis的全范围

条形图中的条形图没有try 赋予它们的 colored颜色