title: How to perform Post Hoc test, Tukey in agricolae package?
body: <p>I did this with</p>
  
  <pre><code>library(lsmeans)
  </code></pre>
  
  <p>and</p>
  
  <pre><code>library(multcomp)
  lm(Chlorophyll ~ Treatment + Stage + Treatment:Stage, "")
  </code></pre>
  
  <p>but I'm interested how to perform a post-hoc test after TW ANOVA in R. In <code>agricolae</code> package?</p>
  
  <pre><code>structure( list(Treatment = structure(c(3L, 3L, 3L, 3L, 3L, 3L,
      3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 
      1L, 1L, 1L, 1L, 1L), .Label = c("Control", "Nitrogen", "Salt"
      ), class = "factor"), Stage = structure(c(1L, 1L, 1L, 2L, 2L, 
      2L, 3L, 3L, 3L, 1L, 1L, 1L, 2L, 2L, 2L, 3L, 3L, 3L, 1L, 1L, 1L, 
      2L, 2L, 2L, 3L, 3L, 3L), .Label = c("Green", "Pink", "Red"), class = "factor"), 
     Chlorophyll = c(0.2, 0.3, 0.4, 0.5, 0.3, 0.2, 0.5, 0.6, 0.7, 
     0.4, 0.6, 0.9, 0.2, 0.3, 0.5, 0.4, 0.2, 0.3, 0.5, 0.6, 0.8, 
     0.5, 0.4, 0.6, 0.2, 0.3, 0.1)), .Names = c("Treatment", "Stage", 
     "Chlorophyll"), class = "data.frame", row.names = c(NA, -27L)
      )
  </code></pre>
  
  <p>After mean separation, How Can we plot the graphs for each stage with the group numbers/letters?</p>
  
  <p><a href="https://i.stack.imgur.com/QaOvo.png" rel="nofollow noreferrer"><img src="https://i.stack.imgur.com/QaOvo.png" alt="plot"></a></p>