Resultados
ClassificaçãoDoença
dados_pacientes <- c("moderado", "leve", "leve", "severo", "leve", "moderado",
"moderado", "moderado", "leve", "leve", "severo", "leve",
"moderado", "moderado", "leve", "severo", "moderado",
"moderado", "moderado", "leve")
tabela <- table(dados_pacientes)
porcentagens <- round(100 * tabela / sum(tabela), 1)
rotulos <- paste(porcentagens, "%", sep="")
cores <- c("green", "blue", "salmon")
pie(tabela,
labels = rotulos,
col = cores,
main = "Classificação dos Estágios da Doença, ou seja n=20")
legend("topright",
legend = names(tabela),
fill = cores,
title = "Estágio",
cex = 0.8)

Teorema
# Carrega o dataset "flu" (Gripe Espanhola na Suíça, 1918)
flu <- read.csv('C:/Users/Merti/Downloads/flu.csv')
table(flu$age)
##
## 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
## 6409 1307 723 504 445 368 368 300 248 271 233 270 238 281 335 396
## 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31
## 544 620 696 840 923 944 986 1062 1076 1117 1105 1144 1166 1100 1075 1027
## 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47
## 992 997 832 877 872 785 733 757 751 689 727 738 598 622 629 617
## 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63
## 697 574 591 661 656 745 794 791 792 814 799 849 848 947 909 839
## 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79
## 891 994 966 1039 1135 1154 1054 1190 1229 1221 1171 1189 1112 1030 1039 922
## 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95
## 753 727 677 579 501 392 354 299 224 144 133 100 56 43 30 16
## 96 97 98 99 100
## 11 6 6 3 1
# Histograma e curva de densidade
hist(flu$age,
freq = FALSE,
col = "lavender",
main = "Idades das Mortes",
xlab = "Idade",
breaks = 30)
lines(density(flu$age, na.rm=TRUE), col = "red", lwd = 2)

# Média por Amostragem - 200 médias com n = 35
n_amostras <- 200
tam_amostra <- 35
medias_amostrais <- numeric(n_amostras)
set.seed(72)
for(i in 1:n_amostras) {
amostra <- sample(flu$age, size = tam_amostra, replace = FALSE)
medias_amostrais[i] <- mean(amostra)
}
# Histograma e curva de densidade das médias das amostras
hist(medias_amostrais,
probability = TRUE,
col = "green",
main = "Médias Amostras n=35",
xlab = "Média da Idade",
breaks = 20)
lines(density(medias_amostrais), col = "black", lwd = 2)

VADeaths
head(VADeaths)
## Rural Male Rural Female Urban Male Urban Female
## 50-54 11.7 8.7 15.4 8.4
## 55-59 18.1 11.7 24.3 13.6
## 60-64 26.9 20.3 37.0 19.3
## 65-69 41.0 30.9 54.6 35.1
## 70-74 66.0 54.3 71.1 50.0
# Definindo as cores dos grupos
cores <- c("pink", "orange", "lavender", "purple", "cyan")
barplot(height = VADeaths,
beside = TRUE,
col = cores,
legend.text = TRUE,
main = "Taxas de Mortalidade na Virgínia (1940)",
xlab = "Grupos Demográficos",
ylab = "Taxa de Mortalidade (por 1000)",
ylim = c(0, 100),
args.legend = list(title = "Faixa Etária", x = "topright")
)
