title: Multivariate linear mixed model using glmer in R package lme4 - inconsistent error between updates
body: <p>I am trying to run a multivariate linear mixed model and need to use a remote workstation to reduce compute time.</p>
  
  <p>When I run glmer() from lme4 on my personal computer (R version 2.15.1, lme4 version 0.999999-0, 64 bit Unix) my model runs appropriately.</p>
  
  <p>When I switch to the remove workstation (R version 3.0.2, lme4 version 1.0-6, 64 bit Linux) my model does not run and I receive an error and a warning message:</p>
  
  <pre><code>Error in eval(expr, envir, enclos) : updateMu: Size mismatch
  In addition: Warning message:
    calling glmer() with family=gaussian (identity link) as a shortcut to lmer() is deprecated; please call lmer() directly
  </code></pre>
  
  <p>Interestingly, if I run lmer() on either of the dependent variables, the model will run without error or warning.</p>
  
  <p>Here is sample code:</p>
  
  <pre><code>library(lme4)
  
  DV.1 &lt;- rnorm(100,2,0.3)
  DV.2 &lt;- rnorm(100,5,2)
  FE.1 &lt;- sample(c(1:3),100,replace=TRUE)
  FE.2 &lt;- sample(c(1:2),100,replace=TRUE)
  RE.1 &lt;- rep(c(1:10),10)
  
  df &lt;- as.data.frame(cbind(DV.1,DV.2,FE.1,FE.2,RE.1))
  
  glmer(cbind(DV.1,DV.2)~FE.1+FE.2+(1|RE.1),family=gaussian,data=df)
  </code></pre>
  
  <p>Here is output from my personal computer:</p>
  
  <pre><code>&gt; glmer(cbind(DV.1,DV.2)~FE.1+FE.2+(1|RE.1),family=gaussian,data=df)
  Linear mixed model fit by REML 
  Formula: cbind(DV.1, DV.2) ~ FE.1 + FE.2 + (1 | RE.1) 
     Data: df 
     AIC   BIC logLik deviance REMLdev
   76.94 89.97 -33.47    53.66   66.94
  Random effects:
   Groups   Name        Variance Std.Dev.
   RE.1     (Intercept) 0.00000  0.00000 
   Residual             0.10322  0.32128 
  Number of obs: 100, groups: RE.1, 10
  
  Fixed effects:
              Estimate Std. Error t value
  (Intercept)  2.02330    0.13377  15.125
  FE.1        -0.04980    0.04082  -1.220
  FE.2         0.06238    0.06451   0.967
  
  Correlation of Fixed Effects:
       (Intr) FE.1  
  FE.1 -0.626       
  FE.2 -0.719 -0.037
  </code></pre>
  
  <p>Here is output from the remote workstation with a newer version of lme4:</p>
  
  <pre><code>&gt; glmer(cbind(DV.1,DV.2)~FE.1+FE.2+(1|RE.1),family=gaussian,data=df)
  Error in eval(expr, envir, enclos) : updateMu: Size mismatch
  In addition: Warning message:
  In glmer(cbind(DV.1, DV.2) ~ FE.1 + FE.2 + (1 | RE.1), family = gaussian,  :
    calling glmer() with family=gaussian (identity link) as a shortcut to lmer() is deprecated; please call lmer() directly
  &gt; 
  </code></pre>
  
  <p>I am wondering if the model is formed incorrectly, or if there is a complication with one version of lme4.
  Any thoughts or suggestions are greatly appreciated.
  Thanks!!</p>