# Cabeçalho ----
# 2º Censo Distrital da População em Situação de Rua

# IPEDF - Instituto de Pesquisa e Estatística do Distrito Federal

# Acesso a serviços e animais de estimação

# Contato - Email: politicas.sociais@ipe.df.gov.br


# Carregar pacotes e bases ----

# Pacotes
library(readxl)
library(dplyr)
library(survey)
library(srvyr)
library(janitor)
library(tidyr)
library(data.table)
library(openxlsx)


# carregar a base censitária
dados_censitária <- read_xlsx("S:/AREA_TECNICA/2025/Pesquisas Dipos/1. Pastas por projeto/População em situação de rua/Censo 2/Dados/Base censitária/Censo 2 - base censitária.xlsx")

# carregar a base amostral
amostra <- read_xlsx("S:/AREA_TECNICA/2025/Pesquisas Dipos/1. Pastas por projeto/População em situação de rua/Censo 2/Dados/Base amostral/Censo 2 - base amostral - com pesos.xlsx")

# Criar variáveis categóricas
amostra <- amostra %>%
  mutate(
    fx.etaria = cut(`4.1`, c(-Inf, 17, 64, +Inf), labels = c("00-18", "18-59", "60+")))

# Criar variável de acesso aos restaurantes comunitários ao menos 1 vez
amostra <- amostra %>%
  mutate(
    acessou_1_vez_ou_mais = case_when(
      `16.1.14` == "Nenhuma vez" ~ 0L,
      !is.na(`16.1.14`) ~ 1L,
      TRUE ~ NA_integer_
    )
  )

# Criar objeto de amostra complexa com pesos
amostra_pes <- survey::svydesign(id = ~ID,
                                 weights = ~fator,
                                 nest=TRUE,
                                 data=amostra)

amostra_pes <- srvyr::as_survey(amostra_pes)


# Criando variável de faixa etária na base censitária
dados_censitária <- dados_censitária %>%
  mutate(
    faixa_etaria = case_when(
      !is.na(`4.4`) & `4.4` <= 11 ~ "Até 11",
      !is.na(`4.4`) & `4.4` >= 12 & `4.4` <= 17 ~ "12 a 17",
      !is.na(`4.4`) & `4.4` >= 18 & `4.4` <= 30 ~ "18 a 30",
      !is.na(`4.4`) & `4.4` >= 31 & `4.4` <= 49 ~ "31 a 49",
      !is.na(`4.4`) & `4.4` >= 50 & `4.4` <= 59 ~ "50 a 59",
      !is.na(`4.4`) & `4.4` >= 60 & `4.4` <= 69 ~ "60 a 69",
      !is.na(`4.4`) & `4.4` >= 70 & `4.4` <= 79 ~ "70 a 79",
      !is.na(`4.4`) & `4.4` >= 80 & `4.4` <= 89 ~ "80 a 89",
      !is.na(`4.4`) & `4.4` >= 90 & `4.4` <= 99 ~ "90 a 99",
      # Se a idade está faltando, usa diretamente a faixa de 4.4_Observação (caso válida)
      `4.4_Observação` %in% c("Até 11", "12 a 17", "18 a 30", 
                              "31 a 49", "50 a 59", "60 a 69",
                              "70 a 79", "80 a 89", "90 a 99") ~ `4.4_Observação`,
      TRUE ~ "Não foi possível identificar/estimar"
    )
  )

# Animais de estimação ----
# Você tem algum animal?
Tem_animais <- dados_censitária %>%
  filter(tem_animal != "Não se aplica", Responsável != "Criança") %>%
  count(tem_animal) %>%
  filter(tem_animal != "Não se aplica") %>%
  mutate(proporcao = n / sum(n)) %>%  # Calcula a proporção
  arrange(desc(n))  # Ordena do maior para o menor


# Por tipos de animais:
Quais_animais_gatos <- dados_censitária %>%
  filter(`13.1.1` != "Não se aplica", Responsável != "Criança") %>%
  count(`13.1.1`) %>%
  mutate(proporcao = n / sum(n)) %>%  # Calcula a proporção
  arrange(desc(n))  # Ordena do maior para o menor

Quais_animais_cachorros <- dados_censitária %>%
  filter(`13.1.2` != "Não se aplica", Responsável != "Criança") %>%
  count(`13.1.2`) %>%
  mutate(proporcao = n / sum(n)) %>%  # Calcula a proporção
  arrange(desc(n))  # Ordena do maior para o menor

Quais_animais_cavalo <- dados_censitária %>%
  filter(`13.1.3` != "Não se aplica", Responsável != "Criança") %>%
  count(`13.1.3`) %>%
  mutate(proporcao = n / sum(n)) %>%  # Calcula a proporção
  arrange(desc(n))  # Ordena do maior para o menor

Quais_animais_outro <- dados_censitária %>%
  filter(`13.1.4` != "Não se aplica", Responsável != "Criança") %>%
  count(`13.1.4`) %>%
  mutate(proporcao = n / sum(n)) %>%  # Calcula a proporção
  arrange(desc(n))  # Ordena do maior para o menor

Quais_animais_não <- dados_censitária %>%
  filter(`13.1.5` != "Não se aplica", Responsável != "Criança") %>%
  count(`13.1.5`) %>%
  mutate(proporcao = n / sum(n)) %>%  # Calcula a proporção
  arrange(desc(n))  # Ordena do maior para o menor

Quais_animais_nr <- dados_censitária %>%
  filter(`13.1.6` != "Não se aplica", Responsável != "Criança") %>%
  count(`13.1.6`) %>%
  mutate(proporcao = n / sum(n)) %>%  # Calcula a proporção
  arrange(desc(n))  # Ordena do maior para o menor


# Quantos gatos?
gatos <- dados_censitária %>%
  filter(`13.1.1_Gato` != "Não se aplica", Responsável != "Criança") %>% 
  tabyl(`13.1.1_Gato`) %>%  # Conta a frequência de cada categoria
  mutate(proporcao = n / sum(n)) %>%  # Calcula a proporção
  arrange(desc(n))  # Ordena do maior para o menor

# Quantos gatos são vacinados?
  gatos_vacinados <- dados_censitária %>%
    filter(`13.1.1.1` != "Não se aplica", Responsável != "Criança") %>% 
    tabyl(`13.1.1.1`) %>%  # Conta a frequência de cada categoria
    mutate(proporcao = n / sum(n)) %>%  # Calcula a proporção
    arrange(desc(n))  # Ordena do maior para o menor
  
# Quantos cachorros?
  cachorros <- dados_censitária %>%
    filter(`13.1.2_Cachorro` != "Não se aplica", Responsável != "Criança") %>% 
    tabyl(`13.1.2_Cachorro`) %>%  # Conta a frequência de cada categoria
    mutate(proporcao = n / sum(n)) %>%  # Calcula a proporção
    arrange(desc(n))  # Ordena do maior para o menor

# Quantos cachorros são vacinados?
  cachorros_vacinados <- dados_censitária %>%
    filter(`13.1.2.1` != "Não se aplica", Responsável != "Criança") %>% 
    tabyl(`13.1.2.1`) %>%  # Conta a frequência de cada categoria
    mutate(proporcao = n / sum(n)) %>%  # Calcula a proporção
    arrange(desc(n))  # Ordena do maior para o menor
  
# Por estar com seu/s animal/is, você alguma vez já teve dificuldade em acessar seviços 
# ou lugares? 
Dificuldade_centro_pop <- dados_censitária %>%
  filter(tem_animal == 1, `13.2.1` != "Não se aplica", Responsável != "Criança") %>% 
  tabyl(`13.2.1`) %>%  # Conta a frequência de cada categoria
  arrange(desc(n))  # Ordena do maior para o menor
  
Dificuldade_serv_acolh <- dados_censitária %>%
  filter(tem_animal == 1, `13.2.2` != "Não se aplica", Responsável != "Criança") %>% 
  tabyl(`13.2.2`) %>%  # Conta a frequência de cada categoria
  arrange(desc(n))  # Ordena do maior para o menor 
  
Dificuldade_creas <- dados_censitária %>%
  filter(tem_animal == 1, `13.2.3` != "Não se aplica", Responsável != "Criança") %>% 
  tabyl(`13.2.3`) %>%  # Conta a frequência de cada categoria
  arrange(desc(n))  # Ordena do maior para o menor 

Dificuldade_estabelecimentos <- dados_censitária %>%
  filter(tem_animal == 1, `13.2.4` != "Não se aplica", Responsável != "Criança") %>% 
  tabyl(`13.2.4`) %>%  # Conta a frequência de cada categoria
  arrange(desc(n))  # Ordena do maior para o menor 

Dificuldade_UBS <- dados_censitária %>%
  filter(tem_animal == 1, `13.2.5` != "Não se aplica", Responsável != "Criança") %>% 
  tabyl(`13.2.5`) %>%  # Conta a frequência de cada categoria
  arrange(desc(n))  # Ordena do maior para o menor 

Dificuldade_rc <- dados_censitária %>%
  filter(tem_animal == 1, `13.2.6` != "Não se aplica", Responsável != "Criança") %>% 
  tabyl(`13.2.6`) %>%  # Conta a frequência de cada categoria
  arrange(desc(n))  # Ordena do maior para o menor 

Dificuldade_outros <- dados_censitária %>%
  filter(tem_animal == 1, `13.2.7` != "Não se aplica", Responsável != "Criança") %>% 
  tabyl(`13.2.7`) %>%  # Conta a frequência de cada categoria
  arrange(desc(n))  # Ordena do maior para o menor

Sem_dificuldades <- dados_censitária %>%
  filter(tem_animal == 1, `13.2.8` != "Não se aplica", Responsável != "Criança") %>% 
  tabyl(`13.2.8`) %>%  # Conta a frequência de cada categoria
  arrange(desc(n))  # Ordena do maior para o menor

nr_dificuldades <- dados_censitária %>%
  filter(tem_animal == 1, `13.2.9` != "Não se aplica", Responsável != "Criança") %>% 
  tabyl(`13.2.9`) %>%  # Conta a frequência de cada categoria
  arrange(desc(n))  # Ordena do maior para o menor

 

# Serviços públicos e privados acessados ----

# Você recebe algum benefício do Governo?
recebe_benefício <- dados_censitária %>%
  filter(`8.1` != "Não se aplica", Responsável != "Criança") %>% 
  tabyl(`8.1`) %>%  # Conta a frequência de cada categoria
  arrange(desc(n))  # Ordena do maior para o menor

# recebe benefício por sexo
recebe_benefício_sexo <- dados_censitária %>%
  filter(`8.1` != "Não se aplica", Responsável != "Criança") %>%
  count(`4.1.1`, `8.1`) %>%
  group_by(`4.1.1`) %>%
  mutate(
    proporcao = n / sum(n)
  ) %>%
  arrange(`4.1.1`, desc(n))


# Qual(is) dos benefícios abaixo?
BPC <- dados_censitária %>%
  filter(`8.2.1` != "Não se aplica", Responsável != "Criança") %>% 
  tabyl(`8.2.1`) %>%  # Conta a frequência de cada categoria
  arrange(desc(n))  # Ordena do maior para o menor

Seguro_desemprego <- dados_censitária %>%
  filter(`8.2.2` != "Não se aplica", Responsável != "Criança") %>% 
  tabyl(`8.2.2`) %>%  # Conta a frequência de cada categoria
  arrange(desc(n))  # Ordena do maior para o menor

bolsa_família <- dados_censitária %>%
  filter(`8.2.3` != "Não se aplica", Responsável != "Criança") %>% 
  tabyl(`8.2.3`) %>%  # Conta a frequência de cada categoria
  arrange(desc(n))  # Ordena do maior para o menor

DF_social <- dados_censitária %>%
  filter(`8.2.4` != "Não se aplica", Responsável != "Criança") %>% 
  tabyl(`8.2.4`) %>%  # Conta a frequência de cada categoria
  arrange(desc(n))  # Ordena do maior para o menor

Prato_cheio <- dados_censitária %>%
  filter(`8.2.5` != "Não se aplica", Responsável != "Criança") %>% 
  tabyl(`8.2.5`) %>%  # Conta a frequência de cada categoria
  arrange(desc(n))  # Ordena do maior para o menor

Temporário <- dados_censitária %>%
  filter(`8.2.6` != "Não se aplica", Responsável != "Criança") %>% 
  tabyl(`8.2.6`) %>%  # Conta a frequência de cada categoria
  arrange(desc(n))  # Ordena do maior para o menor 

Excepcional <- dados_censitária %>%
  filter(`8.2.7` != "Não se aplica", Responsável != "Criança") %>% 
  tabyl(`8.2.7`) %>%  # Conta a frequência de cada categoria
  arrange(desc(n))  # Ordena do maior para o menor 

Aposentadoria <- dados_censitária %>%
  filter(`8.2.8` != "Não se aplica", Responsável != "Criança") %>% 
  tabyl(`8.2.8`) %>%  # Conta a frequência de cada categoria
  arrange(desc(n))  # Ordena do maior para o menor 

Pensão_por_morte <- dados_censitária %>%
  filter(`8.2.9` != "Não se aplica", Responsável != "Criança") %>% 
  tabyl(`8.2.9`) %>%  # Conta a frequência de cada categoria
  arrange(desc(n))  # Ordena do maior para o menor 

Outros <- dados_censitária %>%
  filter(`8.2.10` != "Não se aplica", Responsável != "Criança") %>% 
  tabyl(`8.2.10`) %>%  # Conta a frequência de cada categoria
  arrange(desc(n))  # Ordena do maior para o menor 

não_sabe <- dados_censitária %>%
  filter(`8.2.11` != "Não se aplica", Responsável != "Criança") %>% 
  tabyl(`8.2.11`) %>%  # Conta a frequência de cada categoria
  arrange(desc(n))  # Ordena do maior para o menor 

não_respondeu <- dados_censitária %>%
  filter(`8.2.12` != "Não se aplica", Responsável != "Criança") %>% 
  tabyl(`8.2.12`) %>%  # Conta a frequência de cada categoria
  arrange(desc(n))  # Ordena do maior para o menor 

outro_qual <- dados_censitária %>%
  filter(`8.2_Outro` != "Não se aplica", Responsável != "Criança") %>% 
  tabyl(`8.2_Outro`) %>%  # Conta a frequência de cada categoria
  arrange(desc(n))  # Ordena do maior para o menor



# Agora eu vou listar alguns serviços e locais e você poderia me dizer quantas vezes acessou 
# nos últimos 6 meses?
Centro_Pop <- amostra_pes %>%
  group_by(`16.1.1`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))


CRAS <- amostra_pes %>%
  group_by(`16.1.2`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))


CREAS <- amostra_pes %>%
  group_by(`16.1.3`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))


Seas_ipês <- amostra_pes %>%
  group_by(`16.1.4`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))


Serv_acolh <- amostra_pes %>%
  group_by(`16.1.5`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))


CAPS <- amostra_pes %>%
  group_by(`16.1.6`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))


Delegacia <- amostra_pes %>%
  group_by(`16.1.7`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))


Defensoria <- amostra_pes %>%
  group_by(`16.1.8`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))


ag_trabalhador <- amostra_pes %>%
  group_by(`16.1.9`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))


Escolas <- amostra_pes %>%
  group_by(`16.1.10`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))


UBS <- amostra_pes %>%
  group_by(`16.1.11`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))


UPA <- amostra_pes %>%
  group_by(`16.1.12`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))


Consultório_rua <- amostra_pes %>%
  group_by(`16.1.13`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))


RC <- amostra_pes %>%
  group_by(`16.1.14`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

# Acesso a Restaurante Comunitário ao menos 1 vez
RC2 <- amostra_pes %>%
  summarise(
    prop_acesso = survey_mean(acessou_1_vez_ou_mais, vartype = "cv", na.rm = TRUE)
  )


CT <- amostra_pes %>%
  group_by(`16.1.15`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))


MP <- amostra_pes %>%
  group_by(`16.1.16`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))



# Nos últimos 6 meses, você foi atendido por alguma instituição social ou religiosa?
Atendim_instituição_social <- amostra_pes %>%
  group_by(`16.2`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))



# Você já foi impedido de entrar em algum local ou serviço?
Impedido <- amostra_pes %>%
  group_by(`16.5`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))


Impedido_sexo <- amostra_pes %>%
  group_by(`4.2`, `16.5`) %>%
  summarise(
    n_ponderado = survey_total(vartype = "cv", na.rm = TRUE),
    .groups = "drop"
  ) %>%
  group_by(`4.2`) %>%
  mutate(
    n_ponderado    = round(as.numeric(n_ponderado)),
    proporcao      = n_ponderado / sum(n_ponderado),
    n_ponderado_cv = n_ponderado_cv,
    cv             = n_ponderado_cv * 100
  ) %>%
  arrange(`4.2`, desc(n_ponderado)) %>%
  ungroup()


# Em quais desses locais você já foi impedido de entrar?
Shopping <- amostra_pes %>%
  filter(`16.5` == "Sim") %>% 
  group_by(`16.5.1.1`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

restaurantes <- amostra_pes %>%
  filter(`16.5` == "Sim") %>% 
  group_by(`16.5.1.2`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

banco <- amostra_pes %>%
  filter(`16.5` == "Sim") %>% 
  group_by(`16.5.1.3`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

serv_saúde <- amostra_pes %>%
  filter(`16.5` == "Sim") %>% 
  group_by(`16.5.1.4`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

serv_assistência <- amostra_pes %>%
  filter(`16.5` == "Sim") %>% 
  group_by(`16.5.1.5`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))


outros_órgãos_públicos <- amostra_pes %>%
  filter(`16.5` == "Sim") %>% 
  group_by(`16.5.1.6`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))


outro_ambiente_público <- amostra_pes %>%
  filter(`16.5` == "Sim") %>% 
  group_by(`16.5.1.7`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

NR <- amostra_pes %>%
  filter(`16.5` == "Sim") %>% 
  group_by(`16.5.1.8`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))


outro_qual <- amostra_pes %>%
  filter(`16.5` == "Sim", `16.5.1_Outro` != "Não se aplica") %>% 
  group_by(`16.5.1_Outro`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))



# Avaliação serviços de acolhimento ----

# Você já utilizou os Serviços de Acolhimento?
usou_serv_acolh <- amostra_pes %>%
  filter(`17.1` != "Não se aplica") %>% 
  group_by(`17.1`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

# Uso do serv de acolhimento por sexo:
usou_serv_acolh_sexo <- amostra_pes %>%
  filter(`17.1` != "Não se aplica") %>% 
  group_by(`17.1`, `4.2`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE),
    .groups = "drop"
  ) %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado))  # total arredondado
  ) %>%
  group_by(`4.2`) %>%
  mutate(
    proporcao = n_ponderado / sum(n_ponderado)    # proporção por sexo
  ) %>%
  ungroup() %>%
  arrange(desc(n_ponderado))



# Por que você não utilizou os Serviços de Acolhimento?
falta_vaga <- amostra_pes %>%
  filter(`17.1` == "Não") %>% 
  group_by(`17.2.1`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

guardar_material_trabalho <- amostra_pes %>%
  filter(`17.1` == "Não") %>% 
  group_by(`17.2.2`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

não_tem_interesse <- amostra_pes %>%
  filter(`17.1` == "Não") %>% 
  group_by(`17.2.3`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

regras <- amostra_pes %>%
  filter(`17.1` == "Não") %>% 
  group_by(`17.2.4`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

dific_relacionamento <- amostra_pes %>%
  filter(`17.1` == "Não") %>% 
  group_by(`17.2.5`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

ñ_animais <- amostra_pes %>%
  filter(`17.1` == "Não") %>% 
  group_by(`17.2.6`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

infraestrutura <- amostra_pes %>%
  filter(`17.1` == "Não") %>% 
  group_by(`17.2.7`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

comida <- amostra_pes %>%
  filter(`17.1` == "Não") %>% 
  group_by(`17.2.8`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

Drogas <- amostra_pes %>%
  filter(`17.1` == "Não") %>% 
  group_by(`17.2.9`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

Furtos <- amostra_pes %>%
  filter(`17.1` == "Não") %>% 
  group_by(`17.2.10`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

Atividades <- amostra_pes %>%
  filter(`17.1` == "Não") %>% 
  group_by(`17.2.11`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

Atendimento <- amostra_pes %>%
  filter(`17.1` == "Não") %>% 
  group_by(`17.2.12`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

Localização_ruim <- amostra_pes %>%
  filter(`17.1` == "Não") %>% 
  group_by(`17.2.13`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

Outro <- amostra_pes %>%
  filter(`17.1` == "Não") %>% 
  group_by(`17.2.14`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

# Entre ótimo e péssimo, como você avalia os Serviços de Acolhimento?
avalia_serv_acolh <- amostra_pes %>%
  filter(`17.3` != "Não se aplica") %>% 
  group_by(`17.3`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))


# O que você acha que o Serviço de Acolhimento pode lhe proporcionar?
Liberdade <- amostra_pes %>%
  filter(`17.4.1` != "Não se aplica") %>% 
  group_by(`17.4.1`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

Autonomia <- amostra_pes %>%
  filter(`17.4.2` != "Não se aplica") %>% 
  group_by(`17.4.2`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

Privacidade <- amostra_pes %>%
  filter(`17.4.3` != "Não se aplica") %>% 
  group_by(`17.4.3`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

Segurança <- amostra_pes %>%
  filter(`17.4.4` != "Não se aplica") %>% 
  group_by(`17.4.4`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

Familiares <- amostra_pes %>%
  filter(`17.4.5` != "Não se aplica") %>% 
  group_by(`17.4.5`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

Deixar_rua <- amostra_pes %>%
  filter(`17.4.6` != "Não se aplica") %>% 
  group_by(`17.4.6`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

Outro <- amostra_pes %>%
  filter(`17.4.7` != "Não se aplica") %>% 
  group_by(`17.4.7`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

Nada <- amostra_pes %>%
  filter(`17.4.8` != "Não se aplica") %>% 
  group_by(`17.4.8`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))


# Na sua opinião, existem problemas do Serviço de Acolhimento?
Problemas_serv <- amostra_pes %>%
  filter(`17.5` != "Não se aplica") %>% 
  group_by(`17.5`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))


# Quais probelmas existem?
liberdade <- amostra_pes %>%
  filter(`17.5.1.1` != "Não se aplica") %>% 
  group_by(`17.5.1.1`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

dific_relacionamento_abrigados <- amostra_pes %>%
  filter(`17.5.1.2` != "Não se aplica") %>% 
  group_by(`17.5.1.2`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

dific_relacionamento_funcionários <- amostra_pes %>%
  filter(`17.5.1.3` != "Não se aplica") %>% 
  group_by(`17.5.1.3`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))


dific_relacionamento_funcionários <- amostra_pes %>%
  filter(`17.5.1.3` != "Não se aplica") %>% 
  group_by(`17.5.1.3`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

animais <- amostra_pes %>%
  filter(`17.5.1.4` != "Não se aplica") %>% 
  group_by(`17.5.1.4`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

infraestrutura <- amostra_pes %>%
  filter(`17.5.1.5` != "Não se aplica") %>% 
  group_by(`17.5.1.5`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

comida <- amostra_pes %>%
  filter(`17.5.1.6` != "Não se aplica") %>% 
  group_by(`17.5.1.6`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

uso_drogas <- amostra_pes %>%
  filter(`17.5.1.7` != "Não se aplica") %>% 
  group_by(`17.5.1.7`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

limpeza <- amostra_pes %>%
  filter(`17.5.1.8` != "Não se aplica") %>% 
  group_by(`17.5.1.8`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

limpeza <- amostra_pes %>%
  filter(`17.5.1.8` != "Não se aplica") %>% 
  group_by(`17.5.1.8`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

roubos <- amostra_pes %>%
  filter(`17.5.1.9` != "Não se aplica") %>% 
  group_by(`17.5.1.9`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

falta_fazer <- amostra_pes %>%
  filter(`17.5.1.10` != "Não se aplica") %>% 
  group_by(`17.5.1.10`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

atendimento_ruim <- amostra_pes %>%
  filter(`17.5.1.11` != "Não se aplica") %>% 
  group_by(`17.5.1.11`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

localização <- amostra_pes %>%
  filter(`17.5.1.12` != "Não se aplica") %>% 
  group_by(`17.5.1.12`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

material_trab <- amostra_pes %>%
  filter(`17.5.1.13` != "Não se aplica") %>% 
  group_by(`17.5.1.13`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

Outro <- amostra_pes %>%
  filter(`17.5.1.14` != "Não se aplica") %>% 
  group_by(`17.5.1.14`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))


# Acesso à internet ----
# Você tem celular?
Tem_cel <- amostra_pes %>%
  group_by(`13.1`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

# Você aceessa a internet nesse celular ou em outros aparelhos?
Acesso_net <- amostra_pes %>%
  filter(`13.1` == "Sim") %>% 
  group_by(`13.2.1`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

Acesso_net2 <- amostra_pes %>%
  filter(`13.1` == "Sim") %>% 
  group_by(`13.2.2`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

Acesso_net3 <- amostra_pes %>%
  filter(`13.1` == "Sim") %>% 
  group_by(`13.2.3`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

# Em qual local você usa internet?
dados <- amostra_pes %>%
  filter(`13.3.1` != "Não se aplica") %>% 
  group_by(`13.3.1`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

casa <- amostra_pes %>%
  filter(`13.3.2` != "Não se aplica") %>% 
  group_by(`13.3.2`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

serv_públicos <- amostra_pes %>%
  filter(`13.3.3` != "Não se aplica") %>% 
  group_by(`13.3.3`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

comércio <- amostra_pes %>%
  filter(`13.3.4` != "Não se aplica") %>% 
  group_by(`13.3.4`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

inst_sociais <- amostra_pes %>%
  filter(`13.3.5` != "Não se aplica") %>% 
  group_by(`13.3.5`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))


# Acesso a documentos ----
# Quais desses documentos você tem ou teve?

rg <- amostra_pes %>%
  group_by(`16.7.1`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

CTrab <- amostra_pes %>%
  group_by(`16.7.2`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))


CPF <- amostra_pes %>%
  group_by(`16.7.3`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))


título_eleitor <- amostra_pes %>%
  group_by(`16.7.4`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

SUS <- amostra_pes %>%
  group_by(`16.7.5`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

certidão_nascimento <- amostra_pes %>%
  group_by(`16.7.6`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

RANI <- amostra_pes %>%
  group_by(`16.7.7`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

# Qual dos documentos abaixo a criança/adolescente possui?
SUS <- dados_censitária %>%
  filter(`11.15.1` != "Não se aplica", Responsável != "Entrevistado") %>% 
  tabyl(`11.15.1`) %>%  # Conta a frequência de cada categoria
  arrange(desc(n))  # Ordena do maior para o menor

RG <- dados_censitária %>%
  filter(`11.15.2` != "Não se aplica", Responsável != "Entrevistado") %>% 
  tabyl(`11.15.2`) %>%  # Conta a frequência de cada categoria
  arrange(desc(n))  # Ordena do maior para o menor

CPF <- dados_censitária %>%
  filter(`11.15.3` != "Não se aplica", Responsável != "Entrevistado") %>% 
  tabyl(`11.15.3`) %>%  # Conta a frequência de cada categoria
  arrange(desc(n))  # Ordena do maior para o menor

Certidão_nasc <- dados_censitária %>%
  filter(`11.15.4` != "Não se aplica", Responsável != "Entrevistado") %>% 
  tabyl(`11.15.4`) %>%  # Conta a frequência de cada categoria
  arrange(desc(n))  # Ordena do maior para o menor

Carteira_vacina <- dados_censitária %>%
  filter(`11.15.5` != "Não se aplica", Responsável != "Entrevistado") %>% 
  tabyl(`11.15.5`) %>%  # Conta a frequência de cada categoria
  arrange(desc(n))  # Ordena do maior para o menor

RANI <- dados_censitária %>%
  filter(`11.15.6` != "Não se aplica", Responsável != "Entrevistado") %>% 
  tabyl(`11.15.6`) %>%  # Conta a frequência de cada categoria
  arrange(desc(n))  # Ordena do maior para o menor


# Desagregar por faixa etária das crianças e adolescentes
SUS_fx_etária <- dados_censitária %>%
  filter(`11.15.1` != "Não se aplica", Responsável != "Entrevistado") %>% 
  count(faixa_etaria, `11.15.1`) %>%  # Conta a frequência de cada categoria
  group_by(faixa_etaria) %>% 
  mutate(proporcao = n / sum(n)) %>% 
  arrange(faixa_etaria)  # Ordena do maior para o menor

RG_fx_etária <- dados_censitária %>%
  filter(`11.15.2` != "Não se aplica", Responsável != "Entrevistado") %>% 
  count(faixa_etaria, `11.15.2`) %>%  # Conta a frequência de cada categoria
  group_by(faixa_etaria) %>% 
  mutate(proporcao = n / sum(n)) %>% 
  arrange(faixa_etaria)  # Ordena do maior para o menor

CPF_fx_etária <- dados_censitária %>%
  filter(`11.15.3` != "Não se aplica", Responsável != "Entrevistado") %>% 
  count(faixa_etaria, `11.15.3`) %>%  # Conta a frequência de cada categoria
  group_by(faixa_etaria) %>% 
  mutate(proporcao = n / sum(n)) %>% 
  arrange(faixa_etaria)  # Ordena do maior para o menor

Certidão_fx_etária <- dados_censitária %>%
  filter(`11.15.4` != "Não se aplica", Responsável != "Entrevistado") %>% 
  count(faixa_etaria, `11.15.4`) %>%  # Conta a frequência de cada categoria
  group_by(faixa_etaria) %>% 
  mutate(proporcao = n / sum(n)) %>% 
  arrange(faixa_etaria)  # Ordena do maior para o menor

Vacina_fx_etária <- dados_censitária %>%
  filter(`11.15.5` != "Não se aplica", Responsável != "Entrevistado") %>% 
  count(faixa_etaria, `11.15.5`) %>%  # Conta a frequência de cada categoria
  group_by(faixa_etaria) %>% 
  mutate(proporcao = n / sum(n)) %>% 
  arrange(faixa_etaria)  # Ordena do maior para o menor

Rani_fx_etária <- dados_censitária %>%
  filter(`11.15.6` != "Não se aplica", Responsável != "Entrevistado") %>% 
  count(faixa_etaria, `11.15.6`) %>%  # Conta a frequência de cada categoria
  group_by(faixa_etaria) %>% 
  mutate(proporcao = n / sum(n)) %>% 
  arrange(faixa_etaria)  # Ordena do maior para o menor


# Sentimentos e percepções sobre a situação de rua ----
# Desde que está em situação de rua, por quais dessas situações você já passou? 
Discriminação <- amostra_pes %>%
  group_by(`16.6.1`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

Vigilância <- amostra_pes %>%
  group_by(`16.6.2`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

Violência <- amostra_pes %>%
  group_by(`16.6.3`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

Violência_sexual <- amostra_pes %>%
  group_by(`16.6.4`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

tentativa_feminicídio <- amostra_pes %>%
  group_by(`16.6.5`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

roubado <- amostra_pes %>%
  group_by(`16.6.6`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

internado <- amostra_pes %>%
  group_by(`16.6.7`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

preso <- amostra_pes %>%
  group_by(`16.6.8`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

negligenciado <- amostra_pes %>%
  group_by(`16.6.9`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

negligenciado <- amostra_pes %>%
  group_by(`16.6.9`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

outra_situação <- amostra_pes %>%
  group_by(`16.6.10`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))


# Você tem interesse em sair da situação de rua?
sair_rua <- amostra_pes %>%
  group_by(`16.3`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))


# Do que você precisa para sair dessa situação?
emprego <- amostra_pes %>%
  group_by(`16.4.1`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

moradia <- amostra_pes %>%
  group_by(`16.4.2`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

ac_médico <- amostra_pes %>%
  group_by(`16.4.3`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

vínculos_família <- amostra_pes %>%
  group_by(`16.4.4`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

não_quer <- amostra_pes %>%
  group_by(`16.4.5`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

outras_necessidades <- amostra_pes %>%
  group_by(`16.4.6`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))

outro_qual <- amostra_pes %>%
  filter(`16.4_Outro` != "Não se aplica") %>% 
  group_by(`16.4_Outro`) %>%  
  summarise(
    n_ponderado  = survey_total(vartype = "cv", na.rm = TRUE)
  ) %>%
  ungroup() %>%
  mutate(
    n_ponderado = round(as.numeric(n_ponderado)),      # total arredondado
    proporcao   = n_ponderado / sum(n_ponderado),      # proporção simples
  ) %>%
  arrange(desc(n_ponderado))