関数の動作方法を変更して理解を深めようとしています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
。
基本的に、私は2つの質問があります:
プロットは信頼区間を適切に出力しますが、どのタイプですか?
survfit()
それらは、オブジェクトの作成に使用した関数内で指定したexp_fit
ものと同じですか? そうでない場合、同じタイプの信頼区間が使用されるように変更する必要があるものは何ですか?凡例のフォント サイズは非常に小さいです。どうすれば大きくできますか?で作業しようとしました
theme()
が、残念ながらうまくいきませんでした。
あなたの答えに感謝します。前もって感謝します。