1

R Studio では、インターセプト ラインを追加すると、生存プロットからリスク テーブルが削除されます。この問題を克服する方法はありますか? リスクテーブルを示す私の最初のコードは次のとおりです。

    library(readxl)
    library(ggplot2)
    library(survival)
    library(survminer)
    library(gmodels) 
    library(tidyverse)

    demo2<-read_excel('data.xlsx',sheet=1)

    #tte1=time to event 1
    tte1<-read_excel("data.xlsx", sheet=1)
    tte1$group<-factor(tte1$group,levels=c("C","I")) #turn group into a factor, and sets order of         
    levels

    #Kaplan-Meier plot for event
    fit1 <- survfit(Surv(Time, Event) ~ group,data = tte1)
    summary(fit1)
    p<-ggsurvplot(fit1, data = tte1, 
                  risk.table = TRUE, #adds risk table
                  risk.table.title = "Number at risk",
                  risk.table.subtitle = " ",
                  cumcensor.title= " ",
                  conf.int = TRUE,   #adds shading
                  surv.scale="percent",
                  censor= FALSE,
                  conf.int.style= 'ribbon',
                  legend.title= '',
                  legend.labs = c("Control", "Intervention"),
                  break.time.by=60,  #adds intervals
                  tables.theme = theme_cleantable(), # theme for tables
                  tables.y.text = FALSE,
                  xlim=c(0,290),
                  surv.plot.height= 0.75,
                  risk.table.height=0.1
    )
    p+xlab("Time (days)")

しかし、このコードの最終行を変更してインターセプト行を追加すると、リスク テーブルが突然消えてしまいます。インセプト行を追加するために、コードの最後の行を から に変更しp+xlab("Time (days)")ました

    p$plot+geom_vline(xintercept = 186, linetype="dashed",
                      color = "black", size=1.0, label= "median")+ 
      geom_text(aes(x=225, y=.20, label="Median follow-up time"), size=5, family="sans", 
    fontface="plain")

リスク テーブルを失わずにインターセプト ラインを追加する方法を知っている人はいますか?

ありがとう！

4

1 に答える 1