volcanoPlot: Making a volcano plot

Description Usage Arguments Value See Also Examples

View source: R/stattest.R

Description

Makes a volcano plot of the data.

Usage

1
volcanoPlot(header, dataset, col4plotx, col4ploty, col4anno, plotTitle, sigLevel, showPlot)

Arguments

header

the header of a dataset file generated with generateDatasetFile

dataset

an R data frame generated with generateDatasetFile

col4plotx

a character vector specifying the name of the column containing the intensity values, usually SigIntensity

col4ploty

a character vector specifying the name of the dataset column containing the corresponding p-values

col4anno

a character string specifying the name of the dataset column to be used to define the replicate, e.g. "GeneName" or "Internal_GeneID".

plotTitle

the plot title

sigLevel

the significance level for the p-value, indicating where a horizontal green line will be drawn

showPlot

0 or 1. 1 will open one or several plot windows in the R GUI, 0 will only save the plot(s) without opening windows

Value

Plots the intensity values against the negative decadic logarithm of the p-values. A green horizontal line is drawn at the specified significance level.

The plot is saved in a pdf and a png file named after the experiment name specified in the header concatenated with the plotTitle.

The function returns the plot name.

See Also

Ttest

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
data(exampleHeader, package="RNAither")
data(pValVec1, package="RNAither")
data(scoredDataset1, package="RNAither")

##for details on the generation of pValVec1 and scoredDataset1, see the example of the Ttest function linked above.

scoredHits1 <- hitselectionPval(scoredDataset1, pValVec1, "SigIntensity", "pValue.ttest_l", 0.05, 
"GeneName", "pvalue_testfile1.txt")

hitDataset1 <- scoredHits1[[1]]
hitvector1 <- scoredHits1[[2]]

volcano_name <- volcanoPlot(header, hitDataset1, "SigIntensity", "pValue.ttest_l", "GeneName", 
"Volcano Plot", 0.05, 1)

Example output

Loading required package: topGO
Loading required package: BiocGenerics
Loading required package: parallel

Attaching package: 'BiocGenerics'

The following objects are masked from 'package:parallel':

    clusterApply, clusterApplyLB, clusterCall, clusterEvalQ,
    clusterExport, clusterMap, parApply, parCapply, parLapply,
    parLapplyLB, parRapply, parSapply, parSapplyLB

The following objects are masked from 'package:stats':

    IQR, mad, sd, var, xtabs

The following objects are masked from 'package:base':

    Filter, Find, Map, Position, Reduce, anyDuplicated, append,
    as.data.frame, cbind, colMeans, colSums, colnames, do.call,
    duplicated, eval, evalq, get, grep, grepl, intersect, is.unsorted,
    lapply, lengths, mapply, match, mget, order, paste, pmax, pmax.int,
    pmin, pmin.int, rank, rbind, rowMeans, rowSums, rownames, sapply,
    setdiff, sort, table, tapply, union, unique, unsplit, which,
    which.max, which.min

Loading required package: graph
Loading required package: Biobase
Welcome to Bioconductor

    Vignettes contain introductory material; view with
    'browseVignettes()'. To cite Bioconductor, see
    'citation("Biobase")', and for packages 'citation("pkgname")'.

Loading required package: GO.db
Loading required package: AnnotationDbi
Loading required package: stats4
Loading required package: IRanges
Loading required package: S4Vectors

Attaching package: 'S4Vectors'

The following object is masked from 'package:base':

    expand.grid


Loading required package: SparseM

Attaching package: 'SparseM'

The following object is masked from 'package:base':

    backsolve


groupGOTerms: 	GOBPTerm, GOMFTerm, GOCCTerm environments built.

Attaching package: 'topGO'

The following object is masked from 'package:IRanges':

    members

Loading required package: RankProd
Loading required package: Rmpfr
Loading required package: gmp

Attaching package: 'gmp'

The following objects are masked from 'package:base':

    %*%, apply, crossprod, matrix, tcrossprod

C code of R package 'Rmpfr': GMP using 64 bits per limb


Attaching package: 'Rmpfr'

The following objects are masked from 'package:IRanges':

    cbind, pmax, pmin, rbind

The following objects are masked from 'package:S4Vectors':

    cbind, rbind

The following objects are masked from 'package:BiocGenerics':

    cbind, pmax, pmin, rbind

The following objects are masked from 'package:stats':

    dbinom, dnorm, dpois, pnorm

The following objects are masked from 'package:base':

    cbind, pmax, pmin, rbind

Loading required package: prada
Loading required package: RColorBrewer
Loading required package: grid

Attaching package: 'grid'

The following object is masked from 'package:topGO':

    depth

Loading required package: rrcov
Loading required package: robustbase

Attaching package: 'robustbase'

The following object is masked from 'package:Biobase':

    rowMedians

Scalable Robust Estimators with High Breakdown Point (version 1.4-4)


Attaching package: 'RNAither'

The following object is masked from 'package:IRanges':

    trim

RNAither documentation built on Nov. 8, 2020, 8:06 p.m.