关于 R 代码 cph 错误的问题:未找到对象 'n'”

a problem about error in cph of R code: object 'n' not found"

我正在尝试使用以下 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)

我们可以将 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