generateFakeWholeGenomeGeneticMap {sim1000G}

Generates a fake genetic map that spans the whole genome.

Description

Generates a fake genetic map that spans the whole genome.

Usage

generateFakeWholeGenomeGeneticMap(vcf)

Arguments

vcf

A vcf file read by function readVCF.

Examples

library("sim1000G")

examples_dir = system.file("examples", package = "sim1000G")
vcf_file = sprintf("%s/region.vcf.gz", examples_dir)
vcf = readVCF( vcf_file, maxNumberOfVariants = 100 ,
               min_maf = 0.12 ,max_maf = NA)

# For realistic data use the function
# downloadGeneticMap
generateFakeWholeGenomeGeneticMap(vcf)

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