selSpots             package:maigesPack             R Documentation

_S_e_l_e_c_t _s_p_o_t_s _t_o _u_s_e _i_n _n_o_r_m_a_l_i_s_a_t_i_o_n

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

     Function to select spots to be used in microarray normalisation.

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

     selSpots(obj=NULL, sigNoise=1, rmFlag=NULL, gLabelsID=c("Name"),
              remove=list(c("BLANK","DAP","LYS","PHE","Q_GENE","THR","TRP")),
              badSpots=NULL, badLabel=NULL)

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

     obj: object of class 'maigesRaw'.

sigNoise: positive real number indicating the cutoff to remove spots
          with signal-to-noise ratio below it.

  rmFlag: vector of flag symbols to be removed (for normalisation).
          These flags are stored in the slot 'Flags' from 'maigesRaw'
          class.

gLabelsID: character vector indicating the gene labels to be searched
          to exclude that ones specified in 'remove' argument.

  remove: list of same length as 'GlabelsID' containing character
          vector indicating the symbols of spots to be removed,
          according to the 'GlabelsID' argument.

badSpots: index of bad spots (numeric or logical) identifying bad
          spots. May be the gene labels, with label ID specified by the
          argument 'badLabel'.

badLabel: character string specifying the gene label ID for remove
          'badSpots'.

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

     This function takes the object of class 'maigesRaw' and actualise
     the slot 'UseSpots' according with the arguments passed to the
     function. This slot is read by the normalisation functions
     'normLoc', 'normOLIN', 'normRepLoess', 'normScaleLimma' and
     'normScaleMarray' to use only the spots that passed the criteria
     specified here.

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

     This function returns another object of class 'maigesRaw' with the
     'UseSpots' slot actualised.

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

     Gustavo H. Esteves <gesteves@vision.ime.usp.br>

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

     ## Loading the dataset
     data(gastro)

     ## Filtering all spots with signal2noise ratio (Sf/Sb or Rf/Rb) greater
     ## or equal to 1 and that have 'Name' label as 'BLANK', 'DAP', ..., 'TRP'.
     gastro.raw2 = selSpots(gastro.raw, sigNoise=1, rmFlag=NULL, gLabelsID="Name",
       remove=list(c("BLANK","DAP","LYS","PHE","Q_GENE","THR","TRP")))

     ## To see the number of spots that suvived the filtering above do
     apply(gastro.raw2@UseSpots, 2, sum)

     ## To do the same filtering as above, also filtering flags marcke as 1
     ## and 4 do
     gastro.raw2 = selSpots(gastro.raw, sigNoise=1, rmFlag=c(1,4), gLabelsID="Name",
       remove=list(c("BLANK","DAP","LYS","PHE","Q_GENE","THR","TRP")))

     apply(gastro.raw2@UseSpots, 2, sum)

