Kaplan-Meier 生存估计的核或贝塞尔曲线平滑

Question

我是一名土木工程研究生，从事桥梁生存建模工作。我正在使用带有 Kaplan-Meier 方法的 SAS 来分析数据。我想使用内核平滑或贝塞尔平滑用平滑曲线覆盖生存概率图。但是由于这些方法在统计包中不可用，我将生存概率导出到 Excel 以便手动生成（我这样做可能是错误的）。经过几次尝试，我无法弄清楚如何在 Excel 中执行此操作。

我正在附上一个带有数据和我正在工作的 Kaplan Meier 曲线的 excel 文件图像，如果有人可以帮助我编写内核平滑或贝塞尔平滑步骤并覆盖曲线并将文件重新发布给我，我将不胜感激.或者，更有效地为我提供执行此操作的 excel 插件或执行此操作的 SAS 代码。

Answer 1

如果提供生存数据的平滑描述的参数模型不合适，您可以使用平滑样条曲线或黄土曲线，例如在this page 上进行比较。为了平滑，将多个Bézier curves组合起来形成一种样条，黄土是kernel smoother的一种形式。这些很容易在标准统计软件包中实现；例如，在 SAS 中，TRANSREG 提供样条曲线，LOESS 实现了局部加权回归。*

也就是说，临床研究中的基线生存曲线通常以您所显示的阶跃函数类型呈现，没有进行平滑处理。此外，您的基础数据似乎只能以 1 年的间隔提供信息。如果是这样，那么呈现平滑曲线可能会在某种程度上歪曲数据的平滑度。如果您确实对每座桥梁都有单独的事件和审查时间，那么您的平滑曲线应该基于每个事件时间的生存概率，以提供最接近的数据表示。

---

*我会谨慎使用 Excel 进行任何严肃的统计工作。我发现很容易出现无法在 Excel 中检测到的非强制错误。