PASUA1<-svydesign(id=~SDMVPSU,strata =~SDMVSTRA, weights =~WTMEC18YR, nest = TRUE, data = PASUA)#设置survey参数
PASUA2<-as.svrepdesign(PASUA1)
fit<-svycoxph(Surv(PERMTH_EXM,ALLMORT) ~ rcs(LBXGLU,4)+SEXCAT+RACECAT+EDUCAT+SMKCAT+ALCCAT+ACTCAT,design=PASUA1)
Pre0 <-rms::Predict(fit,LBXGLU,type="predictions",ref.zero=F,conf.int = 0.95,digits=2);
Error in rms::Predict(fit, LBXGLU, type = "predictions", ref.zero = F, :
predictors(s) not in model: LBXGLU