BUG:
 * fixed default method to Newton: solves the problem of high iterations.
 * new value of z: allows estimate alternative specific variables + individual specific variables (new example 4).