VADeaths
barplot(VADeaths, main = "Mortes a cada mil.",
xlab = "Grupos", ylab = "Mortes", col = rainbow(5),
beside = F, legend.text = T)

ClassificaçãoDoença
x <- c(8, 9, 3)
labels <- c("leve", "moderado", "severo")
pct <- round(x / sum(x) * 100)
lbls <- paste(pct, "%", sep = "")
pie(x, labels = lbls, col = rainbow(4))
legend("topright", legend = labels, cex = 0.8, fill = rainbow(length(x)))

Teorema
flu <- read.csv("flu.csv")
density_flu <- density(flu$age)
hist(flu$age, probability = TRUE, col = "grey")
lines(density_flu, col = "red")

n <- 200
media <- 35
xbar <- rep(NA, n)
for(i in 1:n){
amostra <- sample(flu$age, size = media)
xbar[i] <- mean(amostra)
}
density_xbar <- density(xbar)
hist(xbar, probability = TRUE, col = "grey")
lines(density_xbar, col = "red")
