I am trying to build a model where both dependent and independent variables are time series some of which may not be stationary (this can apply to both types of variables).
My first port of call was (S)ARIMAX. But I have a few questions:
(1) Can ARIMAX capture the lagged impact of the independent time series?
I get the impression that in ARIMAX the independent time series is expected to have a simultaneous impact on the dependent one - which is not necessarily true in my case. How do I handle such a lagged impact?
Also, I am not really sure what are the lags. So my other question is:
(2) How can I estimate the order of such lagged impact?
Also, I am not sure what to do about the stationarity. It seems that either I need to treat (e.g. by differencing, power-transforming) all variables (even those which are stationary) or none at all. So my other question is:
(3) Do I need to treat the non-stationarity of some of the variables?
And my last question:
(4) Can you point to publicly accessible worked examples of such problems?
I would like to see some publicly available worked example of such problems (i.e. independent variables are time series with lagged impact on the dependent time series) with python or R code.
References to books on Amazon (or to some other non-open sources) won't be that helpful.