| 所在主题: | |
| 文件名: 代码.txt | |
| 资料下载链接地址: https://bbs.pinggu.org/a-3731565.html | |
| 附件大小: | |
|
各位大佬,想请教一个问题,这个问题困扰我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、本站为非盈利性质的学术交流网站,鼓励和保护原创作品,拒绝未经版权人许可的上传行为。本站如接到版权人发出的合格侵权通知,将积极的采取必要措施;同时,本站也将在技术手段和能力范围内,履行版权保护的注意义务。 (如有侵权,欢迎举报) |
|
京ICP备16021002号-2 京B2-20170662号
京公网安备 11010802022788号
论坛法律顾问:王进律师
知识产权保护声明
免责及隐私声明