搜索
人大经济论坛 附件下载

附件下载

所在主题:
文件名:  代码.txt
资料下载链接地址: https://bbs.pinggu.org/a-3731565.html
附件大小:
455 Bytes   举报本内容
各位大佬,想请教一个问题,这个问题困扰我3天了,问题如下library(survival)
library(rms)
dt<-cancer
dd<-datadist(dt)
options(datadist='dd')
fit <- cph(Surv(dt$time,dt$status==1) ~ rcs(age,4), x=TRUE, y=TRUE,data=dt)
pr1<-rms::Predict(fit,age=dt$age,fun=exp,type="predictions",ref.zero=T,conf.int = 0.95,digits=2)

这组代码在外部运行得很好,求出了预测数据pr1,我想写个function在已知数据和fit的情况下求pr1



f<-function(data,fit,x){
fit <- fit
x<-dt[,x]
dt<-data
pr1 <-rms::Predict(fit,x=dt$x,fun=exp,type="predictions",ref.zero=T,conf.int = 0.95,digits=2)
pr1
}

f(dt,fit,"age")



我调式了一下,主要错在x=dt$x,这一句,在外部运行的是时候是age=dt$age,而在function中运行变成了"age"=dt$age多了一个双引号,导致失败。我想请教一下各位大佬有什么办法可以实现我function中的代码,得出pr1吗?先谢谢了,论坛编辑代码太难了,请看附件附有代码。






    熟悉论坛请点击新手指南
下载说明
1、论坛支持迅雷和网际快车等p2p多线程软件下载,请在上面选择下载通道单击右健下载即可。
2、论坛会定期自动批量更新下载地址,所以请不要浪费时间盗链论坛资源,盗链地址会很快失效。
3、本站为非盈利性质的学术交流网站,鼓励和保护原创作品,拒绝未经版权人许可的上传行为。本站如接到版权人发出的合格侵权通知,将积极的采取必要措施;同时,本站也将在技术手段和能力范围内,履行版权保护的注意义务。
(如有侵权,欢迎举报)
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

GMT+8, 2026-2-8 11:56