extractConsequenceEnrichment: Analyze data for enrichment of specific consequences

Description Usage Arguments Details Value Author(s) References See Also Examples

View source: R/analyze_switch_consequences.R

Description

This functions analyzes for enrichment of specific consequnces by for each set of opposing consequnces (fx. domain gain vs loss), by analyzing the fraction of events belonging to one of them.

Usage

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
extractConsequenceEnrichment(
    switchAnalyzeRlist,
    consequencesToAnalyze = 'all',
    alpha=0.05,
    dIFcutoff = 0.1,
    countGenes = TRUE,
    analysisOppositeConsequence=FALSE,
    plot=TRUE,
    localTheme = theme_bw(base_size = 12),
    minEventsForPlotting = 10,
    returnResult=TRUE,
    returnSummary=TRUE
)

Arguments

switchAnalyzeRlist

A switchAnalyzeRlist object where analyzeSwitchConsequences() have been run to identify consequences of isoform switches

consequencesToAnalyze

A string indicating which consequnces should be considered. See detail section of analyzeSwitchConsequences for describtion . Default is all consequences analyzed with analyzeSwitchConsequences.

alpha

The cutoff which the (callibrated) fdr correct p-values must be smaller than for calling significant switches. Default is 0.05.

dIFcutoff

The cutoff which the changes in (absolute) isoform usage must be larger than before an isoform is considered switching. This cutoff can remove cases where isoforms with (very) low dIF values are deemed significant and thereby included in the downstream analysis. This cutoff is analogous to having a cutoff on log2 fold change in a normal differential expression analysis of genes to ensure the genes have a certain effect size. Default is 0.1 (10%).

countGenes

A logic indicating whether it is the number of genes (if TRUE) or isoform switches (if FALSE) which primarly result in gain/loss that are counted. Default is TRUE.

analysisOppositeConsequence

A logic indicating whether reverse the analysis meaing if "Domain gains"" are analyze using default parameters setting analysisOppositeConsequence=TRUE will case the analysis to be performed on "Domain loss". The main effect is for the visual appearence of plot which will be mirrored (around the 0.5 fraction). Default is FALSE.

plot

A logic indicting whether the analysis should be plotted. Default is TRUE.

localTheme

General ggplo2 theme with which the plot is made, see ?ggplot2::theme for more info. Default is theme_bw(base_size = 14).

minEventsForPlotting

The minimum number of events (total gain/loss) must be pressent before the result is visualized. Default is 10.

returnResult

A logic indicating whether the analysis should be returned as a data.frame. Default is TRUE.

returnSummary

A logic indicating whether to return the statistical summary (if TRUE) or the underlying data (if FALSE). Default is TRUE.

Details

The significance test is performed with R's build in prop.test() with default parameters and resulting p-values are corrected via p.adjust() using FDR (Benjamini-Hochberg).

Value

If returnResult=TRUE a data.frame with the statisitcal summary for each oposing consequences in each comparison. If plot=TRUE a plot summarizing the proportions is also created of switches with specific consequences is created.

Author(s)

Kristoffer Vitting-Seerup

References

See Also

analyzeSwitchConsequences
extractSwitchSummary
extractConsequenceEnrichmentComparison
extractConsequenceGenomeWide

Examples

1
2
3
4
### Load exampled data
data("exampleSwitchListAnalyzed")

extractConsequenceEnrichment( exampleSwitchListAnalyzed)

kvittingseerup/IsoformSwitchAnalyzeR documentation built on July 20, 2019, 8:54 a.m.