library(ggrcs)
library(rms)
library(ggplot2)
library(scales)
library(cowplot)
data<-na.omit(lung)
str(data)
dd<-datadist(data)
options(datadist='dd')
fit<- cph(Surv(time,status==1) ~ rcs(wt.loss,4)+sex, x=TRUE, y=TRUE,data=data)
ggrcs(data=data,fit=fit,x="wt.loss")
R语言实现限制性立方样条拟合(1)
最新推荐文章于 2025-02-28 13:23:29 发布