Generates a uniform genetic map by approximating 1 cm / Mbp. Only used for examples.
1 |
1 2 3 4 5 6 7 8 9 10 11 12 13 | 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 readGeneticMap
generateUniformGeneticMap()
pdf(file=tempfile())
plotRegionalGeneticMap(seq(1e6,100e6,by=1e6/2))
dev.off()
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.