Valor Adicionado do Ceará

Introdução

Este relatório apresenta a participação dos principais setores econômicos no Valor Adicionado Bruto (VAB) dos municípios do Ceará no ano de 2021.Os dados foram obtidos diretamente do sistema SIDRA/IBGE utilizando a linguagem R. O trabalho mostra:participação da agropecuária;participação da indústria;participação dos serviços;mapas temáticos dos municípios;gráfico comparativo; tabela interativa com os resultados municipal.

Pacotes Necessários

Os pacotes abaixo são utilizados para coleta, manipulação, visualização e apresentação dos dados.

library(sidrar)
library(tidyverse)
library(geobr)
library(sf)
library(patchwork)
library(DT)

Coleta de Dados

Nesta etapa são coletados os dados do Valor Adicionado Total e dos setores econômicos para todos os municípios do Ceará.

va <- get_sidra(
  x        = 5938,
  variable = 498,
  period   = "2021",
  geo      = "City",
  geo.filter = list("State" = 23)
) |>
  select(Município, `Município (Código)`, Valor) |>
  rename(va = Valor, code_muni = `Município (Código)`, MUN = Município) |>
  mutate(code_muni = as.double(code_muni))


agro <- get_sidra(
  x        = 5938,
  variable = 513,
  period   = "2021",
  geo      = "City",
  geo.filter = list("State" = 23)
) |>
  select(Município, `Município (Código)`, Valor) |>
  rename(va_agro = Valor, code_muni = `Município (Código)`, MUN = Município) |>
  mutate(code_muni = as.double(code_muni))



ind <- get_sidra(
  x        = 5938,
  variable = 517,
  period   = "2021",
  geo      = "City",
  geo.filter = list("State" = 23)
) |>
  select(Município, `Município (Código)`, Valor) |>
  rename(va_ind = Valor, code_muni = `Município (Código)`, MUN = Município) |>
  mutate(code_muni = as.double(code_muni))

serv1 <- get_sidra(
  x        = 5938,
  variable = 6575,
  period   = "2021",
  geo      = "City",
  geo.filter = list("State" = 23)
) |>
  select(Município, `Município (Código)`, Valor) |>
  rename(va_serv1 = Valor, code_muni = `Município (Código)`, MUN = Município) |>
  mutate(code_muni = as.double(code_muni))

serv2 <- get_sidra(
  x        = 5938,
  variable = 525,
  period   = "2021",
  geo      = "City",
  geo.filter = list("State" = 23)
) |>
  select(Município, `Município (Código)`, Valor) |>
  rename(va_serv2 = Valor, code_muni = `Município (Código)`, MUN = Município) |>
  mutate(code_muni = as.double(code_muni))

serv <- inner_join(serv1, serv2, by = c("code_muni", "MUN")) |>
  mutate(va_serv = va_serv1 + va_serv2) |>
  select(code_muni, MUN, va_serv )

rm(serv1, serv2)

participacao <- va |>
  inner_join(agro, by = c("code_muni", "MUN")) |>
  inner_join(ind,  c("code_muni", "MUN")) |>
  inner_join(serv, c("code_muni", "MUN")) |>
  mutate(
    part_agro = va_agro / va,
    part_ind  = va_ind  / va,
    part_serv = va_serv / va
  ) |>
  select(code_muni, MUN, part_agro, part_ind, part_serv)

Mapas Os mapas mostram a participação setorial em cada município do Ceará.

ce_muni <- read_municipality(code_muni = 23, year = 2020)

ce_map <- ce_muni |>
  left_join(participacao, by = "code_muni")

p1 <- ggplot(ce_map, aes(fill = part_agro)) +
  geom_sf() +
  scale_fill_distiller(palette = "Greens", direction = 1, labels = scales::percent) +
  labs(title = "Agropecuária", fill = NULL) +
  theme_void()

p2 <- ggplot(ce_map, aes(fill = part_ind)) +
  geom_sf() +
  scale_fill_distiller(palette = "Oranges", direction = 1, labels = scales::percent) +
  labs(title = "Indústria", fill = NULL) +
  theme_void()

p3 <- ggplot(ce_map, aes(fill = part_serv)) +
  geom_sf() +
  scale_fill_distiller(palette = "Blues", direction = 1, labels = scales::percent) +
  labs(title = "Serviços", fill = NULL) +
  theme_void()

p1 + p2 + p3 +
  plot_annotation(title = "Participação setorial no Valor Adicionado — Ceará (2021)")

Tabela Municipios

participacao |>
  mutate(
    part_agro = scales::percent(part_agro, accuracy = 0.1),
    part_ind  = scales::percent(part_ind,  accuracy = 0.1),
    part_serv = scales::percent(part_serv, accuracy = 0.1)
  ) |>
  select(
    Município = MUN,
    Agropecuária = part_agro,
    Indústria = part_ind,
    Serviços = part_serv
  ) |>
  datatable(
    caption = "Participação setorial no Valor Adicionado dos municípios do Ceará — 2021",
    rownames = FALSE,
    filter = "top",
    options = list(
      pageLength = 10,
      language = list(url = "//cdn.datatables.net/plug-ins/1.13.6/i18n/pt-BR.json")
    )
  )

Gráfico de Participação Média dos Setores

O gráfico abaixo mostra a participação média dos setores econômicos no Ceará.

grafico_setores <- participacao |>
  summarise(
    Agropecuária = mean(part_agro, na.rm = TRUE),
    Indústria = mean(part_ind, na.rm = TRUE),
    Serviços = mean(part_serv, na.rm = TRUE)
  ) |>
  pivot_longer(cols = everything(),
               names_to = "Setor",
               values_to = "Participacao")

ggplot(grafico_setores,
       aes(x = Setor,
           y = Participacao,
           fill = Setor)) +
  geom_col() +
  scale_y_continuous(labels = scales::percent) +
  labs(
    title = "Participação média dos setores econômicos do Ceará",
    x = "Setores",
    y = "Participação"
  ) +
  theme_minimal()