我在R中进行了Kaplan Meier分析，研究了疲劳测试中纤维的存活率。我尚未为限制均值预定义上限。 R如何计算或确定上限以计算限制均值？ 我正在使用以下代码：
fit = survfit(Surv(cyclicdata[,1], cyclicdata[,2]) ~ cyclicdata[,3])
print(fit, print.rmean=TRUE,rmean="common")
答案 0 :(得分：1)
根据过去几年的计算经验，我认为默认情况下，受限均值是根据每个组中最长寿命的均值来计算的。
例如，我遇到了两个生活如下的群体：
第1组的生活：
22 23 25 26 30 32 32 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 132134135151154157169180
第2组的生活：
5 7 30 41 44 56 59 64 67 79 86 101 110 120 120 123 125138150163163164167167199201214235236237237242245270 272274282283284287296300300310314321322325340342 345 355 371 375 376 398 414 419 422 428 442 444 449 474 511 516 549 552560563581608608618628637638675685702702782782817885 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