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)
}