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, 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 |
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. |
Adds elements to the vcf
INFO
field
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)