library(asbio)
data(shad)

tapply(shad$length, shad$site, mean)
tapply(shad$length, shad$site, median)
tapply(shad$length, shad$site, sd)

library(MASS)
data(topo)

qqnorm(topo$z); qqline(topo$z)

averager <- function(a) {
  X <- runif(n=10000, min = 0, max = a)
  X <- matrix(X, nrow=10, ncol=1000)
  return(X)
}

av = averager(1)

plot(density(apply(av, 2, max)))

print(paste0("(", as.character(max(av[,1])), " ",as.character(max(av[,1]) / 0.05^(1/10)), ")"))

counter = 0;
for (i in 1:1000){
  if (max(av[,i]) <= 1 && (max(av[,i])/(0.05^(1/10))) >= 1) {
    counter = counter + 1
  }
}
print(counter)