使用pathview時(shí)遇到一個(gè)error投蝉,初始代碼如下:
pathview(gene.data = df_gene_KO$KEGG_ko,
pathway.id = mapID,
species = "ko",
limit = list(gene=1, cpd=1),
out.suffix = "suffix")
部分pathway可以正常生成,如ko00970
但是部分pathway無法生成征堪,報(bào)錯(cuò)瘩缆,如
ko01230
仔細(xì)對比發(fā)現(xiàn),這兩幅圖的區(qū)別在于:
ko00970
圖中基因Kxxxxx是矩形方框佃蚜,所以可以添加顏色庸娱;ko01230
中基因Kxxxxx是線條,無法添加顏色
解決辦法:
不要將顏色添加在同一圖層爽锥,修改參數(shù)same.layer = T
涌韩,代碼如下:
pathview(gene.data = df_gene_KO$KEGG_ko,
pathway.id = mapID,
species = "ko",
limit = list(gene=1, cpd=1),
out.suffix = "suffix",
same.layer = T)
生成圖形: