问题描述
我正在使用 gggenes 绘制基因箭头图,但是基因彼此相距很远,这使得该图不可读。不知道有没有什么方法可以删除不含基因的区域。
数据示例:
Genome start end state id gene
koreense -68765 -69505 reverse WP_165607793.1 dsrN
koreense -40289 -40585 reverse WP_069393961.1 sat
koreense -23366 -24769 reverse WP_084228866.1 dsrA
koreense -14754 -14957 reverse WP_085302191.1 sat
koreense -7192 -7686 reverse WP_069391389.1 dsrA
koreense -4858 -5154 reverse WP_069393961.1 sat
koreense -3072 -3554 reverse WP_085304442.1 sat
koreense -1258 -1503 reverse WP_085303636.1 sat
koreense 198 1013 forward WP_085304930.1 sat
koreense 29324 29524 forward WP_085304899.1 dsrB
koreense 33707 35119 forward WP_085304857.1 dsrN
koreense 152477 152875 forward WP_069392762.1 sat
我的代码是:
ggplot(gene,aes(xmin = start,xmax = end,y = Genome,fill = gene)) +
geom_gene_arrow() +
facet_wrap(~ Genome,scales = "free",ncol = 1) +
scale_fill_brewer(palette = "Set3") +
theme_genes() +
theme(axis.line = element_line())
示例数据:
enter code here
我得到的情节是:
我期望的是使用“//”来打破 x 轴,如:
提前致谢!
解决方法
暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!
如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@)