gggenes icon indicating copy to clipboard operation
gggenes copied to clipboard

Gene length issue

Open chantisakee opened this issue 3 months ago • 1 comments

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. 01 02

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 03

`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

chantisakee avatar Apr 03 '24 03:04 chantisakee