profile.merMod() now allows  which = "beta_" or "theta_"