1. My example in post #379 above, we have the special situation that the ec1 form is overparameterized because the underlying model is an ARDL(1,0,2) model with 0 lags for the ln_inc regressor. By forcing ln_inc to enter with its first lag in the long-run relationship, we have to artifically construct a corresponding short-run term that would not be present in the ec form. The coefficient of this short-run term is not freely estimated but determined by a nonlinear constraint. This overparameterization is not preserved when the results are stored with option regstore() because the covariance matrix of the coefficients would be singular. As you can see in the output, the regressor ln_consump (which has a positive number of lags) still appears with its first lag in the newey output when you use ardl with the ec1 option. You can avoid the problem for the ln_inc regressor by forcing it to have exactly 1 lag with the lags(. 1 .) option of the ardl command. Alternatively, there is the undocumented minlag1 option that enforces a minimum of 1 lag for each regressor.
2. The parameterization of the long-run coefficients cannot be preserved for the newey command because it is a nonlinear transformation of the underlying ARDL coefficients. That remark holds both for the ec and the ec1 option. To obtain the Newey West standard errors for the long-run coefficients, you would need to use the nlcom command directly after the newey command, as on slide 42 of my 2018 London Stata Conference presentation:
Code:
nlcom (_b[ln_inc] / (1 - _b[L.ln_inv])) (_b[L.ln_consump] / (1 - _b[L.ln_inv]))
Leave a comment: