VADeaths

# Carregando o dataset
data(VADeaths)

# Cores para cada grupo (Rural Male, Rural Female, Urban Male, Urban Female)
cores <- c(
  "Rural Male"   = "#E63946",
  "Rural Female" = "#457B9D",
  "Urban Male"   = "#2A9D8F",
  "Urban Female" = "#E9C46A"
)

# Transpondo a matriz para que as faixas etárias sejam as categorias do eixo X
# e os grupos (Rural/Urban x Male/Female) sejam as barras lado a lado
dados <- t(VADeaths)

# Criando o gráfico de barras agrupadas
barplot(
  dados,
  beside      = TRUE,
  col         = cores,
  main        = "Taxa de Mortalidade na Virgínia (1940) por Faixa Etária e Grupo",
  xlab        = "Faixa Etária",
  ylab        = "Taxa de Mortalidade (por 1000 habitantes)",
  ylim        = c(0, max(dados) * 1.3),
  legend.text = rownames(dados),
  args.legend = list(
    title  = "Grupo",
    x      = "topleft",
    fill   = cores,
    bty    = "n",
    inset  = 0.01
  ),
  border      = "white",
  las         = 1
)