如何用R的lubridate计算旅行者在指定时段内的停留天数？

你可以通过以下步骤计算旅行者在指定时段内的总停留天数：

1. 定义所有目标观测时段：将每年5月1日至9月1日的时间段转换为lubridate区间对象
2. 计算单个时段的重叠天数：通过求旅行者停留区间与目标时段的交集，得到重叠的时长并转换为天数
3. 求和所有重叠天数：将各时段的重叠天数累加得到总天数

完整代码示例 #

library(lubridate)
library(purrr) # 用于批量处理时段

# 旅行者的起止日期
traveler1_startdate <- mdy("3/11/2020")
traveler1_enddate <- mdy("7/1/2023")

# 创建旅行者的完整停留区间
travel_interval <- interval(traveler1_startdate, traveler1_enddate)

# 定义需要统计的目标年份，生成对应的5.1-9.1区间
target_years <- 2020:2023
target_intervals <- map(target_years, function(year) {
  start <- ymd(paste(year, "05", "01", sep = "-"))
  end <- ymd(paste(year, "09", "01", sep = "-"))
  interval(start, end)
})

# 计算每个目标时段的重叠天数，空交集返回0
overlap_days <- map_dbl(target_intervals, function(interval) {
  overlap <- intersect(travel_interval, interval)
  if (is.na(overlap)) 0 else time_length(overlap, "days")
})

# 总停留天数
total_overlap <- sum(overlap_days)
total_overlap

代码说明 #

• interval()：创建时间区间对象，用于表示连续的时间段
• intersect()：计算两个区间的交集，得到旅行者在该目标时段内实际停留的时间段
• time_length(..., "days")：将交集区间转换为天数，若两个区间无交集（返回NA），则记为0天
• map()/map_dbl()：批量处理多个目标时段，避免重复代码

运行上述代码后，total_overlap即为该旅行者在所有指定时段内的总停留天数。

内容的提问来源于stack exchange，提问作者SugarRay