R做列线图相关问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

各位大虾，想请教一个问题，我想做一个列线图
命令如下： （但是我想问一下h <- .02*exp(.024*(age-50)+.011*(GTV-60)+0.835*(LDH=='1'))
中的0.02,代表什么？其它的是B系数）
library(rms)
n <- 1000
set.seed(731)
age <-rnorm(n,50,15)
GTV<-rnorm(n,40,20)
label(age) <-"Age"
LDH<- factor(sample(c('0','1'), n, rep=TRUE, prob=c(.5, .5)))
cens <- 15*runif(n)
h <- .02*exp(.024*(age-50)+.011*(GTV-60)+0.835*(LDH=='1'))
dt <- -log(runif(n))/h
label(dt) <- 'Time'
e <- ifelse(dt <= cens,1,0)
dt <- pmin(dt, cens)
units(dt) <- "Months"
ddist <- datadist(age, GTV,LDH)
options(datadist='ddist')
S <- Surv(dt,e)
f <- cph(S ~ rcs(age,4)+rcs(GTV,4)+LDH, x=TRUE, y=TRUE)
f<-psm(Surv(dt,e)~ rcs(age,4)+rcs(GTV,4)+LDH,dist='lognormal')
surv<-Survival(f)
nom <- nomogram(f, fun=list(function(x) surv(36, x),
                            function(x) surv(60, x)),
            funlabel=c("3-year Survival Probability",
                       "5-year Survival Probability"))
plot(nom, xfrac=.7)
另外bootstrap怎么做？不同模型的C-index怎么求？

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：library factor

先顶再看好习惯