我正在尝试对数据使用statsmodels进行一些Newey-West OLS来估计我的参数,以下是我这样做的代码:
from __future__ import print_function, division
import xlrd as xl
import numpy as np
import scipy as sp
import pandas as pd
import statsmodels.formula.api as smf
import statsmodels.api as sm
file_loc = "/Python/dataset_3.xlsx"
workbook = xl.open_workbook(file_loc)
sheet = workbook.sheet_by_index(0)
tot = sheet.nrows
data = [[sheet.cell_value(r, c) for c in range(sheet.ncols)] for r in
range(sheet.nrows)]
rv1 = []
rv5 = []
rv22 = []
rv1fcast = []
T = []
price = []
time = []
retnor = []
for i in range(1, tot):
t = data[i][0]
ret = data[i][1]
ret5 = data[i][2]
ret22 = data[i][3]
ret1_1 = data[i][4]
retn = data[i][5]
t = xl.xldate_as_tuple(t, 0)
rv1.append(ret)
rv5.append(ret5)
rv22.append(ret22)
rv1fcast.append(ret1_1)
retnor.append(retn)
T.append(t)
df = pd.DataFrame({'RVFCAST':rv1fcast, 'RV1':rv1, 'RV5':rv5, 'RV22':rv22,})
df = df[df.RV1.notnull()]
model = smf.OLS(formula = 'df.RVFCAST ~ df.RV1 + df.RV5 + df.RV22', data = df)
当我查看数组或数据框时,一切看起来都很好,但返回的只是:TypeError:init()至少接受2个参数(给定1个)
我尝试了多种方法,但看不到缺少的内容。
当我运行它时,以下错误消息显示:
---------------------------------------------------------------------------
TypeError Traceback (most recent call last)
/Python/harrv.py in <module>()
41 df = df[df.RV1.notnull()]
42
---> 43 model = smf.OLS(formula = 'df.RVFCAST ~ df.RV1 + df.RV5 + df.RV22', data = df)
44
45 #mdl = model.get_robustcov_results(cov_type='HAC',maxlags=1)
TypeError: __init__() takes at least 2 arguments (1 given)
打印rv1可为您提供:
Out[318]:
[0.015538008996147568,
0.008881670570720125,
0.010421778063375802,
.....
0.003151044550868834,
0.0029676428110974166,
0.005236329928710288,
0.004838460533164701,
'']
而另一个rv给出了类似的浮点数。 df只是按照pd.dataframe的方式组装它们,根据文档的支持(http://statsmodels.sourceforge.net/devel/example_formulas.html)。
最佳答案
问题在于statsmodels.formula.api中的公式函数是小写字母。大写的OLS与主要的statsmodels.api中的相同。将来会从Formula.api命名空间中删除大写模型,从而避免这种混淆。
这意味着,您需要使用小写的ol,例如model = smf.ols(formula = 'df.RVFCAST ~ df.RV1 + df.RV5 + df.RV22', data = df)
请注意,小写的公式函数只是模型的from_formula
方法的别名。smf.ols
是sm.OLS.from_formula
的快捷方式
关于python - 多个OLS估计TypeError,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/29744909/