如何优化survplot生存曲线的可视化效果?
低死亡率生存数据可视化优化求助
我正在开展一项为期10天的观测实验,对比NTS、NTC、TS、TC四种不同条件。已完成logrank检验,组间差异具有统计学显著性,但实验死亡率极低,当前用ggsurvplot绘制的生存曲线组间重叠严重,希望优化可视化方法,实现组间曲线清晰分离的效果。
现有R脚本
morta <- read.csv("mortality.csv", head=TRUE, sep = ";") surv_obj <- Surv(morta$Time, morta$Statut) survdiff(formula = surv_obj ~ Condition, data = morta)
Logrank检验结果
N Observed Expected (O-E)^2/E (O-E)^2/V Condition=NTC 1200 1069 1087 0.2821 26.46 Condition=NTS 1200 1156 1145 0.1120 10.69 Condition=TC 1200 1126 1130 0.0134 1.27 Condition=TS 1200 1137 1127 0.0900 8.54 Chisq= 35.4 on 3 degrees of freedom, p= 1e-07
数据集信息
- ID:个体标识(无关)
- Statut:0代表死亡,1代表存活
- Time:死亡天数,实验结束存活则记为10
- Condition:分组信息(NTS、NTC、TS、TC)
- Colony:1-6的群体标签
内容的提问来源于stack exchange,提问作者lddc
相关产品推荐
相关产品推荐

