四川生物/药物信息学分析数据科学口碑推荐

ssGSEA基本原理

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

数据要求

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

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

下游分析

免疫细胞浸润分数相关性（corralation）分析 我们团队具备完整的数据分析、数据库构建、软件开发团队。四川生物/药物信息学分析数据科学口碑推荐

Inmmune gene

免疫学研究是目前科研领域争相研究的热点，**免疫细胞浸润是其中一种。**免疫细胞浸润是指免疫细胞从血液中移向**组织发挥作用。我们从**组织中分离出浸润免疫细胞含量，计算基因与浸润免疫细胞含量的相关性，筛选出影响免疫浸润的候选基因。

基本原理：

从基因矩阵数据中提取免疫细胞含量，生成免疫细胞含量矩阵；

计算目标基因与浸润免疫细胞含量的相关性，筛选与浸润免疫细胞含量高度相关的基因。

术语解读：

相关性系数（pearson,spearman, kendall）反应两个变量之间变化趋势的方向以及程度。相关系数范围为-1到+1。0表示两个变量不相关，正值表示正相关，负值表示负相关，值越大表示相关性越强。

数据要求：

**数据表达矩阵 上海组学实验数据科学欢迎咨询参考国内外数据资源，根据需求制定构建方案。

    RoastROAST是一种差异表达分析方法，有助于提高统计能力、组织和解释结果以及在不同实验中的关联表达模式，一般适用于microarray、RNA-seq的表达矩阵，用limma给全部基因做差异表达分析，不需要筛差异表达基因。基本原理：ROAST是一种假设驱动的测试，对结果基因集做富集分析，富集分析考虑基因集中基因的方向性(上调或下调)和强度(log2倍变化)，判断上/下调基因是否***富于集目标基因集；ROAST使用rotation,一种MonteCarlotechnology的多元回归方法，适用于样本数量较少的情况；roast检验一个geneset，对于复杂矩阵，使用mroast做multipleroasttests。富集分析结果用barcodeplot展示，使上/下调基因在目标基因集中的分布可视化。数据要求：表达矩阵。

    STEM基因表达趋势分析基因调控网络是一个连续且复杂的动态系统。当生物体按照一定顺序发生变化或者受到外界环境刺激（如受到不同浓度的化学药物诱导）时，基因表达变化也会呈现趋势特征。趋势分析就是发现基因表达的趋势特征，将相同变化特征的基因集中在一种变化趋势中，从而找到实验变化过程中相当有有代表性的基因群。STEM（ShortTime-seriesExpressionMiner），中文名短时间序列表达挖掘器。该软件主要用于分析短时间实验数据，也可用于多组小样本数据。推荐3至8组数据。一般可应用的研究方向有：多个时间点的时间序列数据，例如多个发育时期、处理后多个时间点取样。基本原理STEM采用了一种新的聚类算法来分析时间序列基因表达趋势。聚类算法首先选择一组不同的、有代表性的时间表达模式（temporalexpressionprofiles）作为模型（modelprofiles）。模型是**于数据选择的，并从理论上保证了所选择的模型剖面具有代表性。然后，根据每个标准化过后的基因表达模式，分配给模型中相关系数比较高的时间表达模式。由于模型的选择是**于数据的，因此该算法可以通过排列测试，确定哪些时间表达模式在统计意义上***富集基因。对每一个基因都分配时间表达模式完成后。 自有服务器机房，可随时调用各计算平台算力，且团队成员有多年科研经历。

    Nomogram列线图（nomogram，诺莫图）是在平面直角坐标系中，用一簇互不相交的线段表示多个临床指标或者生物学特征，用以预测一定的临床结局或者某类事件发生的概率的图。列线图使预测模型的结果更具有可读性，可个性化地计算特定**患者生存率,在临床实践中有较大的价值。一般可应用的研究方向有：将回归的结果进行可视化呈现，对个体样本给出其发病风险或比例风险；根据多个临床指标或生物学特征，判断个体样本的疾病分类或特征。基本原理：列线图的理论于1884年提出，**早用于工程学。它能够将复杂的计算公式以图形的方式，快速、直观、精确的展现出来。列线图通过构建多因素回归模型（例如Cox回归、Logistic回归等），根据模型中各个影响因素对结局变量的影响程度的高低，即回归系数的大小，给每个影响因素的每个取值水平进行赋分。将各个评分相加得到总评分，通过总评分与结局事件发生概率之间的函数转换关系，从而计算出该个体结局事件的预测概率。校准曲线（calibrationcurve）为实际发生率和预测发生率的散点图，常于用于化工行业溶液配制。在这里通过观察预测值与实际值相差情况，判断基于回归模型构建列线图的有效性。 糖尿病药物基因组学分析找到新的作用靶点。天津成果发表指导数据科学经验丰富

甲状腺疾病的靶向药物研究。四川生物/药物信息学分析数据科学口碑推荐

    GeneInteraction基因互作：基因相互作用指miRNA、lncRNA、circRNA或其它RNA介导DNA转录，从而影响mRNA的表达过程。通俗意义上来说，基因互作关系指基于序列预测的靶基因对。miRNA通过与靶mRNA的结合，或促使mRNA降解，或阻碍其翻译，从而***目的基因的表达。竞争性内源RNA网络是靶基因预测的研究深入，简称ceRNA网络。通过进行ceRNA网络的分析，我们能从一个更为宏观的角度来解释转录体如何构建基因表达调控网络，从而进一步挖掘基因在其中的调控机制。基本原理：miRNA主要通过与靶基因的非翻译区（UTR）结合而发挥其作用，对miRNA和mRNA、lncRNA、circRNA结合进行的预测称为靶基因预测。靶基因预测使用软件根据miRNA和靶基因间的结合的规律预测结合基因对。在生物体内，miRNA可以通过与proteincoding特异性结合，影响相关基因的表达，从而参与调控细胞内的各项功能。ceRNA具有miRNA结合位点，能后竞争性地结合miRNA，***miRNA对靶基因的调控。例如lncRNA与miRNA竞争性结合，影响miRNA调控mRNA的过程，**终导致的mRNA表达失调。我们使用基于序列预测的软件对差异分析得到的miRNA与mRNA，lncRNA，circRNA进行靶点预测和ceRNA网络分析。 四川生物/药物信息学分析数据科学口碑推荐