使用R中的幸存者包(ggsurvplot)更改截尾事件的形状

时间:2017-04-06 08:39:41

标签: r ggplot2 survival-analysis

我正在使用survminer包在R中创建一个Kaplan-Meier-Curve,它运行得很好:

plotUICC=survfit(Surv(as.numeric(myData$eventtime),event=myData$eventtype)~myData$UICC, data = myData)

ggsurvplot(plotUICC,data=myData, risk.table = FALSE,pval = TRUE,conf.int = FALSE,legend.labs = c("UICC I", "UICC II"),legend = c(0.25, 0.26),legend.title = "",xlab = "Time in Months"))

Kaplan Meier Curve created with the Cod above

我想将“+”更改为审查事件并使用简单的“|”代替。不幸的是,没有任何帮助(ggsurvplot)和谷歌都可以帮助我。

这些是装载的包装:

    other attached packages:
[1] ggthemes_3.4.0  survminer_0.3.1 ggpubr_0.1.2    ggplot2_2.2.1   survival_2.40-1
[6] readxl_0.1.1    gridExtra_2.2.1

1 个答案:

答案 0 :(得分:2)

最小的工作示例如下:

library(survminer)
library(survival)
fit <- survfit(Surv(time, status) ~ sex, data = lung)
p1 <- ggsurvplot(fit, data = lung)

理想情况下,人们希望改变审查者的形状,例如:像这样:

p2 <- ggsurvplot(fit, data = lung, censor.shape=124)

我发现原始的幸存者包无法做到这一点。在fork幸存者中,我实现了这种行为,以便生成以下图:

enter image description here

censor.shape的可接受值可在此处找到: http://sape.inf.usi.ch/quick-reference/ggplot2/shape

我还没有对此进行过彻底的测试,但这个扩展对我来说已经非常有用了。 &#34; +&#34;审查人员困扰了我一段时间......

编辑:此建议现已合并,ggsurvplot的development version现在接受censor.shape和censor.size参数:

ggsurvplot(fit, data = lung, censor.shape="|", censor.size = 4)