We will incorporate a version of this in R-devel and hence 2.2.0.
A bit more work needs to be done on printing: whereas the existing three
columns are in comparable units, the p values are not and so get printed
to silly accuracies (as can be seen in your example).
You also need to explain what
Dear R-developeRs,
Attached follows a patch against svn 34959 that adds the
printing of p-values to the TukeyHSD.aov function in stats package. I
also updated the corresponding documentation file and added a 'see also'
reference to the simint function of the multcomp package.