source("http://www.matsuka.info/univ/course_folder/cutil.R")
dat<-read.csv("http://www.matsuka.info/data_folder/dktb316.txt",
   col.names = c("dump","method","subj","words"))
dat.aov<-aov(words~method+subj+Error(subj:method),data=dat)

dat<-read.csv("http://www.matsuka.info/data_folder/dktb3211.txt")
interaction.plot(dat$duration,
                 dat$method,    　　　　
                 dat$result,    
                 pch=c(20,20), 
                 col=c("skyblue","orange"),
                 ylab="score",
                 xlab="Duration",
                 lwd=3,　
                 type='b',
                 cex=2, 
                 trace.label="Method")
dat.aov<-aov(result~method*duration+Error(s+s:method),dat)

dat<-read.csv("http://www.matsuka.info/data_folder/dktb3218.txt")

dat.aov<-aov(result~method+duration+method:duration + 
  Error(s+method:s+duration:s+method:duration:s), data=dat)