library(data.table)
strsdata = fread("distance_travelled.csv")

library(ggplot2)
require(ggpubr)
library(ggpubr)
my_comparisons <- list( c("Treatment", "Control-First"), c("Treatment", "Trial"), c("Control", "Treatment"), c("Trial","Control"))

p <- ggboxplot(strsdata, x = "state", y = "ddiff",
               color = "state", palette = "jco",
               add = "jitter", xlab = "Distance Travelled", ylab = "Group")
means <- aggregate(ddiff ~  state, strsdata, median)
means$ddiff = round(means$ddiff, 2)
# Change method
p + stat_compare_means(comparisons = my_comparisons, label = "p.signif")+ 
  geom_text(data = means, aes(label = ddiff, y = 20))


strsdata = fread("radius_of_gyration.csv")
p <- ggboxplot(strsdata, x = "state", y = "freq_rog",
               color = "state", palette = "jco",
               add = "jitter", xlab = "Radius of Gyration", ylab = "Group")
means <- aggregate(freq_rog ~  state, strsdata, median)
means$freq_rog = round(means$freq_rog, 2)
# Change method
p + stat_compare_means(comparisons = my_comparisons, label = "p.signif")+ 
  geom_text(data = means, aes(label = freq_rog, y = 20))