Capítulo 6 Práctica 2

6.1 Preparamos R

Cargamos el paquete de la interfaz gráfica:

install.packages("Rcmdr")
library('Rcmdr')

Sino, puedes cargar una de sus dependencias (aunque tampoco hace falta):

# install.packages("RcmdrMisc")
library(RcmdrMisc)

Y los datos de la sesión:

load('../../data/Datos1.RData')
load('../../data/Framingham.RData')

Especificamos el directorio de trabajo:

setwd("C:/Users/someuser")
getwd() # comprobamos el dir

Visualizamos los datos cargados:

head(Datos1)
##     Sexo Edad Peso Altura
## 1  Mujer   24   58    156
## 2  Mujer   35   75    170
## 3 Hombre   26   62    175
## 4  Mujer   62   61    169
## 5 Hombre   50   93    180
## 6 Hombre   31   67    171
head(Framingham)
##   OBS SEX CHD AGE SBP SBP10 DBP CHOL FRW CIG YRS_CHD YRS_DTH DEATH CAUSE
## 1   1   1   0  45  90   112  50  216  76   5      18      18     0     0
## 2   2   1   0  49 100   108  64  237  97   0      18      18     0     0
## 3   3   1   0  47 100   112  70  215  86  50      18      18     0     0
## 4   4   1   0  48 108   120  70  340  93   0      18      18     0     0
## 5   5   1   0  49 108   118  75  149  95   0      18      18     0     0
## 6   6   1   0  47 108   110  68  165  88   0      18      18     0     0
nrow(Framingham) # n de filas
## [1] 1405

También podríamos obtener más detalles con:

summary(Framingham)

6.2 Cambios de tipo de variable

HDentro de Datos1 hay que especificar que variables son de tipo factores:

Datos1 <- within(Datos1, {
  Sexo <- as.factor(Sexo)
})

6.3 Ordenando la información

  1. Guardamos el orden original de información por si es de interés posteriormente:
Datos1$posOriginal <- 1:nrow(Datos1)
  1. Ordenamos los datos en base a alguna de las vataibles:
Datos1 <- with(Datos1, Datos1[order(Sexo, decreasing = T), ]) 

6.4 Sobreescribir variables

En la base de datos de Framingham, sustituimos las columnas de tipo cualitativo, asignando los nuevos valores:

Framingham <- within(Framingham, {
  SEX <- factor(SEX, labels=c('Hombre', 'Mujer'))
  CHD <- factor(CHD, labels=c('Sin evidencias de CHD', 'Casos prevalentes', 'Casos incidentes 2', 'Casos incidentes 3', 'Casos incidentes 4', 'Casos incidentes 5', 'Casos incidentes 6', 'Casos incidentes 7', 'Casos incidentes 8', 'Casos incidentes 9', 'Casos incidentes 10'))
  DEATH <- factor(DEATH, labels=c('Vivo en el examen 10', 'Fallecido en el examen 2', 'Fallecido en el examen 3', 'Fallecido en el examen 4', 'Fallecido en el examen 5', 'Fallecido en el examen 6', 'Fallecido en el examen 7', 'Fallecido en el examen 8', 'Fallecido en el examen 9', 'Fallecido en el examen 10'))
  CAUSE <- factor(CAUSE, labels=c('Vivo en el examen 10', 'CHD súbita', 'CHD no súbita', 'ACV', 'Otras CHD', 'Cáncer', 'Otras enfermedades'))
})

6.5 Creación de variables dentro de tablas

6.5.1 Operando otra varaible

Creamos una nueva variable AGE.MONTH operando otra variable:

Framingham$AGE.MONTH <- with(Framingham, AGE*12)

Repetimos, calculando el IMC en Data1:

Datos1$IMC <- with(Datos1, Peso/((Altura/100)^2))

6.5.2 Al comparar con valores de otras variables

Variable en función de otra variable:

Framingham$Género[Framingham$SEX=="Hombre"] <- "Masculino"
Framingham$Género[Framingham$SEX=="Mujer"] <- "Femenino"

También se puede hacer especificando varais condicones:

Framingham$AGE.CAT[Framingham$AGE<=49] <- "Menos de 50"
Framingham$AGE.CAT[Framingham$AGE>49 && Framingham$AGE <=54] <- "49 a 54"
## Warning in Framingham$AGE > 49 && Framingham$AGE <= 54: 'length(x) = 1405 > 1'
## in coercion to 'logical(1)'
Framingham$AGE.CAT[Framingham$AGE>55 && Framingham$AGE <=59] <- "54 a 59"
## Warning in Framingham$AGE > 55 && Framingham$AGE <= 59: 'length(x) = 1405 > 1'
## in coercion to 'logical(1)'
Framingham$AGE.CAT[Framingham$AGE>=60] <- "Mayores de 60"

Asignar valores a CHD2 en función del valor de CHD implica un paso previo, cambiar el tipo de variable:

# Casos y no casos de enfermedad coronaria
Framingham$CHD2[as.character(Framingham$CHD)=="Sin evidencias de CHD"] <- "No caso"
Framingham$CHD2[as.character(Framingham$CHD)!="Sin evidencias de CHD"] <- "Caso"

Para el resto de vars:

# Clasificación según el hábito tabaquero
Framingham$CIG2[Framingham$CIG==0] <- "No fuma"
Framingham$CIG2[Framingham$CIG>0] <- "Fuma"

# SBP2 Clasi según presión arterial sistólica recode
Framingham$SBP2[Framingham$SBP<140] <- "Normotenso"
Framingham$SBP2[Framingham$SBP>140] <- "Hipercolesterolemia"

# DBP2
Framingham$DBP2[Framingham$DBP<90] <- "Normotenso"
Framingham$DBP2[Framingham$DBP>=90] <- "Hipotenso"

# CHOL2
Framingham$CHOL2[Framingham$CHOL<200] <- "Normocolestoremia"
Framingham$CHOL2[Framingham$CHOL>=200] <- "Hipercolesterolemia"

# /!\

6.5.3 Agrupando datos

Crear una variable agrupando el valor de otros datos:

# Framingham$AGE.CAT3 <- with(Framingham, binVariable(AGE, bins=3, method='intervals', labels = c('[45, 50)', '[50. 56)', '[56 62)')))
# ^ Requires dependency of R GUI package

#                  start point (not included)   just to specify 1st [
#                              ↓                ↓
cut(Framingham$AGE, breaks = c(44,50,56,62), labels = c('[45, 50)', '[50, 56)', '[56, 62)'))
##    [1] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50)
##    [9] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50)
##   [17] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50)
##   [25] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50)
##   [33] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50)
##   [41] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50)
##   [49] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50)
##   [57] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50)
##   [65] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50)
##   [73] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50)
##   [81] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50)
##   [89] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50)
##   [97] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50)
##  [105] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50)
##  [113] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50)
##  [121] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50)
##  [129] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50)
##  [137] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50)
##  [145] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50)
##  [153] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50)
##  [161] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50)
##  [169] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50)
##  [177] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50)
##  [185] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50)
##  [193] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50)
##  [201] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50)
##  [209] [45, 50) [45, 50) [45, 50) [45, 50) [50, 56) [50, 56) [50, 56) [50, 56)
##  [217] [45, 50) [50, 56) [50, 56) [50, 56) [50, 56) [50, 56) [50, 56) [50, 56)
##  [225] [45, 50) [50, 56) [50, 56) [45, 50) [50, 56) [50, 56) [50, 56) [50, 56)
##  [233] [50, 56) [50, 56) [50, 56) [45, 50) [50, 56) [50, 56) [45, 50) [50, 56)
##  [241] [50, 56) [50, 56) [50, 56) [50, 56) [50, 56) [45, 50) [50, 56) [45, 50)
##  [249] [45, 50) [50, 56) [50, 56) [50, 56) [45, 50) [45, 50) [50, 56) [50, 56)
##  [257] [50, 56) [50, 56) [50, 56) [50, 56) [50, 56) [50, 56) [50, 56) [50, 56)
##  [265] [50, 56) [50, 56) [45, 50) [45, 50) [50, 56) [50, 56) [45, 50) [50, 56)
##  [273] [50, 56) [50, 56) [50, 56) [50, 56) [50, 56) [50, 56) [50, 56) [50, 56)
##  [281] [45, 50) [50, 56) [45, 50) [50, 56) [50, 56) [50, 56) [45, 50) [45, 50)
##  [289] [50, 56) [50, 56) [50, 56) [50, 56) [50, 56) [45, 50) [50, 56) [50, 56)
##  [297] [50, 56) [50, 56) [50, 56) [50, 56) [50, 56) [45, 50) [50, 56) [50, 56)
##  [305] [50, 56) [50, 56) [50, 56) [45, 50) [50, 56) [50, 56) [50, 56) [50, 56)
##  [313] [50, 56) [45, 50) [45, 50) [50, 56) [45, 50) [50, 56) [50, 56) [50, 56)
##  [321] [50, 56) [50, 56) [50, 56) [50, 56) [50, 56) [45, 50) [50, 56) [50, 56)
##  [329] [50, 56) [50, 56) [50, 56) [50, 56) [50, 56) [50, 56) [50, 56) [50, 56)
##  [337] [45, 50) [50, 56) [50, 56) [50, 56) [50, 56) [50, 56) [50, 56) [50, 56)
##  [345] [45, 50) [50, 56) [45, 50) [50, 56) [45, 50) [50, 56) [50, 56) [45, 50)
##  [353] [50, 56) [50, 56) [50, 56) [50, 56) [50, 56) [50, 56) [45, 50) [50, 56)
##  [361] [50, 56) [45, 50) [50, 56) [50, 56) [50, 56) [50, 56) [50, 56) [50, 56)
##  [369] [50, 56) [45, 50) [50, 56) [50, 56) [50, 56) [50, 56) [50, 56) [50, 56)
##  [377] [50, 56) [50, 56) [50, 56) [50, 56) [50, 56) [45, 50) [50, 56) [50, 56)
##  [385] [50, 56) [50, 56) [50, 56) [50, 56) [45, 50) [50, 56) [50, 56) [45, 50)
##  [393] [50, 56) [50, 56) [50, 56) [45, 50) [45, 50) [50, 56) [50, 56) [45, 50)
##  [401] [50, 56) [50, 56) [50, 56) [50, 56) [50, 56) [50, 56) [50, 56) [50, 56)
##  [409] [45, 50) [45, 50) [50, 56) [50, 56) [45, 50) [50, 56) [50, 56) [50, 56)
##  [417] [50, 56) [50, 56) [50, 56) [50, 56) [50, 56) [50, 56) [50, 56) [56, 62)
##  [425] [50, 56) [56, 62) [56, 62) [50, 56) [50, 56) [50, 56) [56, 62) [56, 62)
##  [433] [50, 56) [56, 62) [56, 62) [50, 56) [50, 56) [56, 62) [56, 62) [56, 62)
##  [441] [56, 62) [50, 56) [56, 62) [50, 56) [50, 56) [56, 62) [56, 62) [56, 62)
##  [449] [50, 56) [56, 62) [56, 62) [56, 62) [56, 62) [56, 62) [50, 56) [50, 56)
##  [457] [50, 56) [56, 62) [56, 62) [56, 62) [50, 56) [50, 56) [50, 56) [56, 62)
##  [465] [56, 62) [50, 56) [56, 62) [56, 62) [50, 56) [56, 62) [50, 56) [50, 56)
##  [473] [56, 62) [56, 62) [50, 56) [56, 62) [56, 62) [50, 56) [50, 56) [56, 62)
##  [481] [50, 56) [56, 62) [56, 62) [50, 56) [56, 62) [50, 56) [50, 56) [50, 56)
##  [489] [56, 62) [56, 62) [56, 62) [56, 62) [50, 56) [50, 56) [50, 56) [56, 62)
##  [497] [50, 56) [50, 56) [56, 62) [56, 62) [56, 62) [56, 62) [56, 62) [56, 62)
##  [505] [50, 56) [50, 56) [56, 62) [56, 62) [56, 62) [50, 56) [56, 62) [56, 62)
##  [513] [56, 62) [56, 62) [50, 56) [56, 62) [50, 56) [56, 62) [56, 62) [56, 62)
##  [521] [50, 56) [50, 56) [56, 62) [56, 62) [50, 56) [56, 62) [56, 62) [50, 56)
##  [529] [56, 62) [56, 62) [50, 56) [50, 56) [50, 56) [50, 56) [56, 62) [50, 56)
##  [537] [56, 62) [50, 56) [56, 62) [56, 62) [56, 62) [56, 62) [50, 56) [56, 62)
##  [545] [50, 56) [56, 62) [50, 56) [56, 62) [50, 56) [56, 62) [50, 56) [56, 62)
##  [553] [56, 62) [56, 62) [56, 62) [56, 62) [56, 62) [56, 62) [50, 56) [50, 56)
##  [561] [56, 62) [50, 56) [56, 62) [56, 62) [56, 62) [56, 62) [50, 56) [56, 62)
##  [569] [56, 62) [50, 56) [56, 62) [50, 56) [56, 62) [56, 62) [50, 56) [50, 56)
##  [577] [56, 62) [56, 62) [56, 62) [56, 62) [56, 62) [50, 56) [56, 62) [56, 62)
##  [585] [50, 56) [56, 62) [56, 62) [50, 56) [50, 56) [56, 62) [56, 62) [50, 56)
##  [593] [56, 62) [50, 56) [56, 62) [56, 62) [56, 62) [50, 56) [56, 62) [56, 62)
##  [601] [50, 56) [50, 56) [56, 62) [56, 62) [56, 62) [50, 56) [56, 62) [56, 62)
##  [609] [56, 62) [56, 62) [56, 62) [56, 62) [56, 62) [56, 62) [56, 62) [56, 62)
##  [617] [56, 62) [56, 62) [56, 62) [56, 62) [56, 62) [56, 62) [56, 62) [56, 62)
##  [625] [56, 62) [56, 62) [56, 62) [56, 62) [56, 62) [56, 62) [56, 62) [56, 62)
##  [633] [56, 62) [56, 62) [56, 62) [56, 62) [56, 62) [56, 62) [56, 62) [56, 62)
##  [641] [56, 62) [56, 62) [56, 62) [56, 62) [56, 62) [56, 62) [56, 62) [56, 62)
##  [649] [56, 62) [56, 62) [56, 62) [56, 62) [56, 62) [56, 62) [56, 62) [56, 62)
##  [657] [56, 62) [56, 62) [56, 62) [56, 62) [56, 62) [56, 62) [56, 62) [56, 62)
##  [665] [56, 62) [56, 62) [56, 62) [56, 62) [56, 62) [45, 50) [45, 50) [45, 50)
##  [673] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50)
##  [681] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50)
##  [689] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50)
##  [697] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50)
##  [705] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50)
##  [713] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50)
##  [721] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50)
##  [729] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50)
##  [737] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50)
##  [745] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50)
##  [753] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50)
##  [761] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50)
##  [769] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50)
##  [777] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50)
##  [785] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50)
##  [793] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50)
##  [801] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50)
##  [809] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50)
##  [817] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50)
##  [825] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50)
##  [833] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50)
##  [841] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50)
##  [849] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50)
##  [857] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50)
##  [865] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50)
##  [873] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50)
##  [881] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50)
##  [889] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50)
##  [897] [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [45, 50) [50, 56)
##  [905] [45, 50) [50, 56) [50, 56) [50, 56) [50, 56) [50, 56) [50, 56) [45, 50)
##  [913] [50, 56) [50, 56) [45, 50) [50, 56) [50, 56) [50, 56) [45, 50) [50, 56)
##  [921] [50, 56) [50, 56) [50, 56) [45, 50) [50, 56) [50, 56) [50, 56) [45, 50)
##  [929] [50, 56) [50, 56) [50, 56) [50, 56) [50, 56) [50, 56) [50, 56) [50, 56)
##  [937] [50, 56) [50, 56) [50, 56) [50, 56) [50, 56) [45, 50) [50, 56) [50, 56)
##  [945] [50, 56) [50, 56) [45, 50) [50, 56) [45, 50) [50, 56) [50, 56) [50, 56)
##  [953] [50, 56) [50, 56) [50, 56) [45, 50) [50, 56) [50, 56) [50, 56) [45, 50)
##  [961] [50, 56) [50, 56) [50, 56) [50, 56) [50, 56) [50, 56) [50, 56) [50, 56)
##  [969] [45, 50) [50, 56) [50, 56) [45, 50) [50, 56) [50, 56) [45, 50) [50, 56)
##  [977] [45, 50) [45, 50) [50, 56) [50, 56) [50, 56) [50, 56) [45, 50) [45, 50)
##  [985] [50, 56) [50, 56) [50, 56) [50, 56) [50, 56) [50, 56) [45, 50) [50, 56)
##  [993] [50, 56) [50, 56) [45, 50) [50, 56) [50, 56) [45, 50) [50, 56) [50, 56)
## [1001] [50, 56) [45, 50) [45, 50) [50, 56) [50, 56) [50, 56) [50, 56) [50, 56)
## [1009] [50, 56) [50, 56) [45, 50) [50, 56) [45, 50) [50, 56) [50, 56) [50, 56)
## [1017] [50, 56) [45, 50) [50, 56) [50, 56) [50, 56) [50, 56) [45, 50) [50, 56)
## [1025] [50, 56) [50, 56) [50, 56) [45, 50) [50, 56) [45, 50) [50, 56) [50, 56)
## [1033] [50, 56) [50, 56) [45, 50) [50, 56) [50, 56) [50, 56) [50, 56) [50, 56)
## [1041] [50, 56) [50, 56) [50, 56) [45, 50) [50, 56) [45, 50) [50, 56) [45, 50)
## [1049] [50, 56) [45, 50) [50, 56) [50, 56) [50, 56) [50, 56) [50, 56) [50, 56)
## [1057] [50, 56) [50, 56) [50, 56) [50, 56) [50, 56) [50, 56) [50, 56) [45, 50)
## [1065] [50, 56) [45, 50) [45, 50) [50, 56) [50, 56) [45, 50) [45, 50) [50, 56)
## [1073] [45, 50) [50, 56) [50, 56) [50, 56) [50, 56) [50, 56) [45, 50) [45, 50)
## [1081] [50, 56) [50, 56) [50, 56) [50, 56) [50, 56) [45, 50) [50, 56) [50, 56)
## [1089] [45, 50) [50, 56) [50, 56) [50, 56) [50, 56) [50, 56) [45, 50) [50, 56)
## [1097] [50, 56) [45, 50) [50, 56) [50, 56) [45, 50) [45, 50) [50, 56) [50, 56)
## [1105] [50, 56) [50, 56) [45, 50) [50, 56) [50, 56) [50, 56) [50, 56) [50, 56)
## [1113] [50, 56) [50, 56) [50, 56) [50, 56) [50, 56) [45, 50) [50, 56) [50, 56)
## [1121] [50, 56) [50, 56) [50, 56) [50, 56) [50, 56) [50, 56) [50, 56) [50, 56)
## [1129] [50, 56) [45, 50) [45, 50) [45, 50) [50, 56) [50, 56) [50, 56) [50, 56)
## [1137] [50, 56) [50, 56) [56, 62) [50, 56) [50, 56) [56, 62) [50, 56) [56, 62)
## [1145] [50, 56) [56, 62) [56, 62) [56, 62) [56, 62) [56, 62) [50, 56) [50, 56)
## [1153] [56, 62) [50, 56) [50, 56) [56, 62) [56, 62) [50, 56) [50, 56) [56, 62)
## [1161] [56, 62) [56, 62) [56, 62) [56, 62) [50, 56) [56, 62) [50, 56) [50, 56)
## [1169] [50, 56) [50, 56) [56, 62) [50, 56) [50, 56) [56, 62) [50, 56) [50, 56)
## [1177] [56, 62) [50, 56) [56, 62) [56, 62) [50, 56) [56, 62) [50, 56) [50, 56)
## [1185] [50, 56) [50, 56) [56, 62) [56, 62) [56, 62) [50, 56) [50, 56) [50, 56)
## [1193] [50, 56) [50, 56) [56, 62) [56, 62) [50, 56) [50, 56) [50, 56) [56, 62)
## [1201] [56, 62) [56, 62) [56, 62) [56, 62) [56, 62) [56, 62) [50, 56) [50, 56)
## [1209] [56, 62) [50, 56) [56, 62) [50, 56) [56, 62) [56, 62) [56, 62) [56, 62)
## [1217] [56, 62) [50, 56) [56, 62) [50, 56) [56, 62) [50, 56) [50, 56) [56, 62)
## [1225] [50, 56) [50, 56) [56, 62) [56, 62) [50, 56) [50, 56) [56, 62) [56, 62)
## [1233] [56, 62) [50, 56) [50, 56) [50, 56) [56, 62) [56, 62) [50, 56) [56, 62)
## [1241] [56, 62) [50, 56) [50, 56) [56, 62) [50, 56) [50, 56) [50, 56) [50, 56)
## [1249] [50, 56) [50, 56) [50, 56) [56, 62) [56, 62) [50, 56) [56, 62) [50, 56)
## [1257] [50, 56) [56, 62) [56, 62) [56, 62) [50, 56) [56, 62) [56, 62) [56, 62)
## [1265] [56, 62) [56, 62) [50, 56) [50, 56) [50, 56) [50, 56) [56, 62) [50, 56)
## [1273] [50, 56) [50, 56) [56, 62) [50, 56) [56, 62) [56, 62) [56, 62) [50, 56)
## [1281] [56, 62) [50, 56) [50, 56) [50, 56) [56, 62) [56, 62) [56, 62) [56, 62)
## [1289] [50, 56) [50, 56) [50, 56) [50, 56) [50, 56) [56, 62) [50, 56) [56, 62)
## [1297] [56, 62) [50, 56) [56, 62) [50, 56) [56, 62) [56, 62) [50, 56) [50, 56)
## [1305] [50, 56) [56, 62) [50, 56) [56, 62) [50, 56) [56, 62) [50, 56) [50, 56)
## [1313] [56, 62) [50, 56) [56, 62) [56, 62) [56, 62) [56, 62) [50, 56) [50, 56)
## [1321] [56, 62) [56, 62) [56, 62) [50, 56) [56, 62) [50, 56) [56, 62) [56, 62)
## [1329] [50, 56) [50, 56) [56, 62) [56, 62) [56, 62) [50, 56) [56, 62) [56, 62)
## [1337] [56, 62) [56, 62) [56, 62) [56, 62) [56, 62) [56, 62) [56, 62) [56, 62)
## [1345] [56, 62) [56, 62) [56, 62) [56, 62) [56, 62) [50, 56) [50, 56) [56, 62)
## [1353] [56, 62) [56, 62) [56, 62) [56, 62) [56, 62) [56, 62) [56, 62) [56, 62)
## [1361] [56, 62) [56, 62) [56, 62) [56, 62) [56, 62) [56, 62) [56, 62) [56, 62)
## [1369] [56, 62) [56, 62) [56, 62) [56, 62) [56, 62) [56, 62) [56, 62) [56, 62)
## [1377] [56, 62) [56, 62) [56, 62) [56, 62) [56, 62) [56, 62) [56, 62) [56, 62)
## [1385] [56, 62) [56, 62) [56, 62) [56, 62) [56, 62) [56, 62) [56, 62) [56, 62)
## [1393] [56, 62) [56, 62) [56, 62) [56, 62) [56, 62) [56, 62) [56, 62) [56, 62)
## [1401] [56, 62) [56, 62) [56, 62) [56, 62) [56, 62)
## Levels: [45, 50) [50, 56) [56, 62)
min(Framingham$AGE)
## [1] 45
max(Framingham$AGE)
## [1] 62

6.6 Guardar información

Ya sea guardar archivo original, una parte, o exportar la información:

save(Datos1, file="Datos2.RData")
FraminghamMen <- subset(Framingham, SEX=="Hombre")
save(FraminghamMen, file="Hombres_Framingham.RData")
write.csv(FraminghamMen, "Hombres_Framingham.csv", row.names=F, quote=F)

6.7 Eliminar datos

6.7.1 Eliminar filas que no cumplan una condición

nrow(Datos1) # antes
## [1] 15
Datos1 <- subset(Datos1, Peso<80 & Sexo=='Mujer') # mujeres de menos de 80kg
nrow(Datos1) # después
## [1] 7

6.7.2 Eliminar filas específicas

Datos1 <- Datos1[-c(2,4,5,8:12), ] # - = complement

6.7.3 Eliminar columnas

Datos1$Altura <- NULL # by var name
Datos1 <- Datos1[,-c(3)] # by index

6.8 Cambios de tipo de variable

Las variables de tipo factor están pensadas para almacenar datos cualitativos. Si en nuestros datos, las variables de tipo chr son cualitativas, podemos convertirlas facilmente a varaiables de tipo factor:

fctr.cols <- sapply(Framingham, is.factor) # get vars of type vector
Framingham[, fctr.cols] <- sapply(Framingham[, fctr.cols], as.character)

Limpiamos el cuarto:

## null device 
##           1