R语言 glmnet包用于cox遇到的问题

4关注
0粉丝

初中生

95%

还不是VIP/贵宾

-

0%

威望: 0 级
论坛币: 1 个
通用积分: 0
学术水平: 0 点
热心指数: 0 点
信用等级: 0 点
经验: 43 点
帖子: 2
精华: 0
在线时间: 39 小时
注册时间: 2017-10-31
最后登录: 2022-4-24

楼主

amor_liu

发表于 2018-4-23 23:40:51 |AI写论文

是否 +2 论坛币

k人参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

variable.xlsx (27.71 KB)
library(readxl)

lass<- read_excel("E:/R/prostate/2010+/variable.xlsx")
lass<-cbind(time,status,variable)
time<-lass$survival.months
status<-lass$OSD
variable<-lass[,3:14]
#
variable$subgroup<-relevel(factor(variable$subgroup),ref = "nscNECa")
variable$age<-relevel(factor(variable$age),ref = "70-")
variable$N<-relevel(factor(variable$N),ref = "N0")
variable$M<-relevel(factor(variable$M),ref = "M0")
variable$AJCCstage<-relevel(factor(variable$AJCCstage),ref = "I-III")
variable$csstage<-relevel(factor(variable$csstage),ref = "Localized")
variable$surg.LN<-relevel(factor(variable$surg.LN),ref = "Removed")
variable$radiation<-relevel(factor(variable$radiation),ref = "Yes")
variable$chemotherapy<-relevel(factor(variable$chemotherapy),ref = "Yes")
variable$brain<-relevel(factor(variable$brain),ref = "No")
variable$liver<-relevel(factor(variable$liver),ref = "No")
variable$lung<-relevel(factor(variable$lung),ref = "No")
#
as.matrix(variable)
time<-as.double(time)
status<-as.double(status)
surviv<-Surv(time,status)
#fix(variable)
fit <- glmnet(variable, surviv, family = "cox", alpha = 1)
> fit <- glmnet(variable, surviv, family = "cox", alpha = 1)Error in coxnet(x, is.sparse, ix, jx, y, weights, offset, alpha, nobs, : negative event times encountered; not permitted for Cox family