Я пытаюсь понять, каковы остаточные значения «ld», возникающие в результате выполнения функции невязки на модели survreg? Например library(survival) mod <- survreg(Surv(time, status -1) ~ age , data = lung) residuals( mod , "ldcase")...

(Переходите к концу, чтобы узнать об обновлениях и новом коде и данных) Я пытаюсь построить кривую выживаемости с данными интервального цензора. Мне удалось успешно построить кривую выживаемости, используя данные интервального цензора (показанные...

У меня есть набор данных с сгруппированными наблюдениями в строке. Тем не менее, я хотел бы расширить наблюдение за каждой строкой с одного наблюдения на повторение до заданного числа (в данном случае «20» наблюдений каждое). На прикрепленном...

Я пытаюсь написать набор функций, в которых первая функция соответствует модели cox (через coxph в пакете survival в R ), а вторая функция получает расчетную выживаемость для нового набора данных с учетом подходящего объекта модели из первой...

Когда я запускаю cox.zph на модели Кокса, я получаю доход, отличный от того, что я вижу везде. Я попытался запустить следующий код: library(survival) #version 3.1-7 library(survminer) #version 0.4.6 res.cox <- coxph(Surv(time, status) ~...

Я хотел бы получить прогноз вероятности события из модели Коксма для заданных значений моих независимых переменных и на заданном временном шаге. Мои данные структурированы так: # Generate data set.seed(123) mydata <- data.frame(Site =...

Я пытаюсь проанализировать повторяющийся набор данных событий и изо всех сил пытаюсь соответствовать модели. Подмножество моих данных: outdat <- structure(list(yr = c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,...

Это может показаться глупым вопросом, но мне было интересно, почему медиана из median и медиана из survfit ( пакет выживания ) различаются Я попытался смоделировать учебник на sciencing.com : Составьте список времени выживания всех...

Я надеюсь, что кто-нибудь из вас сможет пролить свет на следующее. Я пытался воспроизвести модель Cox PH из Stata в R. Как вы можете видеть ниже, я получаю одинаковые результаты для моделей Cox PH без tvcs в обеих программах: Стата Кокс PH модель...

Я получаю прогнозы от модели clogit с использованием пакета выживания в R. Я использовал следующую статью о стеке, чтобы понять, как Как получить подходящие значения из модели clogit Однако в этом примере ответ включает формулу подогнанного...

Используя объект survfit, я могу извлечь квантили следующим образом: library('survival') fit <- survfit(Surv(time, status) ~ 1, data = aml) quantile(fit)[[1]] #25 50 75 #12 27 43 Есть ли функция для выполнения обратного? Передать вектор...