kooperberg               package:limma               R Documentation

_K_o_o_p_e_r_b_e_r_g _M_o_d_e_l-_B_a_s_e_d _B_a_c_k_g_r_o_u_n_d _C_o_r_r_e_c_t_i_o_n

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

     This function uses a Bayesian model to background correct data
     from a series of microarray experiments. It currently works only
     with GenePix data.

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

     kooperberg(names, fg="mean", bg="median", a=TRUE, layout)

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

   names: character vector giving the names of data.frames containing
          GenePix data

      fg: character string giving foreground estimator. Choices are
          '"mean"' or '"median"'.

      bg: character string giving foreground estimator. Choices are
          '"mean"' or '"median"'.

       a: logical.  If 'TRUE', the 'a' parameters in the model
          (equation 3 and 4) are estimated for each slide.  If 'FALSE'
          the 'a' parameters are set to unity.

  layout: list containing print layout with components 'ngrid.r',
          'ngrid.c', 'nspot.r' and 'nspot.c'

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

     This function is for use with GenePix data and is designed to cope
     with the problem of large numbers of negative intensities and
     hence missing values on the log-intensity scale. It avoids missing
     values in most cases and at the same time dampens down the
     variability of log-ratios for low intensity spots. See Kooperberg
     et al (2003) for more details.

     'kooperberg' serially extracts the foreground and background
     intensities, standard deviations and number of pixels from GenePix
     data frames. This information is used to compute empirical
     estimates of the model parameters as described in equation 2 of
     Kooperberg et al (2003).

     The foreground and background estimates extracted from the GenePix
     files may be based on means or medians of pixel values. Setting
     'fg="mean"' uses the GenePix column 'F635.Mean' for red foreground
     and the GenePix column 'F532.Mean' for green foreground. Setting
     'fg="median"' uses columns 'F635.Median' and 'F532.Median'.
     Similarly for the background, 'bg="mean"' uses columns 'B635.Mean'
     and 'B532.Mean' while 'bg="median"' uses columns 'B635.Median' and
     'B532.Median'.

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

     An 'RGList' containing the components 

       R: matrix containing the background adjusted intensities for the
          red channel for each spot for each array

       G: matrix containing the background adjusted intensities for the
          green channel for each spot for each array

 printer: list containing printer layout, as input in 'layout'

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

     Matthew Ritchie

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

     Kooperberg, C., Fazzio, T. G., Delrow, J. J., and Tsukiyama, T.
     (2002) Improved background correction for spotted DNA microarrays.
     _Journal of Computational Biology_ *9*, 55-66.

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

     4.Normalization gives an overview of normalization and background
     correction functions defined in the LIMMA package.

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

     #  This is example code for reading and background correcting GenePix data
     #  given GenePix Results (gpr) files in the working directory (data not
     #  provided).
     ## Not run: 
     genepixFiles <- dir(pattern="*\\.gpr$") # get the names of the GenePix image analysis output files in the current directory
     read.series(genepixFiles, suffix=NULL, skip=26, sep="\t") # read the files into data.frames, assuming headers of 26 lines
     layout <- getLayout(get(genepixFiles[1]))
     RGmodel <- kooperberg(genepixFiles, layout=layout) # model-based background correction
     MA <- normalizeWithinArrays(RGModel) # normalize the data
     ## End(Not run)

