爾云間? 一個(gè)專門做科研的團(tuán)隊(duì)

小云在前期工作得到了100個(gè)如下文件

接下來我們要對(duì)BLUE和Prediction進(jìn)行建模，并提取Prediction的回歸系數(shù)

如果靠我們單獨(dú)一個(gè)個(gè)算是很難短時(shí)間完成的，所以小云決定用R語言來做，只需要寫一個(gè)簡單的循環(huán)就可以解決

代碼如下：

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列就是我們需要的回歸系數(shù)

生信果

生信硬核知識(shí)解答

和小果一起學(xué)生信