R语言批量读取文件提取回归系数

尔云间  一个专门做科研的团队

小云在前期工作得到了100个如下文件

接下来我们要对BLUE和Prediction进行建模，并提取Prediction的回归系数

如果靠我们单独一个个算是很难短时间完成的，所以小云决定用R语言来做，只需要写一个简单的循环就可以解决

代码如下：

t<-matrix(0,100,1)#建立空矩阵

file <- list.files(pattern=".sPred.result.csv")#建立文件名列表

t=100

for(i in 1:t){

 

  pred = fread(file[i])

  mod_M6 = lm(pred$BLUE ~ pred$Prediction)#建立线性回归模型

  P<-summary(mod_M6)

  c<-P$coefficients

  c<-as.matrix(c)

  t[i]<-c[2,1]

}#End of for (rep in 1:t)

write.csv(t,"regsBLUP.csv")

 

x列就是我们需要的回归系数

生信果

生信硬核知识解答

和小果一起学生信