懒麻蛇

非负矩阵分解(NMF)的几个相关运用

2022年的存货，也是之前NMF推送中留的坑。

>非负矩阵分解NMF介绍<

本文列出了几项在生物信息学和神经影像学领域中应用NMF分析的研究。

应该是最早将NMF运用到微阵列(Microarray)数据上的文章。作者展示了NMF从癌症相关的微阵列数据中检测生物信息的能力。R包NMF中很多内容都受到了这篇文章的影响，比如默认算法Brunet 2004，甚至输入的矩阵形式采用的都是feature x sample的形式。NMF示意图如下，其中W称之为metagenes，H称之为metagene expression profile。(注意A是feature x sample的矩阵)

rank(k)增加意味着对样本的进一步划分。下图显示了等级k=2、3、4、5时产生的consensusmap。清晰的块状对角线模式证明了2、3和4类模型的稳健性，而5类模型展示出分散性(dispersion)的增加，反映为cophenetic系数的下降。

作者在最后一个实例中使用了四种类型的中枢神经系统胚胎性肿瘤的数据，10个典型的髓母细胞瘤，10个恶性胶质瘤，10个横纹肌瘤，和4个正常样本。分别使用了层次聚类(HC)，自组织映射(SOM)和NMF。

图a. HC虽然展示出了四个类型的分类，但是检查分支后发现，正常样本和恶性胶质瘤在一个分支，髓母细胞瘤和横纹肌瘤在另一个分支。总之层次聚类结果很差。

图b. SOM的三中心的聚类是最稳健的，具有最高的cophenetic系数。与HC的情况一样，在这种情况下，正常和恶性胶质瘤样本始终聚集在一起。而四中心聚类并不稳定，cophenetic系数下降。因此无法使用SOM方法正确识别四种样本。

图c. NMF能更好的识别该数据集中的样本种类，只有两个被错误分类的样本。

什么是SOM自组织映射？

每次迭代更新所有的中心点(神经元)。

Brunet, J. P., Tamayo, P., Golub, T. R., & Mesirov, J. P. (2004). Metagenes and molecular pattern discovery using matrix factorization. Proceedings of the national academy of sciences, 101(12), 4164-4169.

2021年Molecular Psychiatry的文章，对DNA甲基化数据做NMF，下图便是NMF选择rank=2的依据。再对系数做Kmeans确定两个分类。甲基化分析表明，在确定的亚型之间的T细胞比例有明显的差异。由于这一发现，他们还测试了该分区是否可以通过T细胞受体β复合物得到验证，就是额外做一个TCR来验证NMF的分类和相应的假设。Dice相似性系数（DSC）被用来评估甲基化亚型和TCR亚型之间的重叠情况。

Luo, C., Pi, X., Hu, N., Wang, X., Xiao, Y., Li, S., ... & Lui, S. (2021). Subtypes of schizophrenia identified by multi-omic measures associated with dysregulated immune function. Molecular psychiatry, 26(11), 6926-6936.

较早将NMF运用到MRI数据上的研究，文章很多技术细节都很清楚，值得一读。结果表明NMF可以识别出个体间的共同变化(co-vary)区域，并且它们一定程度上遵循了大脑的解剖结构和功能单元。文章使用分别使用了和的数据集系统比较了PCA/ICA/NMF在Sparsity/Incoherence/Regression accuracy/Reproducibility几方面的差异。

例如稀疏性的对比，

各种比较之后表明，相较于PCA和ICA，NMF有更好的特异性(specificity)和普适性(generalizability)。

需要注意的是，文章并没有使用标准的NMF，而是使用了projective NMF，和标准NMF的主要区别就是pNMF得到的系数矩阵更稀疏，并且可以直接用来作为数据分组的信息。那么作为专门提高稀疏性的方法，pNMF和PCA/ICA对比稀疏性，结果难道不是显而易见的吗。自此之后的三年内，UPenn出了一系列基于pNMF的文章，列表见文末。

Sotiras, A., Resnick, S. M., & Davatzikos, C. (2015). Finding imaging patterns of structural covariance via non-negative matrix factorization. Neuroimage, 108, 1-16.

使用NMF做聚类，在三个阿尔茨海默病数据集中确定了结构性MRI的萎缩亚型，并在临床、认知和生物学特征方面对亚型进行了比较。发现四种亚型：(i) 内侧-颞部为主的萎缩，记忆和语言功能最差，年龄较大，CSF tau水平最低，血管病变量最高；(ii) 顶枕部萎缩，执行/注意力和视觉空间功能差，CSF tau高；(iii) 轻度萎缩，认知表现最好，年龄小，但CSF tau水平最高；(iv) 弥漫性皮质萎缩，临床、认知和生物学特征居中。

之后将有前驱症状的(prodromal)被试根据这四种亚型进行分类，并测试他们是否在特定领域表现出认知能力下降。结果发现有前驱症状的被试已经表现出了类似的亚型特征。对比前驱病症被试的不同亚型，内侧-颞部亚型在记忆和语言方面随着时间的推移下降最快；顶枕部亚型在执行/注意力领域下降最快；弥漫性亚型在视觉空间功能方面下降最快；而轻度亚型在所有领域都表现出中间性的下降。

灰质分割成了1024个同等大小的感兴趣区域，获取灰质体积。之后用R做了标准NMF。为什么选的4类？他们在三个数据集中对每个聚类top100的特征计算了Dice系数来评估各数据集的聚类结果的重合度。此外，为了刻画每个亚型的特征，对亚型内的被试计算top100特征的平均值。对于新数据(前驱症状被试)和那个亚型的特征更相似则归为该亚型。top100计算Dice系数和对新数据做预测值得借鉴。

Ten Kate, M., Dicks, E., Visser, P. J., van der Flier, W. M., Teunissen, C. E., Barkhof, F., ... & Alzheimer’s Disease Neuroimaging Initiative. (2018). Atrophy subtypes in prodromal Alzheimer’s disease are associated with cognitive decline. Brain, 141(12), 3443-3456.

2022年1月31日上线的文章，NMF分解生成的系数被用于后续分析，Normative model建模和高斯混合模型聚类属常规操作。分析中的一些步骤存疑，这里仅展示其运用NMF的思路。

对灰质体积(图A)做spare-NMF(也就是在loss函数中加入L1正则)，相比标准NMF，sparse-NMF可以提高分解后矩阵的稀疏性(就不需要像paper4那样选top100特征)。根据Paper3以及后续的研究，pNMF也是提高稀疏性的方法。

对ABIDE2_TD做完sparse-NMF后得到一个特征矩阵(图B)，也就是几个Factors，之后使用该特征矩阵对另外两个数据集做分解(图C/D)。用分解得到系数建模Normative model，之后对deviation用高斯混合模型聚类。

Shan, X., Uddin, L. Q., Xiao, J., He, C., Ling, Z., Li, L., ... & Duan, X. (2022). Mapping the heterogeneous brain structural phenotype of autism spectrum disorder using the normative model. Biological Psychiatry, 91(11), 967-976.

2020年BP-CNN上的文章，做的是MEG。如图A-D，将数据分段，在每个窗口通过计算单试次锁相构建了功能连接矩阵。此处的数据A的维度为 b × (b − 1)/2 × N 其中b是节点数量，N是时间窗口的数量。

如图E所示，对数据A做NMF分解，得到Network和Activations的大小分别为[b×(b-1)/2×k]和(k×N)，其中k是功能网络的数量，这里由误差曲线的拐点确定(如F所示)。然后(G)计算Energy和Entropy并做组间比较。初看以为是对3维数据做NMF感觉很特别，其实从实现上来看应该还是将2维的FC转成1维向量做的，这里用的NMF是sparse-NMF。

Phalen, H., Coffman, B. A., Ghuman, A., Sejdić, E., & Salisbury, D. F. (2020). Non-negative matrix factorization reveals resting-state cortical alpha network abnormalities in the first-episode schizophrenia spectrum. Biological Psychiatry: Cognitive Neuroscience and Neuroimaging, 5(10), 961-970.

使用了sNMF对临床（287变量）和结构数据（204个灰质、白质和脑脊液区域）进行聚类。临床数据研究发现一个基于阳性症状、阴性症状、抑郁和功能的四维临床模型，该模型与所有验证指标(包括早期症状、纵向和精神分裂症多基因风险评分)均显示出关联。结构数据分析揭示了一个与分布式脑体积减少、阴性症状增加、执行IQ降低和精神分裂症多基因风险评分增加相关的亚组。

文章使用了验证集，并且有纵向数据和多基因风险评分。

Dwyer, D. B., Buciuman, M. O., Ruef, A., Kambeitz, J., Dong, M. S., Stinson, C., ... & PRONIA Consortium. (2022). Clinical, brain, and multilevel clustering in early psychosis and affective stages. JAMA psychiatry, 79(7), 677-689.

【这个文章躺在这里的时候还是2022年，目前已发表】NMF找到空间模型后做组间比较，考察的是饮酒程度与这些空间模式皮层厚度减小的交互作用。结果发现在六种非负矩阵分解（NMF）衍生的皮层厚度模式中，无饮酒或低度饮酒者的皮层厚度下降速率在所有年龄组中相似。在中度饮酒和重度饮酒者中，存在不同程度的交互作用。

Github有代码，但是没找到NMF相关的代码，对于模型结果可视化用作补充材料，R代码可以借鉴。

https://htmlpreview.github.io/?https://github.com/sundelinustc/NCANDA_NMF/blob/main/SDL_R_NMF.html

分析的思路借鉴的是2017年PNAS的文章，推测应该用了pNMF。

Sun, D., Adduru, V. R., Phillips, R. D., Bouchard, H. C., Sotiras, A., Michael, A. M., ... & Morey, R. A. (2023). Adolescent alcohol use is linked to disruptions in age-appropriate cortical thinning: an unsupervised machine learning approach. Neuropsychopharmacology, 48(2), 317-326.

Sotiras, Aristeidis, Jon B. Toledo, Raquel E. Gur, Ruben C. Gur, Theodore D. Satterthwaite, and Christos Davatzikos. Patterns of coordinated cortical remodeling during adolescence and their associations with functional specialization and evolutionary expansion. Proceedings of the National Academy of Sciences 114, no. 13 (2017): 3527-3532.

非负矩阵分解(NMF)的几个相关运用

你可能感兴趣的:(矩阵,机器学习,线性代数,人工智能)