我试图复制图1中的kaplan-meier表。图为:
这是我写的代码:
# Python code to create the above Kaplan Meier curve
from lifelines import KaplanMeierFitter
import pandas as pd
df = pd.DataFrame({
'T':[0,0,0,0,0,0,2.5,2.5,2.5,2.5,2.5,4,4,4,4,4,5,5,5,6,6],
'E':[0,0,0,0,0,0,1,0,0,0,0,1,1,0,0,0,1,0,0,0,0],
})
## create a kmf object
kmf = KaplanMeierFitter()
## Fit the data into the model
kmf.fit(df['T'], df['E'],label='Kaplan Meier Estimate')
## Create an estimate
kmf.plot(ci_show=False)
我的输出曲线不同(见比例):
当我打印生存函数时,它是不同的:
Kaplan Meier Estimate
timeline
0.0 1.0000
2.5 0.9375
4.0 0.7500
5.0 0.6000
6.0 0.6000
我假设我没有正确地将数据转换成Dataframe(可能?)。我试图搞乱Dataframe,将1事件添加到时间帧的开始和结束,但这并不重要。有人能告诉我如何复制我正在研究的例子吗?
暂无答案!
目前还没有任何答案,快来回答吧!