Dear list, for an ANOVA analysis I would like to obtain the interaction plots that contains also the confidence intervals (95% Tukey HSD). My anova model is built in a way similar to this: myObj.aov <- aov(response~(factor1+factor2+factor2+factor4+factor5) ^2,data=myDataFrame) To obtain the interaction plots I use the code below for each factor of my model: library(alr3) interaction.plot(myDataFrame$factor1,myDataFrame$factor2,myDataFrame $response,fixed=TRUE,xpd=TRUE,main="Interaction Factor1 & Factor 2",xlab="Factor 1",ylab="response") How could I plot also the confidence intervals? Thanks for any help you might provide. -------------------------------------------------------------------- Max MANFRIN Tel.: +32 (0)2 650 3168 IRIDIA - CoDE, CP 194/6 Fax.: +32 (0)2 650 2715 Universit? Libre de Bruxelles Av. F. D. Roosevelt, 50 1050 Brussels Email: mmanfrin at ulb.ac.be BELGIUM WWW: http://iridia.ulb.ac.be/~mmanfrin gpg DSA: 0x7E67B4C4 gpg fingerprint: C2E9 1689 CADD 7CAE 2FAB A355 8FD9 9DD1 7E67 B4C4 -------------------------------------------------------------------- -------------- next part -------------- A non-text attachment was scrubbed... Name: PGP.sig Type: application/pgp-signature Size: 194 bytes Desc: This is a digitally signed message part Url : https://stat.ethz.ch/pipermail/r-help/attachments/20070402/47537e9e/attachment.bin