Question: How To Use Plink To Compute The Ld Between Two Sets Of Snps?
1
gravatar for Fayue1015
6.7 years ago by
Fayue1015190
Singapore
Fayue1015190 wrote:

I have genotype data and I actually want to know the linkage relationship between two genes? So here is how I did, I map the snps which are reside to genes within +-10kb, then for two genes, each has an associated set of SNPs, how do I use plink to compute the LD between two sets of SNPs? Can I make an conclusion like this? If snp set A mapped to gene A has an average LD>0.8 with snp set B mapped to gene B, I would say gene A and gene B are linked? I find plink can only compute all pairs of LD for a given snp list. Thanks for any advice

ld snp • 3.1k views
ADD COMMENTlink modified 4.5 years ago by chrchang5234.6k • written 6.7 years ago by Fayue1015190
1
gravatar for chrchang523
4.5 years ago by
chrchang5234.6k
United States
chrchang5234.6k wrote:

This isn't directly supported by plink.  Here are two ways to reduce the filtering workload, though:

1. Add "--extract [SNP set {A+B}] --ld-snp-list [SNP set A]".  It'll only be necessary to remove A-A pairs from the final report in that case.

2. Add "--extract [SNP set {A+B}]" and request matrix output.  The values you care about will then all be in the lower-left (or upper-right) part of the matrix.

ADD COMMENTlink written 4.5 years ago by chrchang5234.6k
0
gravatar for Davy
6.7 years ago by
Davy360
United States
Davy360 wrote:

If you combine your two sets of SNPs into one list, and compute all pairwise LD, then simply remove all the LD values from pairs resulting from two SNPs in the same set.

You'd be best of trying something like this however with an LD mapping tool, like Haploview, and LocusZoom. How many SNPs in set A and set B are we talking about?

ADD COMMENTlink written 6.7 years ago by Davy360

set A and set B may contain less 100 snps, I will try Haploview and LocusZoom I do not want to use plink to do all pairwise than remove LD, too trouble some

ADD REPLYlink written 6.7 years ago by Fayue1015190

Hi, are you sure Haploview or LocusZoom has such functions? It seems do not have such function to me

ADD REPLYlink written 6.7 years ago by Fayue1015190
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: 1916 users visited in the last hour