twilight              package:twilight              R Documentation

_E_s_t_i_m_a_t_i_o_n _o_f _l_o_c_a_l _f_a_l_s_e _d_i_s_c_o_v_e_r_y _r_a_t_e

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

     The function performs the successive exclusion procedure (SEP) as
     described in Scheid and Spang (2004).

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

     twilight(xin, lambda = NULL, B = 0, boot.ci = 0.95, clus = NULL, verbose = TRUE)

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

     xin: Numerical vector of u-values, p-values or a 'twilight'
          object. 

  lambda: Numerical value denoting the penalty factor. If not
          specified, the function searchs for an appropriate
          regularization parameter. 

       B: Numerical value specifying the number of bootstrap samples.
          If not specified, no bootstrap estimates are calculated. 

 boot.ci: Numerical value denoting the probability value for bootstrap
          confidence intervals of local false discovery rate and prior
          pi0. 

    clus: A list or numerical value to be passed to
          'makeCluster(clus,...)' in 'library(snow)'. If specified,
          bootstrapping is performed in parallel. No checks on 'clus'
          are performed. Please make sure that 'makeCluster(clus)'
          works properly in your environment. 

 verbose: Logical value for message printing. 

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

     See vignette.

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

     Returns a 'twilight' object consisting of a 'data.frame' named
     'result' with variables 

  pvalue: Sorted input vector.

  qvalue: q-values computed as described in Storey and Tibshirani
          (2003) with new estimate 'pi0'.

   index: Index of the original ordering.

     fdr: Local false discovery rate averaged over 10 runs of SEP.

mean.fdr: Bootstrap estimate of local false discovery rate.

lower.fdr: Lower 'boot.ci'-bootstrap confidence bound.

upper.fdr: Upper 'boot.ci'-bootstrap confidence bound.

     Values are sorted by 'pvalue'.

_N_o_t_e:

     Additional output consists of

       'lambda'    Regularization parameter.
       'pi0'       SEP estimate of prior probability.
       'boot.pi0'  Bootstrap estimate and 'boot.ci'-bootstrap confidence bounds.
       'boot.ci'   Passes 'boot.ci' for plotting.
       'effect'    Histogram of effect size distributions averaged over 10 runs of SEP.

     If 'xin' is of class 'twilight', the remaining slots are filled
     with corresponding input values. If 'xin' is not of class
     'twilight', these slots remain free.

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

     Stefanie Scheid <URL: http://www.molgen.mpg.de/~scheid>

_R_e_f_e_r_e_n_c_e_s:

     Scheid S and Spang R (2004): A stochastic downhill search
     algorithm for estimating the local false discovery rate, _IEEE
     TCBB_ *1(3)*, 98-108.

     Storey JD and Tibshirani R (2003): Statistical significance for
     genomewide studies, _PNAS_ *100(16)*, 9440-9445.

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

     'twilight.pval', 'plot.twilight'

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

     data(exfdr) # twilight object created by twilight with B=1000 bootstraps
     plot(exfdr,which="fdr")

