generateFakeWholeGenomeGeneticMap: Generates a fake genetic map that spans the whole genome.
In adimitromanolakis/sim1000G: Genotype Simulations for Rare or Common Variants Using Haplotypes from 1000 Genomes
Description
Usage
Arguments
Examples
View source: R/recombination.R
Description
Generates a fake genetic map that spans the whole genome.
Usage
1
Arguments
vcf
A vcf file read by function readVCF.
Examples
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()
adimitromanolakis/sim1000G documentation built on Sept. 29, 2021, 3:44 a.m.
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Description Usage Arguments Examples
View source: R/recombination.R
Description
Generates a fake genetic map that spans the whole genome.
Usage
1 |
Arguments
vcf |
A vcf file read by function readVCF. |
Examples
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()
|
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Embedding an R snippet on your website
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.