mlm2lm                package:LMGene                R Documentation

_L_i_n_e_a_r _M_o_d_e_l _c_o_n_v_e_r_t_i_n_g _f_u_n_c_t_i_o_n

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

     This function rule out the specified 'lm' class data out of the
     given 'c("mlm", "lm")' class data.

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

     mlm2lm(lmobj, i)

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

   lmobj: An object of class 'c("mlm", "lm")' 

       i: A specific number that indicates a 'lm' in 'lmobj'  

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

     In case of multiple response from 'lm' function, this function can
     used.

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

 lmobj2 : Selected 'lm' class data

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

     David Rocke and Geun-Cheol Lee

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

     <URL: http://www.idav.ucdavis.edu/~dmrocke/>

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

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

     #library
     library(Biobase)
     library(LMGene)

     #data
     data(sample.eS)
     Smpd0 = sample.eS
     # model information 
     for(i in 1:length(Smpd0@phenoData@varLabels)){
       assign(paste('x', i, sep=''),as.factor(Smpd0@phenoData@pData[,i]))
     }
       
     fchar=''
     for(i in 1:length(Smpd0@phenoData@varLabels)){
       fchar=paste(fchar, paste('x', i, sep=''), ifelse(i<length(Smpd0@phenoData@varLabels), '+', ''), sep='')
     }
     fchar2 <- paste("y ~",fchar)
     #
     # run regression and anovas
     y <- t(as.matrix(Smpd0@exprs))
     formobj <- as.formula(fchar2)
     tmp <- lm(formobj)
     class(tmp)

     tmp2 <- mlm2lm(tmp,i)
     class(tmp2)

