R перекрытие периодов времени

Question

С пакетом "lubridate" в R, я могу узнать, перекрываются ли два периода времени. но есть ли эффективный способ подсчитать, сколько дней они перекрыли. (Например, сколько дней женщины курили во время беременности. Период беременности и период курения могут полностью, частично или не перекрываться)

Вот пример с тремя женщинами:

preg_start<-as.Date(c("2011-01-01","2012-01-01","2013-01-01"))
preg_end<-preg_start+270 # end after 9 months
smoke_start<-as.Date(c("2011-02-01","2012-08-01","2014-01-01"))
smoke_end<-smoke_start+100 # all three smoked 100 days

data<-data.frame(cbind(preg_start,preg_end,smoke_start,smoke_end))

Я хочу добавить переменную, в которой говорится, что первая женщина курила 100 дней во время беременности, вторая курила 30 дней, а третья не курила во время беременности.

Answer 1

Использование interval создать временные интервалы для беременности и курения. Затем рассчитайте intersect из этих интервалов. Из этого вы можете рассчитать period в днях.

library("lubridate")
preg_start<-as.Date(c("2011-01-01","2012-01-01","2013-01-01"))
preg_end<-preg_start+270 # end after 9 months
smoke_start<-as.Date(c("2011-02-01","2012-08-01","2014-01-01"))
smoke_end<-smoke_start+100 # all three smoked 100 days

smoke <- new_interval(smoke_start, smoke_end, tzone="UTC")
preg <- new_interval(preg_start, preg_end, tzone="UTC")
day(as.period(intersect(smoke, preg), "days"))

Я получаю 100, 57 и 0 дней курения во время беременности.