目标编码器的内部交叉拟合#

TargetEncoder 将类别特征的每个类别替换为该类别目标变量的收缩均值。当类别特征与目标之间存在强相关性时,此方法非常有用。为了防止过拟合,TargetEncoder.fit_transform 使用内部交叉拟合方案对训练数据进行编码,以供下游模型使用。该方案涉及将数据分成 k 折,并使用从其他 k-1 折学习到的编码对每一折进行编码。在本示例中,我们演示了交叉拟合过程在防止过拟合方面的重要性。

# Authors: The scikit-learn developers
# SPDX-License-Identifier: BSD-3-Clause

创建合成数据集#

对于本示例,我们构建了一个包含三个类别特征的数据集

  • 一个具有中等基数的信息性特征(“informative”)

  • 一个具有中等基数的非信息性特征(“shuffled”)

  • 一个具有高基数的非信息性特征(“near_unique”)

首先,我们生成信息性特征

import numpy as np

from sklearn.preprocessing import KBinsDiscretizer

n_samples = 50_000

rng = np.random.RandomState(42)
y = rng.randn(n_samples)
noise = 0.5 * rng.randn(n_samples)
n_categories = 100

kbins = KBinsDiscretizer(
    n_bins=n_categories,
    encode="ordinal",
    strategy="uniform",
    random_state=rng,
    subsample=None,
)
X_informative = kbins.fit_transform((y + noise).reshape(-1, 1))

# Remove the linear relationship between y and the bin index by permuting the
# values of X_informative:
permuted_categories = rng.permutation(n_categories)
X_informative = permuted_categories[X_informative.astype(np.int32)]

具有中等基数的非信息性特征是通过置换信息性特征并消除与目标的关系而生成的

X_shuffled = rng.permutation(X_informative)

具有高基数的非信息性特征的生成使其独立于目标变量。我们将展示,在没有交叉拟合的情况下进行目标编码会导致下游回归器出现灾难性的过拟合。这些高基数特征基本上是样本的唯一标识符,通常应从机器学习数据集中删除。在本示例中,我们生成它们是为了展示 TargetEncoder 的默认交叉拟合行为如何自动缓解过拟合问题。

X_near_unique_categories = rng.choice(
    int(0.9 * n_samples), size=n_samples, replace=True
).reshape(-1, 1)

最后,我们组装数据集并执行训练-测试集划分

import pandas as pd

from sklearn.model_selection import train_test_split

X = pd.DataFrame(
    np.concatenate(
        [X_informative, X_shuffled, X_near_unique_categories],
        axis=1,
    ),
    columns=["informative", "shuffled", "near_unique"],
)
X_train, X_test, y_train, y_test = train_test_split(X, y, random_state=0)

训练 Ridge 回归器#

在本节中,我们将训练一个使用和不使用编码的 Ridge 回归器,并探讨目标编码器在有无内部交叉拟合情况下的影响。首先,我们看到在原始特征上训练的 Ridge 模型性能较低。这是因为我们打乱了信息性特征的顺序,这意味着 X_informative 在原始状态下不具有信息性

import sklearn
from sklearn.linear_model import Ridge

# Configure transformers to always output DataFrames
sklearn.set_config(transform_output="pandas")

ridge = Ridge(alpha=1e-6, solver="lsqr", fit_intercept=False)

raw_model = ridge.fit(X_train, y_train)
print("Raw Model score on training set: ", raw_model.score(X_train, y_train))
print("Raw Model score on test set: ", raw_model.score(X_test, y_test))
Raw Model score on training set:  0.0049896314219659565
Raw Model score on test set:  0.004577621581492997

接下来,我们使用目标编码器和 Ridge 模型创建一个流水线。该流水线使用 TargetEncoder.fit_transform,它采用了交叉拟合。我们看到模型很好地拟合了数据并泛化到测试集

from sklearn.pipeline import make_pipeline
from sklearn.preprocessing import TargetEncoder

model_with_cf = make_pipeline(TargetEncoder(random_state=0), ridge)
model_with_cf.fit(X_train, y_train)
print("Model with CF on train set: ", model_with_cf.score(X_train, y_train))
print("Model with CF on test set: ", model_with_cf.score(X_test, y_test))
Model with CF on train set:  0.8000184677460305
Model with CF on test set:  0.7927845601690917

线性模型的系数表明,大部分权重都在列索引为 0 的特征上,这是信息性特征

import matplotlib.pyplot as plt
import pandas as pd

plt.rcParams["figure.constrained_layout.use"] = True

coefs_cf = pd.Series(
    model_with_cf[-1].coef_, index=model_with_cf[-1].feature_names_in_
).sort_values()
ax = coefs_cf.plot(kind="barh")
_ = ax.set(
    title="Target encoded with cross fitting",
    xlabel="Ridge coefficient",
    ylabel="Feature",
)
Target encoded with cross fitting

尽管 TargetEncoder.fit_transform 使用内部交叉拟合方案来学习训练集的编码,但 TargetEncoder.transform 本身不使用。它使用完整的训练集来学习编码并转换类别特征。因此,我们可以使用 TargetEncoder.fit,然后是 TargetEncoder.transform 来禁用交叉拟合。然后将此编码传递给 Ridge 模型。

target_encoder = TargetEncoder(random_state=0)
target_encoder.fit(X_train, y_train)
X_train_no_cf_encoding = target_encoder.transform(X_train)
X_test_no_cf_encoding = target_encoder.transform(X_test)

model_no_cf = ridge.fit(X_train_no_cf_encoding, y_train)

我们评估了在编码时未使用交叉拟合的模型,发现它过拟合了。

print(
    "Model without CF on training set: ",
    model_no_cf.score(X_train_no_cf_encoding, y_train),
)
print(
    "Model without CF on test set: ",
    model_no_cf.score(
        X_test_no_cf_encoding,
        y_test,
    ),
)
Model without CF on training set:  0.858486250088675
Model without CF on test set:  0.6338211367102258

Ridge 模型过拟合是因为它对非信息性的极高基数(“near_unique”)和中等基数(“shuffled”)特征赋予了比模型使用交叉拟合编码特征时更多的权重。

coefs_no_cf = pd.Series(
    model_no_cf.coef_, index=model_no_cf.feature_names_in_
).sort_values()
ax = coefs_no_cf.plot(kind="barh")
_ = ax.set(
    title="Target encoded without cross fitting",
    xlabel="Ridge coefficient",
    ylabel="Feature",
)
Target encoded without cross fitting

结论#

本示例演示了 TargetEncoder 内部交叉拟合的重要性。在将训练数据传递给机器学习模型之前,使用 TargetEncoder.fit_transform 对其进行编码非常重要。当 TargetEncoderPipeline 的一部分并且该流水线被拟合时,流水线将正确调用 TargetEncoder.fit_transform 并在编码训练数据时使用交叉拟合

脚本总运行时间: (0 分钟 0.285 秒)

相关示例

比较目标编码器与其他编码器

比较目标编码器与其他编码器

强异常值数据集上的 HuberRegressor 与 Ridge 比较

强异常值数据集上的 HuberRegressor 与 Ridge 比较

普通最小二乘和岭回归

普通最小二乘和岭回归

绘制岭系数作为正则化的函数

绘制岭系数作为正则化的函数

由 Sphinx-Gallery 生成的画廊