[英]implementing a non-linear relationship within a cox proportional hazards model using the survival package in R
[英]In Survival Analysis with R, what is the purpose of the `surv`function in the Cox Proportional Hazards Model?
我目前正在查看一份说明要使用Cox比例危害模型的文档,该模型是您的公式部分的响应变量
coxph(formula, data=, weights, subset,
na.action, init, control,
ties=c("efron","breslow","exact"),
singular.ok=TRUE, robust=FALSE,
model=FALSE, x=FALSE, y=TRUE, tt, method, ...)
必须在公式部分中surv()。
有人可以告诉我surv()函数的作用吗? 我了解它说这是一个生存对象,但是我不确定这是否是必需的。 谢谢!
在这种情况下,您只需要阅读文档并运行其中的示例即可。 第一个例子是? coxph
? coxph
显示以下内容:
# Create the simplest test data set
test1 <- list(time=c(4,3,1,1,2,2,3),
status=c(1,1,1,0,1,1,0),
x=c(0,2,1,1,1,0,0),
sex=c(0,0,0,0,1,1,1))
# Fit a stratified model
coxph(Surv(time, status) ~ x + strata(sex), test1)
显然,您需要使公式的左侧/响应部分成为Surv
的输出( Surv
的输出也具有清晰的文档;请参见?Surv
)。 如果您查看该对象:
> str(Surv(test1$time,test1$status))
Surv [1:7, 1:2] 4 3 1 1+ 2 2 3+
- attr(*, "dimnames")=List of 2
..$ : NULL
..$ : chr [1:2] "time" "status"
- attr(*, "type")= chr "right"
并查看它如何反映time
和status
列中包含的信息:
> with(test1, cbind.data.frame(time, status, Surv(time,status)))
time status Surv(time, status)
1 4 1 4
2 3 1 3
3 1 1 1
4 1 0 1+
5 2 1 2
6 2 1 2
7 3 0 3+
然后,要回答有关是否必要的问题,您可以尝试在没有它的情况下运行coxph
并查看会发生什么:
> coxph(time ~ x + strata(sex), test1)
Error in coxph(time ~ x + strata(sex), test1) :
Response must be a survival object
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.