data <- c(40, 40, 5, 41, 41, 4, 41, 5, 42, 42, 43, 43, 45, 45, 45,
          45, 46, 47, 5, 48, 5, 48, 50, 50, 50, 51, 52, 52, 53, 55)

data.frame(row = data)

mean(data)

median(data)

factorx <- factor(cut(data, breaks=seq(0,max(data),5)))
factorx
xout <- as.data.frame(table(factorx))
xout

library("ggplot2")

ggplot(data.frame(row = data), aes(x = row)) + geom_freqpoly(binwidth = 5)


#last bsp
alpha <- 3
beta <- 1/30

data <- rgamma(40000,alpha,beta)
mean(data)
sd(data)

hist(data)

sum(data >= 80)/ 40000

#b
n <- 100
tmp <- sample(data,2000*n,replace = TRUE)
m <- matrix(tmp,nrow = 2000,byrow = TRUE)
yo <- rowMeans(m)
hist(yo)

qqnorm(yo)
qqline(yo)