LocalFDR {anapuce}R Documentation

Local FDR

Description

Estimation of local false discovery rates

Usage

LocalFDR(dataf = dataf, filegraph = "Graphics", filedon = "FDRFile", method = NULL, lambda0 = 0.5, smoothing = "1", seuil = c(0.01, 0.05, 0.1, 0.2), mm = c(3, 5, 15, NA))

Arguments

dataf input data file with two columns (1 = gene name, 2= ordered raw pvalues).
filegraph name of the graphics file.
filedon name of the output data file.
method how to estimate m0 ? , choice = adaptive,conservative,bootstrap,smoother or NULL.
lambda0 value used in calculating m0.
smoothing "1" for the initial published method, "2" for the PAVA method isotonic (monotonely increasing nonparametric) least squares regression - see P. Broberg (2005).
seuil threshold
mm

Value

- A file of graphics (.ps)
- A data file with 4 columns : gene name, raw pvalues and two columns corresponding to smoothed FDR values.

Author(s)

J. Aubert

References

Aubert J, Bar-Hen A, Daudin J.J, Robin S (2004) Determination of the differentially expressed genes in microarrays experiments using local FDR, BMC Bioinformatics, 5:125.

Aubert J , Bar-Hen A, Daudin J.J, Robin S (2005) Correction: Determination of the differentially expressed genes in microarray experiments using local FDR, BMC Bioinformatics, 6:42.

Per Broberg (2005) A comparative review of estimates of the proportion unchanged genes and the false discovery rate, BMC Bioinformatics 6:199


[Package anapuce version 2.0 Index]