GSEA全名为GeneSetEnrichmentAnalysis(基因集富集分析)。用以分析特定基因集(如关注的GO条目或KEGGPathway)在两个生物学状态(如**与对照,高龄与低龄)中是否存在差异。能够研究基因变化的生物学意义。SubtypeGSEA是在GSEA的基础上对不同亚型样本中重要通路的富集情况进行组间比较,能直观比较不同亚型中相同通路富集情况。基本原理GSEA主要分为基因集进行排序、计算富集分数(EnrichmentScore,ES)、估计富集分数的***性水平并进行多重假设检验三个步骤。**步对输入的所有基因集L进行排序,通常来说初始输入的基因数据为表达矩阵,排序的过程相当于特定两组中(case-control、upper-lower等等)基因差异表达分析的过程。根据所有基因在两组样本的差异度量不同(共有六种差异度量,默认是signal2noise,GSEA官网有提供公式,也可以选择较为普遍的foldchange),对基因进行排序,并且Z-score标准化。第二步是GSEA的**步骤,通过分析预先定义基因集S在**步获得的基因序列上的分布计算富集指数EnrichmentScore,并绘制分布趋势图Enrichmentplot。每个基因在基因集S的EnrichmentScore取决于这个基因是否属于基因集S及其差异度量(如foldchange)。 云生物提供数据科学服务。上海生物/药物信息学分析数据科学口碑推荐
ssGSEA基本原理
对于一个基因表达矩阵,ssGSEA首先对样本的所有基因的表达水平进行排序获得其在所有基因中的秩次rank。然后对于输入的基因集,从基因集中寻找表达数据里存在的基因并计数,并将这些基因的表达水平求和。接着基于上述求值,计算通路中每个基因的富集分数,并进一步打乱基因顺序重新计算富集分数,重复一千次,***根据基因富集分数的分布计算p值整合基因集**终富集分数。
数据要求
1、特定感兴趣的基因集(通常为免疫细胞表面marker genes),列出基因集中基因
2、基因表达矩阵,为经过log2标准化的芯片数据或者RNA-seq count数数据(基因名形式与基因集对应)
下游分析
免疫细胞浸润分数相关性(corralation)分析 湖北算法还原与开发数据科学售后分析自有服务器机房,可随时调用各计算平台算力,且团队成员有多年科研经历。
survivalCurve生存分析生存分析(survivalCurve)旨在更好地分析对不同因素对患者预后的影响,从而找到影响患者疾病的关键因素。生存曲线(Kaplan-Meier曲线)是生存分析的基本步骤,展示分类样本的生存曲线,从而揭示不同因素对疾病预后的影响。一般可应用的研究方向有:患者的生存期跟基因变异的关系、药物处理导致模式动物生存期变化。基本原理Kaplan-Meier法,直接用概率乘法定理估计生存率,故称乘积极限法(product-limitmethod),是一种非参数法。相比其他方法,KM曲线能更好的处理删失数据。先将样本生存时间从小到大排列。若遇到非删失值和删失值相同时,非截删失****。在生存时间后列出与时间相应的死亡人数,期初病例数(即生存期为某时间时尚存活的病例数)。然后计算活过每个时间点的生存率。以生存时间为横坐标,生存率为纵坐标所作的曲线,即为Kaplan-Meier曲线。术语解释风险比(HazardRatio,HR):Kaplan-Meier方法中计算的风险比HR为两分组对生存期影响的比例,用来描述该基因高表达对生存期的危险程度。该方法中的假设检验为两组中样本的生存期是否存在差异,即该因素是否会导致生存期的改变。删失(censored):在生存分析中。
单细胞测序数据挖掘:GEO目前收录的单细胞研究样本已经超过2万例,单细胞测序几乎成为生物医学领域CNS***文章的标配。实验费用高昂,阻断了CNS梦,既然其他数据可以挖,单细胞测序数据照样可以挖。已知公共数据库中单细胞测序数据涉及各种疾病类型,包括**、免疫细胞、炎症类甚至神经、肌肉、骨骼等,样本丰富、数据庞大,你不挖就是失去了一座金山。我们提供各种设计单细胞测序、各种测序、芯片、多组学的公共数据库挖掘、培训、模型构建、临床统计、算法还原服务;你能想到,我能做到;你提供参考文献、思路和目的,我们提供结果;如果没有思路,我们提供付费科研设计服务。示例如下:利用公共数据库的1539个单细胞样本,构建自己的生物学故事。 利用甲基化数据分析样本的拷贝数变异。
pancancer泛**图谱泛*研究是通过整合不同**类型、不同组织起源的**表达数据,查找**之间的共性或者差异的过程。通常使用**数据信息较为***的TCGA数据,通过分裂小提琴图展示某个基因在TCGA**和正常组织中的表达差异。分裂小提琴图(ViolinPlot)结合了箱形图和密度图的特征,主要用来显示数据的分布形状,它一般应用于对比某一基因在TCGA**组织和正常组织基因表达量TPM值或其它表达量数据。基本原理:小提琴图(ViolinPlot)使用一组数据中的最小值、**四分位数、中位数、第三四分位数和**值来反映数据分布的中心位置和散布范围,将多组数据的小提琴图画在同一坐标上,可以清晰地显示各组数据的分布差异。分裂小提琴图在小提琴图的基础上又加入了分组对比项,便于观察多**类型在某一基因上的表达分布情况,或者某一基因在某一**上,其疾病与正常的对比表达差异情况。 构建新的临床预测模型。上海生物/药物信息学分析数据科学口碑推荐
云生物数据分析需要多久?上海生物/药物信息学分析数据科学口碑推荐
术语解读:
TME: Tumormicroenvironment
TMEscore: TMEsignature score(使用PCA算法计算得到,高意味着对病毒和干扰素免疫***和应答敏感。)
PCA:Principal component analysis
CIBERSORT:Cell type identification by estimating relative subset of known RNA transcripts
CYT:Cytolytic activity
EMT:Epithelial-mesenchymal-transition
CR: Completeresponse
PR: Partialresponse
PD:Progressive disease
TMB: Tumormutational burden
数据要求:
各细胞之间的相关关系、pvalue、聚类/分类结果、跟预后的关系表。 上海生物/药物信息学分析数据科学口碑推荐