library(openxlsx)#导入打开文件的包

project1 <- read.xlsx("AAA")#导入文件

shapt1 <- c()#数据框
for(i in c(12:18)){#循环
  print(i)
  fit <- lm(project1[[i]]~project1[[3]],data = project1)#线性回归检验
  sumres <- summary(fit)
  P <- pf(sumres$fstatistic[1L], sumres$fstatistic[2L], sumres$fstatistic[3L], lower.tail = FALSE)#提取P值
  Met1 <- fit$method#检验方法
  result1<-c()#设置内部表格
  result1<-data.frame(rbind(result1,cbind(P,Met1)))
  result1$variate.name1<-names(project1)[3]#加载新列
  result1$variate.name <- names(project1)[i]#加载新列
  shapt1<-rbind(shapt1,result1)
}