我正在尝试找到一种方法来检测危险函数的局部最大值点(使用生命线库绘制)。
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')
感谢任何会帮助我的人。
斯特凡诺