remove(list=ls())
setwd("/Volumes/MichaelsHD/Dropbox/Angewandte Statistik SoSe11/Uebungen/Blatt10/R")
#setwd("/Volumes/Daten/Dropbox/Angewandte Statistik SoSe11/Uebungen/Blatt10/R")
set.seed(23514)

#1a
mensa <- read.table("./mensa.dat",header=T)
essen <- mensa$essen
anzahl <- mensa$anzahl
alle.aov <- aov(anzahl~essen)
alle.aov
summary(alle.aov)
anova(lm(anzahl~essen))

#1b
anzlf <- sum(essen=="LF")
anzgg <- sum(essen=="GG")

#teilstichprobe erstellen
subess <- factor(essen[(anzlf+anzgg+1):length(essen)],level=c("BIO","VEG"))
subanz <- anzahl[(anzlf+anzgg+1):length(essen)]

anovasub <- anova(lm(subanz ~ subess))
anovasub
aovsub <- aov(subanz ~ subess)
aovsub
summary(aovsub)
#Die Nullhypothese wird nicht verworfen

subfless <- factor(essen[1:(anzlf+anzgg)],level=c("GG","LF"))
subflanz <- c(anzahl[essen=="LF"],anzahl[essen=="GG"])
anovafleisch <- anova(lm(subflanz ~ subfless))
anovafleisch
aovfleisch <- aov(subflanz ~ subfless)
aovfleisch
summary(aovfleisch)
#Die Nullhypothese wird knapp nicht verworfen

#Boxplot der Daten
plot(anzahl~essen)