plot {SSPA} | R Documentation |
The function plot
plots an empirical cumulative distribution function of p-values of the pilot data. The test statistics that are given as input to PilotData-class
are transformed to two-sided p-values either using the normal or Student t distribution.
plot(x, y, ...)
x |
object of class PilotData-class |
y |
unused item |
... |
additional arguments given to plot or par |
Empirical cumulative distribution function of p-values. The line at angle of 45 degrees represents the theoretical CDF of a uniform distribution as expected when all genes are non-differentially expressed. An accumulation of p-values near zero indicates a certain number of differentially expressed genes.
Maarten van Iterson
library(multtest) data(golub) teststat <- mt.teststat(golub, golub.cl) table(golub.cl) pd <- pilotData(name="golub", testStatistics=teststat, sampleSizeA=11, sampleSizeB=27) hist(pd) plot(pd)