Question: New to statistical genetics--need advice for calculating LD starting with genomic regions
0
gravatar for DaveC
3.5 years ago by
DaveC0
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

    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.1k views
ADD COMMENTlink modified 11 months ago by RamRS23k • written 3.5 years ago by DaveC0
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: 1921 users visited in the last hour