则则147 2022-11-03 22:30 采纳率: 100%
浏览 54
已结题

GOplot报错value = character(0)

在用GOplot画富集分析图的时候,出现错误

Error in $<-.data.frame(*tmp*, "genes", value = character(0)) :
替换数据里有0行,但数据有811

deg_genelist=read.table("deg_genelist.txt")  #读入差异分析结果
GO_Enrichment=read.table("data/clusterProfiler_GO_result.txt",
                         sep="\t",header = T,check.names = F)  #读入GO富集分析结果
circ=circle_dat(GO_Enrichment,deg_genelist)  #合并上述两个结果
运行结果及报错内容
Error in `$<-.data.frame`(`*tmp*`, "genes", value = character(0)) : 
  替换数据里有0行,但数据有811
我的解答思路和尝试过的方法

有人说是geneID里面基因间不可用/分隔,得用, 分隔,我试了一下
GO_Enrichment$geneID=gsub("/",", ",GO_Enrichment$geneID)
但是还是报错

而且同样的数据我用clusterFiler画图是可以画出来的

我想要达到的结果

想请问这个报错该怎么解决?

  • 写回答

1条回答 默认 最新

  • Bioinfo Guy R语言领域新星创作者 2022-11-04 09:35
    关注

    检查一下deg_genelist和GO_Enrichment是否读入成功了,我看你差异结果导入又没有读取表头,其次GO富集结果要这5列数据都齐

    circle_dat(terms, genes)
    #terms:A data frame with columns for 'category', 'ID', 'term', adjusted p-value ('adj_pval') and 'genes'
    #genes:A data frame with columns for 'ID', 'logFC'
    
    本回答被题主选为最佳回答 , 对您是否有帮助呢?
    评论

报告相同问题?

问题事件

  • 系统已结题 1月31日
  • 已采纳回答 1月23日
  • 创建了问题 11月3日

悬赏问题

  • ¥15 51寻迹小车定点寻迹
  • ¥15 爬虫爬取网站的一些信息
  • ¥15 关于vue2中methods使用call修改this指向的问题
  • ¥15 idea自动补全键位冲突
  • ¥15 请教一下写代码,代码好难
  • ¥15 iis10中如何阻止别人网站重定向到我的网站
  • ¥15 滑块验证码移动速度不一致问题
  • ¥15 Utunbu中vscode下cern root工作台中写的程序root的头文件无法包含
  • ¥15 麒麟V10桌面版SP1如何配置bonding
  • ¥15 Marscode IDE 如何预览新建的 HTML 文件