[major] changed utility model from U = beta * X to U = beta0 + beta * X;
updating algorithm for beta is changed: because beta0 and beta are highly correlated, we instead estimate beta and beta0/beta
case with two covariates (L=3), not complete yet
update comparison and inventory script to the new utility setup
[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
rename file
added case M2.L2
[EM algorithm] added first EM implementation for case M1.L1 (performance seems quite dependent on initial values)
remove unused code
[EM] allow the expectation of nopurchase to be non-integer
[EM] added EM code for Scenario 0 and 3M
[EM] add eps1.mu into the model
tried to estimate eps1.mu but does not seem to be identifiable; for now leave it there
change to X ~ N(2, 0.5); add M1.L2 case
calculate optimal inventory by quantiles
add profit simulation and efficient frontier
new unified MH scaling parameters; new traffic noise setting
[EM] EM code for the censored cases
updated inventory calculation and efficient frontier plot
[EM] added Scenario 0c, 1c, 3Mc for the case M1.L2
added EM models into the frontier comparison