关于Error in `[<-`(`*tmp*`, g, p, value = ifelse(M[g] %in% sub2$genes, 1,: subscript out of bounds报错

我在使用GOplot包的函数chord_dat时发生了报错

我在跑chord <- chord_dat(circ, genelist, go$Term)命令时,报错如下:

Error in `[<-`(`*tmp*`, g, p, value = ifelse(M[g] %in% sub2$genes, 1, :

  subscript out of bounds

我去检查了go和genelist的数据结构发现,genelist里的gene用的是gene名,而go里的基因用的是基因ID,不一样了,所以跑不出结果,所以我把genelist的gene换成了基因ID,就能跑出来了。

我发现在运行kk计算的时候,会随机在基因名和基因ID之间跳转,有时候会生成geneID,有时候又是基因名,所以要根据go生成的数据里生成的是基因名还是基因ID,来调整genelist的gene的形式的选择。

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

推荐阅读更多精彩内容