trying mlogit for psi instead of logit link function, suggested by DNK.