stat_compare_means compare groups not located on x-axis

Question

I\'m stuck on how to plot the mean p-values for grouped boxplots.

This is a portion of my data:

ID<-c(\'E5b\',\'R6\',"S22","E5b",&qu