title: P values from HSD.test
body: <p>I am using the <code>HSD.test</code> function from the <code>agricolae</code> package in <code>R</code> to conduct a TukeyHSD post hoc test. The function works fine except I'm not sure where the p values are hidden. The letters in <code>groups</code> denote significance, but where are the actual p values? Thank you</p>
  
  <pre><code>library(agricolae)
  data(sweetpotato)
  model&lt;-aov(yield~virus, data=sweetpotato)
  out &lt;- HSD.test(model,"virus", group=TRUE,console=TRUE,
  main="Yield of sweetpotato\nDealt with different virus")
  
  
  
  Study: Yield of sweetpotato
  Dealt with different virus
  
  HSD Test for yield 
  
  Mean Square Error:  22.48917 
  
  virus,  means
  
        yield      std r  Min  Max
  cc 24.40000 3.609709 3 21.7 28.5
  fc 12.86667 2.159475 3 10.6 14.9
  ff 36.33333 7.333030 3 28.0 41.8
  oo 36.90000 4.300000 3 32.1 40.4
  
  Alpha: 0.05 ; DF Error: 8 
  Critical Value of Studentized Range: 4.52881 
  
  Minimun Significant Difference: 12.39967 
  
  Treatments with the same letter are not significantly different.
  
        yield groups
  oo 36.90000      a
  ff 36.33333     ab
  cc 24.40000     bc
  fc 12.86667      c
  </code></pre>