Если у меня есть объект survdiff в R, есть ли быстрый способ определить, какая группа имеет наилучшее выживание (если есть значительная разница)?
По сути, я перебираю строки фрейма данных и каждый раз получаю значение survdiff p, но на самом деле меня интересует значение p только в том случае, если конкретная группа показывает наилучшее выживание.
Я нашел ответ.
Для объекта survdiff y = (формула = Surv(время, событие) ~ tclass), y$obs/y$exp, наблюдаемые/ожидаемые события, будут самыми низкими для группы с лучшим выживанием, самыми высокими для группы с лучшим выживанием. лучшее выживание и так далее. Эти значения упорядочены по именам, хранящимся в именах (y $ n), поэтому вы можете проверить, как называется группа с лучшим выживанием:
имена(y$n)[y$obs/y$exp==min(y$obs/y$exp)
