Given calcDinuclFreq
or calcDinuclFreqRef
results, this function
generates a violin plot of dinucleotide frequency
plotDinuclFreq(DNFDataTable)
A ggplot object plotting distribution of dinucleotide content in query regions
DNFDataTable = data.table::data.table(GC = rnorm(400, mean=0.5, sd=0.1),
CG = rnorm(400, mean=0.5, sd=0.5),
AT = rnorm(400, mean=0.5, sd=1),
TA = rnorm(400, mean=0.5, sd=1.5))
DNFPlot = plotDinuclFreq(DNFDataTable)
if (FALSE) { # \dontrun{
query = system.file("extdata", "vistaEnhancers.bed.gz", package="GenomicDistributions")
GRquery = rtracklayer::import(query)
refAssembly = 'hg19'
DNF = calcDinuclFreqRef(GRquery, refAssembly)
DNFPlot2 = plotDinuclFreq(DNF)
} # }