Таблица потери риска при анализе выживания при обрезании хвоста - PullRequest
Новая
       9

Таблица потери риска при анализе выживания при обрезании хвоста

0 голосов
/ 11 января 2019

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

# --functioning code without tail removed
library(survival)
dataset <- read.csv("C:\\Users\\Michael McCarthy\\Desktop\\testing\\new_format.csv", 
                    header = TRUE)
mysurv <- Surv(time=dataset$length_of_time_active_days, event=dataset$censored)
myfit <- survfit(mysurv~dataset$hardware_type)
ggsurvplot(myfit, size=1, data=dataset, 
           font.main = 10, font.x =  10, font.y = 10, font.tickslab = 10,
           legend = "bottom", legend.title = "Control Units",
           conf.int = FALSE, ggtheme = theme_bw(),
           risk.table = TRUE, risk.table.y.text.col = TRUE, tables.height=.35, 
           surv.plot.width=1, fun = function(y) y*100,
           ylab='Survival %', xlab= 'Months with Customers',
           risk.table.y.text = FALSE)

# -- code where tail is removed
library(survival)
dataset <- read.csv("C:\\Users\\Michael McCarthy\\Desktop\\For deleting\\new_format.csv",
                    header = TRUE)
mysurv <- Surv(time=dataset$length_of_time_active_days, event=dataset$censored)
myfit <- survfit(mysurv~dataset$hardware_type)
aa <- subset(surv_summary(myfit,data = dataset), n.risk>=1000)
ggsurvplot(aa, size=1, data=dataset, 
           font.main = 10, font.x = 10, font.y = 10, font.tickslab = 10,
           legend = "bottom", legend.title = "Control Units",
           conf.int = FALSE, ggtheme = theme_bw(),
           risk.table = TRUE, risk.table.y.text.col = TRUE, tables.height=.35,
           surv.plot.width=1, fun = function(y) y*100,
           ylab='Survival %', xlab= 'Months with Customers',
           risk.table.y.text = FALSE)
...