当前位置: 首页>>代码示例 >>用法及示例精选 >>正文


R print.survfit 打印生存曲线的简短摘要


R语言 print.survfit 位于 survival 包(package)。

说明

打印观察数、事件数、限制平均生存期及其标准误差,以及带有中位数置信限的中位生存期。

用法

## S3 method for class 'survfit'
print(x, scale=1, digits = max(options()$digits - 4,3),
    print.rmean=getOption("survfit.print.rmean"),
    rmean = getOption('survfit.rmean'),...)

参数

x

调用 survfit 函数的结果。

scale

用于重新调整生存时间的数值,例如,如果 survfit 的输入数据以天为单位,scale=365 会将打印输出缩放为年。

digits

要打印的位数

print.rmean , rmean

用于计算和显示限制平均值的选项。

...

为了未来的结果

细节

均值及其方差基于截断估计量。也就是说,如果最后一次观察不是死亡,则生存曲线估计不会变为零且平均值未定义。有四种可能的方法可以解决此问题,由 rmean 选项选择。第一个是将上限设置为常量,例如 rmean=365 。在这种情况下,报告的平均值将是前 365 天中每组将经历的预期天数。如果兴趣集中在固定时期,这很有用。其他选项是 "none" (无估计)、"common""individual""common" 选项使用对象中所有曲线的最大时间作为 auc 计算的通用上限。对于 "individual" 选项,平均值计算为每条曲线下的面积,范围从 0 到该曲线的最大观察时间。由于终点是随机的,因此不同曲线的值不可比较,并且打印的标准误差被低估,因为它们没有考虑这种随机变化。提供此选项主要是为了向后兼容,因为此估计是早期版本的代码中的默认(唯一)估计。请注意,SAS(自版本 9.3 起)使用直到每条曲线的最后一个事件时间的积分;我们认为这是最糟糕的选择,并且不提供该计算的选项。

中位数及其置信区间是通过在生存曲线及其置信带图上在 0.5 处画一条水平线来定义的。如果该线不与曲线相交,则中位数未定义。该线与 CI 下带的交点定义了中位数区间的下限,上带也类似。如果任何交点都不是点,则我们使用交点间隔的中心,例如,如果生存曲线在某个间隔内恰好等于 0.5。当数据未经审查时,这与中位数的通常定义一致。

x,设置了不可见标志以防止打印。 (R 中所有打印函数的默认设置是返回传递给它们的对象;print.survfit 符合此模式。如果要捕获这些打印结果以进行进一步处理,请参阅 summary.survfittable 组件。)

副作用

打印观察数、事件数、中位生存期及其置信区间,以及可选的限制平均生存期 (rmean) 及其标准误差。如果有多条曲线,则每条曲线都有一行输出。

参考

Miller, Rupert G., Jr. (1981). Survival Analysis. New York:Wiley, p 71.

也可以看看

summary.survfit , quantile.survfit

相关用法


注:本文由纯净天空筛选整理自R-devel大神的英文原创作品 Print a Short Summary of a Survival Curve。非经特殊声明,原始代码版权归原作者所有,本译文未经允许或授权,请勿转载或复制。