Find SNP proxies in R

get_snp_proxies.ipynb

Let's consider a SNP x. What would it take to find all genome-wide proxies with r2 > 0.8?

We would need to compute the r2 between our chosen SNP x and every other SNP in the genome. If there are 20 million SNPs in the genome, then we need to compute 20 million r2 values and filter for the subset of SNPs that have r2 > 0.8.

The fastest way to compute these values is with the tomahawk software.

Since you do not want to compute these values, you might consider searching for a table that someone else already created and shared publicly. I am not aware of a publicly available LD table with genome-wide values for each SNP, but I'd like to make one when I have some time.

If you're interested to find additional R packages for genetics, this is a good reference:

https://CRAN.R-project.org/view=Genetics

Also check out Bioconductor.

Hi Slowkow,

Thank you for explaning this to me, I didn't think of it that way! I really appreciate your help.

Regards,

Randa