[R] How to embed residual dot plots of BHH2::anovaPlot(...) into plotMeans(...)?

[Xiaoxu LI]

Dear R fans,

I am preparing a lecture discussing the paper of Wilkinson and APA
Task Force on Statistical Inference (1999, American Psychologist, 54,
594-604.) I'd like to demo their Figure 3B (p. 602) with R. That is to
add plotMeans(...) with dots of BHH2::dotPlot(...) . Is there any
elegant script to do the job?



Thanks

---------------
LI, Xiaoxu