K邻近算法（KNN，K-nearest Neighbors Algorithm）

文章目录

- 前言
- 应用场景
- 欧几里得距离（欧氏距离）
- 两类、单一属性（1D）
- 两类、两种属性（2D）
- 两类、两种以上属性（>3D）
Examples in R
- 再来一个
- 补充一下什么是变量
什么是变量？
什么是数值、分类变量？

前言

之前看到一篇文章，方法部分提到了这个K邻近算法，正好自己不是很熟悉其原理，所以学习整理一下😑

K邻近算法（K-nearest neighbors algorithm）是一种常用的机器学习算法，用于分类和回归问题。它基于一个简单的假设：与未知样本最相似的K个已标记样本的类别可以用来预测该未知样本的类别。

在K邻近算法中，输入数据点被表示为n维空间中的向量，并且每个数据点都有一个对应的类别标签。算法的工作原理如下：

计算距离：根据给定的距离度量方法（例如欧氏距离、曼哈顿距离等），计算未知样本与训练集中所有已标记样本之间的距离。
选择最近邻：选择距离未知样本最近的K个已标记样本。
预测类别：对于分类问题，通过多数表决的方式确定未知样本的类别。即，选择K个最近邻中出现次数最多的类别作为预测结果。对于回归问题，可以使用这些最近邻的平均值作为预测结果。

K值的选择很重要，较小的K值会使模型更加敏感和复杂，可能会过拟合。较大的K值则可能导致模型过于简单，无法捕捉到数据集中的复杂模式。K邻近算法简单易懂，且适用于多种问题。然而，它也有一些限制，如对于高维数据和样本不平衡问题的处理效果可能较差。此外，由于该算法需要计算距离矩阵，对于大型数据集来说会变得计算密集。

应用场景

假设有一定数量的物体，每个物体都有其独特的属性，比如，我们有有椅子床和桌子（对象），并知道其对应的长宽高（属性）。如果有人给我们一个具有已知属性的对象，让我们猜测（预测）该对象属于什么，就是说，已知数据的维度，要求预测它是椅子、床还是桌子，就可以使用knn算法。
属性是对象的属性，每个对象都可以看作一个范畴。我们可以检查新对象的属性与任何已知类别的关系。当属性已知，可以画在图上，图形表示能让我们更易于理解并计算新对象与已知类别之间的欧几里得距离，确定新对象最接近哪个类别。

欧几里得距离（欧氏距离）

欧几里得是个人名，不要被这个名字吓倒，它只是简单地表示平面上两点之间的距离。通过简单地使用公式，可以计算两点之间的距离，不管你有多少属性，比如高度、宽度、宽度、重量等等。公式为√（x2−x1）²+（y2−y1）²+（z2−z1）²……（n2-n1）²

两类、单一属性（1D）

我们有两个类别：男性和女性，其各自的高度在下表中给出。

性别	身高（cm）
男	178
男	179
女	163
女	168
男	181
女	170
男	183
女	171

在这里插入图片描述

此时出现一个新成员，并且需要确定它是男是女。已知其身高，在1D平面上画图，看新对象的属性更接近哪里。理想情况下，我们可以计算欧几里德距离，以确定最接近新对象属性的值。如果我们在身高图上画180cm这个点（新对象的属性），它更接近男性的身高。所以推测新对象更可能是男性。

两类、两种属性（2D）

现在再加入一个新的属性：体重，它也可以描述男性和女性的特征，如下表所示。

性别	身高（cm）	体重（kg）
男	178	72
男	179	81
女	163	54
女	168	57
男	181	97
女	170	59
男	184	77
女	171	58

在这里插入图片描述

现在，我们创建一个二维平面，并按照相同的步骤计算新对象与旧对象的距离，它与其中一个类别的距离越近，新对象属于该类别的可能性就越大。可以看到，新成员（属性：身高169cm、体重68kg）与女性更加接近。所以推测其为女性。

两类、两种以上属性（>3D）

通常情况下，有很多属性与分类对象相关联，不能简单画在二维或一维平面上。假设有5个属性：性别、身高、体重、瞳孔的颜色、头发长度和音高，只需使用欧几里德距离公式来计算新对象与给定的对象之间的距离。
$\sqrt{(sex1-sex2)^2 +(height1-height2)^2 + (weight1-weight2)^2+(eye\_color1-eye\_color2)^2+(hair\_length1-hair\_length2)^2+(voice\_pitch1-voice\_pitch2)^2}$
其中1表示已知类别数据点，2表示要确定其类别的新数据点。相比较小，新的对象与已知某类数据的距离（和）越小，则属于该类的可能性越大。

Examples in R

在使用R进行KNN之前，需要说明几点：

KNN，K邻近法中的neighbor数量要有K定义，如果K＝５，那么将将会检测最近的五个neighbor以确定所属类别。如果这５个neighbor中的大部分同属于一个类别，那么可以认定新对象很大可能属于该类。
对象不同的属性会有不同的标度单位，比如：KG、CM等，所以在使用数据之前需要对每个变量进行标准化（特征缩放），比如Min－Max法。
KNN算法更适用于数值变量，但不是说它不能处理分类变量，但如果混合了分类变量和数值变量，那就需要其他方法。如果所有值都是数值，KNN是最好的method。
把数据分成训练集和测试集时，数据应该已经标准化了：数据标准化 → 拆分。
KNN算法不适用于R的有序Factors
K-mean算法和KNN算法是不同的，K-mean用于聚类，是一种无监督的学习算法，而KNN是一种用于分类问题的监督学习算法。

df <- data(iris) 
head(iris, 2) 
##   Sepal.Length Sepal.Width Petal.Length Petal.Width Species
## 1          5.1         3.5          1.4         0.2  setosa
## 2          4.9         3.0          1.4         0.2  setosa# Min-Max方法，标准化
nor <- function(x){(x -min(x))/(max(x)-min(x))}# 生成随机数，包含数据集总行数的80%
ran <- sample(1:nrow(iris), 0.8 * nrow(iris)) # 对数据集的前4列标准化
iris_norm <- as.data.frame(lapply(iris[,c(1:4)], nor))summary(iris_norm)
# Sepal.Length     Sepal.Width      Petal.Length     Petal.Width     
# Min.   :0.0000   Min.   :0.0000   Min.   :0.0000   Min.   :0.00000  
# 1st Qu.:0.2222   1st Qu.:0.3333   1st Qu.:0.1017   1st Qu.:0.08333  
# Median :0.4167   Median :0.4167   Median :0.5678   Median :0.50000  
# Mean   :0.4287   Mean   :0.4406   Mean   :0.4675   Mean   :0.45806  
# 3rd Qu.:0.5833   3rd Qu.:0.5417   3rd Qu.:0.6949   3rd Qu.:0.70833  
# Max.   :1.0000   Max.   :1.0000   Max.   :1.0000   Max.   :1.00000  # 提取训练数据集
iris_train <- iris_norm[ran,] 
# 提取测试数据集
iris_test <- iris_norm[-ran,] 
#提取train数据集的第5列，用作knn函数中的'cl'参数。
iris_target_category <- iris[ran,5]
# 提取test数据集的第5列，衡量预测准确性
iris_test_category <- iris[-ran,5]
library(class)
## 执行KNN，预测
pr <- knn(iris_train,iris_test,cl=iris_target_category,k=13)
(tab <- table(pr,iris_test_category))
#            iris_test_category
#pr           setosa versicolor virginica
#  setosa         12          0         0
#  versicolor      0          8         1
#  virginica       0          0         9# 这个函数将正确的预测除以总预测数，告诉我们模型有多精确。
accuracy <- function(x){sum(diag(x)/(sum(rowSums(x)))) * 100}
accuracy(tab)
# [1] 96.66667

在iris数据集中，使用k近邻法，最终96.7%的准确率。首先，标准化数据集，然后，将标准化值分离为训练和测试数据集。

想象一下，将训练数据集的值绘制在一个图上，然后在运行带有所有必要参数的knn函数后，将测试数据集的值引入到图中，并计算出与图中每个已知数据点的欧氏距离。

我们将花种类的预测值存储在“pr”中，可以将预测值与原始测试数据集的值进行比较。就可以算出模型的准确性，如果有新的50个值，要求预测这50个值的类别，就可以用这个模型。

再来一个

data(diamonds,package = "ggplot2")
dia <- data.frame(diamonds)
head(dia,2)
#   carat     cut color clarity depth table price    x    y    z
# 1  0.23   Ideal     E     SI2  61.5    55   326 3.95 3.98 2.43
# 2  0.21 Premium     E     SI1  59.8    61   326 3.89 3.84 2.31
ran <- sample(1:nrow(dia),0.9 * nrow(dia))
dia_nor <- as.data.frame(lapply(dia[,c(1,5,6,7,8,9,10)], nor))
dia_train <- dia_nor[ran,]
dia_test <- dia_nor[-ran,]
##the 2nd column of training dataset because that is what we need to predict about testing dataset
##also convert ordered factor to normal factor
dia_target <- as.factor(dia[ran,2])##the actual values of 2nd couln of testing dataset to compaire it with values that will be predicted
##also convert ordered factor to normal factor
test_target <- as.factor(dia[-ran,2])
pr <- knn(dia_train,dia_test,cl=dia_target,k=20)
tb <- table(pr,test_target)
accuracy(tb)
## [1] 71.09752