Как рассчитывается ограниченный средний верхний предел в анализе выживаемости в R - PullRequest
0 голосов
/ 17 декабря 2018

Я провел анализ Каплана Мейера в R, изучая выживаемость волокон в тесте на усталость.Я не предопределил верхний предел для среднего значения.Как R вычисляет или определяет верхний предел для вычисления ограниченного среднего?Я использую следующий код:

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

1 Ответ

0 голосов
/ 26 июля 2019

Исходя из опыта этих расчетов за последние несколько лет, я считаю, что ограниченное среднее значение по умолчанию рассчитывается как среднее для самых длинных жизней в каждой группе.

Например, я столкнулся с этим, когда две группы живут следующим образом:

Жизнь в группе 1:

22 23 25 26 30 32 32 34 3738 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 123 125 138 150 163 163 164 167 199 201 214 235 236 237 242 245270 272 274 282 283 284 287 296 300 300 310 314 321 322 325 340 342 345 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 885886 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

...