title: Modifying a Tukey HSD 95% family-wise CL plot in Rstudio
body: <p>I wanted to know how I could change the axis of a plot of Tukey's HSD so I could shorten the words to make each comparison fit and not look ridiculous. </p>
  
  <p>It would help a lot if you could help me with codes for changing axis labels, font size and color. </p>
  
  <p>There are quite a few comparisons in the plot, so I would like to make the significant ones (the comparisons whose intervals are not on the "0" line) stand out more by changing their color.</p>
  
  <pre><code>`Gastropods = read.csv(file = "MaroubraZones.csv", header = TRUE)
  boxplot(Abundance ~ Zone*Species,data = Gastropods, names = c("A.high", "A.mid", "A.low", "C.high", "C.mid", "C.low", "N.high", "N.mid", "N.low"))
  Gastropods.ANOVA = aov(Abundance ~ Zone * Species, data = Gastropods)
  hist(Gastropods.ANOVA$residuals)
  plot(Gastropods.ANOVA)
  
  
  Gastropods$LOGAbundance = log10(Gastropods$Abundance + 1)
  
  Gastropods$SQRTAbundance = sqrt(Gastropods$Abundance + 1)
  
  summary(Gastropods.ANOVA)
  summary(Gastropods$SQRTAbundance.ANOVA)
  
  interaction.plot(Gastropods$Zone, Gastropods$Species, Gastropods$Abundance, main= "Gastropod Interaction Plot", xlab = "Gastropod Zone", ylab= "Mean of Gastropod Abundance",legend = FALSE))
  
  interaction.plot(Gastropods$Zone, Gastropods$Species, Gastropods$Abundance, main= "Gastropod Interaction Plot", xlab = "Gastropod Zone", ylab= "Mean of Gastropod Abundance", legend = FALSE)                 
  
  
  TukeyHSD(Gastropods.ANOVA)
  tuk&lt;-TukeyHSD(Gastropods.ANOVA)
  plot(tuk)`
  </code></pre>
  
  <p>As you can see the axis's are terrible and I want to highlight the significant values outside of the zero interval.</p>
  
  <p><img src="https://i.stack.imgur.com/CPL9W.gif" alt="enter image description here"></p>