文章预览
学员反馈说他在复现一个数据挖掘文章:《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
………………………………