plotRegionalGeneticMap {sim1000G}

Generates a plot of the genetic map for a specified region.

Description

The plot shows the centimorgan vs base-pair positions. The position of markers that have been read is also depicted as vertical lines

Usage

plotRegionalGeneticMap(bp)

Arguments

bp

Vector of base-pair positions to generate a plot for library("sim1000G")

examples_dir = system.file("examples", package = "sim1000G") vcf_file = sprintf(" vcf = readVCF( vcf_file, maxNumberOfVariants = 100, min_maf = 0.12)

# For realistic data use the function readGeneticMap generateUniformGeneticMap()

pdf(file=tempfile()) plotRegionalGeneticMap(seq(1e6,100e6,by=1e6/2)) dev.off()