d = read.table("exmp11-08.txt",header=TRUE,sep=",")
d2 = data.frame(fac=as.factor(c(rep("Perm",8),rep("Die",8),rep("Pla",6))))
d2$y = c(d$Perm,d$Die,d$Pla[1:6])  #Tar ikke med NA
fit = aov(y~fac,data=d2)
summary(fit)
library(car)
TukeyHSD(fit)
plot(TukeyHSD(fit))
abline(v=0)