title: Tukey HSD post-hoc test on a prop.test(..) in R
body: <p>I'm doing a prop.test() on 3 groups:</p>
  
  <pre><code>data = data.frame(a = c(85*0.3412,85*(1-0.3412)),b = c(203*0.2217,203*(1-0.2217)), c = c(231,197))
  prop.test(t(data))
  </code></pre>
  
  <p>The p.value is very low (1.041e-13). I want to know which pair of groups show a significant difference of proportion.</p>
  
  <p>I thought I could run the following code ....</p>
  
  <pre><code>library(agricolae)    
  HSD.test(prop.test(t(data)))
  </code></pre>
  
  <p>... but it is not working because HSD.test can only deal with aov or lm object.</p>
  
  <p><strong>What solution do I have left ?</strong> No matter if the solution implies a Bonferroni (or other) correction instead of a Tukey correction.</p>