PCA主成分分析测序技术的发展使得现在能够从宏观角度分析基因表达，但是也在一定程度上增加了数据分析难度。许多基因之间可能存在相关性，如果分别对每个基因进行分析，分析往往是孤立的，盲目减少指标会损失很多有用的信息。PCA(PrincipalComponentAnalysis)，即主成分分析方法，是一种使用*****的数据降维算法。一般可应用的研究方向有：一组基因在多个分组中的差异情况，多个基因在该样本中的差异情况。基本原理PCA的主要思想是将n维特征映射到k维上，这k维是全新的正交特征也被称为主成分，是在原有n维特征的基础上重新构造出来的k维特征。PCA的工作就是从原始的空间中顺序地找一组相互正交的坐标轴，新的坐标轴的选择与数据本身是密切相关的。其中，**个新坐标轴选择是原始数据中方差**的方向，第二个新坐标轴选取是与**个坐标轴正交的平面中使得方差**的，第三个轴是与第1，2个轴正交的平面中方差**的。依次类推，可以得到n个这样的坐标轴。通过这种方式获得的新的坐标轴，我们发现，大部分方差都包含在前面k个坐标轴中，后面的坐标轴所含的方差几乎为0。于是，我们可以忽略余下的坐标轴，只保留前面k个含有绝大部分方差的坐标轴。事实上。 两个实验组的差异基因比较。成果发表指导数据科学售后服务

    Adonis（置换多元方差分析，分析不同分组或环境因子对样品差异的解释度）：ADONIS置换多元方差分析（Permutationalmultivariateanalysisofvariance，PERMANOVA），又称非参数多因素方差分析（nonparametricmultivariateanalysisofvariance）、或者ADONIS分析。使用PERMANOVA可分析不同分组因素对样品差异的解释度，并使用置换检验进行***性统计。基本原理：置换多元方差分析（PERMANOVA，Adonis）是一种基于F统计的方差分析，依据距离矩阵对总方差进行分解的非参数多元方差分析方法。基本步骤是基于OTU丰度表，计算样本间样本间Bray-curtis距离，然后adonis分析生成结果，绘图展示。术语解读：OTU：operationaltaxonomicunits，分类单元Df：自由度，其值=所比较的分组数量-1；SumsOfSqs：即Sumsofsquares，总方差，又称离差平方和；MeanSqs：即Meansquares，均方（差）；FModel：F检验值；R2：即Variation(R2)，方差贡献，表示不同分组对样品差异的解释度，即分组方差与总方差的比值，R2越大表示分组对差异的解释度越高；Pr(>F)：***性p值，小于***。数据要求：OTU丰度表或者样本距离矩阵。 云南公共数据库挖掘数据科学早期肝疾病的预后基因panel研究。

    GSEA术语解读Enrichmentscore（ES）ES是GSEA**初的结果，反应关注的基因集S在原始基因数据序列L的顶部或底部富集的程度。ES原理：扫描排序序列，当出现一个基因集S中的基因时，增加ES值，反之减少ES值，一个基因的ES值权重与差异表达度相关。ES是个动态值，**终ES是动态扫描过程中获得的**ES值。如果**终ES为正，表示某一功能基因集S富集在排序序列顶部。ES为负，表示某一基因集S富集在排序序列底部。NES由于ES是根据分析的排序序列中的基因是否在一个基因集S中出现来计算的，但各个基因集S中包含的基因数目不同，且不同功能基因集S与原始数据之间的相关性也不同，因此比较数据中基因在不同基因集S中的富集程度要对ES进行标准化处理，也就是计算NES。NES=某一基因集S的ES/数据集所有随机组合得到的ES平均值，NES是主要的统计量。nominalp-value（普通P值）描述的是针对某一功能基因集S得到的富集得分的统计***性，通常p越小富集性越好。FDR（多重假设检验矫正P值）NES确定后，需要判断其中可能包含的错误阳性发现率。FDR=25%意味着对此NES的判断4次可能错1次。GSEA结果中，高亮显示FDR<25%的富集基因集S。因为从这些功能基因集S中**可能产生有意义的假设。大多数情况下。

    ROC机器学习受试者工作特征曲线（receiveroperatingcharacteristiccurve，简称ROC曲线），又称为感受性曲线（sensitivitycurve），是用来验证一个分类器（二分）模型的性能的。一般应用于直观展示敏感性和特异性连续变量的综合指标，如比较多个biomarker或临床参数的诊断表现、比较多个算法的分类效果。基本原理ROC曲线工作原理是，向模型中输入已知正负类的一组数据，对比模型对该组数据的预测，衡量这个模型的性能。术语解读:1、TP（TruePositive,真正,TP）被模型预测为正的正样本（原来为正预测为正）2、TN（TrueNegative,真负,TN）被模型预测为负的负样本（原来为负预测为负）3、FP（FalsePositive,假正,FP）被模型预测为正的负样本（原来为负预测为正）4、FN（FalseNegative,假负,FN）被模型预测为负的正样本（原来为正预测为负）5、真正类率(TruePostiveRate)TPR:TP/(TP+FN),**分类器预测的正类中实际正实例占所有正实例的比例。Sensitivity6、假正类率(FalsePostiveRate)FPR:FP/(FP+TN)，**分类器预测的负类中预测为正实例(实际为负实例)占所有负实例的比例。1-Specificity7、真负类率(TrueNegativeRate)TNR:TN/(FP+TN)。 实验室致病类病原微生物数据分析平台。

industryTemplate云生物立足于上海，提供相关数据科研咨询与服务。湖北生物/药物信息学分析数据科学服务

提供语言润色、图表调整、格式修改等工作模块。成果发表指导数据科学售后服务

ssGSEA基本原理

对于一个基因表达矩阵，ssGSEA首先对样本的所有基因的表达水平进行排序获得其在所有基因中的秩次rank。然后对于输入的基因集，从基因集中寻找表达数据里存在的基因并计数，并将这些基因的表达水平求和。接着基于上述求值，计算通路中每个基因的富集分数，并进一步打乱基因顺序重新计算富集分数，重复一千次，***根据基因富集分数的分布计算p值整合基因集**终富集分数。

数据要求

1、特定感兴趣的基因集（通常为免疫细胞表面marker genes），列出基因集中基因

2、基因表达矩阵，为经过log2标准化的芯片数据或者RNA-seq count数数据（基因名形式与基因集对应）

下游分析

免疫细胞浸润分数相关性（corralation）分析 成果发表指导数据科学售后服务