diagplotVolcano {metaseqR2} | R Documentation

## (Interactive) volcano plots of differentially expressed genes

### Description

This function plots a volcano plot or returns a JSON string which is used to render aninteractive in case of HTML reporting.

### Usage

        diagplotVolcano(f, p, con = NULL, fcut = 1, pcut = 0.05,
            altNames = NULL, output = "x11", path = NULL, ...)

### Arguments

`f` | the fold changes which are to be plotted on the x-axis.
`p` | the p-values whose -log10 transformation is going to be plotted on the y-axis.
`con` | an optional string depicting a name (e.g. the contrast name) to appear in the title of the volcano diagplot.
`fcut` | a fold change cutoff so as to draw two vertical lines indicating the cutoff threshold for biological significance.
`pcut` | a p-value cutoff so as to draw a horizontal line indicating the cutoff threshold for statistical significance.
`altNames` | an optional vector of names, e.g. HUGO gene symbols, alternative or complementary to the unique names of `f` or `p` (one of them must be named!). It is used only in JSON output.
`output` | one or more R plotting device to direct the plot result to. Supported mechanisms: `"x11"` (default), `"png"`, `"jpg"`, `"bmp"`, `"pdf"`, `"ps"` or `"json"`. The latter is currently available for the creation of interactive volcano plots only when reporting the output, through the highcharts javascript library.
`path` | the path to create output files.
`...` | further arguments to be passed to plot devices, such as parameter from `par`.

### Value

The filenames of the plots produced in a named list with names the `which.plot` argument. If `output="x11"`, no output filenames are produced.

### Author(s)

Panagiotis Moulos

### Examples

    dataMatrix <- metaseqR2:::exampleCountData(5000)
    sampleList <- list(A=c("A1","A2"),B=c("B1","B2","B3"))
    contrast <- "A_vs_B"
    M <- normalizeEdger(dataMatrix,sampleList)
    p <- statEdger(M,sampleList,contrast)
    ma <- apply(M[,sampleList$A],1,mean)
    mb <- apply(M[,sampleList$B],1,mean)
    f <- log2(ifelse(mb==0,1,mb)/ifelse(ma==0,1,ma))
    diagplotVolcano(f,p[[1]],con=contrast)
    #j <- diagplotVolcano(f,p[[1]],con=contrast,output="json")

[Package metaseqR2 version 1.0.11 Index]