找到生命线危险函数的局部最大值

Question

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

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')

感谢任何愿意帮助我的人。

斯特凡诺

Answer 1

类似na.smoothed_hazard_(bandwidth)。