hist {SSPA}R Documentation

Histogram of p-values of the pilot data

Description

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.

Usage

hist(x, ...)

Arguments

x object of PilotData
... additional arguments to hist

Details

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.

Value

an object of class "histogram", see hist.

Author(s)

Maarten van Iterson

See Also

hist and plot

Examples

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)

[Package SSPA version 1.0.0 Index]