我在使用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的形式的选择。