Hello, Thank you for the script and the protocol. I tried the...

Alice S.
Alice S.
Nov 17, 2021
protocol Protocol: A Simple Protocol for Informative Visualization of Enriched Gene Ontology Terms
Hello,
Thank you for the script and the protocol. I tried the script with my own data and it worked out well. I had to remove lines 53 and 54 though. I tried to figure out how I could change the script in order to get (1, 2, 3, 4, 5 and 6) on the x-axis instead of (1, 2, 4 and 6). I
Read more Read less more less
answer Answer
recommend recommend recommend recommend Recommend
follow follow follow follow Follow
share share Share
1 answer
Sort by: most helpful / most recent

Titouan Bonnot Author Answered Nov 19, 2021

INRAE - Dijon, Agroecology

Hello Alice S.,
I’m glad to hear that our protocol worked well with your own data.
You can try to add the following code at the end of l.40 of the R script :
scale_y_continuous(breaks = seq(from = 0 , to = 6, by = 1))+
so l.40 of the script will look like :
scale_x_discrete(limits= GO_all$GO_biological_process)+scale_y_continuous(breaks = seq(from = 0 , to = 6, by = 1))+

Please note that the axes are flipped with our script, the x-axis corresponds to the GO terms and the y-axis corresponds to the fold enrichment.
I hope that helps!

Titouan
0 helpful
0 unhelpful
comment 0 comments down up
share Share

My answer

Write your answer...

References (optional)

add Add more

post Post an Answer