
196
13 章 生存分析
も急峻になりますが、それが最も普通の妊娠期間です。
13.2
ハザード関数
生存関数からハザード関数(
hazard function
)
†
が導かれます。妊娠期間については、
ハザード関数が時間
t
に対して、
t
まで妊娠が続いたが
t
で終わってしまった割合を
対応させます。詳しくは次の数式で表されます。
λ
(t)
=
S(t)
−
S(t
+
1)
S(t)
分子が
t
で終わった生存期間の割合で、
PMF(t)
でもあります。
SurvivalFunction
には、ハザード関数を計算する
MakeHazard
があります。
# class SurvivalFunction
def MakeHazard(self, label=''):
ss = self.ss
lams = {}
for i, t in enumerate(self.ts[:-1]):
hazard = (ss[i] - ss[i+1]) / ss[i]
lams[t] = hazard
return HazardFunction(lams, label=label)
HazardFunction
は、
pandas
の
Series
のラッパーです。
class HazardFunction(object):
def __init__(self, d, label= ''):
self.series = pandas.Series(d)
self.label = label
d
はディクショナリなど、他の
Series
を含めて
Series
を初期化できる型であれば ...