r - R中如何计算生存分析中的限制平均上限

Question

我在 R 中进行了 Kaplan Meier 分析，研究了疲劳测试中纤维的存活率。我没有预定义受限均值的上限。R如何计算或决定上限以计算受限均值？我正在使用以下代码：

fit = survfit(Surv(cyclicdata[,1], cyclicdata[,2]) ~ cyclicdata[,3])
print(fit, print.rmean=TRUE,rmean="common")

标签： rsurvival-analysis

Answer 1

根据过去几年的计算经验，我认为受限平均值默认为每个组最长寿命的平均值。

例如，我遇到了两个群体，他们的生活如下：

第 1 组生活：

22 23 25 26 30 32 32 34 37 38 40 43 45 48 48 54 56 59 60 62 70 72 73 73 76 77 78 78 82 86 86 92 92 92 95 98 99 102 104 106 107 112 114 114 115 119 120 123 132 134 135 151 154 157 169 180

第 2 组生活：

5 7 30 41 44 56 59 64 67 79 86 101 110 120 120 120 120 120 123 123 125 163 163 163 163 163 164 164 167 199 2012 214 235 235 236 237 237 242 242 245 272 272 274 274 274 274 282 284 284 284 284 287 296 300 371 375 376 398 414 419 422 428 442 444 449 474 511 516 549 552 560 563 581 608 618 628 637 638 675 685 702 782 782 8417 89 95 958

当我运行我的生存适应并打印输出时，我得到了这个：

* restricted mean with upper limit = 566

这等于：

> mean(max(c(Group1$Lives,Group2$Lives)))

[1] 565.5

或者

> (951+180)/2

[1] 565.5