之前教程提到過Metal是可以做Meta分析涵亏，除了Metal，PLINK也可以進(jìn)行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）的效應(yīng)值是beta，而非OR连锯；
# qt 指的是輸出結(jié)果的效應(yīng)值用beta展示（默認(rèn)輸出效應(yīng)值用OR表示）归苍；
# --meta-analysis-snp-field 指的是輸入文件SNP列用SNP表示；
輸入文件gwas1.plink文件如下所示：
Pasted image 20231126222211.png
輸出結(jié)果如下所示：
Pasted image 20231126222351.png
其中运怖，不同列所代表的意思如下：
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(R) Random-effects meta-analysis p-value
BETA'/'OR' Fixed-effects BETA/OR estimate
BETA(R)'/'OR(R)' Random-effects BETA/OR estimate (DerSimonian and Laird)
Q p-value for Cochran's Q statistic
I2 heterogeneity index (0-100 scale)
致謝橙子牛奶糖（陳文燕）霜医，請(qǐng)用參考模版：We thank the blogger (orange_milk_sugar, Wenyan Chen) for XXX
感謝小可愛們多年來的陪伴， 我與你們一起成長(zhǎng)~
