How to add random effect into the model