#Innlesing av data, bruk av colClasses for aa si fra hva slags type
#de ulike kolonnene er
d = read.table("exmp11-16.txt",header=TRUE,sep=",",
  colClasses=c("numeric","factor","factor"))
#Liste opp navnene paa de ulike variable
names(d)
#Lage interaksjonsplott, se help(interaction.plot)
interaction.plot(d$Variety,d$Density,d$Yield)

#Utfoere variansanalyse, bruker symbolet : for interaksjon
fit = aov(Yield~Density+Variety+Density:Variety,data=d)
#Utskrift av anova tabell
summary(fit)
#Ulike residualplot
plot(fit)

#Se paa individuelle forskjeller ved Tukey's metode
TukeyHSD(fit,which="Variety")
plot(TukeyHSD(fit,which="Variety"));abline(v=0)