Codificación y Transformación

Como ya hemos visto brevemente en el capítulo de tipos de variables, normalmente los datos de un archivo necesitarán ser cuidadosamente preparados antes de poder aplicar, con las mínimas garantías, cualquier análisis estadístico.

La preparación del archivo de datos incluye desde la simple detección y corrección de los posibles errores cometidos al introducir los datos, hasta sofisticadas transformaciones, pasando por la recodificación de los códigos utilizados para los valores de una o más variables, o la creación de nuevas variables a partir de otras ya existentes.

Para ilustrar esta parte, vamos a leer la siguiente base de datos:

data <- read.csv("https://pos.it/r4ds-students-csv")
str(data)

Ahora tendremos que comprobar que los datos se han leído correctamente y no tienen errores

Detección inicial de errores más comunes

1. El nombre de las columnas debe tener un nombre consistente. Las dos palabaras que corresponden a la columnas: Student.ID, Full.Name y favourite.food están separadas por un . y mealPlan se diferencia por una mayúscula. Además todas están en minúsculas y AGE en mayúsculas. Para evitar futuros errores, los nombres deberán ser consistentes.

colnames(data)<-
    c("Student_id","Full_name","Favourite_food","Meal_plan","Age")
colnames(data)           

[1] "Student_id"     "Full_name"      "Favourite_food" "Meal_plan"      "Age"

1. Clasificar las variables en cuantitativas y cualitativas y declararlas de forma correcta

data$Meal_plan<-as.factor(data$Meal_plan)s
tr(data)
data$Age<-as.numeric(data$Age) ###Cuidado cunando salta este warning porque hay algo raro

1. Asegurarnos que los NA se han leído correctamente: Columna Favourite_food tiene un N/A que no está detectado como dato "missing" y un vacío en Age

data$Favourite_food<-car::recode(data$Favourite_food,"'N/A'= NA")
data[data == ""] <- NA ##Esto reemplaza todos los valores vacíos a NA
data

  Student_id        Full_name     Favourite_food     Meal_plan            Age
  1              Sunil Huffmann  Strawberry yoghurt  Lunch only            4            
  2              Barclay Lynn       French fries     Lunch only            5          
  3              Jayendra Lyne          <NA>         Breakfast and lunch   7          
  4              Leon Rossini         Anchovies      Lunch only            <NA>          
  5              Chidiegwu Dunkel    Pizza           Breakfast and lunch   Five          
  6              Güvenç Attila         Ice cream     Lunch only            6

1. Asegurarnos que no hay valores incorrectos: En la variable Age, una de las observaciones tiene "five" en vez del número, esto nos causará problemas a la hora de usar dicha variable como variable numérica como hemos visto al declararla con as.numeric().

data$Age<-car::recode(data$Age,"'five'='5'")
data
data$Age<-as.numeric(data$Age)

Transformación de variables

En ocasiones nos puede interesar incorporar nuevas variables a nuestra base de datos, como ya vimos con la creación de la variable edad

Ejemplo: Para este ejemplo vamos a utilizar el dataset TiposDatos.csv

278B

TiposDatos.csv

Aprovechamos para hacer un brese repaso de los tipos de variables y control de errores

datos <- read.csv("TiposDatos.csv")
str(datos)

## 1. Declaramos de forma consistente el nombre de las variables
colnames(datos)[c(2,6)]<-c("ID","Estado_civil")

## 2. Transformamos las variables que no están bien leídas
datos$Fecha_de_nacimiento<-as.Date(datos$Fecha_de_nacimiento,format="%d/%m/%y")
datos$Estado_civil<-as.factor(datos$Estado_civil)
str(datos)

Cuando una base de datos contiene la fecha de nacimiento, es mucho más útil generar la edad, más fácil de interpretar, analizar como variable en si misma y manejar numéricamente:

##Calcular la edad
datos$Edad<-as.numeric(floor((Sys.Date()-datos$Fecha_de_nacimiento)/365.25))
str(datos)

Podemos querer crear una variable en la base de datos TiposDatos.csv que se llame índice de masa corporal que se obtiene de la siguiente forma:

$$IMC=peso(kg)/altura(m)^2$$

##Crear IMC
datos$IMC  <- datos$Peso / (datos$Altura/100)^2

Categorizar una variable cuantitativa en una cualitativa

Esta opción suele ser muy habitual ya que muchas veces queremos trabajar con las categorías de una variable en vez de la variable cuantitativa.

Ejemplo: La variable cuantitativa IMC que acabamos de crear se suele usar en cuatro grupos:

• bajo peso: IMC <= 18.50

• normal: 18.5 < IMC <=25

• sobrepeso: 25 < IMC <= 30

• obesidad: IMC >30

En realidad buscamos los siguientes tramos de IMC: (minimo, 18.5], (18.5, 25], (25, 30], (30, máximo]

datos$IMC_grupos <- cut(datos$IMC, breaks = c(-Inf,18.5,25,30,Inf), 
                labels = c("bajo peso", "normal", "sobrepreso", "obesidad"))

¿Cuántas personas tenemos en cada uno de los grupos?

table(datos$IMC_grupos)

bajo peso     normal     sobrepreso   obesidad          
    0          3              1          1 

Recodificación de variables

Nos puede interesar trabajar de forma diferente sobre una de las variables que ya tenemos

Ejemplo: Queremos ver la probabilidad de estar casado o de tener un determinado peso si la persona tiene más de 40 años o menos y queremos poder predecir dicha variable, para aplicar en este caso una regresión logística necesitaremos tener una variable codificada como: 0 (<40), 1 (>=40). Esto con el comando if_else() del paquete dplyr es muy fácil:

##Recodificación de variables
library(dplyr)
table(datos$Edad)
datos$Edad2cat<- if_else(datos$Edad < 40, 0, 1)

¿Qué tipo de variable es Edad2Cat?

Es una variable cualitativa que toma dos categorías: 0 (<40), 1 (>=40)

Es una buena práctica recodificar siempre en distintas variables para no perder información.

Ahora es tu turno: Realiza los siguientes ejercicios:

Ejercicio: Codificación y Transformación