Question: Get all SNPs that are in high linkage disequilibrium (based on D') for a set of SNPs
1
gravatar for Mr Locuace
5 months ago by
Mr Locuace90
Chile
Mr Locuace90 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 • 423 views
ADD COMMENTlink modified 5 months ago by chrchang5235.4k • written 5 months ago by Mr Locuace90
2

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

ADD REPLYlink written 5 months ago by zx87547.7k
4
gravatar for chrchang523
5 months ago by
chrchang5235.4k
United States
chrchang5235.4k 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 5 months ago by chrchang5235.4k
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: 1462 users visited in the last hour