GWAS:plink进行meta分析

之前教程提到过Metal是可以做Meta分析,除了Metal,PLINK也可以进行Meta分析。

命令如下所示:

plink --meta-analysis gwas1.plink  gwas2.plink gwas3.plink + logscale qt --meta-analysis-snp-field SNP --meta-analysis-chr-field CHR --meta-analysis-bp-field BP --meta-analysis-a1-field A1 --meta-analysis-a2-field A2  --meta-analysis-se-field SE --out allqt
# logscale,指的是输入文件(比如gwas1.plink)的效应值是beta,而非OR;
# qt 指的是输出结果的效应值用beta展示(默认输出效应值用OR表示);
# --meta-analysis-snp-field 指的是输入文件SNP列用SNP表示;

输入文件gwas1.plink文件如下所示:

输出结果如下所示:

其中,不同列所代表的意思如下:

CHR Chromosome code.

BP Base-pair coordinate.

SNP Variant identifier

A1 Allele 1.

A2 Allele 2.

N Number of valid studies for variant

P Fixed-effects meta-analysis p-value

P® Random-effects meta-analysis p-value

BETA'/'OR' Fixed-effects BETA/OR estimate

BETA®'/'OR®' Random-effects BETA/OR estimate (DerSimonian and Laird)

Q p-value for Cochran's Q statistic

I2 heterogeneity index (0-100 scale)


致谢橙子牛奶糖(陈文燕),请用参考模版:We thank the blogger (orange_milk_sugar, Wenyan Chen) for XXX

感谢小可爱们多年来的陪伴, 我与你们一起成长~

相关推荐
tRNA做科研20 天前
最新保姆级Linux下安装与使用conda:从下载配置到使用全流程
linux·服务器·conda·生物信息学·生物信息·计算生物学·基因组
dundunmm1 个月前
论文阅读:SIMBA: single-cell embedding along with features
论文阅读·人工智能·数据挖掘·embedding·生物信息·多组学细胞数据·单组学
Red Red3 个月前
GEO数据库提取疾病样本和正常样本|GEO数据库区分疾病和正常样本|直接用|生物信息|生信
开发语言·数据库·笔记·学习·r语言·c#·生物信息
tRNA做科研3 个月前
Bio-Linux-shell详解-2-基本Shell命令快速掌握
linux·运维·服务器·生物信息·计算生物学
Red Red3 个月前
GEO数据的下载和处理|GEO数据转换为Gene symbol|GEO注释文件提取symbol|查看样本标签|查看GEO数据疾病或正常|生物信息基础
数据库·笔记·学习·r语言·生物信息·geo数据库
图灵生信4 个月前
python脚本:输入基因名,通过爬虫的方式获取染色体上的location。
生物信息
Crayon小鱼干5 个月前
怎样使用 Juicer tools 的 dump 命令将.hic文件转换为交互矩阵matrix计数文件 (Windows)
生物信息·juicer tools·.hic·hic数据·schic·contact matrix·交互矩阵
qq_214782616 个月前
新工具:轻松可视化基因组,部分功能超IGV~
机器学习·信息可视化·数据分析·生物信息
徐洲更hoptop6 个月前
使用ShinyCell展示你的单细胞数据
生物信息
qq_214782617 个月前
Python散点图矩阵代码模版
数据库·python·mysql·信息可视化·数据挖掘·数据分析·生物信息