EMupdateGets              package:lme4              R Documentation

_U_p_d_a_t_e _a_n _o_b_j_e_c_t _i_n _t_h_e _E_M _a_l_g_o_r_i_t_h_m

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

     A generic replacement function used in the EM and ECME
     optimization algorithms for linear mixed-effects models.

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

     EMupdate(x, nlev) <- value

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

       x: an `lmeStruct' or `pdMat' object

    nlev: integer: the number of levels of the grouping factor
          corresponding to the random-effects structure

   value: a matrix with the triangular factor from an
          orthogonal-triangular decomposition of the modal values of
          the random-effects and their precision matrices

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

     an object of the same class as `x'

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

     Douglas Bates bates@stat.wisc.edu and Saikat DebRoy
     saikat@stat.wisc.edu

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

     library(lme4)
     data(Oxboys, package = "nlme")
     m3 <- pdLogChol(~ age)
     m3
     as(m3, 'pdmatrix') <- crossprod(model.matrix(formula(m3), Oxboys))
     show(m3)
     EMupdate(m3, length(levels(Oxboys$Subject))) <- diag(2)
     show(m3)

