rcModels           package:preprocessCore           R Documentation

_F_i_t _r_o_w-_c_o_l_u_m_n _m_o_d_e_l _t_o _a _m_a_t_r_i_x

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

     These functions fit row-column effect models to matrices

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

     rcModelPLM(y,row.effects=NULL)
     rcModelWPLM(y, w,row.effects=NULL)
     rcModelMedianPolish(y)

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

       y: A numeric matrix

       w: A matrix or vector of weights. These should be non-negative.

row.effects: If these are supplied then the fitting procedure uses
          these (and analyzes individual columns separately)

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

     These functions fit row-column models to the specified input
     matrix. Specifically the model 

                       y_ij = r_i + c_j + e_ij

     with r_i and c_j as row and column effects respectively. Note that
     this functions treat the row effect as the parameter to be
     constrained using sum to zero (for 'rcModelPLM' and 'rcModelWPLM')
     or median of zero (for 'rcModelMedianPolish').

     The 'rcModelPLM' and 'rcModelWPLM' functions use a robust linear
     model procedure for fitting the model.

     The function 'rcModelMedianPolish' uses the median polish
     algorithm.

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

     A list with following items: 

Estimates: The parameter estimates. Stored in column effect then row
          effect order

 Weights: The final weights used

Residuals: The residuals

StdErrors: Standard error estimates. Stored in column effect then row
          effect order

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

     B. M. Bolstad bmb@bmbolstad.com

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

     col.effects <- c(10,11,10.5,12,9.5)
     row.effects <- c(seq(-0.5,-0.1,by=0.1),seq(0.1,0.5,by=0.1))

     y <- outer(row.effects, col.effects,"+")
     w <- runif(50)

     rcModelPLM(y)
     rcModelWPLM(y, w)
     rcModelMedianPolish(y)

     y <- y + rnorm(50)

     rcModelPLM(y)
     rcModelWPLM(y, w)
     rcModelMedianPolish(y)


     rcModelPLM(y,row.effects=row.effects)
     rcModelWPLM(y,w,row.effects=row.effects)

