#Innlesing av data, bruk av colClasses for å si fra hva slags type
#de ulike kollonnene er
d = read.table("exmp11-19.txt",header=TRUE,sep=",",
  colClasses=c("numeric","factor","factor"))
#Liste opp navnene på de ulike variable
names(d)
#Lage interaksjonsplott, to versjoner
par(mfrow=c(1,2))
interaction.plot(d$Source,d$Material,d$Vibrat)
interaction.plot(d$Material,d$Source,d$Vibrat)

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