如何将这个流程用for循环跑1000次,每一次都将样品随机7/3分,然后得到1000个模型,每个模型含有不同的基因集合(our.rf$topvars)这一句代码可以得到每个模型需要的基因 代码如下 不会写for循环 请大佬们指教!谢谢!
A <- read.csv('1035.csv',row.names = 1,header = T)
i <- sample(2, nrow(A), replace = TRUE, prob = c(0.7, 0.3))
train <- A[i==1, ]
test <- A[i==2, ]
disease.rf <- rfsrc(Surv(OS_MONTHS,OS_STATUS)~.,data = train,ntree = 800,mtry = 3,nodesize = 3,splitrule = "logrank")
our.rf <- var.select(object = disease.rf,vdv,method = "vh.vimp",nrep = 50)
our.rf$topvars