setwd("C:/Users/Daniel Westhoff/Documents/Lehre/Angewandte Stochastik II/Übungsblätter/Blatt4")
#2)
#a)
####
s1 = 0
s2 = 0
s3 = 0
for(i in 1:100) {
  # i)
  data = rexp(20, 5)
  # ii)
  s1[i] = (mean(data) - .2)^2
  s2[i] = (sqrt(var(data)) - .2)^2
  s3[i] = (20 * min(data) - .2)^2
}

# iii)
# zuerst alle 3
boxplot(s1, s2, s3)

# 1 und 2 zur besseren Unterscheidung
boxplot(s1, s2)

#b)
lebensdauer <- read.table("lebensdauer.txt", header = TRUE)

means = sapply(lebensdauer, mean)

means[means==max(means)]