mvp结果分析(转载)
摘要:~~~ library(dplyr)library('reshape')#子集的包tmp<-data.frame() #创建空的数据集myfiles <- Sys.glob("*GLM.*csv") #搜索文件夹内所有的gwas结果文件print(paste(myfiles,"has been fo
阅读全文
GWAS流程
摘要:GWAS-3K分析 将现有名称进行统一化处理 第一种办法是从网站下载所有品种的名称信息,从中筛选目标470份样品的名称I找到对应的IRIS名称 第二种办法是利用my_list菜单提交自己的IRIS号,或其他的信息 利用plink提取目标品种的SNP信息 当前目录下,存在bed/bim/fam三个文件
阅读全文