1 Respostas das Questões 1,2 e 3

1.1 VADeaths

# Converter VADeaths para formato longo (tidy)
df <- as.data.frame(VADeaths) %>%
  mutate(GrupoEtario = rownames(VADeaths)) %>%
  pivot_longer(cols = -GrupoEtario,
               names_to = "Categoria",
               values_to = "Valor")
ggplot(df, aes(x = GrupoEtario, y = Valor, fill = Categoria)) +
  geom_bar(stat = "identity", position = "dodge") +
  labs(
    title = "Taxas de Mortalidade - VADeaths",
    x = "Grupos Etários",
    y = "Mortes por 1000 Habitantes",
    fill = "Categorias"
  ) +
  theme_minimal() +
  scale_fill_brewer(palette = "Set2")

1.2 ClassificaçãoDoença

# Classificações do estágio da doença
doencas <- c(
  "moderado", "leve", "leve", "severo", "leve", "moderado", "moderado",
  "moderado", "leve", "leve", "severo", "leve", "moderado", "moderado",
  "leve", "severo", "moderado", "moderado", "moderado", "leve"
)

# Tabela de frequências
freq <- table(doencas)

# Calculando as porcentagens
porcentagens <- round((freq / sum(freq)) * 100, 1)

# Cores para cada categoria
cores <- c("lightblue", "lightgreen", "tomato")

# Gráfico de pizza
pie(
  freq,
  labels = paste0(names(freq), " (", porcentagens, "%)"),
  col = cores,
  main = "Classificação da Doença - Gráfico de Pizza"
)

# Legenda
legend(
  "topright",
  legend = names(freq),
  fill = cores,
  title = "Categorias"
)

1.3 Teorema

# 1) Ler os dados
flu <- read.csv("flu.csv")  

# Extrair somente as idades
idades <- flu$age

# -------------------------------------------------------------
# (1) HISTOGRAMA + CURVA DE DENSIDADE DA POPULAÇÃO
# -------------------------------------------------------------

hist(
  idades,
  breaks = 20,
  probability = TRUE,
  col = "lightblue",
  main = "Histograma das Idades (População)",
  xlab = "Idade"
)

# Curva de densidade da população
lines(density(idades), col = "red", lwd = 2)

# -------------------------------------------------------------
# (2) CRIAR 200 MÉDIAS DE AMOSTRAS COM n = 35
# -------------------------------------------------------------

set.seed(123)

medias <- replicate(
  200,
  mean(sample(idades, size = 35, replace = TRUE))
)


# -------------------------------------------------------------
# (3) HISTOGRAMA + CURVA DE DENSIDADE DAS MÉDIAS
# -------------------------------------------------------------

hist(
  medias,
  breaks = 20,
  probability = TRUE,
  col = "lightgreen",
  main = "Histograma das Médias Amostrais (n = 35)",
  xlab = "Média da Amostra"
)

# Curva de densidade das médias
lines(density(medias), col = "blue", lwd = 2)