PlotSectionEnrichment: Perform enrichment analysis of sites per transcript region...

Description Usage Arguments Value Author(s) Examples

View source: R/plot.R

Description

Perform enrichment analysis of the number of positive sites in txLoc.pos relative to the number of null sites in txLoc.neg per transript region and plot results. Enrichment/depletion is evaluated using (multiple) Fisher's exact test(s). Multiple hypothesis testing correction is applied following the method of Bejamini and Hochberg.

Usage

1
PlotSectionEnrichment(txLoc.pos, txLoc.neg, xAxisLblFmt = 2)

Arguments

txLoc.pos

A txLoc object. These correspond to the positive sites.

txLoc.neg

A txLoc object. These correspond to the negative (null) sites.

xAxisLblFmt

Plot extended axis labels. Default is 2. See ??? for details.

Value

NULL.

Author(s)

Maurits Evers, maurits.evers@anu.edu.au

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
## Not run: 
bedFile <- system.file("extdata",
                       "miCLIP_m6A_Linder2015_hg38.bed",
                       package = "RNAModR")
sites <- ReadBED(bedFile)
pos <- SmartMap(sites, id = "m6A", refGenome = "hg38")
neg <- GenerateNull(pos, method = "nuclAbundance")
PlotSectionEnrichment(
    FilterTxLoc(pos, c("5'UTR", "CDS", "3'UTR")),
    FilterTxLoc(neg, c("5'UTR", "CDS", "3'UTR")))

## End(Not run)

mevers/RNAModR documentation built on Nov. 17, 2019, 9:11 a.m.