Get all SNPs that are in high linkage disequilibrium (based on D') for a set of SNPs
1
2
Entering edit mode
5.5 years ago
Mr Locuace ▴ 180

Hello, I have phased plink files of a population. I would like to obtain all SNPs that are in high linkage disequilibrium (LD) for a set of 10 SNPs. I know how to do it according to a r2 threshold, but I need to set this threshold according to D'

For instance, this script will output all variants in high LD (r2 > 0.6) for a list of 10 SNPs within a range of 1000 Kb.

plink \
    --bfile myfile \
    --r2 --ld-snp-list list_10_snps.txt \
    --ld-window-kb 1000 --ld-window-r2 0.6

Thanks very much

plink SNP LD • 2.4k views
ADD COMMENT
2
Entering edit mode

R Bioconductor SNPRelate package has snpgdsLDpruning function with option of filtering on "dprime". This package can import plink format files.

ADD REPLY
4
Entering edit mode
5.5 years ago

This isn't directly supported by plink, but you can use "--r2 dprime" to include D' in the main LD report, and then follow up with e.g. an awk one-liner to filter on D'.

ADD COMMENT

Login before adding your answer.

Traffic: 2643 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6