LD filter for GWAS with plink
Hi all,
I'm trying to remove SNPs that are in linkage disequilibrium (r2 = 0.8) from my GWAS dataset. For the moment, I've been trying the following command:
plink --vcf myvcf.vcf --const-fid 0 --r2 --ld-window-r2 0.8 --recode vcf --out myvcf.out
It actually works and generates a .ldoutput file with pairwise r2 values for the SNPs. However, I was wondering if there is a method that directly removes those SNPs from the .vcfinstead of having to pass the output list given in the .ldfile to another plink --excludecommand.
Thanks for the help in advance!
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Extracting LD SNPs using PLINK2
written by Sanjar
-
How to Extract SNPs from Linkage Disequilibrium data with r2 > 0.6 threshold?
written by Raju •
-
how to use plink to get all the SNP LD score > 0.8 in matrix format?
written by yuhangchen0802 •
-
How does plink calculates linkage disequilibrium (r2) in unphased data?
written by serpalma.v •
-
Estimating LD score in PLINK
written by GabrielMontenegro
-
Get all SNPs that are in high linkage disequilibrium (based on D') for a set of SNPs
written by Mr Locuace
-
No output estimating LD between SNPs in VCF using PLINK
written by GabrielMontenegro
-
Using vcftools to filter SNPs by a linkage disequilibrium (r2) threshold
written by shyamie •
-
Duplicate ID in bed file
written by janhuang.cn
-
Linkage disequilibrium decay sanity check
written by biogirl
Related posts:
I have not seen any way to directly remove these with PLINK. I do not believe it is too much effort to have to use 2 commands, and it also helps with documentation by producing the output list of SNPs that are in / not in LD.