[R] More basic equivalent of TukeyHSD

[Janet]

This is a simple question, but I'm going on the supposition that the  
only stupid question is the one not asked.

1.  I have many sets of 5 proportions that are different from each  
other (prop.test), and want to know which proportions are different  
from each other.  In other words, I want the equivalent of Tukey's HSD  
test, but for proportions rather than anova.

Here is one of the sets of number of "successes"; the denominator is  
151 for all cases.
                     Parents School Peers Church No one
Menstruation             79     58    55      3     38

As proportions:
 > round(topics.sum[1,], digits=3)
Parents  School   Peers  Church  No one
   0.523   0.384   0.364   0.020   0.252

Intuitively, it looks like there are 3 or 4 groups that are different  
from each other, but is there a function that does this?
parents > school and peers ?> no one > church

2.  I am displaying this data in a bar graph.  How do I show  
graphically which proportions are different?  I'm not interested in  
using error bars since they can be ambiguous, but probably something  
like putting a different symbol over each bar to show which bars are  
the same and which differ.  I've made efforts at this before, but  
they've always looked a bit clumsy.

Thanks,

Janet