筰夢祌夢鑒身外身
2021-11-01 16:48
采纳率: 0%
浏览 332

r语言绘图 限制性立方样条 COX回归

想问一下立方样条曲线的置信区间的问题,我应该怎么调整置信区间的图像成喇叭形,代码奉上,第一张是自己画的图,第二张是文献的图。

library(ggplot2)
library(survival)
library(rms)
setwd("D:/work/work/zhang 9-28/2021.10.1")
data=read.csv("11-全部数据汇总 加入其它预测模型(2).csv")  
dd <- datadist(data) 
head(data)
options(datadist='dd') 
fit<- cph(Surv(recurrence_time,recurrence) ~ rcs(DURATION,4) + CAAPAF,data=data)
cox.zph(fit, "rank")  
ggcoxzph(cox.zph(fit, "rank"))
anova(fit)  
view(dd)
dd$limits$DURATION[1] <- 50
fit=update(fit)
HR<-Predict(fit,DURATION,fun=exp,ref.zero = TRUE) 
P1<-ggplot(HR)
P1
P2<-ggplot()+geom_line(data=HR, aes(DURATION,yhat),linetype="solid",size=1,alpha = 0.7,colour="red")+
 geom_ribbon(data=HR, aes(DURATION,ymin = lower, ymax = upper),alpha = 0.1,fill="red")
P2<-P2+theme_classic()+geom_hline(yintercept=1, linetype=2,size=1)+ labs(title = "RCS", x="EeSD", y="HR (95%CI)") 
P2

另:为啥这句代码更改数值不好用呢?dd$limits$DURATION[1] <- 50

img

img

1条回答 默认 最新

相关推荐 更多相似问题