########################################################################################################################
#Pesquisa: Levantamento do perfil sociodemográfico, necessidades e barreiras de acesso a serviços públicos por pessoas
#com síndrome de Down no Distrito Federal
#Instituto de Pesquisa e Estatística do Distrito Federal - IPEDF Codeplan
#Diretoria de Estudos e Políticas Sociais - DIPOS
#Contato - Email: politicas.sociais@ipe.df.gov.br
########################################################################################################################

##### Pacotes ####
rm(list = ls())
library(readxl)
library(tidyverse)
library(data.table)
library(dplyr)
library(janitor)
library(srvyr)

##### Abrir a base #####
base <- read_excel("Base_de_dados-Sindrome-de-Down.xlsx")
base <- data.table(base)

##### Análise dos dados ####

## Número de respondentes por perfil
perfil_resp <- base %>% 
  select(`B1-Você é:`) %>% 
  # remove linhas com todas as colunas vazias
  remove_empty("rows") %>% 
  # mutate(H2=ifelse(is.na(H2), "", H2)) %>% 
  # mutate(M1=ifelse(is.na(M1), "", M1)) %>% 
  unite(perfil, c(`B1-Você é:`), sep = " ", remove = FALSE) %>% 
  select(-`B1-Você é:`) %>% 
  mutate(perfil = stringr::str_trim(perfil, side = "both")) %>% 
  count(perfil) %>% 
  mutate(Perc = prop.table(n)*100) %>% 
  mutate(Perc = round(Perc, 2))

## Sexo
sexo_ep <- base %>% 
  select(`H2-Qual o sexo que..`) %>% 
  # remove linhas com todas as colunas vazias
  remove_empty("rows") %>%
  mutate(`H2-Qual o sexo que..`=ifelse(is.na(`H2-Qual o sexo que..`), "", `H2-Qual o sexo que..`)) %>% 
  # mutate(M1=ifelse(is.na(M1), "", M1)) %>% 
  unite(sexo, c(`H2-Qual o sexo que..`), sep = " ", remove = FALSE) %>% 
  # select(-'H2-Qual o sexo que..') %>% 
  mutate(sexo = stringr::str_trim(sexo, side = "both")) %>% 
  count(sexo) %>% 
  mutate(Perc = prop.table(n)*100) %>% 
  mutate(Perc = round(Perc, 2))

## Renda mensal do domicílio da pessoa com síndrome de Down
table(base$H11.Somando.a.renda.., useNA = "always")
table(base$`M10-Somando a renda..`, useNA = "always")

rename(base, c())

base[`H11-Somando a renda..`==  "ATÉ 1 SALÁRIO MÍNIMO (ATÉ R$ 1.212,00 INCLUSIVE)", `H11-Somando a renda..` := "Até 2 salários-mínimos"]
base[`H11-Somando a renda..`==  "DE 1 A 2 SALÁRIOS MÍNIMOS (DE R$ 1.212,00 ATÉ R$ 2.424,00 INCLUSIVE)", `H11-Somando a renda..` := "Até 2 salários-mínimos"]
base[`H11-Somando a renda..`==  "DE 2 A 5 SALÁRIOS MÍNIMOS (DE R$ 2.424,00 ATÉ R$ 6.060,00 INCLUSIVE)", `H11-Somando a renda..` := "De 2 a 5 salários-mínimos"]
base[`H11-Somando a renda..`==  "DE 5 A 10 SALÁRIOS MÍNIMOS (DE R$ 6.060,00 ATÉ R$ 12.120,00 INCLUSIVE)", `H11-Somando a renda..` := "Mais de 5 salários-mínimos"]
base[`H11-Somando a renda..`==  "DE 10 A 20 SALÁRIOS MÍNIMOS (DE R$ 12.120,00 ATÉ R$ 24.240,00 INCLUSIVE)", `H11-Somando a renda..` := "Mais de 5 salários-mínimos"]
base[`H11-Somando a renda..`==  "DE 20 A 30 SALÁRIOS MÍNIMOS (DE R$ 24.240,00 ATÉ R$ 36.360,00 INCLUSIVE)", `H11-Somando a renda..` := "Mais de 5 salários-mínimos"]
base[`H11-Somando a renda..`==  "MAIS DE 30 SALÁRIOS MÍNIMOS (MAIS DE R$ 36.360,00)", `H11-Somando a renda..` := "Mais de 5 salários-mínimos"]
base[`H11-Somando a renda..`==  "PREFIRO NÃO RESPONDER", `H11-Somando a renda..` := "Não sabe/não respondeu"]
base[`H11-Somando a renda..`==  "SEI INFORMAR", `H11-Somando a renda..` := "Não sabe/não respondeu"]

base[`M10-Somando a renda..`==  "ATÉ 1 SALÁRIO MÍNIMO (ATÉ R$ 1.212,00 INCLUSIVE)", `M10-Somando a renda..` := "Até 2 salários-mínimos"]
base[`M10-Somando a renda..`==  "DE 1 A 2 SALÁRIOS MÍNIMOS (DE R$ 1.212,00 ATÉ R$ 2.424,00 INCLUSIVE)", `M10-Somando a renda..` := "Até 2 salários-mínimos"]
base[`M10-Somando a renda..`==  "DE 2 A 5 SALÁRIOS MÍNIMOS (DE R$ 2.424,00 ATÉ R$ 6.060,00 INCLUSIVE)", `M10-Somando a renda..` := "De 2 a 5 salários-mínimos"]
base[`M10-Somando a renda..`==  "DE 5 A 10 SALÁRIOS MÍNIMOS (DE R$ 6.060,00 ATÉ R$ 12.120,00 INCLUSIVE)", `M10-Somando a renda..` := "Mais de 5 salários-mínimos"]
base[`M10-Somando a renda..`==  "DE 10 A 20 SALÁRIOS MÍNIMOS (DE R$ 12.120,00 ATÉ R$ 24.240,00 INCLUSIVE)", `M10-Somando a renda..` := "Mais de 5 salários-mínimos"]
base[`M10-Somando a renda..`==  "DE 20 A 30 SALÁRIOS MÍNIMOS (DE R$ 24.240,00 ATÉ R$ 36.360,00 INCLUSIVE)", `M10-Somando a renda..` := "Mais de 5 salários-mínimos"]
base[`M10-Somando a renda..`==  "MAIS DE 30 SALÁRIOS MÍNIMOS (MAIS DE R$ 36.360,00)", `M10-Somando a renda..` := "Mais de 5 salários-mínimos"]
base[`M10-Somando a renda..`==  "PREFIRO NÃO RESPONDER", `M10-Somando a renda..` := "Não sabe/não respondeu"]
base[`M10-Somando a renda..`==  "SEI INFORMAR", `M10-Somando a renda..` := "Não sabe/não respondeu"]

Renda_Total <- base %>% 
  select(`H11-Somando a renda..`,`M10-Somando a renda..`) %>% 
  # remove linhas com todas as colunas vazias
  remove_empty("rows") %>% 
  mutate(`H11-Somando a renda..`=ifelse(is.na(`H11-Somando a renda..`), "", `H11-Somando a renda..`)) %>% 
  mutate(`M10-Somando a renda..`=ifelse(is.na(`M10-Somando a renda..`), "", `M10-Somando a renda..`)) %>% 
  unite(Renda, c(`H11-Somando a renda..`,`M10-Somando a renda..`), sep = " ", remove = FALSE) %>% 
  select(-`H11-Somando a renda..`, -`M10-Somando a renda..`) %>% 
  mutate(Renda = stringr::str_trim(Renda, side = "both")) %>% 
  count(Renda) %>% 
  mutate(Perc = prop.table(n)*100) %>% 
  mutate(Perc = round(Perc, 2)) %>% 
  data.table(.)


## Região administrativa
table(base$`H8-Em qual região ..`, useNA = "always")
table(base$`M7-Em qual região ..`, useNA = "always")

Reg_Adm <- base %>% 
  select(`H8-Em qual região ..`,`M7-Em qual região ..`) %>% 
  # remove linhas com todas as colunas vazias
  remove_empty("rows") %>% 
  mutate(`H8-Em qual região ..`=ifelse(is.na(`H8-Em qual região ..`), "", `H8-Em qual região ..`)) %>% 
  mutate(`M7-Em qual região ..`=ifelse(is.na(`M7-Em qual região ..`), "", `M7-Em qual região ..`)) %>% 
  unite(RA, c(`H8-Em qual região ..`,`M7-Em qual região ..`), sep = " ", remove = FALSE) %>% 
  select(-`H8-Em qual região ..`, -`M7-Em qual região ..`) %>% 
  mutate(RA = stringr::str_trim(RA, side = "both")) %>% 
  count(RA) %>% 
  mutate(Perc = prop.table(n)*100) %>% 
  mutate(Perc = round(Perc, 2))

## Dificuldades para a obtenção do diagnóstico
# base[`D4[SQ001]-Por quais motiv..  [Por falta de co.. ]`==  "ATÉ 1 SALÁRIO MÍNIMO (ATÉ R$ 1.212,00 INCLUSIVE)", `H11-Somando a renda..` := "Até 2 salários-mínimos"]

base$`D4[SQ001]` = ifelse(base$`D4[SQ001]-Por quais motiv..  [Por falta de co.. ]`=="Sim", "Falta de Conhecim. dos Profis. de Saúde Sobre a SD", NA)
base$`D4[SQ002]` = ifelse(base$`D4[SQ002]-Por quais motiv..  [Por falta de ac.. ]`=="Sim", "Por Falta de Acomp. Pré-Natal Adequado", NA)
base$`D4[SQ003]` = ifelse(base$`D4[SQ003]-Por quais motiv..  [Por falta de es.. ]`=="Sim", "Falta de Estrut. Hosp. Adeq. Para A Realiz. do Diag.", NA)
base$`D4[SQ004]` = ifelse(base$`D4[SQ004]-Por quais motiv..  [Não sei informar]`=="Sim", "Não sei", NA)
base$`D4[SQ005]` = ifelse(base$`D4[SQ005]-Por quais motiv..  [Prefiro não res.. ]`=="Sim", "Prefiro não responder", NA)
base$`D4[other]` = ifelse(base$`D4[other]-Por quais motiv..  [Outros]`!="Não" | !is.na(base$`D4[other]-Por quais motiv..  [Outros]`), "Outros", NA)

base$`D4[other]` = ifelse(base$`D4[other]-Por quais motiv..  [Outros]`=="x", NA, base$`D4[other]`)

teste <- base %>% 
  select(`D4[other]-Por quais motiv..  [Outros]`, 'D4[other]')

dific_diagn <- base %>% 
  select(`D4[SQ001]`, `D4[SQ002]`, `D4[SQ003]`, `D4[SQ004]`, `D4[other]`) %>% 
  pivot_longer(
    cols = c(`D4[SQ001]`, `D4[SQ002]`, `D4[SQ003]`, `D4[SQ004]`, `D4[other]`),
    names_to = "sample",
    values_to = "Dificuldade",
    values_drop_na = TRUE
  ) %>% 
  select(-sample) %>% 
  arrange(Dificuldade) %>%
  count(Dificuldade) %>% 
  # 109 respostas válidas neste bloco
  mutate(Porcentagem = (n/109)*100) %>% 
  mutate(Porcentagem = round(Porcentagem, 2))

## Especialidades em saúde acessadas por pessoas com síndrome de Down
especialidades <- base %>% 
  select(`D8[SQ001]-A pessoa com sí..  [Geneticista]`, `D8[SQ002]-A pessoa com sí..  [Cardiologista]`, `D8[SQ003]-A pessoa com sí..  [Neurologista]`, 	
         `D8[SQ004]-A pessoa com sí..  [Clínico geral]`, `D8[SQ007]-A pessoa com sí..  [Oftalmologista]`, `D8[SQ008]-A pessoa com sí..  [Otorrinolaringo.. ]`, 	
         `D8[SQ009]-A pessoa com sí..  [Psicólogo]`, `D8[SQ010]-A pessoa com sí..  [Psiquiatra]`, `D8[SQ011]-A pessoa com sí..  [Fisioterapeuta]`, 	
         `D8[SQ012]-A pessoa com sí..  [Terapeuta ocupa.. ]`, `D8[SQ013]-A pessoa com sí..  [Fonoaudiólogo]`, `D8[SQ014]-A pessoa com sí..  [Nutricionista]`, 	
         `D8[SQ016]-A pessoa com sí..  [Pediatra]`, `D8[SQ019]-A pessoa com sí..  [Endocrinologista]`, `I5[SQ001]-Você faz acompa..  [Geneticista]`, 	
         `I5[SQ002]-Você faz acompa..  [Cardiologista]`, `I5[SQ003]-Você faz acompa..  [Neurologista]`, `I5[SQ004]-Você faz acompa..  [Clínico geral]`, 	
         `I5[SQ007]-Você faz acompa..  [Oftalmologista]`, `I5[SQ008]-Você faz acompa..  [Otorrinolaringo.. ]`, `I5[SQ009]-Você faz acompa..  [Psicólogo]`, 	
         `I5[SQ010]-Você faz acompa..  [Psiquiatra]`, `I5[SQ011]-Você faz acompa..  [Fisioterapeuta]`, `I5[SQ012]-Você faz acompa..  [Terapeuta ocupa.. ]`, 	
         `I5[SQ013]-Você faz acompa..  [Fonoaudiólogo]`, `I5[SQ014]-Você faz acompa..  [Nutricionista]`, `I5[SQ016]-Você faz acompa..  [Pediatra]`, 	
         `I5[SQ019]-Você faz acompa..  [Endocrinologista]`) %>% 
  pivot_longer(
    cols = c(`D8[SQ001]-A pessoa com sí..  [Geneticista]`, `D8[SQ002]-A pessoa com sí..  [Cardiologista]`, `D8[SQ003]-A pessoa com sí..  [Neurologista]`, 	
             `D8[SQ004]-A pessoa com sí..  [Clínico geral]`, `D8[SQ007]-A pessoa com sí..  [Oftalmologista]`, `D8[SQ008]-A pessoa com sí..  [Otorrinolaringo.. ]`, 	
             `D8[SQ009]-A pessoa com sí..  [Psicólogo]`, `D8[SQ010]-A pessoa com sí..  [Psiquiatra]`, `D8[SQ011]-A pessoa com sí..  [Fisioterapeuta]`, 	
             `D8[SQ012]-A pessoa com sí..  [Terapeuta ocupa.. ]`, `D8[SQ013]-A pessoa com sí..  [Fonoaudiólogo]`, `D8[SQ014]-A pessoa com sí..  [Nutricionista]`, 	
             `D8[SQ016]-A pessoa com sí..  [Pediatra]`, `D8[SQ019]-A pessoa com sí..  [Endocrinologista]`, `I5[SQ001]-Você faz acompa..  [Geneticista]`, 	
             `I5[SQ002]-Você faz acompa..  [Cardiologista]`, `I5[SQ003]-Você faz acompa..  [Neurologista]`, `I5[SQ004]-Você faz acompa..  [Clínico geral]`, 	
             `I5[SQ007]-Você faz acompa..  [Oftalmologista]`, `I5[SQ008]-Você faz acompa..  [Otorrinolaringo.. ]`, `I5[SQ009]-Você faz acompa..  [Psicólogo]`, 	
             `I5[SQ010]-Você faz acompa..  [Psiquiatra]`, `I5[SQ011]-Você faz acompa..  [Fisioterapeuta]`, `I5[SQ012]-Você faz acompa..  [Terapeuta ocupa.. ]`, 	
             `I5[SQ013]-Você faz acompa..  [Fonoaudiólogo]`, `I5[SQ014]-Você faz acompa..  [Nutricionista]`, `I5[SQ016]-Você faz acompa..  [Pediatra]`, 	
             `I5[SQ019]-Você faz acompa..  [Endocrinologista]`),
    names_to = "Variavel",
    values_to = "Resposta",
    values_drop_na = TRUE
  )

especialidades$Titulo <- substr(especialidades$Variavel, 31, 70)
especialidades$Titulo <- str_sub(especialidades$Titulo, -30, -2)
table(especialidades$Titulo, useNA = "always")
especialidades$Variavel <- NULL
especialidades$Titulo <- ifelse(especialidades$Titulo=="Otorrinolaringo.. ", "Otorrinolaringologista",
                                ifelse(especialidades$Titulo=="Terapeuta ocupa.. ", "Terapeuta ocupacional", especialidades$Titulo))
table(especialidades$Titulo, useNA = "always")
table(especialidades$Resposta)

especialidades = especialidades

especialidades$Resposta <- substr(especialidades$Resposta, 1, 3)
table(especialidades2$Resposta)

especialidades <- especialidades %>% 
  drop_na(Resposta) %>% 
  # select(-Variavel) %>% 
  arrange(Titulo) %>%
  count(Titulo, Resposta) %>% 
  # mutate(Resposta=ifelse(is.na(Resposta), "", Resposta)) %>% 
  subset(Titulo != "") %>% 
  # 533 respostas neste bloco
  mutate(Porcentagem = (n/542)*100) %>%
  mutate(Porcentagem = round(Porcentagem, 2)) %>% 
  select(-n) %>% 
  group_by(Resposta, Titulo) %>% 
  reshape2::dcast(Titulo ~ Resposta) %>% 
  replace(is.na(.), 0) 
  
## Região administrativa
table(base$`D11-Especifique a r..`, useNA = "always")

Reg_Adm <- base %>% 
  select(`D11-Especifique a r..`) %>% 
  # remove linhas com todas as colunas vazias
  remove_empty("rows") %>% 
  mutate(`D11-Especifique a r..`=ifelse(is.na(`D11-Especifique a r..`), "", `D11-Especifique a r..`)) %>% 
  unite(RA, `D11-Especifique a r..`, sep = " ", remove = FALSE) %>%
  # select(-`D11-Especifique a r..`, -`I8-Em qual região ..`) %>% 
  mutate(RA = stringr::str_trim(RA, side = "both")) %>% 
  count(RA) %>% 
  mutate(Perc = prop.table(n)*100) %>% 
  mutate(Perc = round(Perc, 2))

## Dificuldades enfrentadas por pessoas com síndrome de Down ao acessar os serviços de atendimento na área da saúde

base$`D15[SQ001]` = ifelse(base$`D15[SQ001]-Quais dificulda..  [Dificuldade de .. ]`=="Sim", "Dificuldade de agendamento", NA)
base$`D15[SQ002]` = ifelse(base$`D15[SQ002]-Quais dificulda..  [Falta de profis.. ]`=="Sim", "Falta de profissional especializado na rede pública", NA)
base$`D15[SQ003]` = ifelse(base$`D15[SQ003]-Quais dificulda..  [Falta de profis.. ]`=="Sim", "Falta de profissional especializado na rede particular", NA)
base$`D15[SQ004]` = ifelse(base$`D15[SQ004]-Quais dificulda..  [Profissional se.. ]`=="Sim", "Profissional sem preparo para lidar com pacientes com síndrome de Down", NA)
base$`D15[SQ005]` = ifelse(base$`D15[SQ005]-Quais dificulda..  [Indisponibilida.. ]`=="Sim", "Indisponibilidade de serviço perto de onde mora", NA)
base$`D15[SQ006]` = ifelse(base$`D15[SQ006]-Quais dificulda..  [Alto custo das .. ]`=="Sim", "Alto custo das consultas/acompanhamento nas especialidades que não atendem na rede pública", NA)
base$`D15[SQ007]` = ifelse(base$`D15[SQ007]-Quais dificulda..  [Falta de inform.. ]`=="Sim", "Falta de informação sobre onde buscar atendimento para uma dada especialidade", NA)
base$`D15[SQ008]` = ifelse(base$`D15[SQ008]-Quais dificulda..  [Dificuldade de .. ]`=="Sim", "Dificuldade de deslocamento ao local de atendimento", NA)
base$`D15[SQ009]` = ifelse(base$`D15[SQ009]-Quais dificulda..  [Longa lista de .. ]`=="Sim", "Longa lista de espera para acessar o serviço de saúde", NA)
base$`D15[SQ010]` = ifelse(base$`D15[SQ010]-Quais dificulda..  [Horário de trab.. ]`=="Sim", "Horário de trabalho/estudos incompatível com o horário do serviço de atendimento", NA)
base$`D15[SQ011]` = ifelse(base$`D15[SQ011]-Quais dificulda..  [Discriminação/p.. ]`=="Sim", "Discriminação/preconceito ao buscar atendimento/no atendimento", NA)
base$`D15[SQ012]` = ifelse(base$`D15[SQ012]-Quais dificulda..  [Recusa ao atend.. ]`=="Sim", "Recusa ao atendimento à pessoa com síndrome de Down", NA)
base$`D15[SQ013]` = ifelse(base$`D15[SQ013]-Quais dificulda..  [Falta de UTI (U.. ]`=="Sim", "Falta de UTI (Unidade de Terapia Intensiva) para pessoas com síndrome de Down", NA)
base$`D15[SQ014]` = ifelse(base$`D15[SQ014]-Quais dificulda..  [Precisei acessa.. ]`=="Sim", "Precisei acessar à justiça para a pessoa com síndrome de Down ter acesso a algum tratamento médico", NA)
base$`D15[SQ015]` = ifelse(base$`D15[SQ015]-Quais dificulda..  [Todas as condiç.. ]`=="Sim", "Todas as condições de saúde da pessoa foram subordinadas à síndrome de Down e, por isso, ela não recebeu o tratamento adequado", NA)
base$`D15[SQ016]` = ifelse(base$`D15[SQ016]-Quais dificulda..  [Não enfrentei d.. ]`=="Sim", "Não enfrentei dificuldades", NA)
base$`D15[SQ017]` = ifelse(base$`D15[SQ017]-Quais dificulda..  [Não sei informar]`=="Sim", "Não sei informar", NA)
base$`D15[SQ018]` = ifelse(base$`D15[SQ018]-Quais dificulda..  [Prefiro não res.. ]`=="Sim", "Prefiro não responder", NA)
base$`D15[other]` = ifelse(base$`D15[other]-Quais dificulda..  [Outros]`=="Sim", "Outros", NA)

dific_serv_saude <- base %>% 
  select(`D15[SQ001]`, `D15[SQ002]`, `D15[SQ003]`, `D15[SQ004]`, `D15[SQ005]`, `D15[SQ006]`, `D15[SQ007]`, `D15[SQ008]`, `D15[SQ009]`, `D15[SQ010]`, 
         `D15[SQ011]`, `D15[SQ012]`, `D15[SQ013]`, `D15[SQ014]`, `D15[SQ015]`, `D15[SQ016]`, `D15[SQ017]`, `D15[SQ018]`, `D15[other]`) %>% 
  pivot_longer(
    cols = c(`D15[SQ001]`, `D15[SQ002]`, `D15[SQ003]`, `D15[SQ004]`, `D15[SQ005]`, `D15[SQ006]`, `D15[SQ007]`, `D15[SQ008]`, `D15[SQ009]`, `D15[SQ010]`, 
             `D15[SQ011]`, `D15[SQ012]`, `D15[SQ013]`, `D15[SQ014]`, `D15[SQ015]`, `D15[SQ016]`, `D15[SQ017]`, `D15[SQ018]`, `D15[other]`),
    names_to = "sample",
    values_to = "Dificuldade",
    values_drop_na = TRUE
  ) %>% 
  select(-sample) %>% 
  arrange(Dificuldade) %>%
  count(Dificuldade) %>% 
  # 530 respostas válidas neste bloco
  mutate(Porcentagem = (n/530)*100) %>% 
  mutate(Porcentagem = round(Porcentagem, 2))


## Nível de escolaridade das pessoas com síndrome de Down no Distrito Federal

Nivel_escolaridade <- base %>%
  select(`H6-Qual é a escola..`, `M5-Qual é a sua es..`) %>%
  pivot_longer(
    cols = c(`H6-Qual é a escola..`, `M5-Qual é a sua es..`),
    names_to = "sample",
    values_to = "Escolaridade",
    values_drop_na = TRUE
  ) %>%
  select(-sample) %>%
  arrange(Escolaridade) %>%
  count(Escolaridade) %>% 
  mutate(Perc = prop.table(n)*100) %>% 
  mutate(Perc = round(Perc, 2))

## Respostas às afirmações sobre a estrutura escolar e o preparo dos profissionais escolares para o acolhimento e formação de pessoas com sndrome de Down
base$`E4[SQ001]` = ifelse(base$`E4[SQ001]-Abaixo, são apr..  [Os(As) professo.. ]`=="Sim", "Os(As) professores(as) são qualificados(as) para o ensino às pessoas com SD", NA)
base$`E4[SQ002]` = ifelse(base$`E4[SQ002]-Abaixo, são apr..  [Os demais profi.. ]`=="Sim", "Os demais profissionais possuem preparo para acompanhar pessoas com SD", NA)
base$`E4[SQ003]` = ifelse(base$`E4[SQ003]-Abaixo, são apr..  [A escola/instit.. ]`=="Sim", "A escola/instituição de ensino possui acessibilidade física e recursos necessários para ensinar e incluir pessoas com SD", NA)
base$`E4[SQ004]` = ifelse(base$`E4[SQ004]-Abaixo, são apr..  [A escola/instit.. ]`=="Sim", "A escola/instituição de ensino permite que o(a) aluno(a) com SD tenha um segundo professor/tutor", NA)
base$`E4[SQ005]` = ifelse(base$`E4[SQ005]-Abaixo, são apr..  [A escola oferec.. ]`=="Sim", "A escola oferece Atendimento Educacional Especializado (AEE) no contraturno quando necessário", NA)
base$`E4[SQ006]` = ifelse(base$`E4[SQ006]-Abaixo, são apr..  [Houve cobrança .. ]`=="Sim", "Houve cobrança ou tentativa de cobrança de valor adicional (sobretaxa de matrícula ou mensalidade) para inclusão da pessoa com SD na escola/instituição de ensino", NA)
base$`E4[SQ007]` = ifelse(base$`E4[SQ007]-Abaixo, são apr..  [O currículo da .. ]`=="Sim", "O currículo da escola passou por adequações metodológicas necessárias ao ensino de pessoas com SD", NA)
base$`E4[SQ008]` = ifelse(base$`E4[SQ008]-Abaixo, são apr..  [A escola/instit.. ]`=="Sim", "A escola/instituição de ensino dispõe de um(a) pedagogo(a) para acompanhar o desenvolvimento do aluno com SD", NA)
base$`E4[SQ009]` = ifelse(base$`E4[SQ009]-Abaixo, são apr..  [A escola/instit.. ]`=="Sim", "A escola/instituição de ensino dispõe de um(a) psicopedagogo(a) para acompanhar o desenvolvimento do aluno com SD", NA)

preparo_profissionais <- base %>% 
  select(`E4[SQ001]`, `E4[SQ002]`, `E4[SQ003]`, `E4[SQ004]`, `E4[SQ005]`, `E4[SQ006]`, `E4[SQ007]`, `E4[SQ008]`, `E4[SQ009]`) %>% 
  pivot_longer(
    cols = c(`E4[SQ001]`, `E4[SQ002]`, `E4[SQ003]`, `E4[SQ004]`, `E4[SQ005]`, `E4[SQ006]`, `E4[SQ007]`, `E4[SQ008]`, `E4[SQ009]`),
    names_to = "sample",
    values_to = "Dificuldade",
    values_drop_na = TRUE
  ) %>% 
  select(-sample) %>% 
  arrange(Dificuldade) %>%
  count(Dificuldade) %>% 
  # 365 respostas válidas neste bloco
  mutate(Porcentagem = (n/365)*100) %>% 
  mutate(Porcentagem = round(Porcentagem, 2)) %>% 
  select(-n)

## Respostas às afirmações sobre experiências no ambiente escolar por pessoas com síndrome de Down

exp_amb_escolar <- base %>% 
  select(`E5[SQ002]-Abaixo, são apr..  [A pessoa com sí.. ]`, `E5[SQ003]-Abaixo, são apr..  [A pessoa com sí.. ]`, 
         `E5[SQ004]-Abaixo, são apr..  [A pessoa com sí.. ]`, `E5[SQ005]-Abaixo, são apr..  [A pessoa com sí.. ]`, 
         `E5[SQ006]-Abaixo, são apr..  [A pessoa com sí.. ]`, `E5[SQ007]-Abaixo, são apr..  [A pessoa com sí.. ]`, 
         `E5[SQ008]-Abaixo, são apr..  [A pessoa com sí.. ]`, `E5[SQ009]-Abaixo, são apr..  [A pessoa com sí.. ]`, 
         `E5[SQ010]-Abaixo, são apr..  [O convívio com .. ]`, `E5[SQ011]-Abaixo, são apr..  [O convívio com .. ]`) %>% 
  pivot_longer(
    cols = c(`E5[SQ002]-Abaixo, são apr..  [A pessoa com sí.. ]`, `E5[SQ003]-Abaixo, são apr..  [A pessoa com sí.. ]`, 
             `E5[SQ004]-Abaixo, são apr..  [A pessoa com sí.. ]`, `E5[SQ005]-Abaixo, são apr..  [A pessoa com sí.. ]`, 
             `E5[SQ006]-Abaixo, são apr..  [A pessoa com sí.. ]`, `E5[SQ007]-Abaixo, são apr..  [A pessoa com sí.. ]`, 
             `E5[SQ008]-Abaixo, são apr..  [A pessoa com sí.. ]`, `E5[SQ009]-Abaixo, são apr..  [A pessoa com sí.. ]`, 
             `E5[SQ010]-Abaixo, são apr..  [O convívio com .. ]`, `E5[SQ011]-Abaixo, são apr..  [O convívio com .. ]`),
    names_to = "Variavel",
    values_to = "Resposta",
    values_drop_na = TRUE
  )

exp_amb_escolar$Titulo = ifelse(exp_amb_escolar$Variavel=="E5[SQ002]-Abaixo, são apr..  [A pessoa com sí.. ]", "A pessoa com síndrome de Down tem/teve dificuldade para acompanhar as aulas e ou atividades propostas na escola",
                                ifelse(exp_amb_escolar$Variavel=="E5[SQ003]-Abaixo, são apr..  [A pessoa com sí.. ]", "A pessoa com síndrome de Down percebe/percebeu falta de paciência por parte dos(as) professores(as)",
                                       ifelse(exp_amb_escolar$Variavel=="E5[SQ004]-Abaixo, são apr..  [A pessoa com sí.. ]", "A pessoa com síndrome de Down se sente/sentiu isolada dos colegas da escola",
                                              ifelse(exp_amb_escolar$Variavel=="E5[SQ005]-Abaixo, são apr..  [A pessoa com sí.. ]", "A pessoa com síndrome de Down se sente/sentiu discriminada (recebeu/recebe tratamento diferenciado em relação aos seus colegas que a fez/faz se sentir excluída/inferiorizada) no ambiente escolar",
                                                     ifelse(exp_amb_escolar$Variavel=="E5[SQ006]-Abaixo, são apr..  [A pessoa com sí.. ]", "A pessoa com síndrome de Down sofre/sofreu bullying (agressão física, verbal e/ou psicológica) no ambiente escolar",
                                                            ifelse(exp_amb_escolar$Variavel=="E5[SQ007]-Abaixo, são apr..  [A pessoa com sí.. ]", "A pessoa com síndrome de Down é/foi acolhida pela equipe de profissionais da escola",
                                                                   ifelse(exp_amb_escolar$Variavel=="E5[SQ008]-Abaixo, são apr..  [A pessoa com sí.. ]", "A pessoa com síndrome de Down se sente/sentiu segura na escola",
                                                                          ifelse(exp_amb_escolar$Variavel=="E5[SQ009]-Abaixo, são apr..  [A pessoa com sí.. ]", "A pessoa com síndrome de Down participa/participava das atividades desenvolvidas no contexto escolar (educativas, esportivas, culturais ou sociais)",
                                                                                 ifelse(exp_amb_escolar$Variavel=="E5[SQ010]-Abaixo, são apr..  [O convívio com .. ]", "O convívio com os colegas da escola contribui/contribuiu para a autoestima e o desenvolvimento emocional e cognitivo da pessoa com síndrome de Down",
                                                                                        ifelse(exp_amb_escolar$Variavel=="E5[SQ011]-Abaixo, são apr..  [O convívio com .. ]", "O convívio com os colegas da escola contribui/contribuiu para a socialização e a autonomia da pessoa com síndrome de Down", ""))))))))))

exp_amb_escolar <- exp_amb_escolar %>% 
  select(-Variavel) %>% 
  arrange(Titulo) %>%
  count(Titulo, Resposta) %>% 
  subset(Titulo != "") %>% 
  # 365 respostas neste bloco
  mutate(Porcentagem = (n/365)*100) %>%
  mutate(Porcentagem = round(Porcentagem, 2)) %>% 
  select(-n) %>% 
  group_by(Resposta, Titulo) %>% 
  reshape2::dcast(Resposta ~ Titulo) %>% 
  replace(is.na(.), 0)

## Respostas à pergunta “a pessoa com síndrome de Down trabalha atualmente?”
table(base$`F1-A pessoa com sí..`)

trabalha_atualmente <- base %>% 
  select(`F1-A pessoa com sí..`) %>% 
  # remove linhas com todas as colunas vazias
  remove_empty("rows") %>% 
  mutate(`F1-A pessoa com sí..`=ifelse(is.na(`F1-A pessoa com sí..`), "", `F1-A pessoa com sí..`)) %>% 
  unite(trabalha, c(`F1-A pessoa com sí..`), sep = " ", remove = FALSE) %>% 
  select(-`F1-A pessoa com sí..`) %>% 
  mutate(trabalha = stringr::str_trim(trabalha, side = "both")) %>% 
  count(trabalha) %>% 
  mutate(Perc = prop.table(n)*100) %>% 
  mutate(Perc = round(Perc, 2))

## Respostas à pergunta “por quais motivos a pessoa com síndrome de Down não trabalha/nunca trabalhou?”
base$`F2[SQ001]` = ifelse(base$`F2[SQ001]-Por qual/is mot..  [Não há oferta d.. ]`=="Sim", "Não há oferta de vagas para pessoas com síndrome de Down", NA)
base$`F2[SQ002]` = ifelse(base$`F2[SQ002]-Por qual/is mot..  [Os empregadores.. ]`=="Sim", "Os empregadores preferem não contratar pessoas com síndrome de Down", NA)
base$`F2[SQ003]` = ifelse(base$`F2[SQ003]-Por qual/is mot..  [Os empregadores.. ]`=="Sim", "Os empregadores não oferecem apoio ou treinamento às pessoas com síndrome de Down, dificultando a permanência no trabalho", NA)
base$`F2[SQ004]` = ifelse(base$`F2[SQ004]-Por qual/is mot..  [Os empregadores.. ]`=="Sim", "Os empregadores não estão preparados para lidar com pessoas com síndrome de Down no ambiente de trabalho", NA)
base$`F2[SQ005]` = ifelse(base$`F2[SQ005]-Por qual/is mot..  [O deslocamento .. ]`=="Sim", "O deslocamento para o local de trabalho é/era difícil", NA)
base$`F2[SQ006]` = ifelse(base$`F2[SQ006]-Por qual/is mot..  [A pessoa com sí.. ]`=="Sim", "A pessoa com síndrome de Down passou por situações de discriminação (recebeu/recebe tratamento diferenciado que o/a fez se sentir excluído(a)/inferiorizado(a)) no ambiente de trabalho", NA)
base$`F2[SQ007]` = ifelse(base$`F2[SQ007]-Por qual/is mot..  [A pessoa com sí.. ]`=="Sim", "A pessoa com síndrome de Down sofreu bullying (agressão física, verbal e/ou psicológica) no ambiente de trabalho", NA)
base$`F2[SQ008]` = ifelse(base$`F2[SQ008]-Por qual/is mot..  [A pessoa com sí.. ]`=="Sim", "A pessoa com síndrome de Down não possui autonomia para trabalhar", NA)
base$`F2[SQ009]` = ifelse(base$`F2[SQ009]-Por qual/is mot..  [Os familiares/r.. ]`=="Sim", "Os familiares/responsáveis não se sentem seguros em deixar a pessoa com síndrome de Down trabalhar", NA)
base$`F2[SQ010]` = ifelse(base$`F2[SQ010]-Por qual/is mot..  [A pessoa com sí.. ]`=="Sim", "A pessoa com síndrome de Down não teve acesso à formação escolar/capacitação necessária para trabalhar", NA)
base$`F2[SQ011]` = ifelse(base$`F2[SQ011]-Por qual/is mot..  [Nenhum dos moti.. ]`=="Sim", "Nenhum dos motivos anteriores", NA)
base$`F2[other]` = ifelse(base$`F2[other]-Por qual/is mot..  [Outros]`=="Sim", "Outros", NA)

motivos_nao_trabalha <- base %>% 
  select(`F2[SQ001]`, `F2[SQ002]`, `F2[SQ003]`, `F2[SQ004]`, `F2[SQ005]`, `F2[SQ006]`, `F2[SQ007]`, `F2[SQ008]`, `F2[SQ009]`, `F2[SQ010]`, `F2[SQ011]`, `F2[other]`) %>% 
  pivot_longer(
    cols = c(`F2[SQ001]`, `F2[SQ002]`, `F2[SQ003]`, `F2[SQ004]`, `F2[SQ005]`, `F2[SQ006]`, `F2[SQ007]`, `F2[SQ008]`, `F2[SQ009]`, `F2[SQ010]`, `F2[SQ011]`, `F2[other]`),
    names_to = "sample",
    values_to = "Dificuldade",
    values_drop_na = TRUE
  ) %>% 
  select(-sample) %>% 
  arrange(Dificuldade) %>%
  count(Dificuldade) %>% 
  # 140 respostas válidas neste bloco
  mutate(Porcentagem = (n/140)*100) %>% 
  mutate(Porcentagem = round(Porcentagem, 2)) %>% 
  select(-n)

## Respostas às afirmações sobre dificuldades à inclusão enfrentadas por pessoas com síndrome de Down no Distrito Federal
dific_inclusao <- base %>% 
  select(`G3[SQ001]-Avalie o quanto..  [Preconceito/dis.. ]`, `G3[SQ002]-Avalie o quanto..  [Falta de inform.. ]`, `G3[SQ003]-Avalie o quanto..  [Preconceito/não.. ]`, 
         `G3[SQ004]-Avalie o quanto..  [Falta de vagas .. ]`, `G3[SQ005]-Avalie o quanto..  [Falta de prepar.. ]`, `G3[SQ006]-Avalie o quanto..  [Falta de oportu.. ]`, 
         `G3[SQ007]-Avalie o quanto..  [Falta de profis.. ]`, `G3[SQ008]-Avalie o quanto..  [Alto custo dos .. ]`, `G3[SQ009]-Avalie o quanto..  [Falta de progra.. ]`, 
         `G3[SQ010]-Avalie o quanto..  [Condições da pr.. ]`) %>% 
  pivot_longer(
    cols = c(`G3[SQ001]-Avalie o quanto..  [Preconceito/dis.. ]`, `G3[SQ002]-Avalie o quanto..  [Falta de inform.. ]`, `G3[SQ003]-Avalie o quanto..  [Preconceito/não.. ]`, 
             `G3[SQ004]-Avalie o quanto..  [Falta de vagas .. ]`, `G3[SQ005]-Avalie o quanto..  [Falta de prepar.. ]`, `G3[SQ006]-Avalie o quanto..  [Falta de oportu.. ]`, 
             `G3[SQ007]-Avalie o quanto..  [Falta de profis.. ]`, `G3[SQ008]-Avalie o quanto..  [Alto custo dos .. ]`, `G3[SQ009]-Avalie o quanto..  [Falta de progra.. ]`, 
             `G3[SQ010]-Avalie o quanto..  [Condições da pr.. ]`),
    names_to = "Variavel",
    values_to = "Resposta",
    values_drop_na = TRUE
  )

dific_inclusao$Titulo = ifelse(dific_inclusao$Variavel=='G3[SQ001]-Avalie o quanto..  [Preconceito/dis.. ]', "Preconceito/discriminação da sociedade", 
                               ifelse(dific_inclusao$Variavel=='G3[SQ002]-Avalie o quanto..  [Falta de inform.. ]', "Falta de informação/desconhecimento da sociedade sobre a síndrome de Down", 
                                      ifelse(dific_inclusao$Variavel=='G3[SQ003]-Avalie o quanto..  [Preconceito/não.. ]', "Preconceito/não aceitação por parte da família da pessoa com síndrome de Down", 
                                             ifelse(dific_inclusao$Variavel=='G3[SQ004]-Avalie o quanto..  [Falta de vagas .. ]', "Falta de vagas para pessoas com síndrome de Down em escolas regulares", 
                                                    ifelse(dific_inclusao$Variavel=='G3[SQ005]-Avalie o quanto..  [Falta de prepar.. ]', "Falta de preparo e de qualificação das escolas para receber alunos com síndrome de Down", 
                                                           ifelse(dific_inclusao$Variavel=='G3[SQ006]-Avalie o quanto..  [Falta de oportu.. ]', "Falta de oportunidades para pessoas com síndrome de Down no mercado de trabalho", 
                                                                  ifelse(dific_inclusao$Variavel=='G3[SQ007]-Avalie o quanto..  [Falta de profis.. ]', "Falta de profissionais de saúde na rede pública especializados no cuidado de pessoas com síndrome de Down", 
                                                                         ifelse(dific_inclusao$Variavel=='G3[SQ008]-Avalie o quanto..  [Alto custo dos .. ]', "Alto custo dos tratamentos/acompanhamentos médicos", 
                                                                                ifelse(dific_inclusao$Variavel=='G3[SQ009]-Avalie o quanto..  [Falta de progra.. ]', "Falta de programas de governo/políticas públicas voltados para pessoas com síndrome de Down", 
                                                                                       ifelse(dific_inclusao$Variavel=='G3[SQ010]-Avalie o quanto..  [Condições da pr.. ]', "Condições da própria síndrome (dificuldade de fala, dificuldade de locomoção, dificuldade de aprendizagem, etc)", NA))))))))))

dific_inclusao <- dific_inclusao %>% 
  select(-Variavel) %>% 
  arrange(Titulo) %>%
  count(Titulo, Resposta) %>% 
  subset(Titulo != "") %>% 
  # 430 respostas neste bloco
  mutate(Porcentagem = (n/430)*100) %>%
  mutate(Porcentagem = round(Porcentagem, 2)) %>% 
  select(-n) %>% 
  group_by(Resposta, Titulo) %>% 
  reshape2::dcast(Resposta ~ Titulo) %>% 
  replace(is.na(.), 0)