回归系数(回归系数的含义)没想到

尔云间 一个专门做科研的团队原创 小果 生信果 关注我们小果在前期工作得到了100个如下文件接下来我们要对BLUE和Prediction进行建模

尔云间 一个专门做科研的团队原创 小果 生信果 关注我们

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

推荐阅读多种花样都能懂，富集图的看法相关性图谱之相关性热图学习小果教你快速优雅使用Genecards数据库从UCSC下载TCGA数据简单高效利用Batch Entrez批量获取基因别名关注小果，小果将会持续为你带来更多生信干货哦。

“生信果”，生信入门、R语言、生信图解读与绘制、软件操作、代码复现、生信硬核知识技能、服务器、生物信息学的教程，以及基于R的分析和可视化等原创内容，一起见证小白和大佬的成长。