plot.twilight: Plot function for twilight objects
In twilight: Estimation of local false discovery rate

Description Usage Arguments Details Value Author(s) References See Also Examples

Interface to plotting of twilight objects. Produces one of five possible plots.

1 2	## S3 method for class 'twilight' plot(x, which = NULL, grayscale = FALSE, legend = TRUE, ...)

`x`	An object of class twilight.
`which`	A character string specifying the plot to be made.
`grayscale`	Logical value. Specifying whether plots should be colored or grayscaled. Only necessary for "scores" and "fdr".
`legend`	Logical value. Produces a legend for "scores" and "effectsize". A legend for "fdr" is only available if bootstrap estimates exist.
`...`	Additional graphical arguments.

Option which="scores" plots the expected versus the observed test statistics and draws confidence lines calculated from permutations. This plot is similar to plots in Tusher et al. (2001).

Option which="qvalues" plots q-values versus the number of rejected hypotheses.

Option which="fdr" plots p-values versus 1 - local false discovery rate, that is the conditional probability of being significant given the corresponding p-value, plus bootstrap estimates if computed. Bottom ticks are 1%-quantiles of p-values.

Option which="volcano" results in the volcano plot, that is observed score versus 1 - local false discovery rate. Bottom ticks are 1%-quantiles of scores.

Option which="effectsize" plots the observed fold change equivalent score distribution overlaid by the estimated effect size distribution, that is distribution of scores under the alternative. This plot is only available if function twilight.pval was applied with method="fc" as fold change equivalent scores are computed from log ratios.

Option which="table" tabulates histogram "effectsize".

No value is returned except for "table".

Stefanie Scheid

Scheid S and Spang R (2004): A stochastic downhill search algorithm for estimating the local false discovery rate, IEEE TCBB 1(3), 98–108.

Scheid S and Spang R (2005): twilight; a Bioconductor package for estimating the local false discovery rate, Bioinformatics 21(12), 2921–2922.

Scheid S and Spang R (2006): Permutation filtering: A novel concept for significance analysis of large-scale genomic data, in: Apostolico A, Guerra C, Istrail S, Pevzner P, and Waterman M (Eds.): Research in Computational Molecular Biology: 10th Annual International Conference, Proceedings of RECOMB 2006, Venice, Italy, April 2-5, 2006. Lecture Notes in Computer Science vol. 3909, Springer, Heidelberg, pp. 338-347.

Tusher VG, Tibshirani R and Chu G (2001): Significance analysis of mircroarrays applied to the ionizing response, PNAS 98(9), 5116–5121.

twilight.pval, twilight

1
2
3

### contains a twilight object created by function twilight
data(exfdr)
plot(exfdr)

Loading required package: splines
Loading required package: Biobase
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, basename, cbind, colMeans, colSums, colnames,
    dirname, 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

Welcome to Bioconductor

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