Description Usage Arguments Examples
Generates a fake genetic map that spans the whole genome.
1 |
vcf |
A vcf file read by function readVCF. |
1 2 3 4 5 6 7 8 9 10 11 12 13 14 | 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()
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.