aspettl — May 10, 2013, 10:49 AM
load("tatort2013.rda")
# a)
hist(tatort2013$Schaetzung$personC, breaks=seq(25,55,3), freq=F, main="Histogramm Altersschätzungen Person C", xlab="Alter", ylab="relative Häufigkeiten")
# Hinweis: breaks=10 liefert falsche Klassenanzahl, da 10 nur ein Vorschlag an R ist
# b)
mean <- mean(tatort2013$Schaetzung$personC)
var <- var(tatort2013$Schaetzung$personC)
print(data.frame(mean=mean, var=var))
mean var
1 43.16 27.16
# c)
curve(dnorm(x, mean, sqrt(var)), col="darkgreen", lwd=2, add=T)
# Achtung: dnorm erwartet Standardabweichung, nicht Varianz
# (stärkere Liniendicke nur für bessere Sichtbarkeit)
abline(v=tatort2013$Alter[3], col="red", lwd=3)