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