在Python中使用StatsModels运行GLM

Question

我尝试使用自己的数据复制文档中的示例：

>>> import numpy as np
>>> import statsmodels.api as sm
>>> data = sm.datasets.longley.load()
>>> data.exog = sm.add_constant(data.exog)
>>> ols_resid = sm.OLS(data.endog, data.exog).fit().resid
>>> res_fit = sm.OLS(ols_resid[1:], ols_resid[:-1]).fit()
>>> rho = res_fit.params
>>> from scipy.linalg import toeplitz
>>> order = toeplitz(np.arange(16))
>>> sigma = rho**order
>>> gls_model = sm.GLS(data.endog, data.exog, sigma=sigma)
>>> gls_results = gls_model.fit()
>>> print(gls_results.summary())

我的数据：

^{pr2}$

我不明白错误信息。在我看来，我忠实地复制了文档中的示例。在

Answer 1

设计矩阵exog有dtype=object，没有为其定义数值运算。在

你需要把exog转换成float

exog = exog.astype(np.float64)