I am trying to simulate data from a liner mixed model in R using the sim() function from the arm packag