setMappingBiasVcf {PureCN} | R Documentation |
Function to set mapping bias for each variant in the provided
CollapsedVCF
object. By default, it returns the same value for all
variants, but a mapping bias file can be provided for position-specific
mapping bias calculation.
setMappingBiasVcf(vcf, tumor.id.in.vcf = NULL, mapping.bias.file = NULL, normal.panel.vcf.file = NULL, min.normals = 2, smooth = TRUE, smooth.n = 5)
vcf |
|
tumor.id.in.vcf |
Id of tumor in case multiple samples are stored in VCF. |
mapping.bias.file |
A precomputed mapping bias database
obtained by |
normal.panel.vcf.file |
Deprecated, use |
min.normals |
Deprecated. Minimum number of normals with heterozygous SNP for
calculating position-specific mapping bias. Requires
|
smooth |
Impute mapping bias of variants not found in the panel by
smoothing of neighboring SNPs. Requires |
smooth.n |
Number of neighboring variants used for smoothing. |
A data.frame
with elements
bias |
A |
pon.count |
A |
shape1, shape2 |
Fit of a beta distribution. |
Markus Riester
# This function is typically only called by runAbsoluteCN via # fun.setMappingBiasVcf and args.setMappingBiasVcf. vcf.file <- system.file("extdata", "example.vcf.gz", package="PureCN") vcf <- readVcf(vcf.file, "hg19") vcf.bias <- setMappingBiasVcf(vcf)