I have a model consisting 4 variables, two of which are stationary at level and two that are I1 stationary. There seems to be a cointegration relationship between the two I1 variables after running the Johansen test.
Do I use a VAR model but taking the first differences of these I1 variables and ignore cointegration, or I proceed with a VECM model? If is the latter how do I proceed?