mvp结果分析

library(dplyr)
library(‘reshape‘)#子集的包
tmp<-data.frame() #创建空的数据集
myfiles <- Sys.glob("*GLM.*csv") #搜索文件夹内所有的gwas结果文件
print(paste(myfiles,"has been found")) #显示所有找到的结果文件
for(i in myfiles) #循环开始
{new_trait<-read.csv(file=i,header = TRUE,sep = ",") #读取文件
names(new_trait)[8]<-‘p.value‘
newdata <- subset(new_trait, p.value < 1e-4)#筛选pvalue列中小于1e-6的行，提取至新的文档newdata
newdata$trait<-i#创建新的列，标注该表型名称
tmp<-rbind(tmp,newdata) #与之前的文件合并
}
write.csv(tmp,‘rMVP.GLM.csv‘)

mvp结果分析

原文：https://www.cnblogs.com/xiaosagege/p/14288755.html