-
Notifications
You must be signed in to change notification settings - Fork 50
Description
Hi @wilkox ,
Thank you for this amazing package. It's really helpful :)
I have a problem with the visualization of gene length (group_506) as I have attached the files.
Pic_01 is the plot that i need, and Pic_02 is what i got. I'm quite not sure whether anything i can do to increase the length of group_506.
this is the code i used:
ggplot(input, aes(xmin = start, xmax = end, y = molecule, fill = gene)) + geom_gene_arrow() + facet_wrap(~ molecule, scales = "free", ncol = 1) + scale_fill_brewer(palette = "Set3") + theme_genes()
this is my dataset:
molecule gene start end orientation
1 NC_006142.1 group_506 906635 906880 TRUE
2 NC_006142.1 group_505 841537 841782 TRUE
3 B9991CWPP group_506 908080 908325 TRUE
4 B9991CWPP group_505 842201 842446 TRUE
5 TH1527 group_506 907512 907757 TRUE
6 TH1527 group_505 841633 841878 TRUE
And when i tried to use dummies to align by group_505, Pic_03 is the plot i got
`dummies <- make_alignment_dummies(
input,
aes(xmin = start, xmax = end, y = molecule, id = gene),
on = "group_505"
)
ggplot(input, aes(xmin = start, xmax = end, y =
molecule, fill = gene, label = gene)) +
geom_gene_arrow(arrowhead_height = unit(3, "mm"), arrowhead_width = unit(1, "mm")) +
geom_gene_label(align = "left") +
geom_blank(data = dummies) +
facet_wrap(~ molecule, scales = "free", ncol = 1) +
scale_fill_brewer(palette = "Set3") +
theme_genes()`
Thank you