Is it possible to use the xwith command? Or would this lead to biased data due to fact that even when X and Z are normally distributed, X² will not. Thus, b5 will be overestimated because of the sensitivity of the LMS approach to nonnormal data. Or isn't this a problem, when I constrain the b-weights?