【发布时间】:2016-01-04 12:23:45
【问题描述】:
我正在尝试同时使用 SVR 和 GP 来预测时间序列。
时间序列实际上是一个以pandas.DatetimeIndex 为索引的pandas.Series。
这两种算法都在 scikit-learn 中实现。而 SVR 在修改后实际上接受使用 X 轴进行预测:
train_size = 100
svr = GridSearchCV(SVR(kernel='rbf', gamma=0.001), cv=5,
param_grid={"C": [1e0, 1e1, 1e2, 1e3],
"gamma": np.logspace(-2, 2, 5)})
# Chose train_size items randomly in the data to predict
data_train = ts.sample(n=train_size)
data_train_y = data_train.values
data_train_x = (data_train.index.values).reshape(train_size, 1)
t0 = time.time()
svr.fit(data_train_x, data_train_y)
svr_fit = time.time() - t0
print("SVR complexity and bandwidth selected and model fitted in %.3f s"
% svr_fit)
sv_ratio = svr.best_estimator_.support_.shape[0] / train_size
print("Support vector ratio: %.3f" % sv_ratio)
min_date = min(ts.index.values)
max_date = max(ts.index.values) + np.timedelta64('1','D')
X_plot = pd.date_range(min_date, max_date, freq='H')
X_plot = X_plot.to_series()
X_plot = X_plot.reshape([len(X_plot),1])
t0 = time.time()
y_svr = svr.predict(X_plot)
svr_predict = time.time() - t0
print("SVR prediction for %d inputs in %.3f s"
% (X_plot.shape[0], svr_predict))
高斯过程不想处理:
from sklearn import gaussian_process
gp = gaussian_process.GaussianProcess(theta0=1e-2, thetaL=1e-4, thetaU=1e-1)
gp.fit(data_train_x, data_train_y)
它导致:
TypeError: ufunc add cannot use operands with types dtype('<M8[ns]') and dtype('<M8[ns]')
【问题讨论】:
标签: python pandas scikit-learn time-series gaussian