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

原创 小果 生信果

小云在前期工作得到了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列就是我们需要的回归系数

好了，小伙伴们，今天的主要内容就是这些了，小伙伴们有什么需要分享和讨论的都可以来找小云啊，小云非常高兴和大家交流。