天津成果发表指导数据科学共同合作

ssGSEA基本原理

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

数据要求

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

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

下游分析

免疫细胞浸润分数相关性（corralation）分析 协助构建各类科研、临床数据库。天津成果发表指导数据科学共同合作

    cox风险比例回归模型：产品详情产品评论(0)比例风险回归模型，又称Cox回归模型，是由英国统计学家。模型可以用来描述了不随时间变化的多个特征对于在某一时刻死亡率的影响。它是生存分析中的一个重要的模型。应用场景cox比例风险回归模型，由英国统计学家主要用于**和其他慢性疾病的预后分析，也可用于队列研究的病因探索单因素cox分析主要探索单个基因的**预后影响cox分析可用于转录组，甲基化，miRNA,LncRNA,可变剪切等等基本原理：在这里，是一个与时间有关的基准危险率，其选择具有充分的灵活度，一种可能的选择是采用概率论中的Weibull分布。是模型的参数。由于只要给定数据，就能够通过极大似然估计求出模型的参数，而的选择具有很大的灵活性，所以我们称之为一个半参数模型。对公式进行变形，得到：通过这个公式，我们可以发现，模型中各危险因素对危险率的影响不随时间改变，且与时间无关,同时，对数危险率与各个危险因素呈线性相关。这就是Cox回归中的两个基本假设。参数的极大似然估计：术语解读：1.输入变量，由m个影响因素组成：2.生存函数，输入为X时，在t时刻仍然存活的概率：3.死亡函数，输入为X时，在t时刻已经死亡的概率：4死亡密度函数，输入为X时。 北京生物/药物信息学分析数据科学服务可对接各类公共数据库，切入各类接口，并对公共数据库进行大规模数据挖掘。

sankey

桑基图（sankey）是一种数据流图，每条边**一条数据流，宽度**数据流的大小。一套数据集可能有多重属性，每层属性之间有交叉，就可以用这种图来展示。一般应用场景：分组与基因为多对多关系，展示高频突变基因所处的分组；miRNA和靶基因的关系；人群按性别、年龄、家族史等特征分组，展示不同分组得**的规律。

数据要求：

多个分组及其关系，包括且不限于基因表达、突变。

下游分析：

1.   补充展示部分的已有相关研究

2.   解释展示部分对研究课题的意义

    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中**可能产生有意义的假设。大多数情况下。 指导科研方案纠偏，更好更快发表文章。

    GSVA（基因集变异分析，反映了样本和感兴趣的通路之间的联系）：GSVA全名Genesetvariationanalysis（基因集变异分析），是一种非参数，无监督的算法。与GSEA不同，GSVA不需要预先对样本进行分组，可以计算每个样本中特定基因集的富集分数。换而言之，GSVA转化了基因表达数据，从单个基因作为特征的表达矩阵，转化为特定基因集作为特征的表达矩阵。GSVA对基因富集结果进行了量化，可以更方便地进行后续统计分析。如果用limma包做差异表达分析可以寻找样本间差异表达的基因，同样地，使用limma包对GSVA的结果（依然是一个矩阵）做同样的分析，则可以寻找样本间有***差异的基因集。这些“差异表达”的基因集，相对于基因而言，更加具有生物学意义，更具有可解释性，可以进一步用于**subtype的分型等等与生物学意义结合密切的探究。 利用甲基化数据分析样本的拷贝数变异。北京生物/药物信息学分析数据科学服务

云生物提供数据科学服务。天津成果发表指导数据科学共同合作

industryTemplate天津成果发表指导数据科学共同合作