视频实操SCI作图课(3)：OPLS-DA分析，组间差异的挖掘神器

在上一场小工具讲解中，小姐姐给大家介绍了PLS-DA的原理及用途，而在代谢组学数据分析中，除去PLS-DA以外，OPLS-DA分析也是非常常见的，仅一个字母之差，那二者到底有何差别，我们一起来一探究竟！

■ ■■■■

 视频解说教程

OPLS-DA分析，全称正交偏最小二乘法判别分析（OrthogonalPartialLeast Squares-DiscriminantAnalysis），它结合了正交信号矫正（OSC）和PLS-DA方法，能够将X矩阵分解成与Y相关和不相关的两类信息，通过去除不相关的差异来筛选差异变量。

OPLS-DA不同于PCA，它是一种有监督的判别分析统计方法。运用偏最小二乘回归建立代谢物表达量与样本类别之间的关系模型，来实现对样本类别的预测。OPLS-DA分析需要样本变量矩阵和样本分类矩阵两个文件来确立样本关系，如下所示：

X矩阵，样本-变量矩阵

Y矩阵，样本分类矩阵

OPLS-DA建模时，将X矩阵信息分解成与Y相关和不相关的两类信息，其中与Y相关的变量信息为预测主成分，与Y不相关的变量信息为正交主成分。根据OPLS-DA模型分析代谢组数据，绘制各分组的得分图，进一步展示各个分组之间的差异（Thévenotet al., 2015）。

OPLS-DA分析在实现降维的同时考虑了分组信息，因此它可以用于特征选择以及分类，也就是在代谢组学数据分析中，可以用于筛选不同组之间的差异代谢物。通过OPLS-DA分析，每个代谢物可以得出一个VIP值，即变量重要性投影（VariableImportance inProjection，VIP），VIP值越大，代表该物质对于区分两组所具有的贡献越大，因此我们在挑选差异代谢物时，通常会将VIP值作为其中一项重要的考察指标。

OPLS-DA分析结果中最常用的图就是OPLS-DA得分图，图中横坐标表示预测主成分，因此横坐标方向可以看出组间的差距；纵坐标表示正交主成分，因此纵坐标方向可以看出组内的差距；百分比表示该成分对数据集的解释率。图中的每个点表示一个样品，同一个组的样品使用同一种颜色表示，Group为分组。

■ ■■■■

OPLS-DA得分图

除去得分图以外，OPLS-DA分析还可以得到S-plot图，S-plot图的横坐标表示主成份与代谢物的协方差，纵坐标表示主成份与代谢物的相关系数。S-plot图一般用来挑选与OSC过程中主要成分的相关性比较强的代谢物，从另一方面同时也可以挑选与Y相关性强的代谢物。越靠近两个角的代谢物重要度越强。S-plot图中红色的点表明这些代谢物的VIP值大于等于1，绿色的点表示这些代谢物的VIP值小于等于1。

■ ■■■■

OPLS-DA的S-plot图

并非所有的数据都适合使用OPLS-DA模型进行分析，因此在模型建立之后，我们需要通过模型验证来对模型质量进行评价。

OPLS-DA评价模型的参数有R2X，R2Y和Q2，其中R2X和R2Y分别表示所建模型对X和Y矩阵的解释率，Q2表示模型的预测能力，这三个指标越接近于1时表示模型越稳定可靠，Q2 >0.5时可认为是有效的模型，Q2 >0.9时为出色的模型。

■ ■■■■

OPLS-DA模型验证图

上图为OPLS-DA模型验证图，图中横坐标表示模型R²Y，Q²值，纵坐标是模型分类效果出现的频数，即本模型对数据进行200次随机排列组合实验，若Q² 的p= 0.02，说明在此次Permutation检测中共有4个随机分组模型的预测能力优于本OPLS-DA模型，若R²Y的p= 0.545，说明在此次Permutation检测中共有109个随机分组模型其对Y矩阵的解释率优于本OPLS-DA模型。一般情况下，p< 0.05 时模型最佳。

免

费

注

册

迈维代谢

武汉迈特维尔生物科技有限公司（简称“迈维代谢”）位于武汉国家生物产业基地--光谷生物城，专注于提供领先的代谢组学及多组学技术开发及服务。迈维代谢创新了"广泛靶向代谢组"技术，基于“代谢组+基因组+转录组”的技术路线，近年来以通讯作者身份在Cell、Nature Genetics、Nature Communications、PNAS、National Science Reviewzen等国际学术期刊发表多篇论文，引领基因组时代的代谢生物学研究新方向。

咨询电话：027-62433042

咨询微信：metware888

邮箱：support@metware.cn

网址：www.metware.cn

我就知道你“在看”