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)
}
用R语言写循环函数,来对数据进行线性回归分析,最终提取p值,并将数据放入数据框中
最新推荐文章于 2024-05-13 16:47:16 发布