Suppose I run the GC-test:
grangercausalitytests(np.vstack((df['target'], df['feature'])).T, maxlag=5)
I can pick the lag of the "feature" variable, which most likely Granger-causes the "target" variable.
- But what number of lags does the "target" variable have in this model?
- Further, how do I estimate this ADL model (some autoregressive lags + some lags of the independent variable)? I've seen somewhere, that ADL should be substituted with OLS/FGLS in Python, since there is no package for ADL. Yet I do not understand how to do that