原文链接:http://tecdat.cn/?p=19077
原文出处:拓端数据部落公众号
导入
自组织映射(SOM)是一种工具,通过生成二维示意来可视化高维数据中的模式,在高维构造中显示有意义的模式。通过以下形式应用给定的数据(或数据样本)对 SOM 进行“训练”:
-
定义了网格的大小。
-
网格中的每个单元都在数据空间中调配了一个初始化向量。
-
例如,如果要创立 22 维空间的地图,则会为每个网格单元调配一个 22 维向量。
-
数据被重复输出到模型中进行训练。每次输出训练向量时,都会执行以下过程:
-
辨认具备最靠近训练向量的代表向量的网格单元。
-
随着训练向量的屡次输出,收敛的参数使调整变得越来越小,从而使地图稳固。
该算法赋予 SOM 的要害特色:数据空间中靠近的点在 SOM 中更靠近。因而,SOM 可能是示意数据中的空间聚类的好工具。
Kohonen 映射类型
上面的示例将应用 2015/16 NBA 赛季的球员统计数据。咱们将查看每 36 分钟更新一次的球员统计信息。这些数据可从 http://www.basketball-reference.com/ 取得。咱们曾经清理了数据。
NBA <- read.csv("NBA_cleaned.csv",
sep = ",", header = T, check.names = FALSE)
根本 SOM
在创立 SOM 之前,咱们须要抉择要在其中搜寻模式的变量。
colnames(NBA)
## \[1\] """Player" "Pos" "Age" "Tm" "G" "GS"
## \[8\] "MP" "FG" "FGA" "FG%" "3P" "3PA" "3P%"
## \[15\] "2P" "2PA" "2P%" "FT" "FTA" "FT%" "ORB"
## \[22\] "DRB" "TRB" "AST" "STL" "BLK" "TOV" "PF"
## \[29\] "PTS"
咱们从简略示例开始:
som(scale(NBA\[res1\], grid = somgrid(6, 4, "rectangular")
请留神,咱们标准化了训练数据,并定义了网格大小。规范 SOM 图可为网格单元的代表矢量创立这些饼图示意,其中半径对应于特定维度上的大小。
热图 SOM
咱们能够通过将每个球员调配到具备最靠近该球员状态的代表向量来辨认地图。“计数”类型的 SOM 依据球员数量创立了一个热图。
# 色带
colors <- function(n, alpha = 1) {rev(heat.colors(n, alpha))
}
绘图点
您能够应用“映射”类型的 SOM 将球员绘制为网格上的点。咱们与惯例 SOM 进行可视化比拟。
每个地图单元格的代表性矢量显示在右侧。左侧是依据其状态与这些代表向量的靠近水平绘制的球员图表。
环形 SOM
下一个示例是一种更改几何形态的办法。在为上述示例训练 SOM 时,咱们应用了矩形网格。因为边缘(尤其是拐角处)的单元比外部单元具备更少的街坊,因而偏向于将更多的极其值推到边缘。
par(mfrow = c(1, 2))
plot(NBA.SOM2, type = "mapping", pchs = 20, main = "Mapping Type SOM")
plot(NBA.SOM2, main = "Default SOM Plot")
映射间隔
当用绘制时 type = "dist.neighbours"
,单元格将依据与它们最近的街坊的间隔着色,这使咱们能够直观地看到高维空间中不同因素之间的间隔。
plot(SOM2, type = "dist.neighbours")
有监督 SOM
有监督的 SOM 使咱们能够进行分类。到目前为止,咱们仅将三维数据映射到二维。当咱们解决更高维度的数据时,SOM 的实用性变得更加显著,因而让咱们应用扩大的球员统计信息列表来做这个受监督的示例:
咱们创立有监督的 SOM,并依据球员在球场上的地位对其进行分类。咱们将数据随机分为训练集和测试集。
indices <- sample(nrow(NBA), 200)
training <- scale(NBA\[indices, NBA.measures2\])
testing <- scale(NBA\[-indices, NBA.measures2\], center = attr(training,
"scaled:center"), scale = attr(training, "scaled:scale"))
请留神,当咱们从新标准化测试数据时,咱们须要依据训练数据的形式对其进行标准化。
您能够在训练算法中对训练变量(NBA.training
)与预测变量(NBA$Pos
)进行加权。当初让咱们查看预测的准确性:
##
## Center Point Guard Power Forward Shooting Guard
## Center 16 0 26 1
## Point Guard 0 49 0 12
## Power Forward 10 1 29 5
## Shooting Guard 0 8 4 38
## Small Forward 0 0 15 9
##
## Small Forward
## Center 4
## Point Guard 11
## Power Forward 8
## Shooting Guard 19
## Small Forward 38
可视化预测:
这次,咱们应用 xweight
参数为权重掂量球员统计数据。
应用 type = "codes"
咱们进行绘制,能够 失去规范的可视化球员状态(Codes X
)和球员地位预测(Codes Y
)。
add.cluster.boundaries(NBA.SOM4, NBA.SOM4.hc)
该视图使咱们能够将球员统计数据与地位预测进行比拟。
可视化预测:自定义 SOM
在最初一个示例中,咱们将对该 type = mapping
图进行一些自定义,以便咱们能够同时示意理论球员地位和 SOM 的预测地位。咱们将从可视化开始。
背景色彩绘制的球员点的背景代表其实在地位。
bg.pallet <- c("red", "blue", "yellow", "purple", "green")
# 为所有单元格制作仅背景色彩的矢量
base.color.vector <- bg.pallet\[match(position.predictions, levels(NBA$Pos))\]
# 设置 alpha 以最大的预测置信度标准化
max.conf <- apply(NBA.SOM4$codes$Y, 1, max)
最受欢迎的见解
1.R 语言 k -Shape 算法股票价格工夫序列聚类
2.R 语言中不同类型的聚类办法比拟
3.R 语言对用电负荷工夫序列数据进行 K -medoids 聚类建模和 GAM 回归
4.r 语言鸢尾花 iris 数据集的档次聚类
5.Python Monte Carlo K-Means 聚类实战
6. 用 R 进行网站评论文本开掘聚类
7. 用于 NLP 的 Python:应用 Keras 的多标签文本 LSTM 神经网络
8.R 语言对 MNIST 数据集剖析 摸索手写数字分类数据
9.R 语言基于 Keras 的小数据集深度学习图像分类