[major] allow for choice-specific constant terms in utility
Utility[m] = Beta.const[m] + Beta.coef %*% X[m]
Beta.const is M-dimensional, Beta.coef is L-dimensional, so a total of M+L unknown parameters to estimate