Dear professor:
        I apologize for disturbing you! I'm appreciating your package 'lme4()' in the R ,because it's very useful when we deal with the linear mixed model. However,I get struck in doing my Master's degree paper , I'm hoping that you can help me.
       Y = X*β + B*a + ε   is a linear mixed model , X : dim=n*k .  β: dim=k*1 represent fixed effects .B :dim=n*(n-2) .  a: dim=p*1 represent random effects with a ∼ N(0, τI). For example, n=100,k=15 ,X and B is known,how can I use lmer() to  get  β and τ .How to set parameters in lmer()?
      Thank you very much! Wish you a success in your work!  
	[[alternative HTML version deleted]]

