plotHistWSAFPlotly: WSAF histogram

In DEploid: Deconvolute Mixed Genomes with Unknown Proportions

WSAF histogram

Description

Produce histogram of the allele frequency within sample.

Usage

plotHistWSAFPlotly(obsWSAF, exclusive = TRUE, title = "Histogram 0<WSAF<1")

Arguments

obsWSAF	Observed allele frequency within sample
exclusive	When TRUE 0 < WSAF < 1; otherwise 0 <= WSAF <= 1.
title	Figure title, "Histogram 0<WSAF<1" by default

Value

histogram

Examples

# Example 1
refFile <- system.file("extdata", "PG0390-C.test.ref", package = "DEploid")
altFile <- system.file("extdata", "PG0390-C.test.alt", package = "DEploid")
PG0390Coverage <- extractCoverageFromTxt(refFile, altFile)
obsWSAF <- computeObsWSAF(PG0390Coverage$altCount, PG0390Coverage$refCount)
plotHistWSAFPlotly(obsWSAF)
myhist <- plotHistWSAFPlotly(obsWSAF)

# Example 2
vcfFile <- system.file("extdata", "PG0390-C.test.vcf.gz", package = "DEploid")
PG0390CoverageV <- extractCoverageFromVcf(vcfFile, "PG0390-C")
obsWSAF <- computeObsWSAF(PG0390CoverageV$altCount, PG0390CoverageV$refCount)
plotHistWSAFPlotly(obsWSAF)
myhist <- plotHistWSAFPlotly(obsWSAF)

DEploid documentation built on April 4, 2025, 2:16 a.m.