0

我正在尝试找到一种方法来检测危险函数的局部最大值点(使用生命线库绘制)。

Lifelines 库使您能够直观地绘制危害函数( plot.hazard() ),但是 - 如上所述 - 它只是函数的视觉表示,而不是函数本身。这意味着您显然不能将此函数作为参数传递给方法(即:来自 scipy.signals 库的 find_peaks)。

我想知道是否有人有解决方案来检测这些问题,请在下面找到我如何计算我的危险函数(在我突出显示局部最大值的最右边的图):

fig, (ax1, ax2, ax3) = plt.subplots(1, 3,figsize=(18,8))

km.fit(duration, event_observed = event)
kap = km.plot_survival_function(ax=ax1,legend="")
kap.set_title("Kaplan-Meier",fontsize=15)

na.fit(duration, event_observed = event)

na1 = na.plot_cumulative_hazard(ax=ax2,legend="",color="green")
na1.set_title("Cumulative Hazard",fontsize=15)

bandwidth = 5
na2 = na.plot_hazard(ax=ax3,xlim=(0,200),ylim=(0,0.03),bandwidth=bandwidth,color="orange",legend="")
na2.set_title("Hazard Function",fontsize=15)
ax3.set_xlabel('timeline',fontsize = 10)
plt.grid(linestyle="-.", color='grey')

在此处输入图像描述

感谢任何会帮助我的人。

斯特凡诺

4

1 回答 1

-1

类似的东西na.smoothed_hazard_(bandwidth)

于 2021-08-18T19:36:09.493 回答