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

Question

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

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

Answer 1

你好，我是有问必答小助手，非常抱歉，本次您提出的有问必答问题，技术专家团超时未为您做出解答


本次提问扣除的有问必答次数，将会以问答VIP体验卡（1次有问必答机会、商城购买实体图书享受95折优惠）的形式为您补发到账户。


因为有问必答VIP体验卡有效期仅有1天，您在需要使用的时候【私信】联系我，我会为您补发。