當前位置: 首頁>>代碼示例 >>用法及示例精選 >>正文


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。非經特殊聲明,原始代碼版權歸原作者所有,本譯文未經允許或授權,請勿轉載或複製。