专栏名称: 生信技能树
生物信息学学习资料分析,常见数据格式及公共数据库资料分享。常见分析软件及流程,基因检测及癌症相关动态。
今天看啥  ›  专栏  ›  生信技能树

想富集哪个通路就富集哪个通路?

生信技能树  · 公众号  ·  · 2024-10-07 11:24

文章预览

学员反馈说他在复现一个数据挖掘文章:《Diagnostic and Predictive Values of Ferroptosis-Related Genes in Child Sepsis》,但是同样的差异分析后的上下调基因的富集分析结果里面的铁死亡这个通路并没有统计学显著性,因为没有出现在他可视化的图里面。所以我就引导他去看背后的表。 我给出来的代码如下所示:   kk.up                       organism     =  'hsa' ,                        #universe     = gene_all,                       pvalueCutoff =  0.9 ,                       qvalueCutoff = 0.9 )   head(kk.up)[, 1 : 6 ]   dotplot(kk.up)    kk.up=DOSE::setReadable(kk.up, OrgDb= 'org.Hs.eg.db' ,                        keyType= 'ENTREZID' ) #按需替换  tmp=kk.up@result   tmp   head(  tmp[, 4 : 7 ]) #  write.csv(kk.up@result, file = paste0(path,'/','KEGG_up ………………………………

原文地址:访问原文地址
快照地址: 访问文章快照
总结与预览地址:访问总结与预览