changed some modeling, added simulation estimation of coefficients