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

相关推荐
biomooc12 小时前
R 语言 | 绘图的文字格式(绘制上标、下标、斜体、文字标注等)
开发语言·r语言
Tiger Z2 天前
R 语言科研绘图第 6 期 --- 散点图-基础
r语言·贴图
LvManBa2 天前
R 常用的内置软件包及功能介绍
开发语言·r语言·rstudio
新知图书3 天前
R语言的数据结构-数据框
开发语言·r语言
新知图书4 天前
R语言的字符串操作
开发语言·r语言
新知图书5 天前
R语言函数简介
开发语言·python·r语言
还没想好1235 天前
mlr3机器学习AUC的置信区间提取
人工智能·机器学习·r语言·学习笔记
拓端研究室TRL6 天前
R语言混合模型回归GBTM群组轨迹模型绘图可视化研究
开发语言·人工智能·数据挖掘·回归·r语言
_刘文凯_7 天前
R环境配置 以及Debug方法 (VSCode, conda, 远程R)
vscode·r语言·conda
新知图书8 天前
R语言的数据结构-数组
开发语言·数据结构·r语言