R语言实现GWAS meta分析(1)

1、基于数据集的Meta分析

复制代码
datafilename=c("data1.txt","data2.txt"),
setwd(workdir)
library(Metalgwas)
a=1
name1=c()
for(i in datafilename){  assign(paste("file",a,sep=""),data.table::fread(paste(getwd(),"/","data",a,".txt",sep="")))
    name1=c(name1,paste(workdir,"//","data",a,".txt",sep=""))
    a=a+1}
METAL_gwas(GWASfile =name1, Metal_path = "XXXX",#/metal.exe : 放入下载好的软件
                    save_path = "./", GWAS_name = "dat1", 
                    type = "outcome",Analysis_Scheme = "SE", 
                    heterogeneity = TRUE, Average_FREQ = TRUE, 
                    MinMax_FREQ = TRUE, Filter_N = FALSE, Filter_MAF = 0.01, 
                    Genomic_Control_Correction = FALSE, Sample_Overlap_Correction = FALSE)

需要下载metal:Metal - Meta Analysis Helper for SNP data

得出来结果在文件夹中dat1.txt

相关推荐
Tiger Z1 天前
R 语言科研绘图第 41 期 --- 桑基图-基础
开发语言·r语言·贴图
jerry2011082 天前
R语言之rjava版本不匹配解决方法
开发语言·r语言
Tiger Z5 天前
R 语言科研绘图 --- 饼状图-汇总
开发语言·人工智能·程序人生·r语言·贴图
maizeman1265 天前
R语言——方差分析2
开发语言·r语言·可视化·方差分析·单因素·多重t检验多重比较·tukey多重比较
zhanghongyi_cpp8 天前
R语言操作练习2
r语言
jerry2011088 天前
R语言之环境清理
开发语言·r语言
TDengine (老段)9 天前
TDengine 语言连接器(R语言)
大数据·数据库·物联网·r语言·时序数据库·tdengine·iotdb
maizeman12610 天前
R语言——分布和公式
开发语言·r语言·公式·随机数·分布
Tiger Z10 天前
R 语言科研绘图第 39 期 --- 饼状图-旭日
开发语言·程序人生·r语言·贴图
清同趣科研11 天前
扩增子分析|基于R语言microeco包进行微生物群落网络分析(network网络、Zi-Pi关键物种和subnet子网络图)
r语言