カプランマイヤー推定量について
概要
- ノンパラメトリックの回帰の一種
- 生存する確率を求める類に使用することができる
数式
\[\hat{S} = \begin{cases} 1 \\ \prod \frac{n_i - d_i}{n_i} \end{cases}\]google colabによる再現
import numpy as np
import pandas as pd
y = np.random.random(25)
df = pd.DataFrame({"死亡時期": y})
probs = []
for i in np.arange(0, 1, 0.0001):
prob = len(df.query(f'死亡時期 > {i}')) / len(df)
probs.append(prob)
r=pd.DataFrame(dict(prob=probs))
plt.figure(figsize=(15, 15)) # ここを大きくする
ax = sns.scatterplot(data=r, x=r.index, y=r.prob)
ax.set(title="カプランマイヤー推定量", xlabel="index", ylabel="推定量")
ax
