Question: Get all SNPs that are in high linkage disequilibrium (based on D') for a set of SNPs
1
gravatar for Mr Locuace
14 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 • 828 views
ADD COMMENTlink modified 14 months ago by chrchang5236.7k • written 14 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 14 months ago by zx87549.1k
4
gravatar for chrchang523
14 months ago by
chrchang5236.7k
United States
chrchang5236.7k 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 14 months ago by chrchang5236.7k
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: 1329 users visited in the last hour