论文笔记：Unsupervised Deep Embedding for Clustering Analysis

前言

转载自：https://blog.csdn.net/SCS199411/article/details/90759274

至今已有许多工作致力于距离函数与嵌入方法的研究，用于执行聚类的特征空间无监督学习的的研究工作还比较少。对于所有最简单的图像数据集，使用欧氏距离的聚类完全失效。本文希望找到一种方案，可以联合解决特征空间学习与聚类关系判别。

基于此本文提出了一种深度嵌入聚类方法(Deep embedded clustering)

论文链接：https://arxiv.org/abs/1511.06335

1.介绍

本文定义了一个由数据空间X到低维特征空间Z的参数化非线性映射，在低维空间优化聚类目标。此前的工作在数据空间或是浅层线性嵌入空间上进行操作，本文在聚类目标上采用反向传播的SGD来学习映射，该映射被深度神经网络进行了参数化，被称为DEC。

如何优化DEC？本文希望同时解决聚类分配与潜在特征表示的学习，然而无法通过标签数据训练DNN。因此本文提出使用目前的soft cluster assignment派生出的辅助目标分布来迭代细化聚类。这个过程可以同时改进聚类与特征表示。

2.深度嵌入聚类方法（Deep embedded clustering）

考虑$n$个点的聚类问题，将划分为$k$个类，$j=1,2\dots \dots$，$k$表示$k$个聚类中心。并非直接在聚类空间$X$中聚类，我们首先将数据使用非线性映射:$X\to Z$，$\theta$是学习到的参数，$Z$是潜在特征空间。$Z$的维度一般来说要比$X$小的多。使用DNN来参数化。

DEC算法同时学习特征空间$Z$的$k$个聚类中心和用于映射的参数$\theta$。DEC有两个阶段：

• (1)使用deep autoencoder初始化参数
• (2)参数优化(聚类)，该过程中，本文在计算辅助目标分布于最小化KL之间交替进行。本文从参数优化过程开始阐述，给定$\theta$与参数中心$\mu$的初始估计。

2.1 通过 KL divergence 聚类

给定$\theta$与参数中心$u$的初始估计，本文提出使用分两步迭代进行的非监督算法来改进聚类。

• 第一步：计算嵌入点与聚类中心间的soft assignment（当某点与某聚类中心依概率符合分布时，将其分配给该中心）；
• 第二步：更新映射，通过从"使用辅助目标分布的高自信assignments"中学习来细化聚类中心。这个过程直到某种收敛准则符合而停止。

2.1.1 SOFT ASSIGNMENT

本文使用学生$t$-分布作为kernel来衡量嵌入点$z_i$与中心点${\mu }_j$的相似性：

其中， 对应于的嵌入点，$\alpha$是$t$分布的自由度，可以解释为将样本$i$分配给聚类$j$的概率（soft assignment）。因为在非监督学习中无法交叉验证$\alpha$，所以没必要学习它，统一设置为1。

软分配概念：首先明确基本思想：硬分配的意思就是，样本1就是类型1，样本2就是类型2，软分配的意思是样本1有30%的可能是类型1，70%的可能是类型2

2.1.2 KL DIVERGENCE MINIMIZATION

在辅助目标函数的帮助下，本文从高自信assignments中学习，迭代精炼聚类。本文模型通过将soft assignment与目标分布进行匹配来进行训练。为实现该目的，本文将目标函数设置为soft assignments与辅助分布之间的KL divergence损失：

目标分布$P$的选择对于DEC的表现来说至关重要。一般做法是将每个$p_i$设置为高于置信度阈值的数据点的delta分布（到最近的质心），并忽略其余部分。然而，由于$q_i$是soft assignments，因此随$q_i$使用softer probabilistic targets是一件自然而然的选择。本文希望target distribution具有以下性质：

• （1）强化预测
• （2）更加重视高可信度地分配的数据点
• （3）规范每个质心的损失贡献，以防止大类扭曲隐藏的特征空间

本文实验通过首先将$q_i$提高到第次幂然后按每个簇的频率归一化来计算$p_i$

其中，$f_j=\sum i{q}_{ij}$是软聚类频率。

2.1.3 OPTIMIZATION

本文联合优化聚类中心$u_j$和DNN的参数$\theta$（通过动量SGD）。L的梯度与数据点$z_i$的嵌入特征空间和每个聚类中心$u_j$有关，如下计算：

梯度$\partial L/\partial z$i随后被传递给DNN，用于标准反向传播来计算DNN的参数梯度$\partial L/\partial \theta$为了发现聚类分配的目的，本文设置当少于总量的tol%数据点在连续迭代间改变其所属类时停止优化。

2.2 参数初始化

使用堆叠自动编码器（SAE）初始化DEC的参数$\theta$；
使用k-means对用SAE得到的嵌入表示聚类，得到聚类中心作为DEC模型的初始化质心。

总体来说，第一步是最小化重建误差（初始化），第二步（DEC模型）是最小化KL散度。论文定义了基于质心的概率分布，并最小化其与辅助目标分布的KL偏差，以同时改进聚类分配和特征表示，而不是最小化KL散度以产生忠实于原始数据空间中的距离的嵌入。