Question: Calculate Ld Between All Pairs Of Snps, Both Within And Between Chromosomes
1
gravatar for misconstruction
5.3 years ago by
European Union
misconstruction20 wrote:

Hi Biostar!

After coming across a few papers claiming to find long range LD between SNPs on different chromosomes, I decided it would be fun to do this on my own dataset.

The calculation itself is not terribly difficult, but a bit bothersome due to the huge RAM requirements.

I have a an already filtered dataset with SNP data as a 01-matrix for every chromosome for several populations. I would like to calculate LD (either D' or r2) between every pair of SNP, so obviously a huge amount of comparisons (in the millions at the very least).

The plan is then to apply some stringent threshold to that output to only get SNPs in significant LD (and then compare that across populations).

I'm aware of several programs for calculating LD, but they seem to be mainly focused on LD-decay within each chromosome, and not between. I'm also worried they can't handle this amount of data.

Thanks in advance for any help!

genetics • 3.6k views
ADD COMMENTlink modified 5.2 years ago by chrchang5234.9k • written 5.3 years ago by misconstruction20

Can you get your SNPs into VCF format by any chance?

ADD REPLYlink written 5.3 years ago by Adrian Pelin2.3k
1
gravatar for chrchang523
5.2 years ago by
chrchang5234.9k
United States
chrchang5234.9k wrote:

PLINK 1.9 can handle this. See https://www.cog-genomics.org/plink2/ld#r . Parallelization of the calculation across multiple machines is supported.

ADD COMMENTlink written 5.2 years ago by chrchang5234.9k
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: 1028 users visited in the last hour