Prova - 2va

Respostas

VADeaths

data(VADeaths)
print(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

grafico_barras <- barplot(VADeaths, 
                         beside = TRUE,  
                         col = c("darkred", "darkgreen", "darkorange", "darkblue", "darkmagenta"),
                         main = "Taxas de Mortalidade em Virginia (1940)",
                         xlab = "Grupos",
                         ylab = "Taxas de Mortalidade (por 1000)")

legend("topright", 
       legend = rownames(VADeaths),
       fill = c("darkred", "darkgreen", "darkorange", "darkblue", "darkmagenta"),
       title = "Faixas Etárias")

ClassificaçãoDoença

status_doenca <- c("moderado", "leve", "leve", "severo", "leve", "moderado", "moderado", 
                  "moderado", "leve", "leve", "severo", "leve", "moderado", "moderado", 
                  "leve", "severo", "moderado", "moderado", "moderado", "leve")

frequencia <- table(status_doenca)
porcentagem <- round(100 * frequencia / length(status_doenca), 1)

paleta_cores <- c("leve" = "lightcyan", "moderado" = "mediumaquamarine", "severo" = "indianred")

pizza <- pie(frequencia, 
             main = "Classificação de Pacientes por Estágio da Doença",
             col = paleta_cores,
             labels = paste0(names(frequencia), " (", porcentagem, "%)"))

legend("bottomright", 
       legend = names(frequencia),
       fill = paleta_cores,
       title = "Estágios da Doença")

Teorema

dados_flu <- read.csv("flu.csv")

densidade_original <- density(dados_flu$age)
hist(dados_flu$age, 
     probability = TRUE, 
     col = "lavender", 
     border = "white",
     main = "Distribuição das Idades das Mortes - Gripe Espanhola", 
     xlab = "Idade", 
     ylab = "Densidade")
lines(densidade_original, col = "darkblue", lwd = 2)

set.seed(321) 
num_amostras <- 200
tamanho_amostral <- 35

medias_amostrais <- vector("numeric", num_amostras)

for(indice in seq_len(num_amostras)) {
  amostra_atual <- sample(dados_flu$age, size = tamanho_amostral, replace = TRUE)
  medias_amostrais[indice] <- mean(amostra_atual)
}

densidade_medias <- density(medias_amostrais)
hist(medias_amostrais, 
     probability = TRUE, 
     col = "lavender", 
     border = "white",
     main = "Distribuição das Médias Amostrais - Teorema do Limite Central", 
     xlab = "Média das Idades", 
     ylab = "Densidade")
lines(densidade_medias, col = "darkblue", lwd = 2)