pdf("C:\\Rmisc\\graphs\\caffeine1.pdf") caffeine <- read.fwf("http://www.stat.ufl.edu/~winner/data/caffeine.dat", width=c(1,5,8), col.names=c("subject_c", "dose_c", "endtime")) attach(caffeine) dose <- factor(dose_c) subject_c <- factor(subject_c) caffeine1.rbd <- aov(endtime ~ subject_c + dose_c) summary(caffeine1.rbd) TukeyHSD(caffeine1.rbd, "dose_c") interaction.plot(subject_c,dose_c,endtime) caffeine.res <- residuals(caffeine1.rbd) hist(caffeine.res) friedman.test(endtime ~ dose_c | subject_c) dev.off()