RMSTの算出 - 盆栽日記

生存時間分析でRMSTを出したい。 こちらではsurvRM2パッケージを使うと良いとある。 nshi.jp survRM2パッケージのvignetteはこちら。 https://cran.r-project.org/web/packages/survRM2/vignettes/survRM2-vignette3-2.html こんな感じでRMSTを算出してくれる。 > smp <- pbc %>% filter(status %in% c(0,2), !is.na(trt)) %>% + mutate(status = if_else(status == 2, 1, 0), + trt = trt-1 # 後述のrmst2()のamr引数が0,1しか受け付けないため変換 + ) # サンプルデータ >library(survRM2) > with(smp, rmst2(time, stat