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

相关推荐
Biomamba生信基地20 小时前
R语言入门课| 04 R语言基本函数
开发语言·r语言·生信
十三画者1 天前
【工具】BioPred一个用于精准医疗中生物标志物分析的 R 软件包
数据挖掘·数据分析·r语言·数据可视化
小艳加油1 天前
MaxEnt物种分布建模全流程;R+ArcGIS+MaxEnt模型物种分布模拟、参数优化方法、结果分析制图与论文写作
arcgis·r语言·物种分布模拟·maxent模型
weixin_贾1 天前
MaxEnt模型进阶:基于R语言自动化与GIS空间建模的物种栖息地精准预测
r语言·生物多样性·物种分布
预测模型的开发与应用研究1 天前
R语言实现轨迹分析--traj和lcmm包体会
机器学习·数据分析·r语言
Hefin_H3 天前
怎么把 GitHub 中的 R 包项目手动安装到 conda 环境中?
r语言·github·conda
Faxxtty5 天前
【R语言】无法调用stats.dll的问题解决方案
开发语言·r语言
医学生_R语言5 天前
【R语言可视化】相关系数热图
开发语言·python·r语言
维维180-3121-14555 天前
R语言结构方程模型(SEM)在生态学领域中的应用
r语言·生态·结构方程
数据知道5 天前
【R】一文快速掌握 R 语言指令(R 备忘清单)
开发语言·r语言