インターベンション直後に起きやすいイベント

• 少し具体的にする
• 今、１００万人の集団があるとする
• そのうち、半分が一斉にインターベンションを受ける
• インターベンションを受けた人だけに着目し、イベントが起きたら、「インターベンションの○日後にイベントが起きた」と記録することにする
• このイベントは、１００万人の集団で、１年に１００件発生し、それは、どの日もどの人も同じ確率で起きるという

Nall<-10^6
Y<-365
D<-150 # 年間のイベント数
p<-D/Nall/Y

# Nallのうち0.5がインターベンションを受けるとする
N <-Nall*0.5


# あるインターベンションを受けた人が、インターベンション後、(d-1)〜d日にイベントに遭遇する確率は
d<-1:Y
Prd<-pexp(d,p)-pexp(d-1,p)

plot(d,Prd,ylim=c(0,max(Prd)))
sum(Prd)


# d日目まで、だれにもイベントが誰にも起きない確率は

FreeProb<-1-pexp(d,p,lower.tail=FALSE)^N
showD<-10
plot(1:showD,FreeProb[1:showD])