python - 具有外生变量矩阵的 statsmodels SARIMAX 大小不同

Question

我正在运行 SARIMAX 模型，但遇到了指定外生变量的问题。在第一段代码(下方)中，我指定了一个外生变量 lesdata['LESpost'] 并且模型运行没有问题。但是，当我添加另一个外生变量时，我最终得到一条错误消息(请参阅堆栈跟踪)。

ar = (1,0,1)      #  AR(1 3)
ma = (0)  #  No MA terms
mod1 = sm.tsa.statespace.SARIMAX(lesdata['emadm'], exog= (lesdata['LESpost'],lesdata['QOF']), trend='c', order=(ar,0,ma), mle_regression=True)

Traceback (most recent call last):

  File "<ipython-input-129-d1300aeaeffc>", line 4, in <module>
    mle_regression=True)

  File "C:\Users\danie\Anaconda2\lib\site-packages\statsmodels\tsa\statespace\sarimax.py", line 510, in __init__
    endog, exog=exog, k_states=k_states, k_posdef=k_posdef, **kwargs

  File "C:\Users\danie\Anaconda2\lib\site-packages\statsmodels\tsa\statespace\mlemodel.py", line 84, in __init__
    missing='none')

  File "C:\Users\danie\Anaconda2\lib\site-packages\statsmodels\tsa\base\tsa_model.py", line 43, in __init__
    super(TimeSeriesModel, self).__init__(endog, exog, missing=missing)

  File "C:\Users\danie\Anaconda2\lib\site-packages\statsmodels\base\model.py", line 212, in __init__
    super(LikelihoodModel, self).__init__(endog, exog, **kwargs)

  File "C:\Users\danie\Anaconda2\lib\site-packages\statsmodels\base\model.py", line 63, in __init__
    **kwargs)

  File "C:\Users\danie\Anaconda2\lib\site-packages\statsmodels\base\model.py", line 88, in _handle_data
    data = handle_data(endog, exog, missing, hasconst, **kwargs)

  File "C:\Users\danie\Anaconda2\lib\site-packages\statsmodels\base\data.py", line 630, in handle_data
    **kwargs)

  File "C:\Users\danie\Anaconda2\lib\site-packages\statsmodels\base\data.py", line 80, in __init__
    self._check_integrity()

  File "C:\Users\danie\Anaconda2\lib\site-packages\statsmodels\base\data.py", line 496, in _check_integrity
    super(PandasData, self)._check_integrity()

  File "C:\Users\danie\Anaconda2\lib\site-packages\statsmodels\base\data.py", line 403, in _check_integrity
    raise ValueError("endog and exog matrices are different sizes")

ValueError: endog and exog matrices are different sizes

我在这里明显遗漏了什么吗？变量的长度都是一样的，没有缺失数据。

感谢阅读，希望对您有所帮助!

Answer 1

应用 numpy.asarray 后，二维数据需要在行中有观察值，在列中有变量。

exog = (lesdata['LESpost'],lesdata['QOF'])

将 asarray 应用于此元组会将变量放在行中，这是来自 C 源的 numpy 默认值，这不是 statsmodels 想要的。

DataFrame 已经以适当的方式进行了整形，因此一种选择是使用具有所需列的 DataFrame

exog = lesdata[['LESpost', 'QOF']]

array_likes 列表或元组的另一种选择是使用 numpy.column_stack，例如

exog = np.column_stack((lesdata['LESpost'].values,lesdata['QOF'].values))