pulse <- read.table("pulse.dat",header=T)
pulse$task <- factor(pulse$task)
pulse.fit <- aov(rate ~ task,pulse)