norm1b                 package:nudge                 R Documentation

_F_u_n_c_t_i_o_n _f_o_r _n_o_r_m_a_l_i_z_i_n_g _t_h_e _m_e_a_n _a_n_d _v_a_r_i_a_n_c_e (_o_r _j_u_s_t _t_h_e _v_a_r_i_a_n_c_e) _o_f _s_i_n_g_l_e _r_e_p_l_i_c_a_t_e _l_o_g _r_a_t_i_o_s

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

     This performs a robust loess normalization of the variance of the
     log ratios in a single replicate experiment by regressing the
     absolute (mean normalized) log ratios on the log intensities and
     using the fitted values to scale the (mean normalized) log ratio
     for each gene.

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

     norm1b(logratio, logintensity, span1 = 0.6, span2 = 0.2, mean.norm=TRUE)

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

logratio: A vector or single-column matrix of log (base 2) ratios of
          gene expressions in two samples, if mean.norm is FALSE the
          log ratios should be already mean normalized.

logintensity: A vector or single-column matrix of log (base 2) total
          intensities (defined as the product) of gene expressions in
          the two samples.

   span1: Proportion of data used to fit the loess regression of the
          log ratios on the log intensities for the mean normalization.

   span2: Proportion of data used to fit the loess regression of the
          absolute (mean normalized) log ratios on the log intensities
          for the variance normalization.

mean.norm: A logical value indicating whether or not a mean
          normalization should be performed prior to the variance
          normalization.

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

     A vector or single-column matrix of mean and variance normalized
     log (base 2) ratios of gene expressions in two samples.

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

     N. Dean and A. E. Raftery

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

     N. Dean and A. E. Raftery (2005). Normal uniform mixture
     differential gene expression detection for cDNA microarrays.  BMC
     Bioinformatics. 6, 173-186.

     <URL: http://www.biomedcentral.com/1471-2105/6/173>

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

     'norm1a','norm1c','norm1d','norm2c','norm2d'

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

     data(like)
     lR<-log(like[,1],2)-log(like[,2],2)
     lI<-log(like[,1],2)+log(like[,2],2)

     lRnorm<-norm1b(lR,lI,mean.norm=TRUE)

