randomEffects function

Random effects estimation