Question: Calculating genetic distance from VCF to MAP format
0
gravatar for Cookie-san
6 months ago by
Cookie-san90
Cookie-san90 wrote:

I have a VCF from which I got the .ped and .map files using plink. However, the genetic distance (3rd column) has been set to 0. I want the map file to include the genetic distances. For my analysis, I am assuming a constant recombination rate of 15 kb/Cm.

Does plink have a command where you can provide the recombination rate so the output map file has the genetic distances? I was not able to find how the calculation is done from their manual.

I have thought about calculating them manually, but I am slightly confused with the formula -

If I have a VCF file having the following SNPs...

Chrom_01_v3     613
Chrom_01_v3     714
Chrom_01_v3     779

is it 714/15k or (714 - 613)/15k? In the latter case, since I have SNPs that are quite close to one another, I am getting genetic distances that are too small for my analysis.

My goal is to do IBD through the R package isoRelate which requires the input as .ped and .map format. My end goal is to get a Manhattan Plot.

population genomics • 263 views
ADD COMMENTlink modified 6 months ago by chrchang5237.3k • written 6 months ago by Cookie-san90
0
gravatar for chrchang523
6 months ago by
chrchang5237.3k
United States
chrchang5237.3k wrote:

https://www.cog-genomics.org/plink/1.9/input#cm_map

ADD COMMENTlink written 6 months ago by chrchang5237.3k
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: 1649 users visited in the last hour