mrlm                  package:limma                  R Documentation

_F_i_t _L_i_n_e_a_r _M_o_d_e_l _t_o _M_i_c_r_o_r_r_a_y _D_a_t_a _b_y _R_o_b_u_s_t _R_e_g_r_e_s_s_i_o_n

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

     Fit a linear model genewise to expression data from a series of
     arrays. The fit is by robust M-estimation allowing for a small
     proportion of outliers. This is a utility function for 'lmFit'.

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

     mrlm(M,design=NULL,ndups=1,spacing=1,weights=NULL,...)

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

       M: numeric matrix containing log-ratio or log-expression values
          for a series of microarrays, rows correspond to genes and
          columns to arrays.

  design: numeric design matrix defining the linear model, with rows
          corresponding to arrays and columns to comparisons to be
          estimated. The number of rows must match the number of
          columns of 'M'. Defaults to the unit vector meaning that the
          arrays are treated as replicates.

   ndups: a positive integer giving the number of times each gene is
          printed on an array. 'nrow(M)' must be divisible by 'ndups'.

 spacing: the spacing between the rows of 'M' corresponding to
          duplicate spots, 'spacing=1' for consecutive spots.

 weights: numeric matrix of the same dimension as 'M' containing
          weights. If it is of different dimension to 'M', it will be
          filled out to the same size. 'NULL' is equivalent to equal
          weights.

     ...: any other arguments are passed to 'rlm.default'.

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

     This is a utility function used by the higher level function
     'lmFit'. Most users should not use this function directly but
     should use 'lmFit' instead.

     This function fits a linear model for each gene by calling the
     function 'rlm' from the MASS library.

     Warning: don't use weights with this function unless you
     understand how 'rlm' treats weights. The treatment of weights is
     somewhat different from that of 'lm.series' and 'gls.series'.

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

     A list with components 

coefficients: numeric matrix containing the estimated coefficients for
          each linear model. Same number of rows as 'M', same number of
          columns as 'design'.

stdev.unscaled: numeric matrix conformal with 'coef' containing the
          unscaled standard deviations for the coefficient estimators.
          The standard errors are given by 'stdev.unscaled * sigma'.

   sigma: numeric vector containing the residual standard deviation for
          each gene.

df.residual: numeric vector giving the degrees of freedom corresponding
          to 'sigma'.

      qr: QR decomposition of 'design'.

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

     Gordon Smyth

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

     'rlm'.

     An overview of linear model functions in limma is given by
     06.LinearModels.

