Я пытаюсь изменить и лучше понять, как работает функция ggsurvplot()
(пакет: survminer
):
Прежде всего, некоторый код:
suppressPackageStartupMessages(library(survival))
suppressPackageStartupMessages(library(survminer))
suppressPackageStartupMessages(library(tidyverse))
set.seed(2020)
exp <- sample_n(lung, size = 10, replace = F) %>%
select(time, status)
exp <- sample_n(lung, size = 10, replace = F) %>%
select(time, status)
exp$survival <- Surv(exp$time, exp$status == 2)
exp_fit <- survfit(survival ~ 1, data = exp, conf.type = "log-log")
ggsurvplot(exp_fit,
data = exp,
palette = "black",
censor.shape = "I",
censor.size = 5.5,
xlab = "days",
ylab = "survival probability",
conf.int = T,
surv.scale = "percent",
legend.title = "censored data",
legend.lab = "",
legend ="bottom")
Фрейм данных используется называется lung
. Данные содержатся в пакете survival
.
В основном у меня два вопроса:
На графике печатаются доверительные интервалы, как и должно быть, НО какой тип? Являются ли они такими же, как я указал в функции survfit()
, которую я использовал для создания объекта exp_fit
, и если нет, то что нужно изменить, чтобы использовать тот же тип доверительных интервалов?
Размер шрифта легенды довольно маленький. Как я могу сделать это больше? Я пытался работать с theme()
, но, к сожалению, это не сработало.
Буду признателен за ваши ответы. Заранее спасибо.