volcanoPlot {circRNAprofiler} | R Documentation

## Plot differential circRNA expression results

### Description

The function volcanoPlot() generates a volcano plot with the results of the differential expression analysis.

### Usage

    volcanoPlot(res, log2FC = 1, padj = 0.05, title = "", gene = FALSE,
      setxLim = FALSE, xlim = c(-8, 8), setyLim = FALSE, ylim = c(0,
      5), color = "blue")

### Arguments

`res` | A data frame containing the the differential expression resuls. It can be generated with `getDeseqRes` or `getEdgerRes`.
`log2FC` | An integer specifying the log2FC cut-off. Deafult value is 1.
`padj` | An integer specifying the adjusted P value cut-off. Deafult value is 0.05.
`title` | A character string specifying the title of the plot.
`gene` | A logical specifying whether to add the circRNA host gene names to the plot. Deafult value is FALSE.
`setxLim` | A logical specifying whether to set x scale limits. If TRUE the value in xlim will be used. Deafult value is FALSE.
`xlim` | A numeric vector specifying the lower and upper x axis limits. Deafult values are c(-8 , 8).
`setyLim` | A logical specifying whether to set y scale limits. If TRUE the value in ylim will be used. Deafult value is FALSE.
`ylim` | An integer specifying the lower and upper y axis limits Deafult values are c(0, 5).
`color` | A string specifying the color of the differentially expressed circRNAs. Default value is "blue".

### Value

A ggplot object.

### Examples

    # Load data frame containing detected back-spliced junctions
    data("mergedBSJunctions")

    pathToExperiment <- system.file("extdata", "experiment.txt",
        package ="circRNAprofiler")

    # Filter circRNAs
    filterdCirc <- filterCirc(
        mergedBSJunctions,
        allSamples = FALSE,
        min = 5,
        pathToExperiment)

    # Find differentially expressed circRNAs
     deseqResBvsA <- getDeseqRes(
        filterdCirc,
        condition = "A-B",
        pAdjustMethod = "BH",
        pathToExperiment)

    # Plot
    p <- volcanoPlot(
        deseqResBvsA,
        log2FC = 1,
        padj = 0.05,
        title = "",
        setxLim = TRUE,
        xlim = c(-8 , 7.5),
        setyLim = FALSE,
        ylim = c(0 , 4),
        gene = FALSE)
    p

[Package circRNAprofiler version 1.0.0 Index]