GSEA分析:GSEA全名为GeneSetEnrichmentAnalysis(基因集富集分析)。用以分析特定基因集(如关注的GO条目或KEGGPathway)在两个生物学状态(如**与对照,高龄与低龄)中是否存在差异。能够研究基因变化的生物学意义。普通GO/KEGG富集的思路是先筛选差异基因,然后确定这些差异基因的GO/KEGG注释,然后通过超几何分布计算出哪些通路富集到了,再通过p值或FDR等阈值进行筛选。挑选用于富集的基因有一定的主观性,没有关注到的基因的信息会被忽视,所以有一定的局限性。在这种情况下有了GSEA(GeneSetEnrichmentAnalysis),其思路是发表于2005年的Genesetenrichmentanalysis:aknowledge-basedapproachforinterpretinggenome-wideexpressionprofiles。主要是要有两个概念:预先定义的基因集S(基于先验知识的基因注释信息)和待分析基因集L(一般初始输入是表达矩阵);然后GSEA目的就是为了判断S基因集中的基因是随机分布于L(按差异表达程度对基因进行排序),还是聚集分布在L的顶部或者底部(也就是存在差异性富集)。如果基因集中的基因***富集在L的顶部或者底部,这说明这些基因的表达对定义的分组(预先分组)的差异有***影响(一致性)。在富集分析的理论中。 可对接各类公共数据库,切入各类接口,并对公共数据库进行大规模数据挖掘。湖北组学数据处理数据科学怎么样
GSEA基本原理从方法上来讲,GSEA主要分为基因集进行排序、计算富集分数(EnrichmentScore,ES)、估计富集分数的***性水平并进行多重假设检验三个步骤。**步对输入的所有基因集L进行排序,通常来说初始输入的基因数据为表达矩阵,排序的过程相当于特定两组中(case-control、upper-lower等等)基因差异表达分析的过程。根据所有基因在两组样本的差异度量不同(共有六种差异度量,默认是signal2noise,GSEA官网有提供公式,也可以选择较为普遍的foldchange),对基因进行排序,并且Z-score标准化。第二步是GSEA的**步骤,通过分析预先定义基因集S在**步获得的基因序列上的分布计算富集指数EnrichmentScore,并绘制分布趋势图Enrichmentplot。每个基因在基因集S的EnrichmentScore取决于这个基因是否属于基因集S及其差异度量(如foldchange)。差异度量越大基因的EnrichmentScore权重越大,如果基因在基因集S中则EnrichmentScore取正,反则取负。将基因集L在基因集S里的所有基因的EnrichmentScore一个个加起来,就是Enrichmentplot上的EnrichmentScore趋势,直到EnrichmentScore达到**值,就是基因集S**终的EnrichmentScore。第三步是为了检验第二部获得结果的统计学意义。 重庆数据科学售后服务云生物提供数据科学服务。
t-SNE(t分布随机邻域嵌入)是一种用于探索高维数据的非线性降维算法。它将多维数据映射到适合于人类观察的两个或多个维度。t-SNE非线性降维算法通过基于具有多个特征的数据点的相似性识别观察到的簇来在数据中找到模式。另外t-SNE的输出可以作为其他分类算法的输入特征。因为t-SNE算法定义了数据的局部和全局结构之间的软边界。t-SNE几乎可用于所有高维数据集,广泛应用于图像处理,自然语言处理和语音处理。在生物信息中可广泛应用于基因表达数据、基因甲基化数据、基因突变数据等,能够直观地对不同数据集进行比较。基本原理从方法上来讲,t-SNE本质上是基于流行学习(manifoldlearning)的降维算法,不同于传统的PCA和MMD等方法,t-SNE在高维用normalizedGaussiankernel对数据点对进行相似性建模。相应的,在低维用t分布对数据点对进行相似性(直观上的距离)建模,然后用KL距离来拉近高维和低维空间中的距离分布。
PPImodule蛋白质互作蛋白质-蛋白质相互作用(protein-proteininteraction,PPI)是指两个或两个以上的蛋白质分子通过非共价键形成蛋白质复合体(proteincomplex)的过程。PPImodule是指共表达蛋白模块或蛋白质相互作用模块。蛋白质相互作用形成人体复杂的蛋白质相互作用网络,对蛋白质相互作用网络进行聚类形成模块从而帮助我们理解细胞的功能。我们一般使用PPImodule把基因列表跟蛋白相互作用网络联系起来。例如RNA-seq获得的差异表达基因,看他们在蛋白相互作用网络中,哪些基因处于同一module。基本原理:蛋白质在细胞中的功能取决于它与其他蛋白质、核酸和小分子相互作用关系,对蛋白质相互作用网络进行聚类形成模块,各个蛋白模块发挥不同的功能,我们将基因列表重叠于模块上,查找基因列表所在的功能模块,从而发现基因列表中的基因可能发挥的细胞功能。我们通过PPI数据库找到共表达蛋白中的module,然后从模块中筛选出基因列表的产物蛋白,筛选出的结果就是基因列表***表达的PPImodule。 检测服务及数据分析助力取得2020年国自然面上十项、青年基金十八项。
LASSO回归:更多的变量在拟合时往往可以给出一个看似更好的模型,但是同时也面临过度拟合的危险。此时如果用全新的数据去验证模型(Validation),通常效果很差。一般来说,变量数大于数据点数量很多,或者某一个离散变量有太多独特值时,都有可能过度拟合。LASSO回归复杂度调整的程度由参数λ来控制,λ越大对变量较多的线性模型的惩罚力度就越大,从而**终获得一个变量较少的模型。LASSO回归与Ridge回归同属于一个被称为ElasticNet的广义线性模型家族。这一家族的模型除了相同作用的参数λ之外,还有另一个参数α来控制应对高相关性(highlycorrelated)数据时模型的性状。LASSO回归α=1,Ridge回归α=0,一般ElasticNet模型0<α<1。LASSO过程中我们通常会进行多次交叉验证(crossvalidation)拟合(1000次)进而选取模型,从而对模型的性能有一个更准确的估计。 根据委托方提供的参考文献和要求进行个性化特定分析。诊疗软件开发数据科学活动
微生物多样性分析桌面软件。湖北组学数据处理数据科学怎么样
下游分析针对LASSO获得的基因模型(或称基因Panel)的验证:1.计算风险指数RiskScore2.绘制ROC曲线、DCA曲线、列线图进行验证3.绘制生KM存曲线对基因模型中的基因进行解释和分析:1.基因注释2.靶向药物分析应用示例:文献1:PrognosticandpredictivevalueofamicroRNAsignatureinstageIIcoloncancer:amicroRNAexpressionanalysis.于2013年12月发表在LancetOncol.,影响因子。一个miRNA特征集在stageII结肠*的预后预测作用分析文章对stageII结肠*组织和*旁正常组织的miRNA芯片数据进行了差异表达分析,并通过LASSOCox回归对获得的差异表达miRNA进行筛选,获得了6个miRNA的可以预测预后情况的miRNA特征集。文献2:PrognosticValueofaBCSC-associatedMicroRNASignatureinHormoneReceptor-PositiveHER2-NegativeBreastCancer(于2016年9月发表在EBioMedicine.上,影响因子)文章将符合条件的患者划分为训练集和测试集,首先分析获得了**干细胞相关的miRNA,接着通过LASSO对**干细胞相关的miRNA进行筛选,构建了10个miRNA的预后预测模型,并计算风险指数绘制了生存曲线和ROC曲线。 湖北组学数据处理数据科学怎么样