normRepLoess           package:maigesPack           R Documentation

_B_o_o_t_s_t_r_a_p _o_f _L_O_W_E_S_S _n_o_r_m_a_l_i_s_a_t_i_o_n

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

     This function normalises a microarray object re-doing the LOWESS
     fitting several times, selecting a pre-specified proportion of
     points aleatorily.

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

     normRepLoess(raw, span=0.4, propLoess=0.5, nRep=50, func="none",
                  bkgSub="none", ...)

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

     raw: an object of class 'maigesRaw' to be normalised.

    span: real number in (0,1) representing the proportion of points to
          use in the loess regression.

propLoess: real number in (0,1) representing the proportion of points
          (spots) to be used in each iteration of loess.

    nRep: number of repetitions for loess procedure.

    func: character string giving the function to estimate the final W
          value. You must use 'mean', 'median' or 'none' (default).

  bkgSub: character with background subtraction method, using the
          function 'backgroundcorrect' from _limma_ package.

     ...: additional parameters for function 'loessFit' from _limma_
          package.

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

     The LOWESS fitting for normalising microarray data is a
     computational intensive task, so pay attention to not specify a
     very large argument in 'nRep'. If you do so, your process will
     take so much time to conclude.

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

     The result of this function is an object of class 'maiges'.

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

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

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

     'loessFit'.

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

     ## Loading the dataset
     data(gastro)

     ## Doing the repetition loess with default parameters. Be carefull, this
     ## is very time consuming
     ## Not run: 
     gastro.norm = normRepLoess(gastro.raw2)
     ## End(Not run)

     ## Do the same normalization selecting 60% dos spots with 10
     ## repetitions and estimating the W by the mean value.
     ## Not run: 
     gastro.norm = normRepLoess(gastro.raw2, propLoess=0.6, nRep=10, func="mean")
     ## End(Not run)

