Question

我想使用scikit-learn的GridSearchCV来确定随机森林模型的一些超参数。我的数据是时间相关的，看起来像

import pandas as pd

train = pd.DataFrame({'date': pd.DatetimeIndex(['2012-1-1', '2012-9-30', '2013-4-3', '2014-8-16', '2015-3-20', '2015-6-30']), 
'feature1': [1.2, 3.3, 2.7, 4.0, 8.2, 6.5],
'feature2': [4, 4, 10, 3, 10, 9],
'target': [1,2,1,3,2,2]})

>>> train
        date  feature1  feature2  target
0 2012-01-01       1.2         4       1
1 2012-09-30       3.3         4       2
2 2013-04-03       2.7        10       1
3 2014-08-16       4.0         3       3
4 2015-03-20       8.2        10       2
5 2015-06-30       6.5         9       2

如何实现以下交叉验证折叠技术？

train:(2012, 2013) - test:(2014)
train:(2013, 2014) - test:(2015)

也就是说，我想用2年的历史观察来训练模型，然后在接下来的一年里测试它的准确性。

Answer 1

您只需将带有拆分的iterable传递给GridSearchCV。此拆分应采用以下格式：

[
 (split1_train_idxs, split1_test_idxs),
 (split2_train_idxs, split2_test_idxs),
 (split3_train_idxs, split3_test_idxs),
 ...
]

要获得idx，您可以执行以下操作：

groups = df.groupby(df.date.dt.year).groups
# {2012: [0, 1], 2013: [2], 2014: [3], 2015: [4, 5]}
sorted_groups = [value for (key, value) in sorted(groups.items())] 
# [[0, 1], [2], [3], [4, 5]]

cv = [(sorted_groups[i] + sorted_groups[i+1], sorted_groups[i+2])
      for i in range(len(sorted_groups)-2)]

这看起来像这样：

[([0, 1, 2], [3]),  # idxs of first split as (train, test) tuple
 ([2, 3], [4, 5])]  # idxs of second split as (train, test) tuple

然后你可以这样做：

GridSearchCV(estimator, param_grid, cv=cv, ...)

Answer 2

sklearn中还有一个TimeSeriesSplit函数，该函数在训练/测试集中拆分时间序列数据（即以固定的时间间隔）。请注意，与标准的交叉验证方法不同，连续的训练集是它们之前的训练集的超集，即在每个分组中，测试索引必须比以前更高，因此在交叉验证器中改组是不合适的。

Answer 3

使用GroupShuffleSplit标准sklearn方法。来自文档：

根据第三方提供的组提供随机训练/测试索引以分割数据。该组信息可用于将样本的任意域特定分层编码为整数。

例如，这些组可能是样本收集年份，因此可以对基于时间的分割进行交叉验证。

非常方便您的使用案例。这是它的样子：

cv = GroupShuffleSplit().split(X, y, groups)

并像以前一样将其传递给GridSearchCV：

GridSearchCV(estimator, param_grid, cv=cv, ...)

scikit-learn交叉验证时间序列数据的自定义拆分

3 个答案: