如何在R中计算生存分析中的限制平均上限
我在R中进行了Kaplan Meier分析，观察疲劳测试中纤维的存活率。我还没有预先定义限制平均值的上限。为了计算受限平均值，R如何计算或决定上限？我使用了以下代码：
fit = survfit(Surv(cyclicdata[,1], cyclicdata[,2]) ~ cyclicdata[,3])
print(fit, print.rmean=TRUE,rmean="common")回答 1
Stack Overflow用户
发布于 2019-07-27 01:18:22
根据过去数年的经验，我相信受限制的平均数是按每组最长寿命的平均数计算。
例如，我遇到过这样的两个群体，他们的生活如下：
第一组生活：
22 23 25 26 30 32 32 34 37 38 40 43 45 48 54 56 59 60 62 70 73 73 76 77 78 78 82 86 92 92 95 98 99 102 104 106 107 112 114 115 119 120 123 132 134 151 154 157 169 180
第二组生活：
5 7 30 41 44 56 59 64 67 79 86 101 110 120 123 125 138 150 163 164 167 199 201 214 235 236 237 242 245 270 272 274 283 283 284 287 296 300 300 310 314 321 322 325 340 342 355 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 817 885 886 946 947 951
当我运行我的生存拟合并打印输出时，我得到以下结果：
* restricted mean with upper limit = 566
这等于：
> mean(max(c(Group1$Lives,Group2$Lives)))
[1] 565.5
或
> (951+180)/2
[1] 565.5
https://stackoverflow.com/questions/53813167
复制相似问题
腾讯云开发者
