szpiech / selscan

Haplotype based scans for selection
GNU General Public License v3.0
111 stars 33 forks source link

Question about genetic map #95

Open Youpu-Chen opened 1 year ago

Youpu-Chen commented 1 year ago

Hi, I have a question about how to analyze the 1000 Genomes Project (hg38) dataset. I searched across the Internet, and I did not find a usable genetic map for the dataset, it always report the error below: 135fd9ae6783eac082ede4d77989d2d

Is there any solution to solve this error?

szpiech commented 1 year ago

I've used the hg38 genetic map here: https://bochet.gcc.biostat.washington.edu/beagle/genetic_maps/. You may need to slightly reformat it for use with selscan. I believe this map is a liftOver of a map inferred on an earlier build.

TanXinjiang commented 1 year ago

You can use "predictGMAP"to interpolate the loci not included in the genetic maps from the 1000 Genomes project

Youpu-Chen commented 1 year ago

Thank you! I've already tried predictGMAP out, but I think PLINK might just be useful enough for the analysis I've been conducting right now.

TanXinjiang @.***> 于2023年4月23日周日 16:57写道:

You can use "predictGMAP https://github.com/szpiech/predictGMAP"to interpolate the loci not included in the genetic maps from the 1000 Genomes project

— Reply to this email directly, view it on GitHub https://github.com/szpiech/selscan/issues/95#issuecomment-1519000235, or unsubscribe https://github.com/notifications/unsubscribe-auth/AU7RW2FIY5SHW33SSHN46JTXCTVFZANCNFSM6AAAAAAWLUZATY . You are receiving this because you authored the thread.Message ID: @.***>

Youpu-Chen commented 1 year ago

I've used the hg38 genetic map here: https://bochet.gcc.biostat.washington.edu/beagle/genetic_maps/. You may need to slightly reformat it for use with selscan. I believe this map is a liftOver of a map inferred on an earlier build.

Thank you, I also found that Eagle also provide the hg38 genetic map.