Mixed-effects models in S and S-PLUS