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

附件下载

所在主题：
文件名:  dat.rank.1.txt
资料下载链接地址: https://bbs.pinggu.org/a-3507097.html
附件大小:
289.12 KB   举报本内容
将数据随机分为建模组与验证组，为了获得比较好的建模效果，进行1000次随机分组，种子值分别为1:1000，循环1000次，在迭代到i=140报错，循环跳出，第140行及以后没有数值。特向朋友们请教：
1.我想忽略某一次迭代错误，不要跳出循环，继续下一次循环迭代，如何实现？
2.如何将大括号内的部分代码封装成函数，例如:C_index.test <- function(seed=i) {set.seed(i) 其余略}，最后用for (i in 139:141) C_index.test ()的形式调用该函数？我前面尝试过，没有成功。
多谢朋友们！
  1. library(rms)

  3. dat.rank.1 <- read.table('dat.rank.1.txt', sep='\t', header=T)
  4. units(dat.rank.1$OS) <- "Month"

  6. C_index.results1000 <- matrix(nrow = 145, ncol = 4)


  9. # 迭代到i=140报错，循环跳出，第140行及以后没有数值
  10. for (i in 139:141) {
  11. set.seed(i)
  12. ID_mod <- sample(1:nrow(dat.rank.1), nrow(dat.rank.1)*2/3)
  13. dat.model <- dat.rank.1[ID_mod, ]
  14. dat.model$mod <- 'Training'

  16. dat.validate <- dat.rank.1[-ID_mod, ]
  17. dat.validate$mod <- 'Validation'


  20. ddist <- datadist(dat.model)
  21. options(datadist = "ddist")


  24. f <- cph(Surv(OS, Event) ~ miR.index+T+M, x=T, y=T, surv=T, data=dat.model, time.inc=36)
  25. validate(f, method="boot", B=1000, dxy=T)
  26. rcorrcens(Surv(OS, Event) ~ predict(f), data = dat.model)
  27. C_index <- 1-rcorrcens(Surv(OS, Event) ~ predict(f), data = dat.model)[1]


  30. fev <- cph(Surv(OS, Event) ~ predict(f, newdata=dat.validate), x=T, y=T, surv=T, data=dat.validate, time.inc=36)
  31. validate(fev, method="boot", B=1000, dxy=T)
  32. rcorrcens(Surv(OS, Event) ~ predict(f, newdata=dat.validate), data = dat.validate)
  33. C_index2 <- 1-rcorrcens(Surv(OS, Event) ~ predict(f, newdata=dat.validate), data = dat.validate)[1]

  35. C_index.results1000[i, ] <- c(i, C_index, C_index2, C_index+C_index2)
  36. }

  38. View(C_index.results1000)# 第141行有数值说明后续迭代成功了
复制代码
代码及数据以网盘及附件两种方式提供，选择一种即可。
链接: https://pan.baidu.com/s/18mQ1takpdV24R28g-LyHsA 提取码: apca


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

扫码加我 拉你入群

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

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

京ICP备16021002号-2 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问：王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2025-12-26 02:44