Please use VolcanoPlotW() from package wrGraph. This function does NOT produce a plot any more.
VolcanoPlotW2(
Mvalue,
pValue = NULL,
useComp = 1,
filtFin = NULL,
ProjNa = NULL,
FCthrs = NULL,
FdrList = NULL,
FdrThrs = NULL,
FdrType = NULL,
subTxt = NULL,
grayIncrem = TRUE,
col = NULL,
pch = 16,
compNa = NULL,
batchFig = FALSE,
cexMa = 1.8,
cexLa = 1.1,
limM = NULL,
limp = NULL,
annotColumn = NULL,
annColor = NULL,
cexPt = NULL,
cexSub = NULL,
cexTxLab = 0.7,
namesNBest = NULL,
NbestCol = 1,
sortLeg = "descend",
NaSpecTypeAsContam = TRUE,
useMar = c(6.2, 4, 4, 2),
returnData = FALSE,
callFrom = NULL,
silent = FALSE,
debug = FALSE
)deprecated - returns nothing
(numeric or matrix) data to plot; M-values are typically calculated as difference of log2-abundance values and 'pValue' the mean of log2-abundance values;
M-values and p-values may be given as 2 columsn of a matrix, in this case the argument pValue should remain NULL
(numeric, list or data.frame) if NULL it is assumed that 2nd column of 'Mvalue' contains the p-values to be used
(integer, length=1) choice of which of multiple comparisons to present in Mvalue (if generated using moderTestXgrp())
(matrix or logical) The data may get filtered before plotting: If FALSE no filtering will get applied; if matrix of TRUE/FALSE it will be used as optional custom filter, otherwise (if Mvalue if an MArrayLM-object eg from limma) a default filtering based on the filtFin element will be applied
(character) custom title
(numeric) Fold-Change threshold (display as line) give as Fold-change and NOT log2(FC), default at 1.5, set to NA for omitting
(numeric) FDR data or name of list-element
(numeric) FDR threshold (display as line), default at 0.05, set to NA for omitting
(character) FDR-type to extract if Mvalue is 'MArrayLM'-object (eg produced by from moderTest2grp etc);
if NULL it will search for suitable fields/values in this order : 'FDR','BH',"lfdr" and 'BY'
(character) custom sub-title
(logical) if TRUE, display overlay of points as increased shades of gray
(character) custom color(s) for points of plot (see also par)
(integer) type of symbol(s) to plot (default=16) (see also par)
(character) names of groups compared
(logical) if TRUE figure title and axes legends will be kept shorter for display on fewer splace
(numeric) font-size of title, as expansion factor (see also cex in par)
(numeric) size of axis-labels, as expansion factor (see also cex in par)
(numeric, length=2) range of axis M-values
(numeric, length=2) range of axis FDR / p-values
(character) column names of annotation to be extracted (only if Mvalue is MArrayLM-object containing matrix $annot).
The first entry (typically 'SpecType') is used for different symbols in figure, the second (typically 'GeneName') is used as prefered text for annotating the best points (if namesNBest allows to do so.)
(character or integer) colors for specific groups of annoatation (only if Mvalue is MArrayLM-object containing matrix $annot)
(numeric) size of points, as expansion factor (see also cex in par)
(numeric) size of subtitle, as expansion factor (see also cex in par)
(numeric) size of text-labels for points, as expansion factor (see also cex in par)
(integer or character) number of best points to add names in figure; if 'passThr' all points passing FDR and FC-filtes will be selected;
if the initial object Mvalue contains a list-element called 'annot' the second of the column specified in argument annotColumn will be used as text
(character or integer) colors for text-labels of best points
(character) sorting of 'SpecType' annotation either ascending ('ascend') or descending ('descend'), no sorting if NULL
(logical) consider lines/proteins with NA in Mvalue$annot[,"SpecType"] as contaminants (if a 'SpecType' for contaminants already exits)
(numeric,length=4) custom margings (see also par)
(logical) optional returning data.frame with (ID, Mvalue, pValue, FDRvalue, passFilt)
(character) allow easier tracking of message(s) produced
(logical) suppress messages
(logical) additional messages for debugging
this function was replaced by plotPCAw)
set.seed(2005); mat <- matrix(round(runif(900),2), ncol=9)
Run the code above in your browser using DataLab