Stata：非唯一日期的多日期范围重叠天数计算问题

问题背景 #

基于Stack Overflow回答《Stata: Number of overlapping days within multiple date ranges》，需计算多日期范围内的重叠天数，同时要将日期范围展开为完整日期列表。原解决方案在处理非唯一日期（即不同日期范围存在重叠）时出现报错，原代码如下：

/* Fake Data */
input str9(start  end)      
"01jan2000" "04apr2000"
"30may2000" "15aug2000"
"07jan2000" "31dec2000"
end

foreach var of varlist start end {
    gen d = date(`var', "DMY")
    drop `var'
    gen `var' = d
    format %td `var'
    drop d
}

/* Count Overlapping Days */
rename (start end) date=
gen spell = _n
reshape long date, i(spell) j(range) string
drop range
xtset spell date, delta(1 day)
tsfill
bys date: keep if _N>1

调整后的解决方案 #

原代码依赖xtset构建时间序列，要求每个spell-date组合唯一，当日期范围存在重叠或同一序列内有重复日期时会触发报错。改用expand生成完整日期序列的方式，可兼容非唯一日期场景，调整后代码如下：

/* Fake Data */
input str9(start  end)      
"01jan2000" "04apr2000"
"30may2000" "15aug2000"
"07jan2000" "31dec2000"
end

/* 转换字符串日期为数值型日期 */
foreach var of varlist start end {
    gen d = date(`var', "DMY")
    drop `var'
    gen `var' = d
    format %td `var'
    drop d
}

/* 生成每个日期范围内的所有日期 */
gen spell = _n
gen days_in_spell = end - start + 1
expand days_in_spell
bys spell: gen date = start + _n - 1
format %td date

/* 筛选并统计重叠日期 */
bys date: gen overlap_count = _N
// 保留被至少2个日期范围覆盖的日期
keep if overlap_count > 1

// 若需统计总重叠天数，去重后计数
duplicates drop date, force
display "总重叠天数: " _N

代码说明 #

1. 日期转换：保持原逻辑，将字符串日期转为Stata数值型日期并设置格式。
2. 生成完整日期序列：通过days_in_spell计算每个日期范围的天数，用expand扩展观测，再按spell分组生成该范围内的所有日期。
3. 筛选重叠日期：按date分组统计覆盖次数，保留覆盖次数大于1的日期；若需总重叠天数，去重后用count或display输出结果。

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