hist                  package:SSPA                  R Documentation

_H_i_s_t_o_g_r_a_m _o_f _p-_v_a_l_u_e_s _o_f _t_h_e _p_i_l_o_t _d_a_t_a

_D_e_s_c_r_i_p_t_i_o_n:

     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'.

_U_s_a_g_e:

     hist(x, ...)

_A_r_g_u_m_e_n_t_s:

       x: object of 'PilotData'

     ...: additional arguments to 'hist'

_D_e_t_a_i_l_s:

     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.

_V_a_l_u_e:

     an object of class "histogram", see 'hist'.

_A_u_t_h_o_r(_s):

     Maarten van Iterson

_S_e_e _A_l_s_o:

     'hist' and 'plot'

_E_x_a_m_p_l_e_s:

     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)

