检查日期列表是否在 r 的日期范围内

Question

下面显示了我的数据的三个主要列（squirrel_id = 唯一的个人 ID，byear = 出生年份，dyear = 死亡年份）：

> summary(complete)
  squirrel_id        byear          dyear     
 Min.   :  416   Min.   :1989   Min.   :1989  
 1st Qu.: 4152   1st Qu.:1997   1st Qu.:1998  
 Median : 7870   Median :2003   Median :2004  
 Mean   :10419   Mean   :2004   Mean   :2004  
 3rd Qu.:16126   3rd Qu.:2011   3rd Qu.:2012  
 Max.   :23327   Max.   :2017   Max.   :2017 

我正在尝试将第二条数据（如下所示）合并到上述数据集中。

mast.yr<-c("1993", "1998", "2005", "2010", "2014")

我正在尝试做两件事：

1. 添加一列，说明个人 (squirrel_id) 在任何 mast.yr 年期间是否还活着（dyear- byear = 活着的年限（包括 byear 和 dyear）。
2. 添加另一列，计算每个人 (squirrel_id) 在一生中经历的 mast.yr 年数 (dyear- byear = 生存年数范围（包括 byear 和 dyear）。李>

为了生成第一列，我一直在使用 dplyr 包中的 mutate 函数，但我只能让它分别为 byear 和 dyear 工作，如下所示：

complete <- complete %>%
    mutate (mast = ifelse (byear %in% c("1993", "1998", "2005", "2010", "2014"), 1, 0), 
    mast = ifelse (dyear %in% c("1993", "1998", "2005", "2010", "2014"), 1, 0))) 

但它并没有给出所需的输出，因为它是单独考虑byear 和dyear，而不是作为一个连续的时间段。我已经尝试了here 和here 发布的解决方案，但没有运气。

任何建议将不胜感激！

可以在here 找到我的数据副本。为了将来的重现性，这里有一个示例：

> head(DF)

> squirrel_id  byear dyear
  <dbl>        <int> <dbl>
  6715         2006  2006
  22274        2016  2017
  20445        2014  2017
  19528        2013  2013
  2674         1995  1995
  1419         1992  1993

Answer 1

# put target years in a table
mastDF = data_frame(year = as.integer(mast.yr))

# count based on conditions    
dat %>% 
  mutate(in_mast = count_matches(., mastDF, year >= byear, year <= dyear) > 0) %>%
  as.tbl

# A tibble: 100 x 4
   squirrel_id byear dyear in_mast
         <int> <int> <int>   <lgl>
 1        6715  2006  2006   FALSE
 2       22274  2016  2017   FALSE
 3       20445  2014  2017    TRUE
 4       19528  2013  2013   FALSE
 5        2674  1995  1995   FALSE
 6        1419  1992  1992   FALSE
 7       15014  2004  2004   FALSE
 8       10946  2009  2012    TRUE
 9        4369  1998  1999    TRUE
10        4344  1992  1999    TRUE
# ... with 90 more rows

count_matches 是一个辅助函数：

library(data.table)
count_matches = function(DF, targetDF, ...){
  onexpr = substitute(list(...))
  data.table(targetDF)[data.table(DF), on=eval(onexpr), allow.cart=TRUE, .N, by=.EACHI]$N
}

---

如果您想要计数和计数是否非零，可以通过将其分解为mutate 参数序列来完成：

dat %>% 
  mutate(
    n_mast = count_matches(., mastDF, year >= byear, year <= dyear),
    in_mast = n_mast > 0
) %>% as.tbl

# A tibble: 6 x 5
  squirrel_id byear dyear n_mast in_mast
        <int> <int> <int>  <int>   <lgl>
1        6715  2006  2006      0   FALSE
2       22274  2016  2017      0   FALSE
3       20445  2014  2017      1    TRUE
4       19528  2013  2013      0   FALSE
5        2674  1995  1995      0   FALSE
6        1419  1992  1993      1    TRUE

Answer 2

虽然@Frank 已经提供了一个优雅的解决方案，但sqldf 为非 equii 连接提供了更简单的方法。使用sqldf的解决方案可以是：

mast.yr<-c("1993", "1998", "2005", "2010", "2014")

mastDf <- data.frame(year = as.integer(mast.yr))


library(sqldf)


sqldf("select dat.*, IFNULL(Mast.inMast,0) as n_Mast, IFNULL(Mast.inMast,0) >0 as inMast
       from dat left outer join 
      (select *, count(squirrel_id) as inMast 
      from dat, mastDf 
      where mastDf.year between dat.byear AND dat.dyear
      group by squirrel_id) Mast on
      dat.squirrel_id = Mast.squirrel_id")

#     squirrel_id byear dyear n_Mast inMast
# 1          6715  2006  2006      0      0
# 2         22274  2016  2017      0      0
# 3         20445  2014  2017      1      1
# 4         19528  2013  2013      0      0
# 5          2674  1995  1995      0      0
# 6          1419  1992  1992      0      0
# 7         15014  2004  2004      0      0
# 8         10946  2009  2012      1      1
# 9          4369  1998  1999      1      1
# 10         4344  1992  1999      2      1
#....90 more rows