Question: New to statistical genetics--need advice for calculating LD starting with genomic regions
gravatar for DaveC
4.4 years ago by
DaveC0 wrote:

I'm new to human statistical genetics so am a little lost.

Specifically, I would like to know what the most practical approach would be for me to tackle the following task:

Starting with a set of genomic regions (hg19) in bam format, how can I efficiently generate a set of SNPs (from 1000 Genomes, EUR) that are in LD (for example, r^2> 0.6) with each of these regions?

What I think I need to do are:

  1. find all SNPs in the regions of interest


  2. use software (eg Plink) to calculate pairwise LD's of all SNPs within (for example) 100kb of each of those SNPs

My questions are:

Am I on the right track? Is there and easier way? For instance, I have been considering SNAP and HaploReg but at least the web version of SNAP appears to be out of date.

What is the best (least error prone) way for me to achieve #1 above (ie, finding all SNP IDs in a given region? -- I'm completely baffled by how SNP IDs can change from build-to build

Thank you

snp • 1.3k views
ADD COMMENTlink modified 22 months ago by RamRS27k • written 4.4 years ago by DaveC0
Please log in to add an answer.


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