gravatar for Jautis

2 hours ago by

United States

Hello, I have a recombination map calculated from LDhelmet which has estimates of recombination rates for windows of the genome (in recombination events per generation per bp) and a separate file of SNP locations which I would like to transform from base number to location in centiMorgans. Is there an easy way to do this computationally? The best strategy I've been able to devise is, for each site, taking the sum of previous windows and then the proportion of the current window. However, that seems very inefficient when given hundreds of thousands of SNPs. Any suggestions for a way to calculate this (or program out there) that would be more efficient?

##Recombination rate data
 start     end      rate/bp     length 
 101       400      0.004       300
 401       900      0.0001      500

## SNPs
base pair
 800

# So for this SNP, the position would be 300*0.004 + (400/500)*0.0001

link

written
2 hours ago
by

Jautis290



Source link