logitT                  package:plw                  R Documentation

_l_o_g_i_t-_t _a_n_d _t-_t_e_s_t _b_y _r_o_w

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

     Functions for the logit-t test (Lemon et al. 2003) and the
     ordinary t-test computed for each row of an matrix.

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

     logitTTransform(pm)
     logitTStat(affy.batch,group)
     studenttTTest(x, group)

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

      pm: A matrix of Pm intensities

affy.batch: An AffyBatch object

   group: A group indicator vector, should have values 1 and 2 only.

       x: A matrix

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

     See the definition (R-code) of each function for details.

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

     logitTTransform returns a matrix

     logitTStat returns a vector with the logit-t statistic for each
     probe set.

     studenttTTest returns a vector with t-statistic for each row of x.

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

     Magnus Astrand

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

     Lemon et al. (2003). A high performance test of differential gene
     expression for oligonucleotide arrays. Genome Biol. 2003;
     4(10):R67

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

     # ------------------------------------------
     # Example analyzing the 6 arrays in the 
     # AffySpikeU95Subset data set

     # Loading the data
     data(AffySpikeU95Subset)

     # Vector with groups assignment 
     group<-factor(rep(1:2,each=3))

     # logit-T statistic
     logitT<-logitTStat(AffySpikeU95Subset,
                        as.numeric(group))

     # Computing RMA expression index
     data.rma<-exprs(rma(AffySpikeU95Subset))

     # Ordinary t-test by row/gene
     studentT<-studenttTTest(data.rma, as.numeric(group))

     # Comparing genes ranked top-20
     logitTTop20  <- rank(-abs(logitT)) < 21
     studentTTop20<- rank(-abs(studentT)) < 21
     table(logitTTop20,studentTTop20)

