论文部分内容阅读
在生物信息爆炸的今天,尤其是基因芯片技术的快速发展,使得越来越多的科学家被吸引到生物信息领域当中,而基因聚类也成为该领域当中热门的课题,其有着非常重要的研究价值。基因表达数据为疾病检测提供了一种快速而有效的诊断方法,但其表达数据维数高、样本量小,以及大噪音等特点,这使得选择或提取相关的特征信息成为一项有挑战性的工作。非负矩阵分解(NMF),一种有效的处理海量非负数据的方法,具有实现简单,分解速度快以及分解结果有明确物理意义等优点,已成为高维数据分析和降维的重要研究方向。 本文首先对现有的基于二范数和KL散度的矩阵分解算法背景知识,研究现状,实现意义进行一个概要描述。同时介绍了基因聚类问题在生物信息中的重要性,以及目前基因聚类的研究动向和发展前景,在此基础上提出了β散度的非负矩阵分解算法。由于该算法中参数β选择的不确定性,为我们研究矩阵分解提供了更多的选择,即我们可以通过不同的数据类型选择不同的参数值,这使得该算法使用起来更加的灵活。论文通过梯度下降的方法得出β散度非负矩阵分解迭代规则,对算法的收敛性进行了证明。并应用β散度矩阵分解对肿瘤基因表达数据进行降维,特征的选择,然后将其与K-means相结合用于聚类分析。论文讨论了参数β对实验结果的影响,通过实验得出β为0.5时该算法的分解结果能够较好的表示基因表达数据,其聚类效果较显著。针对β散度分解算法在提取数据特征信息不够显著做了如下两个方面的的改进:一是基于正交约束的β散度非负矩阵分解,在β散度非负矩阵分解目标函数中加入正交约束。二是将L0范数约束引入到目标函数中,在使得目标函数最小化的条件下对数据矩阵进行分解。算法在保证分解数据低维特征的条件下,使得分解误差最小,提高分解矩阵的稀疏表示能力,然后用 K-means对分解的基因表达数据进行聚类研究。最后将改进的算法分别与β散度分解、最大熵分解、传统的矩阵分解等算法的结果进行聚类对比,通过对实验的分析我们得出改进算法能够较好地提高聚类精度。