runWilcox: Computing Multiple Wilcoxon Tests

Description Usage Arguments Value Author(s) See Also Examples

View source: R/anadiff.R

Description

This function provides a convenient way to compute the wilcoxon statistics, for each row of a data frame using the multtest package. It returns the raw and adjusted pvalues for each genes as well as the significance of the genes and a quantile-quantile plot.

Usage

1
runWilcox(data, labels, typeFDR = "FDR-BH", q = 0.05, plot = TRUE)

Arguments

data

A matrix, a data frame, or an ExpressionSet object. Each row of 'data' (or 'exprs(data)', respectively) must correspond to a gene, and each column to a sample.

labels

A vector of integers corresponding to observation (column) class labels. For 2 classes, the labels must be 0 and 1.

typeFDR

The method to apply fo the multiple testing correction.

q

A numeric value specifying the pvalue threshold.

plot

A logical value specifying if drawing plots or not.

Value

A matrix with the probes ID, the statistics, the raw p-value and the adjust p-value

Author(s)

Nicolas Servant, Eleonore Gravier, Pierre Gestraud, Cecile Laurent, Caroline Paccard, Anne Biton, Jonas Mandel, Bernard Asselain, Emmanuel Barillot, Philippe Hupe

See Also

mt.teststat,multiple.correction

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
## load data
data(marty)

## filtering data
##marty <- expFilter(marty, threshold=3.5, graph=FALSE)

##Class label 0/1
marty.type.num <- ifelse(marty.type.cl=="Her2+",0,1)

## run differential analysis
example.subset<-marty[1:100,]
out <- runWilcox(example.subset, labels=marty.type.num, typeFDR="FDR-BH", plot=FALSE)

EMA documentation built on March 26, 2020, 8:40 p.m.