十四、K 最近邻

作者:Chris Albon

译者:飞龙

协议:CC BY-NC-SA 4.0

确定 K 的最佳值

十四、K 最近邻 - 图1

  1. # 加载库
  2. from sklearn.neighbors import KNeighborsClassifier
  3. from sklearn import datasets
  4. from sklearn.preprocessing import StandardScaler
  5. from sklearn.pipeline import Pipeline, FeatureUnion
  6. from sklearn.model_selection import GridSearchCV
  7. # 加载数据
  8. iris = datasets.load_iris()
  9. X = iris.data
  10. y = iris.target
  11. # 创建标准化器
  12. standardizer = StandardScaler()
  13. # 标准化特征
  14. X_std = standardizer.fit_transform(X)
  15. # 拟合 5 个邻居的 KNN 分类器
  16. knn = KNeighborsClassifier(n_neighbors=5, metric='euclidean', n_jobs=-1).fit(X_std, y)
  17. # 创建流水线
  18. pipe = Pipeline([('standardizer', standardizer), ('knn', knn)])
  19. # 创建候选值空间
  20. search_space = [{'knn__n_neighbors': [1, 2, 3, 4, 5, 6, 7, 8, 9, 10]}]
  21. # 创建网格搜索
  22. clf = GridSearchCV(pipe, search_space, cv=5, verbose=0).fit(X_std, y)
  23. # 最佳邻居大小(K)
  24. clf.best_estimator_.get_params()['knn__n_neighbors']
  25. # 6

KNN 分类

K 最近邻分类器(KNN)是一种简单而强大的分类学习器。

KNN 有三个基本部分

  • 十四、K 最近邻 - 图2: 观测的类别(我们试图在测试数据中预测的东西)。
  • 十四、K 最近邻 - 图3: 观察的预测因子/ IV /属性。
  • 十四、K 最近邻 - 图4: 研究者指定的正数。 K 表示最接近特定观测的观测数,它定义了“邻域”。 例如,K = 2意味着每个观测都有一个邻域,包含最接近它的另外两个观测。

想象一下,我们有一个观测,我们知道它的自变量 十四、K 最近邻 - 图5,但不知道它的类别 十四、K 最近邻 - 图6。 KNN 学习器找到最接近 十四、K 最近邻 - 图7 的K个其他观测,并使用他们已知的类别,将类别分配给 十四、K 最近邻 - 图8

  1. import pandas as pd
  2. from sklearn import neighbors
  3. import numpy as np
  4. %matplotlib inline
  5. import seaborn

这里我们创建三个变量,test_1test_2是我们的自变量,outcome是我们的因变量。 我们将使用这些数据来训练我们的学习器。

  1. training_data = pd.DataFrame()
  2. training_data['test_1'] = [0.3051,0.4949,0.6974,0.3769,0.2231,0.341,0.4436,0.5897,0.6308,0.5]
  3. training_data['test_2'] = [0.5846,0.2654,0.2615,0.4538,0.4615,0.8308,0.4962,0.3269,0.5346,0.6731]
  4. training_data['outcome'] = ['win','win','win','win','win','loss','loss','loss','loss','loss']
  5. training_data.head()
test_1test_2outcome
00.30510.5846win
10.49490.2654win
20.69740.2615win
30.37690.4538win
40.22310.4615win

这不是必需的,但因为我们只有三个变量,所以我们可以绘制训练数据集。 X 轴和 Y 轴是自变量,而点的颜色是它们的类别。

  1. seaborn.lmplot('test_1', 'test_2', data=training_data, fit_reg=False,hue="outcome", scatter_kws={"marker": "D","s": 100})
  2. # <seaborn.axisgrid.FacetGrid at 0x11008aeb8>

png

scikit-learn库需要将数据格式化为numpy数组。 这是重新格式化的代码。

  1. X = training_data.as_matrix(columns=['test_1', 'test_2'])
  2. y = np.array(training_data['outcome'])

这是我们的重点。 我们使用“观测的邻域是其三个最近的邻居”的参数来训练 KNN 学习器。 weights ='uniform'可以当做所用的投票系统。 例如,uniform意味着所有邻居对观测的类别进行同等权重的“投票”,而weight ='distance'则告诉学习器根据到我们正在分类的观测的距离,来调整每个观测的“投票”。

  1. clf = neighbors.KNeighborsClassifier(3, weights = 'uniform')
  2. trained_model = clf.fit(X, y)

与训练数据相比,我们训练的模型有多好?

  1. trained_model.score(X, y)
  2. # 0.80000000000000004

我们的模型准确率达 80%!

注:在任何现实世界的例子中,我们都希望将训练的模型与一些保留的测试数据进行比较。 但由于这是一个玩具示例,我使用了训练数据。

现在我们已经训练了我们的模型,我们可以预测班级的任何新观测,十四、K 最近邻 - 图10。 我们现在就这样做吧!

  1. # 使用 'test_1' 第一个和第二个自变量的值
  2. # 创建一个新观测,为 .4 和 .6
  3. x_test = np.array([[.4,.6]])
  4. # 将学习者应用于新的未分类的观测。
  5. trained_model.predict(x_test)
  6. # array(['loss'], dtype=object)

好哇! 我们可以看到学习器预测的新观测的类是“输”。

我们甚至可以查看学习器分配给每个分类的概率:

  1. trained_model.predict_proba(x_test)
  2. # array([[ 0.66666667, 0.33333333]])

根据这个结果,模型预测观测结果是“输”的概率约为 67%,“赢”的概率为 33%。 因为观测有更大的“输”的概率,所以它预测为这个分类。

  • K 的选择对创建的分类器有重大影响。
  • K 越大,决策边界越线性(高偏差和低方差)。
  • 有多种方法可以测量距离,两种流行的方法是简单的欧几里德距离和余弦相似度。

基于半径的 KNN 分类器

十四、K 最近邻 - 图11

  1. # 加载库
  2. from sklearn.neighbors import RadiusNeighborsClassifier
  3. from sklearn.preprocessing import StandardScaler
  4. from sklearn import datasets
  5. # 加载数据
  6. iris = datasets.load_iris()
  7. X = iris.data
  8. y = iris.target
  9. # 创建标准化器
  10. standardizer = StandardScaler()
  11. # 标准化特征
  12. X_std = standardizer.fit_transform(X)

在 scikit-learn 中,RadiusNeighborsClassifierKNeighborsClassifier非常相似,但有两个参数除外。 首先,在RadiusNeighborsClassifier中,我们需要指定固定区域的半径,用于确定观测是否是半径内的邻居。 将半径设置为某个值,最好将其视为任何其他超参数,并在模型选择期间对其进行调整。 第二个有用的参数是outlier_label,它表示半径内没有观测的观测的标签 - 这本身通常可以是识别异常值的有用工具。

  1. # 训练半径邻居分类器
  2. rnn = RadiusNeighborsClassifier(radius=.5, n_jobs=-1).fit(X_std, y)
  3. # 创建两个观测
  4. new_observations = [[ 1, 1, 1, 1]]
  5. # 预测两个观测的类别
  6. rnn.predict(new_observations)
  7. # array([2])