normalizeByReference       package:tilingArray       R Documentation

_P_r_o_b_e-_s_p_e_c_i_f_i_c _n_o_r_m_a_l_i_z_a_t_i_o_n _o_f _h_y_b_r_i_d_i_z_a_t_i_o_n _i_n_t_e_n_s_i_t_i_e_s _f_r_o_m _a_n _o_l_i_g_o_n_u_c_l_e_o_t_i_d_e _m_i_c_r_o_a_r_r_a_y

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

     Adjust the hybridization intensities from an oligonucleotide
     microarray for probe-specific response effect by using one or
     several reference hybridizations.

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

     normalizeByReference(x, reference, pm, background, refSig, nrStrata=10,
       cutoffQuantile=0.05, plotFileNames, verbose=FALSE)

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

       x: ExpressionSet containing the data to be normalized.

reference: ExpressionSet with the same number of features as 'x',
          containing the reference signal, on the raw scale
          (non-logarithmic). This argument can be used to directly
          input the data from a set of replicate DNA hybridizations.
          Alternatively, the argument 'refSig' can be specified.

      pm: Indices specifying the perfect match features in 'reference'
          (see Details). This can be either an integer vector with
          values between 1 and 'nrow(exprs(reference))' or a logical
          vector.

background: Indices specifying a set of background features in 'x' (see
          Details). This can be either an integer vector with values
          between 1 and 'nrow(exprs(x))' or a logical vector.

  refSig: A numeric vector of the same length as 'pm' with estimates of
          probe response effects, on a logarithm-like scale. This
          argument can be specified alternatively to 'reference'.

nrStrata: Integer (length 1), number of strata for the estimation of
          the background function.

cutoffQuantile: Numeric (length 1), the probes whose reference signal
          is below this quantile are thrown out.

plotFileNames: Character vector whose length is the same as the number
          of arrays in 'x'. Optional, if missing, no plots are
          produced.

 verbose: Logical of length 1, if 'TRUE', some messages about progress
          are printed.

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

     The intensities in 'x' are adjusted according to the reference
     values. Typically, the reference values are obtained by
     hybridizing a DNA sample to the array, so that the abundance of
     target is the same for all reference probes, and their signal can
     be used to estimate the probe sequence effect. A reference probe
     is a probe that perfectly matches the target genome exactly once.
     Usually, not all probes on a chip are reference probes, hence the
     subset of those that are is specified by the argument 'pm'.

     The background signal is estimated from the probes indicated by
     the argument 'background'. They need to be a strict subset of the
     'reference' probes. I.e., they need to uniquely match the target
     organism's DNA, but are not expected to match any of its
     transcripts. A robust estimation method is used, so a small
     fraction of 'background' probes that do hit transcripts is not
     harmful.

     A limitation of this normalization method is that it only makes
     sense for the data from reference probes, NA values are returned
     for all other probes.

     The functions 'PMindex' and 'BGindex' can be used to produce the
     'pm' and 'background' arguments from a 'probeAnno' environment
     such as provided in the davidTiling package.

     To summarize, a reference probe (indicated by argument 'pm') is a
     probe that perfectly matches the target genome exactly once, a
     background probe (indicated by argument 'background') is a
     reference probe which we expect not to be transcribed. These
     should not be confused with what is called 'perfect match' and
     'mismatch' probes in Affymetrix annotation.

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

     A copy of 'x' with the normalized intensities.

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

     W. Huber huber@ebi.ac.uk

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

     Huber W, Toedling J, Steinmetz, L. Transcript mapping with
     high-density  oligonucleotide tiling arrays. Bioinformatics 22,
     1963-1970 (2006).

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

     'PMindex', 'BGindex'

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

       ## see vignette assessNorm.Rnw in inst/scripts directory

