I tried to run the above. It seems statsmodels in Python now support it. See the below two links:
I still don’t quite get it. Anyone here can provide me a simplified illustration? I think I might have messed up the datatype for the endog and exog. I usually use scikit learn format.
Any help is appreciated.