New to statistical genetics--need advice for calculating LD starting with genomic regions
0
0
Entering edit mode
9.4 years ago
DaveC • 0

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

    then

  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.8k views
ADD COMMENT

Login before adding your answer.

Traffic: 2997 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6