4

我正在使用 glmnet 运行 coxph 模型。以下是我正在使用的代码:

data <- read.csv("CSM Tiers by Org.csv", stringsAsFactors = FALSE, na.strings = "?")
x <- model.matrix(~ data$CSM + data$Exclude +  data$Queue)
y <- data$Time
status <- data$Event
glmnet(x,Surv(y,status), family = "cox")

但是,我收到以下错误:

Error in coxnet(x, is.sparse, ix, jx, y, weights, offset, alpha, nobs,  : 
negative event times encountered;  not permitted for Cox family

但我没有在我的数据中看到任何负面事件时间。

> sum(y<0)
[1] 0

任何帮助是极大的赞赏!

4

1 回答 1

5

检查总和(y <= 0)。函数中的代码表示,如果 y 也等于 0,它将返回该错误。

代码在这里,您引用的错误在第 5 行: https ://github.com/jeffwong/glmnet/blob/master/R/coxnet.R

于 2015-01-16T18:50:30.757 回答