hist {SSPA} | R Documentation |
The function hist
computes a histogram of p-values of the given pilot data. Two-sided p-values are computed either using the normal or Student t distribution based on the test statistics that are stored in an object of PilotData
.
hist(x, ...)
x |
object of PilotData |
... |
additional arguments to hist |
p-values calculated under the null hypothesis (non-differentially expressed genes) are assumed to be uniformly distributed on [0,1]. p-values calculated under the alternative hypothesis (differentially expressed genes) are assumed to accummelated near zero. The height of the flat part of the histogram is an indication of how many hypothesis are calculated under the null hypothesis.
an object of class "histogram", see hist
.
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)