12.1 Arrumando banco de dados: o pacote janitor

O pacote janitor disponibiliza algumas funções para limpar bases de dados.

Primeiramente, instale e carregue o pacote:

install.packages("janitor") # Instale a versão do CRAN, OU:

devtools::install_github("sfirke/janitor") # Versão de desenvolvimento
library(tidyverse)
library(janitor)

12.1.1 Arrumando o nome das variáveis

Utilizaremos a base com informações de pacientes com arritmia cardíaca. O código para obter essa base de exemplo está descrito abaixo:

# URL da base para baixar
url <- "https://github.com/curso-r/livro-material/raw/master/assets/data/dados_nomes_variaveis_livro.xlsx"

# Arquivo de destino onde a base deverá ser salva
destfile <- "dados_nomes_variaveis_livro.xlsx"

# Fazer o download da base
curl::curl_download(url, destfile)

# Importar a base 
dados_brutos <- readxl::read_excel(destfile)

As variáveis presentes na base são:

dados_brutos %>% names() # Ver o nome das variáveis
##  [1] "ID"                   "Sexo"                 "Nascimento"          
##  [4] "Idade"                "Inclusão"             "Cor"                 
##  [7] "Peso"                 "Altura"               "cintura"             
## [10] "IMC"                  "Superfície corporal"  "Tabagismo"           
## [13] "cg.tabag (cig/dia)"   "Alcool (dose/semana)" "Drogas ilícitas"     
## [16] "Cafeína/dia"          "Refrig/dia"           "Sedentario"          
## [19] "ativ. Fisica"

Os nomes das variáveis contém letras maiúsculas, acentos, parênteses, pontos e barras, o que atrapalha na hora da programação. Para resolver esse problema, usamos a função clean_names().

dados_brutos %>% 
  janitor::clean_names() %>% # Limpar os nomes das variáveis
  names() # Ver o nome das variáveis
##  [1] "id"                  "sexo"                "nascimento"         
##  [4] "idade"               "inclusao"            "cor"                
##  [7] "peso"                "altura"              "cintura"            
## [10] "imc"                 "superficie_corporal" "tabagismo"          
## [13] "cg_tabag_cig_dia"    "alcool_dose_semana"  "drogas_ilicitas"    
## [16] "cafeina_dia"         "refrig_dia"          "sedentario"         
## [19] "ativ_fisica"

Veja que a função removeu os parênteses, pontos e barras e substituiu os espaços por _.

E para substituir na base, precisamos atribuir o resultado em um novo objeto:

dados <- dados_brutos %>% 
  janitor::clean_names() # Limpar os nomes das variáveis

12.1.2 Removendo linhas e colunas vazias

Esse banco de dados também contém outro problema: linhas vazias. Na verdade, elas não eram completamente vazias, pois havia algumas informações de identificação da(o) paciente, mas nenhuma outra variável tinha sido computada.

dados %>% 
  dplyr::slice(3) %>% # Apresentar apenas a linha 3
  knitr::kable()
id sexo nascimento idade inclusao cor peso altura cintura imc superficie_corporal tabagismo cg_tabag_cig_dia alcool_dose_semana drogas_ilicitas cafeina_dia refrig_dia sedentario ativ_fisica
3 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA

Para resolver o problema, é possível utilizar a função remove_empty().

dados_sem_linhas_vazias <- dados %>% 
  as.data.frame() %>% 
  dplyr::select(-id) %>% 
  janitor::remove_empty() %>% 
  tibble::rowid_to_column("id") %>% 
  dplyr::select(id, everything()) %>%
  tibble::as_tibble()

dados_sem_linhas_vazias %>% knitr::kable()
id sexo nascimento idade inclusao cor peso altura cintura imc superficie_corporal tabagismo cg_tabag_cig_dia alcool_dose_semana drogas_ilicitas cafeina_dia refrig_dia sedentario ativ_fisica
1 F 1964-01-31 41 2006-02-17 branca 75 1.63 98 28.22839 1.81 N 0 0 N 100 0 S N
2 M 1959-01-28 45 2005-11-29 negra 71 1.70 NT 24.57000 1.83 N 0 35 N 50 300 N insuficiente
3 M 1957-09-13 50 2008-02-13 NT 80 1.64 NT 29.74420 1.87 N 0 0 N 500 0 S N
4 F 1938-02-06 71 2009-06-25 parda 56 1.51 96 24.56033 1,51 N 0 0 N 50 0 S N

Foi necessário converter para data.frame primeiro porque não é possível definir os nomes das linhas de uma tibble. Se a linha estivesse completamente vazia, bastaria usar diretamente a função remove_empty_rows().

Equivalentemente para colunas, existe a função remove_empty_cols().

Outra forma de realizar este mesmo procedimento é utilizando a função drop_na() do pacote tidyr:

dados_sem_linhas_vazias <- dados %>% 
  tidyr::drop_na(-id) %>%
  dplyr::select(-id) %>% 
  tibble::rowid_to_column("id") %>% 
  dplyr::relocate(id, .before = sexo)

dados_sem_linhas_vazias %>% knitr::kable()
id sexo nascimento idade inclusao cor peso altura cintura imc superficie_corporal tabagismo cg_tabag_cig_dia alcool_dose_semana drogas_ilicitas cafeina_dia refrig_dia sedentario ativ_fisica
1 F 1964-01-31 41 2006-02-17 branca 75 1.63 98 28.22839 1.81 N 0 0 N 100 0 S N
2 M 1959-01-28 45 2005-11-29 negra 71 1.70 NT 24.57000 1.83 N 0 35 N 50 300 N insuficiente
3 M 1957-09-13 50 2008-02-13 NT 80 1.64 NT 29.74420 1.87 N 0 0 N 500 0 S N
4 F 1938-02-06 71 2009-06-25 parda 56 1.51 96 24.56033 1,51 N 0 0 N 50 0 S N

12.1.3 Identificando linhas duplicadas

O pacote janitor possui uma função para identificar entradas duplicadas numa base de dados: get_dupes(). Vamos criar uma base genérica para testá-la.

# Criar a base de exemplo
p_nome <- c("Athos", "Daniel", "Fernando", "Julio", "William")
sobrenome <- c("Damiani", "Falbel", "Corrêa", "Trecenti", "Amorim")

base_exemplo <- tibble::tibble(
  nome = sample(p_nome, 25, replace = TRUE),
  sobrenome = sample(sobrenome, 25, replace = TRUE),
  variavel_importante = rnorm(25)
)
# Dar uma espiada na base de exemplo
dplyr::glimpse(base_exemplo)
## Rows: 25
## Columns: 3
## $ nome                <chr> "William", "Julio", "Athos", "Julio", "Athos", "Fe…
## $ sobrenome           <chr> "Amorim", "Damiani", "Damiani", "Amorim", "Trecent…
## $ variavel_importante <dbl> 1.87048800, -0.62852366, 0.19517563, 0.59525152, 0…
# Ver as duplicatas com a função get_dupes()
janitor::get_dupes(base_exemplo, nome, sobrenome)
## # A tibble: 14 x 4
##    nome     sobrenome dupe_count variavel_importante
##    <chr>    <chr>          <int>               <dbl>
##  1 Athos    Amorim             2             -0.808 
##  2 Athos    Amorim             2             -1.26  
##  3 Fernando Falbel             3             -0.344 
##  4 Fernando Falbel             3             -1.62  
##  5 Fernando Falbel             3              1.07  
##  6 Julio    Damiani            2             -0.629 
##  7 Julio    Damiani            2             -1.37  
##  8 Julio    Falbel             2              1.05  
##  9 Julio    Falbel             2              0.0674
## 10 William  Amorim             2              1.87  
## 11 William  Amorim             2             -0.639 
## 12 William  Corrêa             3              0.636 
## 13 William  Corrêa             3             -1.58  
## 14 William  Corrêa             3             -0.0254

Todas as linhas na tibble resultante representam uma combinação de nome-sobrenome repetida.

12.1.4 Outras funções

Por fim, o janitor também tem funções equivalentes à table() para produzir tabelas de frequência:

  • tabyl() - similar a table(), mas pipe-ável e com mais recursos.
  • adorn_totals() - acrescenta o total das linhas ou colunas.
# Fazer uma tabela de frequência
mtcars %>% janitor::tabyl(cyl)
##  cyl  n percent
##    4 11 0.34375
##    6  7 0.21875
##    8 14 0.43750
# Fazer uma tabela de frequência com valores totais,
# e porcentagem
mtcars %>% 
  janitor::tabyl(cyl) %>% 
  janitor::adorn_totals()
##    cyl  n percent
##      4 11 0.34375
##      6  7 0.21875
##      8 14 0.43750
##  Total 32 1.00000
# Fazer uma tabela de frequência com duas variáveis,
# e valores totais
mtcars %>% 
  janitor::tabyl(cyl, am) %>% 
  janitor::adorn_totals(where = "col")
##  cyl  0 1 Total
##    4  3 8    11
##    6  4 3     7
##    8 12 2    14

Esperamos que essas dicas e o pacote janitor ajudem a agilizar as suas análises!


Curso-R