0

我希望一切都好。很抱歉这个天真的问题，因为我对 R 很陌生。我已将数据从 STATA 移动到 R 并尝试做一个生存模型。我做了以下事情：

library(haven)
data <- read_dta("C:/Users/user/Desktop/data.dta")

View(data)

install.packages(c("survival", "survminer"))

library (caret)

library (glmnet)

library (mlbench)

library (psych)

library("survival")

library("survminer")

time (data$finaltime)

event=data$GSTATUS_DTHCNS_KI


Y=cbind ( time, event)


data$DIAB=as.factor(data$DIAB)

data$sex=as.factor(data$sex)

data$ETHCAT=as.factor(data$ETHCAT)

data$AMIS=as.factor(data$AMIS)

data$BMIS=as.factor(data$BMIS)

data$DRMIS=as.factor(data$DRMIS)

data$STEROIDS_MAINT=as.factor(data$STEROIDS_MAINT)

data$donortype=as.factor(data$donortype)

data$dgf=as.factor(data$dgf)

data$induction=as.factor(data$induction)

data$finalcmv=as.factor(data$finalcmv)

data$finalpra=as.factor(data$finalpra)

X=cbind (data$sex, data$BMI_TCR, data$COLD_ISCH_KI, data$SERUM_CREAT, data$finalpra, data$AGE, data$STEROIDS_MAINT, data$induction, data$DIAB, data$dgf, data$timeondialysis, data$ETHCAT, data$KDPI, data$finalcmv)

这非常有效。但是，当我尝试做 coxph 时，它没有用，我得到了以下信息：

 Coxph= coxph(Surv (time, event)~X, method “Breslow”)

错误：“Coxph= coxph(Surv (time, event)~X, method “”中的意外输入

另外，当我尝试时：

fit <- survfit(Surv(time, event) , data = data)

Surv（时间，事件）中的错误：时间变量不是数字

另外，当我尝试时：

Fit=glmnet (X,Y, family="cox")

错误：Cox 模型需要一个包含“时间”（>0）和“状态”（二进制）列的矩阵作为响应；一个“生存”对象就足够了

感谢您的帮助，为什么我会收到这些错误以及如何修复它们

非常感谢

期待您的回音

4

0 回答 0