5 转录组自主分析的差异基因分析和KEGG富集分析时发生错误,求老师解答。

在进行有参转录组自主学习的《差异表达基因功能富集分析》这节内容的学习时,发现KEGG富集分析中断,报错代码如下:

Warning message:

In bitr(DEG_list, fromType = opt$idtype, toType = opt$totype, OrgDb = opt$ann.db) :

  26.02% of input gene IDs are fail to map...

--> No gene can be mapped....

--> Expected input gene ID:

--> return NULL...

NULL

Error in (function (classes, fdef, mtable)  :

  unable to find an inherited method for function ‘dotplot’ for signature ‘"NULL"’

Calls: dotplot -> <Anonymous>

Execution halted

查询问题时候发现,在上一节内容中,练习生成的normal_vs_tumor.DEG.final.tsv 文件与demo中的基因数量不同,demo中319个基因,而自己练习生成的数据是92个,为了排除bug原因,尝试使用demo中normal_vs_tumor.DEG.final.tsv 进行KEGG富集,发现依然如上报错。使用chatgpt查询代码错误,提示检查输入ID和处理未映射的基因ID。查询ID后发现应该没有错误。目前僵持在这一步,求老师解答!
查询发现可能是2023年kegg数据库更新,导致的No gene can be mapped。具体原因还望老师指点!

请先 登录 后评论

最佳答案 2023-11-11 15:53

可能是kegg数据库更新所致,已提交更新计划

请先 登录 后评论

其它 2 个回答

gongqinglong

使用的代码是

Rscript $scriptdir/enrichKEGG_pip.R --deg.file $workdir/5.deg/normal_vs_tumor.DEG.final.tsv -o KEGG -n normal_vs_tumor --pvalueCutoff 1 --ann.db org.Hs.eg.db --organism hsa --idtype ENSEMBL --totype ENTREZID 


请先 登录 后评论
omicsgene - 生物信息
擅长:重测序,遗传进化,转录组,GWAS

你是用的我们提供的docker镜像吗?在镜像里面分析数据我们都经过测试没问题的;

如果没有使用我们的镜像,是你你自己安装的R包有问题,需要更新这个脚本需要的R包到最新:enrichKEGG_pip.R

请先 登录 后评论