关于 R 代码 cph 错误的问题：找不到对象"n"

Question

我正试图用下面的R代码生成一个诺模图，并在最后一行不断出错。错误消息说"；cph中出错(Surv(futime，fustat(~。，x=T，y=T，surv=T，data=rt，：对象'n'未找到"；。有人能帮忙吗？提前谢谢。

library(rms)
setwd("C:\Users\yue\Desktop\04.diff")                  
rt=read.table("risk.txt",sep="t",header=T,row.names=1,check.names=F)   
rt=rt[c(1:(ncol(rt)-2))] 
dd <- datadist(rt)
options(datadist="dd")
f <- cph(Surv(futime, fustat) ~ ., x=T, y=T, surv=T, data=rt, time.inc=1)

Answer 1

我们可以将Surv对象创建为

s1 <- Surv(rt$futime, rt$fustat)
f <- cph(s1 ~ FOXM1 + IGF2BP3, x=TRUE, y=TRUE, surv=TRUE, data=rt, time.inc=1)

-输出

f
Cox Proportional Hazards Model

cph(formula = s1 ~ FOXM1 + IGF2BP3, data = rt, x = TRUE, y = TRUE, 
surv = TRUE, time.inc = 1)

Model Tests    Discrimination    
Indexes    
Obs       232    LR chi2     49.90    R2       0.231    
Events     43    d.f.            2    Dxy      0.550    
Center 1.0881    Pr(> chi2) 0.0000    g        1.019    
Score chi2  70.81    gr       2.771    
Pr(> chi2) 0.0000                      

Coef   S.E.   Wald Z Pr(>|Z|)
FOXM1   0.0812 0.0174 4.65   <0.0001 
IGF2BP3 0.1426 0.0303 4.70   <0.0001