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

相关推荐
医工交叉实验工坊5 小时前
R 语言 ComplexUpset 包实战:替代 Venn 图的高级集合可视化方案
算法·r语言
青春不败 177-3266-05205 小时前
R语言贝叶斯方法在生态环境领域中的高阶技术应用
r语言·贝叶斯·生态学·科研绘图·结构方程·环境科学·混合效应
love530love1 天前
怎么更新 cargo.exe ?(Rust 工具链)
人工智能·windows·python·rust·r语言
烟锁池塘柳05 天前
【R语言】R语言中 rbind() 与 merge() 的区别详解
sql·r语言
夏婵语冰6 天前
实用R语言机器学习指南:从数据预处理到模型实战(附配套学习资源)
开发语言·学习·r语言
Yolo566Q6 天前
“R语言+遥感”的水环境综合评价方法实践技术应用
开发语言·r语言
天意生信云6 天前
生信分析自学攻略 | R语言数据筛选和修改
数据结构·经验分享·r语言
烟锁池塘柳09 天前
【R语言】R 语言中 gsub 与正则表达式详解(含 POSIX 与 Perl 风格实例)
正则表达式·r语言·perl
梦想的初衷~9 天前
“R+遥感”的水环境综合评价方法实践技术应用
r语言·水环境·水文
天意生信云10 天前
生信分析自学攻略 | R软件和Rstudio的安装
经验分享·r语言