Question: Get all SNPs that are in high linkage disequilibrium (based on D') for a set of SNPs
1
gravatar for Mr Locuace
18 months ago by
Mr Locuace100
Chile
Mr Locuace100 wrote:

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

ld plink snp • 1.0k views
ADD COMMENTlink modified 18 months ago by chrchang5237.1k • written 18 months ago by Mr Locuace100
2

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

ADD REPLYlink written 18 months ago by zx87549.4k
4
gravatar for chrchang523
18 months ago by
chrchang5237.1k
United States
chrchang5237.1k wrote:

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 COMMENTlink written 18 months ago by chrchang5237.1k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1240 users visited in the last hour