R построение нескольких кривых выживаемости в одном и том же участке

Question

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

plot(sr_fit_0, col = 'red' , conf.int=TRUE, xlim=c(0, max_m)) 
par(new=TRUE) 
plot(sr_fit_1, col ='blue', conf.int=TRUE, xlim=c(0, max_m))` 

Но теперь я хочу использовать ggsurv построить кривую выживания, и я не знаю, как у них обоих в том же участке (не сюжетных). Любая помощь приветствуется.

Answer 1

Я создал некоторые данные для жизни ниже для жизни хомяков и песчанок. Вы можете использовать функцию survfit(), аналогичную функции других функций подгонки, и определить столбец кадра данных, который разбивает популяцию. Когда вы создаете сюжет с ggsurv(), я думаю, что он отобразит то, что вы ищете.

## Make some data for varmint life 
set.seed(1); l1 <- rnorm(120, 2.5, 1) 
gerbils <- data.frame(life = l1[l1>0]) 
set.seed(3); l2 <- rnorm(120, 3, 1) 
hamsters <- data.frame(life = l2[l2>0]) 

## Load required packages 
require('survival'); require('GGally') 

## Generate fits for survival curves 
## (Note that Surv(x) creates a Survival Object) 
sf.gerbils <- survfit(Surv(life) ~ 1, data = gerbils) 
sf.hamsters <- survfit(Surv(life) ~ 1, data = hamsters) 
ggsurv(sf.gerbils) #Survival plot for gerbils 
ggsurv(sf.hamsters) #Survival plot for hamsters 

## Combine gerbils and hamsters while adding column for identification 
varmints <- rbind((cbind(gerbils, type = 'gerbil')), 
        (cbind(hamsters, type = 'hamster'))) 

## Generate survival for fit for all varmints as a function of type 
sf.varmints <- survfit(Surv(life) ~ type, data = varmints) 

## Plot the survival curves on one chart 
ggsurv(sf.varmints)