
201
201
13.6 生存曲線を推定する
図
13-2
初婚年齢のハザード関数(上)と未婚状態の生存関数(下)
13.6
生存曲線を推定する
ハザード関数が得られれば、生存曲線を推定できます。時間
t
を過ぎて生存する機
会は、
t
に至るまでずっと生存している機会ですから、ハザード関数の補関数の積と
なります。
[1
−λ
(0)][1
−λ
(1)]...[1
−λ
(t)]
HazardFunction
クラスは、この積を計算する
MakeSurvival
を提供します。
# class HazardFunction:
def MakeSurvival(self):
ts = self.series.index
ss = (1 - self.series).cumprod()
cdf = thinkstats2.Cdf(ts, 1-ss)
sf = SurvivalFunction(cdf)
return sf
年齢(歳)