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")