VADeaths
barplot(VADeaths, main = "Grafico barras VADeaths", xlab = "Grupo", ylab = "Mortes", col = rainbow(5),legend.text = T, beside = F)

ClassificaçãoDoença
examinados <- c(8, 9, 3)
labels <- c("Leve", "Moderado", "Severo")
pct <- round(examinados / sum(examinados) * 100)
lbls <- paste(pct, "%", sep="")
pie(examinados, labels=lbls, main="Estágios da doença", col=rainbow(4))
legend("topleft",
legend = labels,
cex = 0.8,
fill = rainbow(length(examinados)))

Teorema
flu <- read.csv("flu.csv")
densidade <- density(flu$age)
hist(flu$age, probability = T, col = "grey", main = "Histograma da Idade com Densidade", xlab = "Idade")
lines(densidade, col = "blue")

n <- 200
sampleM <- 35
sample_means <- rep(NA, n)
for (i in 1:n) {
amostra <- sample(flu$age, size = sampleM)
sample_means[i] <- mean(amostra)
}
sample_means_density <- density(sample_means)
hist(sample_means, probability = T, col = "grey")
lines(sample_means_density, col = "blue")
