如何使用for循环运行1000次代码？

Question

我有这个代码，每次运行它时都会产生不同的“auc”。我想运行这个代码1000次，以计算1000个存储的AUC的平均值。我如何申请循环来做Iam寻找的东西？如果你们为我做并粘贴新代码，将不胜感激。正如我自4天以来一直试图解决的那样。

#iris is a built-in dataset 
main_df<- iris 
# extract data for "setosa" 
setosa <-main_df[main_df$Species=="setosa" ,] 
# extract data for "virginica" 
virginica<-main_df[main_df$Species=="virginica" ,] 
#merge "setosa" and "virginica" as new dataset 
df<- rbind(setosa,virginica) 

Cross.AUC<-rep(0,1000) # create a vector of zeros, here will be stored the auc values from each of 1000 runs 
for (i in seq(1:1000)) { 
#---------------devide data into two datasets 70:30 train:test ---- 
#---------------train dataset 
#select randomly 70% of setosa, generates a 35-by-5 matrix 
setosa_70<-setosa[sample(nrow(setosa),round(0.7*dim(setosa)[1])),] 
#select randomly 70% of virginica, generates a 35-by-5 matrix 
virginica_70<-virginica[sample(nrow(virginica),round(0.7*dim(virginica)[1])),] 
#merge setosa and virginica 
train<-rbind(setosa_70,virginica_70) 
#convert "setosa" to "0" and "virginica" to 1"" 
train$Species<-ifelse(train$Species=="setosa",0,1) 
#select 1st, 2nd and 5th columns 
train <-subset(train,select = c(1,2,5)) 
#--------------test dataset 
#select randomly 30% of setosa, generates a 15-by-5 matrix 
setosa_30<-setosa[sample(nrow(setosa),round(0.3*dim(setosa)[1])),] 
#select randomly 30% of virginica, generates a 15-by-5 matrix 
virginica_30<-virginica[sample(nrow(virginica),round(0.3*dim(virginica)[1])),] 
#merge setosa and virginica 
test<-rbind(setosa_30,virginica_30) 
#convert "setosa" to "0" and "virginica" to 1"" 
test$Species<-ifelse(test$Species=="setosa",0,1) 
#select 1st, 2nd and 5th columns 
test <-subset(test,select = c(1,2,5)) 
#merge "train" and "test" 
train_test<-rbind(train,test) 
#--Model_1-- 
model <-glm(Species~., family = binomial(link = "logit"),data = train_test) 
# install.packages("ROCR") 
library(ROCR) 
p <- predict(model, newdata=test, type="response") 
pr <- prediction(p, test$Species) 
auc <- performance(pr, measure = "auc") 
auc <- [email protected][[1]] 
AUC[i]<-auc 
} 

Answer 1

将AUC替换为Cross.AUC。您可以在范围内定义一个名为AUC的对象，您可以在for循环中访问该对象。

#iris is a built-in dataset 
main_df<- iris 
# extract data for "setosa" 
setosa <-main_df[main_df$Species=="setosa" ,] 
# extract data for "virginica" 
virginica<-main_df[main_df$Species=="virginica" ,] 
#merge "setosa" and "virginica" as new dataset 
df<- rbind(setosa,virginica) 

Cross.AUC<-rep(0,1000) # create a vector of zeros, here will be stored the auc values from each of 1000 runs 
for (i in seq(1:1000)) { 
    #---------------devide data into two datasets 70:30 train:test ---- 
    #---------------train dataset 
    #select randomly 70% of setosa, generates a 35-by-5 matrix 
    setosa_70<-setosa[sample(nrow(setosa),round(0.7*dim(setosa)[1])),] 
    #select randomly 70% of virginica, generates a 35-by-5 matrix 
    virginica_70<-virginica[sample(nrow(virginica),round(0.7*dim(virginica)[1])),] 
    #merge setosa and virginica 
    train<-rbind(setosa_70,virginica_70) 
    #convert "setosa" to "0" and "virginica" to 1"" 
    train$Species<-ifelse(train$Species=="setosa",0,1) 
    #select 1st, 2nd and 5th columns 
    train <-subset(train,select = c(1,2,5)) 
    #--------------test dataset 
    #select randomly 30% of setosa, generates a 15-by-5 matrix 
    setosa_30<-setosa[sample(nrow(setosa),round(0.3*dim(setosa)[1])),] 
    #select randomly 30% of virginica, generates a 15-by-5 matrix 
    virginica_30<-virginica[sample(nrow(virginica),round(0.3*dim(virginica)[1])),] 
    #merge setosa and virginica 
    test<-rbind(setosa_30,virginica_30) 
    #convert "setosa" to "0" and "virginica" to 1"" 
    test$Species<-ifelse(test$Species=="setosa",0,1) 
    #select 1st, 2nd and 5th columns 
    test <-subset(test,select = c(1,2,5)) 
    #merge "train" and "test" 
    train_test<-rbind(train,test) 
    #--Model_1-- 
    model <-glm(Species~., family = binomial(link = "logit"),data = train_test) 
    # install.packages("ROCR") 
    library(ROCR) 
    p <- predict(model, newdata=test, type="response") 
    pr <- prediction(p, test$Species) 
    auc <- performance(pr, measure = "auc") 
    auc <- [email protected][[1]] 
    Cross.AUC[i]<-auc 
} 

cat('Mean AUC:', mean(Cross.AUC), '\n')