I want to perform the Toda-Yamamoto procedure on a financial timeseries. I looked for python packages that already implement it, but couldn’t find any. Does anyone have a suggestion? Or do people tend to implement it themselves using some existing VAR implementation?