EXAMEN PARCIAL 2
pinta l
2024-01-08
INTRODUCCION
En el presente estudio, se ha llevado a cabo una exhaustiva encuesta que abarca diversas dimensiones de la vida cotidiana de los habitantes de la provincia de Chimborazo. Con el objetivo de comprender de manera integral la realidad sociodemográfica de la región, se han recopilado datos sobre aspectos cruciales que van desde la ubicación geográfica hasta los hábitos alimenticios y la infraestructura doméstica. Estas variables proporcionan una panorámica detallada de la diversidad de situaciones que caracterizan a los hogares encuestados, permitiendo así analizar patrones, identificar necesidades y entender mejor la calidad de vida de la población. A través de la exploración de factores como la educación, la salud, la ocupación y la infraestructura, este estudio se propone arrojar luz sobre las condiciones de vida de los habitantes de Chimborazo, contribuyendo de esta manera a la formulación de políticas y estrategias que respondan de manera efectiva a las necesidades y desafíos específicos de esta comunidad.
DESCRIPCIÓN DE LAS VARIABLES
UBICACIÓN: Relacionada con la ubicación geográfica de los encuestados.
Cantón de la provincia de Chimborazo: Indica en qué cantón de la provincia de Chimborazo vive el encuestado.
Edad del encuestado: Representa la edad de la persona que está siendo encuestada.
Sexo: Indica el género del encuestado.
Estado civil: Refleja el estado civil actual del encuestado, como soltero, casado, divorciado, viudo, etc.
Cantidad de personas en el hogar: Representa la cantidad de personas que conforman el hogar del encuestado.
Parentesco con el representante del hogar: Describe la relación de parentesco que el encuestado tiene con el representante del hogar.
Tipo de vivienda: Indica el tipo de vivienda en la que reside el encuestado, como casa, departamento, etc.
Material de las paredes exteriores: Describe el material predominante de las paredes exteriores de la vivienda.
Material del piso de la vivienda: Indica el material predominante del piso de la vivienda.
Número de cuartos de baño con ducha: Representa la cantidad de cuartos de baño con ducha de uso exclusivo en el hogar.
Tipo de servicio higiénico: Describe el tipo de servicio higiénico con el que cuenta el hogar.
Varios servicios y dispositivos electrónicos: Preguntas sobre la disponibilidad de servicios como internet, computadoras, teléfonos, electrodomésticos y otros dispositivos en el hogar.
Nivel de instrucción del Jefe del hogar: Indica el nivel educativo alcanzado por el jefe del hogar.
Afiliación a seguros de salud: Pregunta si alguien en el hogar está afiliado a seguros de salud, tanto públicos como privados.
Seguro de salud privada: Indaga sobre la posesión de seguros de salud privados con diferentes coberturas.
Ocupación del Representante del hogar: Describe la ocupación del representante del hogar.
NSE_2023_CUALI: nivel de estudio del individuo.
Sector de trabajo: Indica el sector en el cual trabaja el encuestado.
Hábitos alimenticios: Preguntas relacionadas con el consumo de alimentos como arroz, queso industrial, huevos y frutas.
DESAROLLO
CARGA DE BASE DE DATOS:
library(readxl)
Hogar <- read_excel("DATA_PARCIAL2.xlsx")
# View(Hogar)NOMBRE DE LAS VARIABLES DE LA BASE DE DATOS
# nombres de las variables
names(Hogar)## [1] "UBICACIÓN"
## [2] "¿En qué cantón de la provincia de Chimborazo vive?"
## [3] "Edad del encuestado"
## [4] "Sexo"
## [5] "Estado civil"
## [6] "Cantidad de personas que conforman el hogar"
## [7] "¿Qué parentesco tiene con el representante del hogar?"
## [8] "¿Cuál es el tipo de vivienda?"
## [9] "El material predominante de las paredes exteriores de la vivienda es de:"
## [10] "EL material predominante del piso de la vivienda es de:"
## [11] "¿Cuántos cuartos de baño con ducha de uso exclusivo tiene este hogar?"
## [12] "Tipo de servicio higiénico con que cuenta este hogar"
## [13] "¿Cuenta con servicio de internet?"
## [14] "¿Cuenta con computadora de escritorio?"
## [15] "¿Cuenta con computadora portátil?"
## [16] "¿Cuántos celulares activados tiene en este hogar?"
## [17] "¿Tiene su hogar servicio de teléfono convencional?"
## [18] "¿Tiene en su hogar cocina con horno?"
## [19] "¿Tiene refrigeradora?"
## [20] "¿Tiene lavadora?"
## [21] "¿Tiene equipo de sonido?"
## [22] "¿Cuántas TV a color tienen en este hogar?"
## [23] "¿Cuántos vehículos de uso exclusivo tiene este hogar?"
## [24] "¿Alguien en el hogar compra vestimenta en centros comerciales?"
## [25] "¿En el hogar alguien ha usado internet en los últimos 6 meses?"
## [26] "¿En el hogar alguien utiliza correo electrónico que no es del trabajo?"
## [27] "¿En el hogar alguien está registrado en una red social?"
## [28] "Exceptuando los libros de texto o manuales de estudio y lecturas de trabajo¿Alguien del hogar ha leído algún libro completo en los últimos 3 meses?"
## [29] "¿Cuál es el nivel de instrucción del Jefe del hogar?"
## [30] "¿Alguien en el hogar está afiliado o cubierto por el seguro del IESS (general, voluntario o campesino) y/o seguro del ISSFA o ISSPOL?"
## [31] "¿Alguien en el hogar tiene seguro de salud privada con hospitalización, seguro de salud privada sin hospitalización, seguro internacional,seguros municipales y de Consejos Provinciales y/o seguro de vida?"
## [32] "¿Cuál es la ocupación del Representante del hogar?"
## [33] "NSE_2023_CUALI"
## [34] "¿Para qué sector trabaja?"
## [35] "¿Consume usted arroz?"
## [36] "¿Consume usted queso industrial (CON MARCA)"
## [37] "¿Consume usted huevos?"
## [38] "¿Consume frutas?"TIPO DE CLASE DE LAS VARIABLES DE LA BASE DE DATOS
# Comprobar el tipo de clase de cada variable
clases <- sapply(Hogar, class)
# clasesCAMBIAR EL TIPO DE CLASE
Variable “UBICACIÓN”
Hogar$UBICACIÓN <- factor(Hogar$UBICACIÓN)
levels(Hogar$UBICACIÓN) <- c("Velasco", "Yaruquíes", "Colta", "Guano", "Guamote", "Veloz", "Maldonado", "Cacha", "Lizarzaburu", "San Luis", "Penipe", "San Juan", "Pallatanga", "Punín", "Pungalá", "Chambo", "Alausí", "Licán", "Cumanda", "Licto", "Calpi", "Chunchi", "Quimiag", "Cubijíes", "Flores")
class(Hogar$UBICACIÓN)## [1] "factor"Variable “¿En qué cantón de la provincia de Chimborazo vive?”
Hogar$`¿En qué cantón de la provincia de Chimborazo vive?` <- factor(Hogar$`¿En qué cantón de la provincia de Chimborazo vive?`)
levels(Hogar$`¿En qué cantón de la provincia de Chimborazo vive?`) <- c("Alausí","Chambo","Chunchi","Colta","Cumanda","Guamote","Guano","Pallatanga","Penipe","Riobamba")
class(Hogar$`¿En qué cantón de la provincia de Chimborazo vive?`)## [1] "factor"Variable “Edad del encuestado”
Numérico
Variable “Sexo”
Hogar$Sexo <- factor(Hogar$Sexo )
levels(Hogar$Sexo ) <- c("Masculino","Femenino")
class(Hogar$Sexo)## [1] "factor"Variable “Estado civil”
Hogar$`Estado civil` <- factor(Hogar$`Estado civil`)
levels(Hogar$`Estado civil`) <- c("Casado","Divorciado","Otro","Soltero","Viudo")
class(Hogar$`Estado civil`)## [1] "factor"Variable “Cantidad de personas que conforman el hogar”
Numérico
Variable “¿Qué parentesco tiene con el representante del hogar?”
Hogar$`¿Qué parentesco tiene con el representante del hogar?` <- factor(Hogar$`¿Qué parentesco tiene con el representante del hogar?`)
levels(Hogar$`¿Qué parentesco tiene con el representante del hogar?`) <- c("Cónyugue o conviviente","Otros parientes","Otros parientes","Yo soy el representante del Hogar","Hijo o hija","Padres o suegros","Otros no parientes","Yerno o nuera")
class(Hogar$`¿Qué parentesco tiene con el representante del hogar?`)## [1] "factor"Variable “¿Cuál es el tipo de vivienda?”
Hogar$`¿Cuál es el tipo de vivienda?` <- factor(Hogar$`¿Cuál es el tipo de vivienda?`)
levels(Hogar$`¿Cuál es el tipo de vivienda?`) <- c("Casa/Villa","Departamento en casa o edificio","Suite de lujo","Cuarto(s) en casa de inquilinato ","Choza/Covacha/Otro","Rancho"," Mediagua")
class(Hogar$`¿Cuál es el tipo de vivienda?`)## [1] "factor"Variable “El material predominante de las paredes exteriores de la vivienda es de:”
Hogar$`El material predominante de las paredes exteriores de la vivienda es de:` <- factor(Hogar$`El material predominante de las paredes exteriores de la vivienda es de:`)
levels(Hogar$`El material predominante de las paredes exteriores de la vivienda es de:`) <- c("Adobe/Tapia","Caña no revestida/ Otros materiales","Ladrillo o bloque","Hormigón"," Caña revestida o bahareque/Madera")
class(Hogar$`El material predominante de las paredes exteriores de la vivienda es de:`)## [1] "factor"Variable “EL material predominante del piso de la vivienda es de:”
Hogar$`EL material predominante del piso de la vivienda es de:` <- factor(Hogar$`EL material predominante del piso de la vivienda es de:`)
levels(Hogar$`EL material predominante del piso de la vivienda es de:`) <- c("Duela, parquet, tablón o piso flotante"," Cerámica, baldosa, vinil o marmetón","Tabla sin tratar","Tierra/ Caña/ Otros materiales","Ladrillo o cemento")
class(Hogar$`El material predominante de las paredes exteriores de la vivienda es de:`)## [1] "factor"Variable “¿Cuántos cuartos de baño con ducha de uso exclusivo tiene este hogar?”
Hogar$`¿Cuántos cuartos de baño con ducha de uso exclusivo tiene este hogar?` <- factor(Hogar$`¿Cuántos cuartos de baño con ducha de uso exclusivo tiene este hogar?`)
levels(Hogar$`¿Cuántos cuartos de baño con ducha de uso exclusivo tiene este hogar?`) <- c("No tiene cuarto de baño exclusivo con ducha en el hogar","Tiene 2 cuartos de baño exclusivos con ducha","Tiene 1 cuarto de baño exclusivo con ducha ","Tiene 3 o más cuartos de baño exclusivos con ducha")
class(Hogar$`¿Cuántos cuartos de baño con ducha de uso exclusivo tiene este hogar?`)## [1] "factor"Variable “Tipo de servicio higiénico con que cuenta este hogar”
Hogar$`Tipo de servicio higiénico con que cuenta este hogar` <- factor(Hogar$`Tipo de servicio higiénico con que cuenta este hogar`)
levels(Hogar$`Tipo de servicio higiénico con que cuenta este hogar`) <- c("Con descarga directa al mar, río, lago o quebrada ","Conectado a pozo séptico","Tiene 1 cuarto de baño exclusivo con ducha Letrina","Conectado a pozo ciego","Conectado a red pública de alcantarillado","Sin Marca")
class(Hogar$`Tipo de servicio higiénico con que cuenta este hogar`)## [1] "factor"Variable “¿Cuenta con servicio de internet?”
Hogar$`¿Cuenta con servicio de internet?` <- factor(Hogar$`¿Cuenta con servicio de internet?`)
levels(Hogar$`¿Cuenta con servicio de internet?`) <- c("No","Si")
class(Hogar$`¿Cuenta con servicio de internet?`)## [1] "factor"Variable “¿Cuenta con computadora de escritorio?”
Hogar$`¿Cuenta con computadora de escritorio?` <- factor(Hogar$`¿Cuenta con computadora de escritorio?`)
levels(Hogar$`¿Cuenta con computadora de escritorio?`) <- c("No","Si")
class(Hogar$`¿Cuenta con computadora de escritorio?`)## [1] "factor"Variable “¿Cuenta con computadora portátil?”
Hogar$`¿Cuenta con computadora portátil?` <- factor(Hogar$`¿Cuenta con computadora portátil?`)
levels(Hogar$`¿Cuenta con computadora portátil?`) <- c("No","Si")
class(Hogar$`¿Cuenta con computadora portátil?`)## [1] "factor"Variable “¿Cuántos celulares activados tiene en este hogar?”
Hogar$`¿Cuántos celulares activados tiene en este hogar?` <- factor(Hogar$`¿Cuántos celulares activados tiene en este hogar?` )
levels(Hogar$`¿Cuántos celulares activados tiene en este hogar?` ) <- c("No tiene celular nadie en el hogar","Tiene 3 celulares","Tiene 1 celular","Tiene 4 ó más celulares","Tiene 2 celulares")
class(Hogar$`¿Cuántos celulares activados tiene en este hogar?` )## [1] "factor"Variable “¿Tiene su hogar servicio de teléfono convencional?”
Hogar$`¿Tiene su hogar servicio de teléfono convencional?` <- factor(Hogar$`¿Tiene su hogar servicio de teléfono convencional?` )
levels(Hogar$`¿Tiene su hogar servicio de teléfono convencional?` ) <- c("No", "Si")
class(Hogar$`¿Tiene su hogar servicio de teléfono convencional?` )## [1] "factor"Variable “¿Tiene en su hogar cocina con horno?”
Hogar$`¿Tiene en su hogar cocina con horno?` <- factor(Hogar$`¿Tiene en su hogar cocina con horno?` )
levels(Hogar$`¿Tiene en su hogar cocina con horno?`) <- c("No", "Si")
class(Hogar$`¿Tiene en su hogar cocina con horno?`)## [1] "factor"Variable “¿Tiene refrigeradora?”
Hogar$`¿Tiene refrigeradora?` <- factor(Hogar$`¿Tiene refrigeradora?`)
levels(Hogar$`¿Tiene refrigeradora?`) <- c("No", "Si")
class(Hogar$`¿Tiene refrigeradora?`)## [1] "factor"Variable “¿Tiene lavadora?”
Hogar$`¿Tiene lavadora?` <- factor(Hogar$`¿Tiene lavadora?` )
levels(Hogar$`¿Tiene lavadora?`) <- c("No", "Si")
class(Hogar$`¿Tiene lavadora?`)## [1] "factor"Variable “¿Tiene equipo de sonido?”
Hogar$`¿Tiene equipo de sonido?` <- factor(Hogar$`¿Tiene equipo de sonido?` )
levels(Hogar$`¿Tiene equipo de sonido?`) <- c("No", "Si")
class(Hogar$`¿Tiene equipo de sonido?`)## [1] "factor"Variable “¿Cuántas TV a color tienen en este hogar?”
Hogar$`¿Cuántas TV a color tienen en este hogar?` <- factor(Hogar$`¿Cuántas TV a color tienen en este hogar?` )
levels(Hogar$`¿Cuántas TV a color tienen en este hogar?`) <- c("No tiene TV a color en el hogar","Tiene 3 ó más TV a color ","Tiene 1 TV a color ","Tiene 2 TV a color ")
class(Hogar$`¿Cuántas TV a color tienen en este hogar?`)## [1] "factor"Variable “¿Cuántos vehículos de uso exclusivo tiene este hogar?”
Hogar$`¿Cuántos vehículos de uso exclusivo tiene este hogar?` <- factor(Hogar$`¿Cuántos vehículos de uso exclusivo tiene este hogar?` )
levels(Hogar$`¿Cuántos vehículos de uso exclusivo tiene este hogar?`) <- c("No tiene vehículo exclusivo para el hogar","Tiene 3 ó más vehículo exclusivo","Tiene 2 vehículo exclusivo","Tiene un vehículo exclusivo ")
class(Hogar$`¿Cuántos vehículos de uso exclusivo tiene este hogar?`)## [1] "factor"Variable “¿Alguien en el hogar compra vestimenta en centros comerciales?”
Hogar$`¿Alguien en el hogar compra vestimenta en centros comerciales?` <- factor(Hogar$`¿Alguien en el hogar compra vestimenta en centros comerciales?`)
levels(Hogar$`¿Alguien en el hogar compra vestimenta en centros comerciales?`) <- c("No", "Si")
class(Hogar$`¿Alguien en el hogar compra vestimenta en centros comerciales?`)## [1] "factor"Variable “¿En el hogar alguien ha usado internet en los últimos 6 meses?”
Hogar$`¿En el hogar alguien ha usado internet en los últimos 6 meses?` <- factor(Hogar$`¿En el hogar alguien ha usado internet en los últimos 6 meses?`)
levels(Hogar$`¿En el hogar alguien ha usado internet en los últimos 6 meses?`) <- c("No", "Si")
class(Hogar$`¿En el hogar alguien ha usado internet en los últimos 6 meses?`)## [1] "factor"Variable “¿En el hogar alguien utiliza correo electrónico que no es del trabajo?”
Hogar$`¿En el hogar alguien utiliza correo electrónico que no es del trabajo?` <- factor(Hogar$`¿En el hogar alguien utiliza correo electrónico que no es del trabajo?`)
levels(Hogar$`¿En el hogar alguien utiliza correo electrónico que no es del trabajo?`) <- c("No", "Si")
class(Hogar$`¿En el hogar alguien utiliza correo electrónico que no es del trabajo?`)## [1] "factor"Variable “¿En el hogar alguien está registrado en una red social?”
Hogar$`¿En el hogar alguien está registrado en una red social?` <- factor(Hogar$`¿En el hogar alguien está registrado en una red social?`)
levels(Hogar$`¿En el hogar alguien está registrado en una red social?`) <- c("No", "Si")
class(Hogar$`¿En el hogar alguien está registrado en una red social?`)## [1] "factor"Variable “Exceptuando los libros de texto o manuales de estudio y lecturas de trabajo¿Alguien del hogar ha leído algún libro completo en los últimos 3 meses?”
Hogar$`Exceptuando los libros de texto o manuales de estudio y lecturas de trabajo¿Alguien del hogar ha leído algún libro completo en los últimos 3 meses?` <- factor(Hogar$`Exceptuando los libros de texto o manuales de estudio y lecturas de trabajo¿Alguien del hogar ha leído algún libro completo en los últimos 3 meses?`)
levels(Hogar$`Exceptuando los libros de texto o manuales de estudio y lecturas de trabajo¿Alguien del hogar ha leído algún libro completo en los últimos 3 meses?`) <- c("No", "Si")
class(Hogar$`Exceptuando los libros de texto o manuales de estudio y lecturas de trabajo¿Alguien del hogar ha leído algún libro completo en los últimos 3 meses?`)## [1] "factor"Variable “¿Cuál es el nivel de instrucción del Jefe del hogar?”
Hogar$`¿Cuál es el nivel de instrucción del Jefe del hogar?` <- factor(Hogar$`¿Cuál es el nivel de instrucción del Jefe del hogar?` )
levels(Hogar$`¿Cuál es el nivel de instrucción del Jefe del hogar?` ) <- c("4 ó más años de educación superior (completa sin postgrado)","Hasta 3 años de educación superior (incompleta)","Postgrado","Primaria completa ","Primaria incompleta","Secundaria completa ","Secundaria incompleta","Sin estudios ")
class(Hogar$`¿Cuál es el nivel de instrucción del Jefe del hogar?` )## [1] "factor"Variable “¿Alguien en el hogar está afiliado o cubierto por el seguro del IESS (general, voluntario o campesino) y/o seguro del ISSFA o ISSPOL?”
Hogar$`¿Alguien en el hogar está afiliado o cubierto por el seguro del IESS (general, voluntario o campesino) y/o seguro del ISSFA o ISSPOL?` <- factor(Hogar$`¿Alguien en el hogar está afiliado o cubierto por el seguro del IESS (general, voluntario o campesino) y/o seguro del ISSFA o ISSPOL?`)
levels(Hogar$`¿Alguien en el hogar está afiliado o cubierto por el seguro del IESS (general, voluntario o campesino) y/o seguro del ISSFA o ISSPOL?`) <- c("No", "Si")
class(Hogar$`¿Alguien en el hogar está afiliado o cubierto por el seguro del IESS (general, voluntario o campesino) y/o seguro del ISSFA o ISSPOL?`)## [1] "factor"Variable “¿Alguien en el hogar tiene seguro de salud privada con hospitalización, seguro de salud privada sin hospitalización, seguro internacional,seguros municipales y de Consejos Provinciales y/o seguro de vida?”
Hogar$`¿Alguien en el hogar tiene seguro de salud privada con hospitalización, seguro de salud privada sin hospitalización, seguro internacional,seguros municipales y de Consejos Provinciales y/o seguro de vida?` <- factor(Hogar$`¿Alguien en el hogar tiene seguro de salud privada con hospitalización, seguro de salud privada sin hospitalización, seguro internacional,seguros municipales y de Consejos Provinciales y/o seguro de vida?`)
levels(Hogar$`¿Alguien en el hogar tiene seguro de salud privada con hospitalización, seguro de salud privada sin hospitalización, seguro internacional,seguros municipales y de Consejos Provinciales y/o seguro de vida?`) <- c("Si", "No")
class(Hogar$`¿Alguien en el hogar tiene seguro de salud privada con hospitalización, seguro de salud privada sin hospitalización, seguro internacional,seguros municipales y de Consejos Provinciales y/o seguro de vida?`)## [1] "factor"Variable “¿Cuál es la ocupación del Representante del hogar?”
Hogar$`¿Cuál es la ocupación del Representante del hogar?` <- factor(Hogar$`¿Cuál es la ocupación del Representante del hogar?`)
levels(Hogar$`¿Cuál es la ocupación del Representante del hogar?`) <- c("Desempleado "," Empleados de oficina","Fuerzas Armadas","Inactivos "," Oficiales operarios y artesanos"," Operadores de instalaciones y máquinas","Personal directivo de la Administración Pública y de empresas"," Profesionales científicos e intelectuales","Técnicos y profesionales de nivel medio","Trabajador calificados agropecuarios y pesqueros","Trabajador de los servicios y comerciantes","Trabajadores no calificados")
class(Hogar$`¿Cuál es la ocupación del Representante del hogar?`)## [1] "factor"Variable “NSE_2023_CUALI”
Hogar$NSE_2023_CUALI <- factor(Hogar$NSE_2023_CUALI)
levels(Hogar$NSE_2023_CUALI) <- c("Alto", "Bajo", "Medio Alto","Medio Bajo","Medio Típico")
class(Hogar$NSE_2023_CUALI)## [1] "factor"Variable “¿Para qué sector trabaja?”
Hogar$`¿Para qué sector trabaja?`<- factor(Hogar$`¿Para qué sector trabaja?`)
levels(Hogar$`¿Para qué sector trabaja?`) <- c("Sector público", "Sector privado", "Ninguno")
class(Hogar$`¿Para qué sector trabaja?`)## [1] "factor"Variable “¿Consume usted arroz?”
Hogar$`¿Consume usted arroz?`<- factor(Hogar$`¿Consume usted arroz?`)
levels(Hogar$`¿Consume usted arroz?`) <- c("No", "Si")
class(Hogar$`¿Consume usted arroz?`)## [1] "factor"Variable “¿Consume usted queso industrial (CON MARCA)”
Hogar$`¿Consume usted queso industrial (CON MARCA)` <- factor(Hogar$`¿Consume usted queso industrial (CON MARCA)`)
levels(Hogar$`¿Consume usted queso industrial (CON MARCA)`) <- c("No", "Si")
class(Hogar$`¿Consume usted queso industrial (CON MARCA)`)## [1] "factor"Variable “¿Consume usted huevos?”
Hogar$`¿Consume usted huevos?` <- factor(Hogar$`¿Consume usted huevos?`)
levels(Hogar$`¿Consume usted huevos?`) <- c("No", "Si")
class(Hogar$`¿Consume usted huevos?`)## [1] "factor"Variable “¿Consume frutas?”
Hogar$`¿Consume frutas?` <- factor(Hogar$`¿Consume frutas?`)
levels(Hogar$`¿Consume frutas?`) <- c("No", "Si")
class(Hogar$`¿Consume frutas?`)## [1] "factor"LIMPIEZA DE LA BASE DE DATOS
# Eliminar filas con valores faltantes
Hogar_clean <- na.omit(Hogar)
# Verificar el resultado
# head(Hogar_clean)ELIMINAR LOS NA
# Eliminar filas con valores NA
Hogar_clean <- na.omit(Hogar)
# Verificar el resultado
# head(Hogar_clean)ANÁLISIS DESCRIPTIVO
Descripción teórica:
Concepto y Fórmula:
Interpretación:
Hogar1 <- Hogar_clean[, c("Edad del encuestado", "Cantidad de personas que conforman el hogar")]
summary(Hogar1)## Edad del encuestado Cantidad de personas que conforman el hogar
## Min. :15.00 Min. : 1.000
## 1st Qu.:26.00 1st Qu.: 3.000
## Median :37.00 Median : 4.000
## Mean :38.16 Mean : 4.319
## 3rd Qu.:47.00 3rd Qu.: 5.000
## Max. :81.00 Max. :37.000SELECCIÓN DE LAS VARIABLES FACTOR
Hogar1 <- Hogar_clean[, c("Edad del encuestado", "Cantidad de personas que conforman el hogar")]
# Verificar el resultado
head(Hogar1)## # A tibble: 6 × 2
## `Edad del encuestado` `Cantidad de personas que conforman el hogar`
## <dbl> <dbl>
## 1 58 6
## 2 28 3
## 3 40 4
## 4 15 4
## 5 23 6
## 6 48 5MÍNIMO
El valor mínimo se refiere al valor más pequeño encontrado de las variables dentro de esa base de datos. Para calcular el valor mínimo de una variable se, puedes utilizar la función min(). Esta función devuelve el valor mínimo dentro de un vector o columna.
\[\text{Valor Mínimo} = \min(x_1, x_2, ..., x_n)\]
# Obtener el valor mínimo de las variables
apply(Hogar1[, c("Edad del encuestado", "Cantidad de personas que conforman el hogar")], 2, min)## Edad del encuestado
## 15
## Cantidad de personas que conforman el hogar
## 1Interpretación:
Min. (Mínimo): El valor más pequeño en la muestra es 15 de las edades del encuestado.
Min. (Mínimo): El valor más pequeño en la muestra es 1 de la cantidad de personas que conforman el hogar .
PRIMER CUARTIL
El primer cuartil, también conocido como Q1 o percentil 25, es un valor estadístico que indica el punto en el cual el 25% de los datos de una muestra se encuentran por debajo de ese valor. En otras palabras, el primer cuartil divide la muestra en cuatro partes iguales, donde el 25% de los datos se encuentran por debajo del primer cuartil y el 75% restante se encuentra por encima.
\[Q1=Percentil . 25\]
# Calcular el primer cuartil de las variables
first_quartiles <- quantile(Hogar1$`Edad del encuestado`, probs = 0.25)
first_quartiles_s <- quantile(Hogar1$`Cantidad de personas que conforman el hogar`, probs = 0.25)
# Verificar los resultados
C <- data.frame(
first_quartiles,
first_quartiles_s
)
print(C)## first_quartiles first_quartiles_s
## 25% 26 3Interpretación:
1st Qu. (Primer cuartil): El 25% de los datos tienen valores menores o iguales a 26 de la edad del encuestado.
1st Qu. (Primer cuartil): El 25% de los datos tienen valores menores o iguales a 3 de la cantidad de personas que conforman el hogar.
MEDIANA
La mediana es una medida de tendencia central que indica el valor que se encuentra en el centro de un conjunto de datos ordenados. Es decir, la mitad de los valores están por encima de la mediana y la otra mitad están por debajo.
\[\text{Mediana} = \text{valor\medio}\]
apply(Hogar1, 2, median)## Edad del encuestado
## 37
## Cantidad de personas que conforman el hogar
## 4Interpretación:
Median (Mediana): La mediana, que divide los datos en dos mitades, es 37. Esto significa que el 50% de los datos tienen valores menores o iguales a 37 de la edad del encuestado.
Median (Mediana): La mediana es 4, lo que significa que el 50% de los datos tienen valores menores o iguales a 4 de la cantidad de personas que conforman el hogar.
MEDIA
El valor de la media se refiere al promedio de los valores de una variable específica. La media es una medida de tendencia central que se calcula sumando todos los valores y dividiendo el resultado entre la cantidad total de valores.
\[\bar{x} = \frac{1}{N} \sum_{i=1}^{N} x_i\]
data.frame(mean(Hogar1$`Edad del encuestado`),
mean(Hogar1$`Cantidad de personas que conforman el hogar`))## mean.Hogar1..Edad.del.encuestado..
## 1 38.16449
## mean.Hogar1..Cantidad.de.personas.que.conforman.el.hogar..
## 1 4.318538Interpretación:
Mean (Media): La media o promedio es aproximadamente 38.16 de la edad del encuestado.
Mean (Media): La media es aproximadamente 4.319 de las personas que conforman el hogar.
TERCER CUARTIL
El valor del tercer cuartil, también conocido como el percentil 75 (Q3), es un valor que divide los datos ordenados en cuatro partes iguales. Esto significa que el 75% de los datos se encuentran por debajo del tercer cuartil y el 25% restante se encuentra por encima.
\[Q_3=Percentil. 75\]
Hogar11 <- Hogar_clean[, c("Edad del encuestado", "Cantidad de personas que conforman el hogar")]
# Calcular el tercer cuartil para cada variable
q3_values <- apply(Hogar11, 2, quantile, probs = 0.75)
# Imprimir los resultados
print(q3_values)## Edad del encuestado
## 47
## Cantidad de personas que conforman el hogar
## 5Interpretación:
3rd Qu. (Tercer cuartil): El 75% de los datos tienen valores menores o iguales a 47 de la edad del encuestado.
3rd Qu. (Tercer cuartil): El 75% de los datos tienen valores menores o iguales a 5 de la cantidad de personas que conforman el hogar.
VALOR MÁXIMO
El valor máximo se refiere al valor más grande que se encuentra en una determinada variable o columna. Es el valor más alto dentro de esa columna específica de un conjunto de datos.
\[\text{Valor Máximo} = \text{max}(x)\]
Hogar12 <- Hogar_clean[, c("Edad del encuestado", "Cantidad de personas que conforman el hogar")]
# Calcular el valor máximo de cada variable
maximos <- apply(Hogar12, 2, max)
# Imprimir los valores máximos
print(maximos)## Edad del encuestado
## 81
## Cantidad de personas que conforman el hogar
## 37Interpretación:
Max. (Máximo): El valor más grande en la muestra es 81 de la edad del encuestado.
Max. (Máximo): El valor más grande en la muestra es 37 de la cantidad de personas que conforman el hogar.
DESVIACIÓN ESTANDAR
La desviación estándar es una medida de dispersión o variabilidad en un conjunto de datos. Indica cuánto se alejan los valores individuales de la media del conjunto de datos.
En términos más simples, la desviación estándar nos dice qué tan dispersos están los valores alrededor de la media. Si la desviación estándar es baja, significa que los valores tienden a estar cerca de la media. Por otro lado, si la desviación estándar es alta, indica que los valores están más dispersos y alejados de la media.
La desviación estándar se calcula tomando la raíz cuadrada de la varianza. Es una medida ampliamente utilizada en estadística para comprender la dispersión de los datos y comparar la variabilidad entre diferentes conjuntos de datos.
\[\sigma = \sqrt{\frac{\sum_{i=1}^{N}(x_i - \bar{x})^2}{N-1}}\]
Hogar13 <- Hogar_clean[, c("Edad del encuestado", "Cantidad de personas que conforman el hogar")]
desviacion_estandar <- apply(Hogar13, 2, sd)
desviacion_estandar## Edad del encuestado
## 13.983034
## Cantidad de personas que conforman el hogar
## 2.343052Interpretación:
La desviación estándar de 13.98 indica la dispersión promedio de las edades de los encuestados respecto a la media. Cuanto mayor sea la desviación estándar, mayor será la dispersión de las edades en el conjunto de datos.
La desviación estándar de 2.34 indica la dispersión promedio en la cantidad de personas que conforman los hogares. Una desviación estándar más baja indica que los tamaños de los hogares tienden a estar más cerca de la media, mientras que una desviación estándar más alta indica una mayor variabilidad en los tamaños de los hogares.
GRÁFICA DE BARRAS
Relación de si consume huevos
library(ggplot2)
grafico_H = ggplot(Hogar_clean,aes(x=`¿Consume usted huevos?`))+
geom_bar(stat="count", width=0.7, fill="blue")+
theme_minimal()
print(grafico_H)
Interpretación:
Se puede observar que en la base de datos de las personas encuestadas si consumen huevos por lo que se puede decir que les gusta o siempre se les apetece comer teniendo la variable (SI) muy alto.
GRÁFICO DE DISPERSIÓN
Edad del encuestado, cantidad de personas que conforman el hogar y consumo de huevos
library(ggplot2)
ggplot(Hogar_clean, aes(x = `Cantidad de personas que conforman el hogar`, y = `Edad del encuestado`, col = factor(`¿Consume usted huevos?`))) +
geom_point() +
labs(title = "Relación entre edad y cantidad de personas que conforman el hogar y consumo de huevos",
x = "cantidad de personas que conforman el hogar",
y = "edad del encuestado")
Interpretación
Los datos en su mayoria de los hogares o personas encuestadas estan en un rango de 0 a 10 y personas que conforman el hogar estan de 20 a 60 años edad del encuestado. El consumo de huevos es alto por lo cual se podría decir que las personas siles gusta o si consumen.
GRÁFICO DE BOXPLOT
Edad del encuestado y consumo de huevos
ggplot(Hogar_clean, aes(x=`¿Consume usted huevos?`, y=`Edad del encuestado`, color=`¿Consume usted huevos?`))+geom_boxplot()
Interpretación
En resumen, este gráfico de boxplot nos permite comparar las distribuciones de los valores del consumo de huevos y las edades, y se identifica datos atípicos posibles en la variable de consumo de huevo (SI), y la media de las cajas esta dispareja o elevada en (SI).
MATRIZ DE CORRELACIÓN DE LAS VARIABLES TIPO FACTOR
Una matriz de correlación se utiliza para investigar als dependencia entre multibles variables a la vez. El resultado de la tabla contiene los coeficientes de correlacion entre las diferentes variables que conforman el dataframe, nos ayuda a entender la relacion entre ellas. Con este grafico podemos seleccionar las variables relevantes y el eliminar aquellas que son redundantes, que tienen covariancia.
deb.cor <- cor(Hogar1, method = "pearson")
round(deb.cor,2)## Edad del encuestado
## Edad del encuestado 1.00
## Cantidad de personas que conforman el hogar 0.03
## Cantidad de personas que conforman el hogar
## Edad del encuestado 0.03
## Cantidad de personas que conforman el hogar 1.00GRÁFICO DE MATRIZ DE CORRELACIÓN
Como ver la correlacion graficamente para identificar rapidamente las variables que nos interesan:
library(corrplot)## corrplot 0.92 loadedcorrplot(deb.cor)
Interpretación
Correlación entre “Edad del encuestado” y “Edad del encuestado”:
El valor es 1.00, que es la correlación perfecta consigo misma (como es de esperar). Esto indica que no hay variabilidad en la relación entre la “Edad del encuestado” y sí misma.
Correlación entre “Edad del encuestado” y “Cantidad de personas que conforman el hogar”:
El valor es 0.03, que es un valor muy cercano a cero. Indica una correlación muy débil entre la “Edad del encuestado” y la “Cantidad de personas que conforman el hogar”. La relación lineal entre estas dos variables es prácticamente inexistente.
Correlación entre “Cantidad de personas que conforman el hogar” y “Cantidad de personas que conforman el hogar”:
El valor es 1.00, indicando una correlación perfecta consigo misma. Esto significa que no hay variabilidad en la relación entre la “Cantidad de personas que conforman el hogar” y sí misma.
En resumen, la matriz de correlación sugiere que no hay una relación lineal fuerte entre la “Edad del encuestado” y la “Cantidad de personas que conforman el hogar”, ya que el coeficiente de correlación es muy bajo (0.03).
GRÁFICO MULTIVARIANTE
library(ggplot2)
library(ggmosaic)
library(plotly)##
## Attaching package: 'plotly'## The following object is masked from 'package:ggplot2':
##
## last_plot## The following object is masked from 'package:stats':
##
## filter## The following object is masked from 'package:graphics':
##
## layoutlibrary(ggrepel)
library(GGally)## Registered S3 method overwritten by 'GGally':
## method from
## +.gg ggplot2##
## Attaching package: 'GGally'## The following object is masked from 'package:ggmosaic':
##
## happy# Imprimir la estructura del conjunto de datos
# str(Hogar_clean)Datos Númericos
# Seleccionar columnas numéricas
numeric_columns <- Hogar_clean[, sapply(Hogar_clean, is.numeric)]Gráfica
# Crear la matriz de gráficos
p <- ggpairs(numeric_columns) + theme_bw()
for (i in 1:p$nrow) {
for(j in 1:p$ncol){
p[i, j] <- p[i, j]+
scale_fill_manual(values=c("#00AFBB", "#E7B800")) +
scale_color_manual(values=c("#00AFBB", "#E7B800"))
}
}
p
# Crear la matriz de gráficos con paleta de colores
p <- ggpairs(numeric_columns, aes(color = "blue")) +
theme_bw() +
scale_color_manual(values = c("#00AFBB", "#E7B800", "#FC4E07"))
# Imprimir la matriz de gráficos
print(p)
clases <- sapply(Hogar_clean, class)
classes <- sapply(Hogar_clean, class)
for (i in 1:ncol(Hogar_clean))
if(classes[i]== "integer")
Hogar_clean[[i]] =as.numeric(Hogar_clean[[i]])
classes = sapply(Hogar_clean, class)
# classesVARIABLE DE CONSUMO DE HUEVOS
# Crear la matriz de gráficos con paleta de colores
w <- ggpairs(Hogar_clean[,c(which(classes =="numeric"),37)], aes(color=`¿Consume usted huevos?`)) + theme_bw()
# change color manually
# loop through each plot changing relevant scales
for (i in 1:w$nrow) {
for (j in 1:w$ncol) {
w[i,j] <- w[i,j]+
scale_fill_manual(values = c("green", "blue"))+
scale_color_manual(values=c("skyblue","red"))
}
}
w## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.
Interpretación
A partir de la gráfica de mosaico se muestra las correlaciones entre las variables del conjunto de datos Hogares, se puede concluir lo siguiente:
No hay una correlación significativa entre las edades del encuestado y la cantidad de personas que conforman el hogar. El valor de la correlación es 0.035.
La variable (NO) es de 0.026 lo que puede estar no relacionado de madera moderada el valor es muy bajo
La variable (SI) es de 0.168 lo que es un valor muy bajo no relacionado
En resumen, las variables, no parecen estar relacionadas de manera moderada .
VARIABLES FICTICIAS
Las variables ficticias, también conocidas como variables dummy, son variables binarias que se utilizan para representar categorías o niveles de una variable categórica en un modelo de regresión. Estas variables toman el valor de 1 si la observación pertenece a una categoría específica y 0 en caso contrario. Se utilizan para capturar el efecto de las diferentes categorías en el modelo de manera individualizada.
Dividiremos el conjunto completo de individuos en dos partes, un para entrenar el modelo, que contiene el 80% de los individuos y otra para validarlo que contiene el resto. Esto es así porque si valoramos el modelo con las misma observaciones que usamos para ajustarlo, la bondad del ajuste quedará sobrevalorada.
Antes de ajustar cualquier modelo, es conveniente escalar las variables numéricas y tenemos que reacondicionar las variables categóricas convirtiéndolas en variables ficticias y usando la primera o la última categoría como comparación.
Escalamiento de númericas
numeric_columns <- sapply(Hogar_clean, is.numeric)
Hogar_clean[, numeric_columns] <- scale(Hogar_clean[, numeric_columns])
head(Hogar_clean)## # A tibble: 6 × 38
## UBICACIÓN ¿En qué cantón de la pr…¹ `Edad del encuestado` Sexo `Estado civil`
## <fct> <fct> <dbl> <fct> <fct>
## 1 Quimiag Riobamba 1.42 Feme… Casado
## 2 Quimiag Riobamba -0.727 Feme… Casado
## 3 Flores Riobamba 0.131 Feme… Casado
## 4 Quimiag Riobamba -1.66 Feme… Soltero
## 5 Cubijíes Riobamba -1.08 Feme… Soltero
## 6 Cubijíes Riobamba 0.703 Feme… Casado
## # ℹ abbreviated name: ¹`¿En qué cantón de la provincia de Chimborazo vive?`
## # ℹ 33 more variables: `Cantidad de personas que conforman el hogar` <dbl>,
## # `¿Qué parentesco tiene con el representante del hogar?` <fct>,
## # `¿Cuál es el tipo de vivienda?` <fct>,
## # `El material predominante de las paredes exteriores de la vivienda es de:` <fct>,
## # `EL material predominante del piso de la vivienda es de:` <fct>,
## # `¿Cuántos cuartos de baño con ducha de uso exclusivo tiene este hogar?` <fct>, …VARIABLE CONSUMO DE HUEVOS
x=model.matrix(`¿Consume usted huevos?`~., data=Hogar_clean)
head(x)## (Intercept) UBICACIÓNYaruquíes UBICACIÓNColta UBICACIÓNGuano UBICACIÓNGuamote
## 1 1 0 0 0 0
## 2 1 0 0 0 0
## 3 1 0 0 0 0
## 4 1 0 0 0 0
## 5 1 0 0 0 0
## 6 1 0 0 0 0
## UBICACIÓNVeloz UBICACIÓNMaldonado UBICACIÓNCacha UBICACIÓNLizarzaburu
## 1 0 0 0 0
## 2 0 0 0 0
## 3 0 0 0 0
## 4 0 0 0 0
## 5 0 0 0 0
## 6 0 0 0 0
## UBICACIÓNSan Luis UBICACIÓNPenipe UBICACIÓNSan Juan UBICACIÓNPallatanga
## 1 0 0 0 0
## 2 0 0 0 0
## 3 0 0 0 0
## 4 0 0 0 0
## 5 0 0 0 0
## 6 0 0 0 0
## UBICACIÓNPunín UBICACIÓNPungalá UBICACIÓNChambo UBICACIÓNAlausí
## 1 0 0 0 0
## 2 0 0 0 0
## 3 0 0 0 0
## 4 0 0 0 0
## 5 0 0 0 0
## 6 0 0 0 0
## UBICACIÓNLicán UBICACIÓNCumanda UBICACIÓNLicto UBICACIÓNCalpi
## 1 0 0 0 0
## 2 0 0 0 0
## 3 0 0 0 0
## 4 0 0 0 0
## 5 0 0 0 0
## 6 0 0 0 0
## UBICACIÓNChunchi UBICACIÓNQuimiag UBICACIÓNCubijíes UBICACIÓNFlores
## 1 0 1 0 0
## 2 0 1 0 0
## 3 0 0 0 1
## 4 0 1 0 0
## 5 0 0 1 0
## 6 0 0 1 0
## `¿En qué cantón de la provincia de Chimborazo vive?`Chambo
## 1 0
## 2 0
## 3 0
## 4 0
## 5 0
## 6 0
## `¿En qué cantón de la provincia de Chimborazo vive?`Chunchi
## 1 0
## 2 0
## 3 0
## 4 0
## 5 0
## 6 0
## `¿En qué cantón de la provincia de Chimborazo vive?`Colta
## 1 0
## 2 0
## 3 0
## 4 0
## 5 0
## 6 0
## `¿En qué cantón de la provincia de Chimborazo vive?`Cumanda
## 1 0
## 2 0
## 3 0
## 4 0
## 5 0
## 6 0
## `¿En qué cantón de la provincia de Chimborazo vive?`Guamote
## 1 0
## 2 0
## 3 0
## 4 0
## 5 0
## 6 0
## `¿En qué cantón de la provincia de Chimborazo vive?`Guano
## 1 0
## 2 0
## 3 0
## 4 0
## 5 0
## 6 0
## `¿En qué cantón de la provincia de Chimborazo vive?`Pallatanga
## 1 0
## 2 0
## 3 0
## 4 0
## 5 0
## 6 0
## `¿En qué cantón de la provincia de Chimborazo vive?`Penipe
## 1 0
## 2 0
## 3 0
## 4 0
## 5 0
## 6 0
## `¿En qué cantón de la provincia de Chimborazo vive?`Riobamba
## 1 1
## 2 1
## 3 1
## 4 1
## 5 1
## 6 1
## `Edad del encuestado` SexoFemenino `Estado civil`Divorciado
## 1 1.4185412 1 0
## 2 -0.7269160 1 0
## 3 0.1312669 1 0
## 4 -1.6566141 1 0
## 5 -1.0844922 1 0
## 6 0.7033888 1 0
## `Estado civil`Otro `Estado civil`Soltero `Estado civil`Viudo
## 1 0 0 0
## 2 0 0 0
## 3 0 0 0
## 4 0 1 0
## 5 0 1 0
## 6 0 0 0
## `Cantidad de personas que conforman el hogar`
## 1 0.7176377
## 2 -0.5627438
## 3 -0.1359500
## 4 -0.1359500
## 5 0.7176377
## 6 0.2908438
## `¿Qué parentesco tiene con el representante del hogar?`Otros parientes
## 1 0
## 2 0
## 3 0
## 4 0
## 5 0
## 6 0
## `¿Qué parentesco tiene con el representante del hogar?`Yo soy el representante del Hogar
## 1 0
## 2 0
## 3 0
## 4 0
## 5 0
## 6 0
## `¿Qué parentesco tiene con el representante del hogar?`Hijo o hija
## 1 0
## 2 0
## 3 0
## 4 0
## 5 0
## 6 0
## `¿Qué parentesco tiene con el representante del hogar?`Padres o suegros
## 1 0
## 2 0
## 3 0
## 4 1
## 5 1
## 6 0
## `¿Qué parentesco tiene con el representante del hogar?`Otros no parientes
## 1 1
## 2 1
## 3 1
## 4 0
## 5 0
## 6 1
## `¿Qué parentesco tiene con el representante del hogar?`Yerno o nuera
## 1 0
## 2 0
## 3 0
## 4 0
## 5 0
## 6 0
## `¿Cuál es el tipo de vivienda?`Departamento en casa o edificio
## 1 0
## 2 0
## 3 0
## 4 0
## 5 0
## 6 0
## `¿Cuál es el tipo de vivienda?`Suite de lujo
## 1 0
## 2 0
## 3 0
## 4 0
## 5 0
## 6 0
## `¿Cuál es el tipo de vivienda?`Cuarto(s) en casa de inquilinato
## 1 0
## 2 1
## 3 0
## 4 0
## 5 1
## 6 0
## `¿Cuál es el tipo de vivienda?`Choza/Covacha/Otro
## 1 0
## 2 0
## 3 0
## 4 0
## 5 0
## 6 0
## `¿Cuál es el tipo de vivienda?`Rancho
## 1 0
## 2 0
## 3 0
## 4 0
## 5 0
## 6 0
## `¿Cuál es el tipo de vivienda?` Mediagua
## 1 0
## 2 0
## 3 0
## 4 0
## 5 0
## 6 0
## `El material predominante de las paredes exteriores de la vivienda es de:`Caña no revestida/ Otros materiales
## 1 0
## 2 0
## 3 0
## 4 0
## 5 0
## 6 0
## `El material predominante de las paredes exteriores de la vivienda es de:`Ladrillo o bloque
## 1 0
## 2 1
## 3 0
## 4 0
## 5 0
## 6 0
## `El material predominante de las paredes exteriores de la vivienda es de:`Hormigón
## 1 0
## 2 0
## 3 0
## 4 0
## 5 0
## 6 0
## `El material predominante de las paredes exteriores de la vivienda es de:` Caña revestida o bahareque/Madera
## 1 1
## 2 0
## 3 1
## 4 1
## 5 1
## 6 1
## `EL material predominante del piso de la vivienda es de:` Cerámica, baldosa, vinil o marmetón
## 1 0
## 2 0
## 3 0
## 4 1
## 5 0
## 6 1
## `EL material predominante del piso de la vivienda es de:`Tabla sin tratar
## 1 0
## 2 0
## 3 0
## 4 0
## 5 0
## 6 0
## `EL material predominante del piso de la vivienda es de:`Tierra/ Caña/ Otros materiales
## 1 0
## 2 0
## 3 0
## 4 0
## 5 1
## 6 0
## `EL material predominante del piso de la vivienda es de:`Ladrillo o cemento
## 1 0
## 2 0
## 3 0
## 4 0
## 5 0
## 6 0
## `¿Cuántos cuartos de baño con ducha de uso exclusivo tiene este hogar?`Tiene 2 cuartos de baño exclusivos con ducha
## 1 0
## 2 0
## 3 0
## 4 1
## 5 1
## 6 0
## `¿Cuántos cuartos de baño con ducha de uso exclusivo tiene este hogar?`Tiene 1 cuarto de baño exclusivo con ducha
## 1 0
## 2 0
## 3 1
## 4 0
## 5 0
## 6 1
## `¿Cuántos cuartos de baño con ducha de uso exclusivo tiene este hogar?`Tiene 3 o más cuartos de baño exclusivos con ducha
## 1 1
## 2 1
## 3 0
## 4 0
## 5 0
## 6 0
## `Tipo de servicio higiénico con que cuenta este hogar`Conectado a pozo séptico
## 1 0
## 2 0
## 3 0
## 4 0
## 5 0
## 6 0
## `Tipo de servicio higiénico con que cuenta este hogar`Tiene 1 cuarto de baño exclusivo con ducha Letrina
## 1 0
## 2 0
## 3 0
## 4 0
## 5 0
## 6 0
## `Tipo de servicio higiénico con que cuenta este hogar`Conectado a pozo ciego
## 1 1
## 2 1
## 3 1
## 4 1
## 5 0
## 6 1
## `Tipo de servicio higiénico con que cuenta este hogar`Conectado a red pública de alcantarillado
## 1 0
## 2 0
## 3 0
## 4 0
## 5 1
## 6 0
## `Tipo de servicio higiénico con que cuenta este hogar`Sin Marca
## 1 0
## 2 0
## 3 0
## 4 0
## 5 0
## 6 0
## `¿Cuenta con servicio de internet?`Si
## 1 1
## 2 1
## 3 1
## 4 1
## 5 1
## 6 1
## `¿Cuenta con computadora de escritorio?`Si
## 1 1
## 2 1
## 3 0
## 4 1
## 5 1
## 6 1
## `¿Cuenta con computadora portátil?`Si
## 1 1
## 2 1
## 3 1
## 4 1
## 5 1
## 6 1
## `¿Cuántos celulares activados tiene en este hogar?`Tiene 3 celulares
## 1 0
## 2 0
## 3 0
## 4 0
## 5 0
## 6 0
## `¿Cuántos celulares activados tiene en este hogar?`Tiene 1 celular
## 1 0
## 2 1
## 3 0
## 4 0
## 5 1
## 6 0
## `¿Cuántos celulares activados tiene en este hogar?`Tiene 4 ó más celulares
## 1 0
## 2 0
## 3 1
## 4 1
## 5 0
## 6 1
## `¿Cuántos celulares activados tiene en este hogar?`Tiene 2 celulares
## 1 1
## 2 0
## 3 0
## 4 0
## 5 0
## 6 0
## `¿Tiene su hogar servicio de teléfono convencional?`Si
## 1 1
## 2 1
## 3 1
## 4 1
## 5 1
## 6 1
## `¿Tiene en su hogar cocina con horno?`Si `¿Tiene refrigeradora?`Si
## 1 1 1
## 2 0 1
## 3 1 1
## 4 1 1
## 5 1 1
## 6 1 1
## `¿Tiene lavadora?`Si `¿Tiene equipo de sonido?`Si
## 1 1 1
## 2 1 1
## 3 0 1
## 4 0 1
## 5 1 1
## 6 1 1
## `¿Cuántas TV a color tienen en este hogar?`Tiene 3 ó más TV a color
## 1 0
## 2 1
## 3 1
## 4 0
## 5 0
## 6 0
## `¿Cuántas TV a color tienen en este hogar?`Tiene 1 TV a color
## 1 0
## 2 0
## 3 0
## 4 1
## 5 0
## 6 1
## `¿Cuántas TV a color tienen en este hogar?`Tiene 2 TV a color
## 1 1
## 2 0
## 3 0
## 4 0
## 5 1
## 6 0
## `¿Cuántos vehículos de uso exclusivo tiene este hogar?`Tiene 3 ó más vehículo exclusivo
## 1 1
## 2 0
## 3 0
## 4 1
## 5 1
## 6 0
## `¿Cuántos vehículos de uso exclusivo tiene este hogar?`Tiene 2 vehículo exclusivo
## 1 0
## 2 0
## 3 0
## 4 0
## 5 0
## 6 0
## `¿Cuántos vehículos de uso exclusivo tiene este hogar?`Tiene un vehículo exclusivo
## 1 0
## 2 0
## 3 0
## 4 0
## 5 0
## 6 1
## `¿Alguien en el hogar compra vestimenta en centros comerciales?`Si
## 1 1
## 2 1
## 3 1
## 4 1
## 5 1
## 6 1
## `¿En el hogar alguien ha usado internet en los últimos 6 meses?`Si
## 1 1
## 2 1
## 3 1
## 4 1
## 5 0
## 6 1
## `¿En el hogar alguien utiliza correo electrónico que no es del trabajo?`Si
## 1 1
## 2 1
## 3 1
## 4 1
## 5 0
## 6 1
## `¿En el hogar alguien está registrado en una red social?`Si
## 1 1
## 2 1
## 3 1
## 4 1
## 5 1
## 6 1
## `Exceptuando los libros de texto o manuales de estudio y lecturas de trabajo¿Alguien del hogar ha leído algún libro completo en los últimos 3 meses?`Si
## 1 1
## 2 1
## 3 0
## 4 1
## 5 1
## 6 1
## `¿Cuál es el nivel de instrucción del Jefe del hogar?`Hasta 3 años de educación superior (incompleta)
## 1 0
## 2 0
## 3 1
## 4 0
## 5 0
## 6 0
## `¿Cuál es el nivel de instrucción del Jefe del hogar?`Postgrado
## 1 0
## 2 0
## 3 0
## 4 1
## 5 0
## 6 0
## `¿Cuál es el nivel de instrucción del Jefe del hogar?`Primaria completa
## 1 0
## 2 0
## 3 0
## 4 0
## 5 0
## 6 0
## `¿Cuál es el nivel de instrucción del Jefe del hogar?`Primaria incompleta
## 1 0
## 2 0
## 3 0
## 4 0
## 5 0
## 6 0
## `¿Cuál es el nivel de instrucción del Jefe del hogar?`Secundaria completa
## 1 0
## 2 0
## 3 0
## 4 0
## 5 1
## 6 0
## `¿Cuál es el nivel de instrucción del Jefe del hogar?`Secundaria incompleta
## 1 0
## 2 1
## 3 0
## 4 0
## 5 0
## 6 0
## `¿Cuál es el nivel de instrucción del Jefe del hogar?`Sin estudios
## 1 0
## 2 0
## 3 0
## 4 0
## 5 0
## 6 0
## `¿Alguien en el hogar está afiliado o cubierto por el seguro del IESS (general, voluntario o campesino) y/o seguro del ISSFA o ISSPOL?`Si
## 1 1
## 2 1
## 3 1
## 4 1
## 5 0
## 6 1
## `¿Alguien en el hogar tiene seguro de salud privada con hospitalización, seguro de salud privada sin hospitalización, seguro internacional,seguros municipales y de Consejos Provinciales y/o seguro de vida?`No
## 1 0
## 2 1
## 3 0
## 4 0
## 5 1
## 6 1
## `¿Cuál es la ocupación del Representante del hogar?` Empleados de oficina
## 1 0
## 2 0
## 3 0
## 4 0
## 5 0
## 6 0
## `¿Cuál es la ocupación del Representante del hogar?`Fuerzas Armadas
## 1 1
## 2 1
## 3 0
## 4 0
## 5 0
## 6 0
## `¿Cuál es la ocupación del Representante del hogar?`Inactivos
## 1 0
## 2 0
## 3 0
## 4 0
## 5 0
## 6 0
## `¿Cuál es la ocupación del Representante del hogar?` Oficiales operarios y artesanos
## 1 0
## 2 0
## 3 0
## 4 1
## 5 0
## 6 0
## `¿Cuál es la ocupación del Representante del hogar?` Operadores de instalaciones y máquinas
## 1 0
## 2 0
## 3 0
## 4 0
## 5 0
## 6 0
## `¿Cuál es la ocupación del Representante del hogar?`Personal directivo de la Administración Pública y de empresas
## 1 0
## 2 0
## 3 0
## 4 0
## 5 0
## 6 0
## `¿Cuál es la ocupación del Representante del hogar?` Profesionales científicos e intelectuales
## 1 0
## 2 0
## 3 0
## 4 0
## 5 1
## 6 0
## `¿Cuál es la ocupación del Representante del hogar?`Técnicos y profesionales de nivel medio
## 1 0
## 2 0
## 3 0
## 4 0
## 5 0
## 6 0
## `¿Cuál es la ocupación del Representante del hogar?`Trabajador calificados agropecuarios y pesqueros
## 1 0
## 2 0
## 3 0
## 4 0
## 5 0
## 6 0
## `¿Cuál es la ocupación del Representante del hogar?`Trabajador de los servicios y comerciantes
## 1 0
## 2 0
## 3 1
## 4 0
## 5 0
## 6 1
## `¿Cuál es la ocupación del Representante del hogar?`Trabajadores no calificados
## 1 0
## 2 0
## 3 0
## 4 0
## 5 0
## 6 0
## NSE_2023_CUALIBajo NSE_2023_CUALIMedio Alto NSE_2023_CUALIMedio Bajo
## 1 0 0 0
## 2 0 1 0
## 3 0 1 0
## 4 0 0 0
## 5 0 0 0
## 6 0 0 0
## NSE_2023_CUALIMedio Típico `¿Para qué sector trabaja?`Sector privado
## 1 0 0
## 2 0 0
## 3 0 0
## 4 0 1
## 5 1 0
## 6 0 1
## `¿Para qué sector trabaja?`Ninguno `¿Consume usted arroz?`Si
## 1 1 1
## 2 1 1
## 3 1 1
## 4 0 1
## 5 1 0
## 6 0 1
## `¿Consume usted queso industrial (CON MARCA)`Si `¿Consume frutas?`Si
## 1 1 1
## 2 1 1
## 3 1 1
## 4 1 1
## 5 1 1
## 6 1 1Interpretación
Observará que las variables numéricas reescaladas contienen los mismos valores iniciales. Cada una de las variables nominales ha sido convertida en variables binarias. El número de variables binarias es siempre el número de categorías de la variable nominal. El resto de variables también son convertida a Binarias con el número de categorias que cada una almacena. Los parámetros correspondientes a las categorías retenidas comparan a éstas con la categoría eliminada.
Dividiremos el conjunto completo de individuos en dos partes, un para entrenar el modelo, que contiene el 80% de los individuos a Binarias con el número de categorias que cada una almacena. Los parámetros correspondientes a las categorías retenidas comparan a éstas con la categoría eliminada.
Dividiremos el conjunto completo de individuos en dos partes, un para entrenar el modelo, que contiene el 80% de los individuos y otra para validarlo que contiene el resto. Esto es asi porque si valoramos el modelo con las misma observaciones que usamos para ajustarlo, la bondad del ajuste quedará sobrevalorada.
Por ejemplo, vamos ausar el 70% de los individuos para ajustar el modelo y el 30% para probarlo. Para que siempre obtengamos los mismos resultados colocamos la semilla de los números pseudo aleatorios con el número de cedula + la fecha de nacimiento.
tr=round(nrow(Hogar_clean)*0.7)
set.seed(5492036)
muestra=sample.int(nrow(Hogar_clean), tr)
Train.hogar=Hogar_clean[muestra,]
val.hogar=Hogar_clean[-muestra,]Ahora disponemos de un conjunto de entrenamiento Train.hogar y un conjunto de validación Val.hogar”.
MODELOS DE CLASIFICACIÓN
Regresión Logística
\[{p_i}=P\left( {Y = 1/x = {x_i}} \right) = {{{e^{{\beta _0} + {\beta_1}{x_{i1}} + \ldots + {\beta_p}{x_{ip}}}}} \over {1 + {e^{{\beta _0} + {\beta _1} {x_{i1}} + \ldots + {\beta _p} {x_{ip}}}}}}\]
Trataremos de modelar la probabilidad de CONSUMO DE HUEVO positiva en función del resto de las variables
gfit3 <- glm(`¿Consume usted huevos?`~., data=Hogar_clean, family=binomial)## Warning: glm.fit: algorithm did not converge## Warning: glm.fit: fitted probabilities numerically 0 or 1 occurredsummary(gfit3)##
## Call:
## glm(formula = `¿Consume usted huevos?` ~ ., family = binomial,
## data = Hogar_clean)
##
## Coefficients: (9 not defined because of singularities)
## Estimate
## (Intercept) 1.004e+03
## UBICACIÓNYaruquíes -4.522e+02
## UBICACIÓNColta -3.890e+02
## UBICACIÓNGuano -2.150e+02
## UBICACIÓNGuamote -2.952e+02
## UBICACIÓNVeloz -4.931e+02
## UBICACIÓNMaldonado -3.953e+02
## UBICACIÓNCacha -3.036e+02
## UBICACIÓNLizarzaburu -3.220e+02
## UBICACIÓNSan Luis -7.345e+01
## UBICACIÓNPenipe -3.887e+02
## UBICACIÓNSan Juan -3.943e+01
## UBICACIÓNPallatanga 3.296e+02
## UBICACIÓNPunín -1.363e+02
## UBICACIÓNPungalá -2.026e+02
## UBICACIÓNChambo 2.101e+02
## UBICACIÓNAlausí -3.053e+02
## UBICACIÓNLicán -1.111e+02
## UBICACIÓNCumanda -2.909e+02
## UBICACIÓNLicto -5.383e+02
## UBICACIÓNCalpi -1.979e+02
## UBICACIÓNChunchi 8.621e+01
## UBICACIÓNQuimiag -2.736e+02
## UBICACIÓNCubijíes -1.152e+02
## UBICACIÓNFlores -1.645e+01
## `¿En qué cantón de la provincia de Chimborazo vive?`Chambo NA
## `¿En qué cantón de la provincia de Chimborazo vive?`Chunchi NA
## `¿En qué cantón de la provincia de Chimborazo vive?`Colta NA
## `¿En qué cantón de la provincia de Chimborazo vive?`Cumanda NA
## `¿En qué cantón de la provincia de Chimborazo vive?`Guamote NA
## `¿En qué cantón de la provincia de Chimborazo vive?`Guano NA
## `¿En qué cantón de la provincia de Chimborazo vive?`Pallatanga NA
## `¿En qué cantón de la provincia de Chimborazo vive?`Penipe NA
## `¿En qué cantón de la provincia de Chimborazo vive?`Riobamba NA
## `Edad del encuestado` 6.237e+01
## SexoFemenino 3.844e+01
## `Estado civil`Divorciado -1.497e+02
## `Estado civil`Otro -4.365e+01
## `Estado civil`Soltero -4.260e+01
## `Estado civil`Viudo 8.397e+01
## `Cantidad de personas que conforman el hogar` -4.132e+01
## `¿Qué parentesco tiene con el representante del hogar?`Otros parientes 4.377e+01
## `¿Qué parentesco tiene con el representante del hogar?`Yo soy el representante del Hogar -5.804e+01
## `¿Qué parentesco tiene con el representante del hogar?`Hijo o hija -1.969e+01
## `¿Qué parentesco tiene con el representante del hogar?`Padres o suegros -5.276e+02
## `¿Qué parentesco tiene con el representante del hogar?`Otros no parientes -7.759e+00
## `¿Cuál es el tipo de vivienda?`Departamento en casa o edificio 3.859e+00
## `¿Cuál es el tipo de vivienda?`Suite de lujo 2.604e+02
## `¿Cuál es el tipo de vivienda?`Cuarto(s) en casa de inquilinato 1.739e+01
## `¿Cuál es el tipo de vivienda?`Choza/Covacha/Otro 5.206e+01
## `¿Cuál es el tipo de vivienda?`Rancho 2.644e+02
## `¿Cuál es el tipo de vivienda?` Mediagua -1.060e+03
## `El material predominante de las paredes exteriores de la vivienda es de:`Caña no revestida/ Otros materiales -9.303e+02
## `El material predominante de las paredes exteriores de la vivienda es de:`Ladrillo o bloque -4.161e+02
## `El material predominante de las paredes exteriores de la vivienda es de:`Hormigón -3.966e+02
## `El material predominante de las paredes exteriores de la vivienda es de:` Caña revestida o bahareque/Madera -3.716e+02
## `EL material predominante del piso de la vivienda es de:` Cerámica, baldosa, vinil o marmetón 1.627e+01
## `EL material predominante del piso de la vivienda es de:`Tabla sin tratar -2.231e+01
## `EL material predominante del piso de la vivienda es de:`Tierra/ Caña/ Otros materiales -6.983e+01
## `EL material predominante del piso de la vivienda es de:`Ladrillo o cemento 7.934e+00
## `¿Cuántos cuartos de baño con ducha de uso exclusivo tiene este hogar?`Tiene 2 cuartos de baño exclusivos con ducha -1.651e+02
## `¿Cuántos cuartos de baño con ducha de uso exclusivo tiene este hogar?`Tiene 1 cuarto de baño exclusivo con ducha -2.087e+02
## `¿Cuántos cuartos de baño con ducha de uso exclusivo tiene este hogar?`Tiene 3 o más cuartos de baño exclusivos con ducha -2.051e+02
## `Tipo de servicio higiénico con que cuenta este hogar`Conectado a pozo séptico 2.980e+02
## `Tipo de servicio higiénico con que cuenta este hogar`Tiene 1 cuarto de baño exclusivo con ducha Letrina -1.308e+02
## `Tipo de servicio higiénico con que cuenta este hogar`Conectado a pozo ciego -2.585e+02
## `Tipo de servicio higiénico con que cuenta este hogar`Conectado a red pública de alcantarillado -5.354e+02
## `Tipo de servicio higiénico con que cuenta este hogar`Sin Marca 3.321e+00
## `¿Cuenta con servicio de internet?`Si -2.780e+02
## `¿Cuenta con computadora de escritorio?`Si -7.396e+01
## `¿Cuenta con computadora portátil?`Si 7.007e+01
## `¿Cuántos celulares activados tiene en este hogar?`Tiene 3 celulares -3.535e+01
## `¿Cuántos celulares activados tiene en este hogar?`Tiene 1 celular -2.334e+01
## `¿Cuántos celulares activados tiene en este hogar?`Tiene 4 ó más celulares -1.918e+00
## `¿Cuántos celulares activados tiene en este hogar?`Tiene 2 celulares -3.002e+01
## `¿Tiene su hogar servicio de teléfono convencional?`Si -1.366e+01
## `¿Tiene en su hogar cocina con horno?`Si 1.342e+02
## `¿Tiene refrigeradora?`Si -4.113e+01
## `¿Tiene lavadora?`Si 1.357e+02
## `¿Tiene equipo de sonido?`Si 4.921e+00
## `¿Cuántas TV a color tienen en este hogar?`Tiene 3 ó más TV a color 1.960e+02
## `¿Cuántas TV a color tienen en este hogar?`Tiene 1 TV a color 2.119e+02
## `¿Cuántas TV a color tienen en este hogar?`Tiene 2 TV a color 8.642e+01
## `¿Cuántos vehículos de uso exclusivo tiene este hogar?`Tiene 3 ó más vehículo exclusivo 1.729e+01
## `¿Cuántos vehículos de uso exclusivo tiene este hogar?`Tiene 2 vehículo exclusivo -1.215e+02
## `¿Cuántos vehículos de uso exclusivo tiene este hogar?`Tiene un vehículo exclusivo -6.777e+01
## `¿Alguien en el hogar compra vestimenta en centros comerciales?`Si 8.609e+01
## `¿En el hogar alguien ha usado internet en los últimos 6 meses?`Si 1.832e+02
## `¿En el hogar alguien utiliza correo electrónico que no es del trabajo?`Si 1.129e+02
## `¿En el hogar alguien está registrado en una red social?`Si -7.700e+01
## `Exceptuando los libros de texto o manuales de estudio y lecturas de trabajo¿Alguien del hogar ha leído algún libro completo en los últimos 3 meses?`Si -3.584e+01
## `¿Cuál es el nivel de instrucción del Jefe del hogar?`Hasta 3 años de educación superior (incompleta) 3.197e+01
## `¿Cuál es el nivel de instrucción del Jefe del hogar?`Postgrado -8.190e+01
## `¿Cuál es el nivel de instrucción del Jefe del hogar?`Primaria completa -6.147e-01
## `¿Cuál es el nivel de instrucción del Jefe del hogar?`Primaria incompleta -1.682e+02
## `¿Cuál es el nivel de instrucción del Jefe del hogar?`Secundaria completa -1.489e+01
## `¿Cuál es el nivel de instrucción del Jefe del hogar?`Secundaria incompleta 1.544e+02
## `¿Cuál es el nivel de instrucción del Jefe del hogar?`Sin estudios 2.669e+02
## `¿Alguien en el hogar está afiliado o cubierto por el seguro del IESS (general, voluntario o campesino) y/o seguro del ISSFA o ISSPOL?`Si 1.086e+02
## `¿Alguien en el hogar tiene seguro de salud privada con hospitalización, seguro de salud privada sin hospitalización, seguro internacional,seguros municipales y de Consejos Provinciales y/o seguro de vida?`No -8.853e+01
## `¿Cuál es la ocupación del Representante del hogar?` Empleados de oficina -3.072e+02
## `¿Cuál es la ocupación del Representante del hogar?`Fuerzas Armadas -2.762e+02
## `¿Cuál es la ocupación del Representante del hogar?`Inactivos -1.626e+02
## `¿Cuál es la ocupación del Representante del hogar?` Oficiales operarios y artesanos 1.971e+02
## `¿Cuál es la ocupación del Representante del hogar?` Operadores de instalaciones y máquinas 2.469e+01
## `¿Cuál es la ocupación del Representante del hogar?`Personal directivo de la Administración Pública y de empresas -1.777e+02
## `¿Cuál es la ocupación del Representante del hogar?` Profesionales científicos e intelectuales -2.603e+02
## `¿Cuál es la ocupación del Representante del hogar?`Técnicos y profesionales de nivel medio -3.310e+02
## `¿Cuál es la ocupación del Representante del hogar?`Trabajador calificados agropecuarios y pesqueros -1.498e+02
## `¿Cuál es la ocupación del Representante del hogar?`Trabajador de los servicios y comerciantes -1.489e+02
## `¿Cuál es la ocupación del Representante del hogar?`Trabajadores no calificados -2.415e+02
## NSE_2023_CUALIBajo 1.451e+02
## NSE_2023_CUALIMedio Alto -4.065e+00
## NSE_2023_CUALIMedio Bajo 1.724e+02
## NSE_2023_CUALIMedio Típico 2.625e+01
## `¿Para qué sector trabaja?`Sector privado 1.801e+02
## `¿Para qué sector trabaja?`Ninguno 8.328e+01
## `¿Consume usted arroz?`Si -2.166e+02
## `¿Consume usted queso industrial (CON MARCA)`Si -4.472e+01
## `¿Consume frutas?`Si 2.113e+02
## Std. Error
## (Intercept) 9.216e+05
## UBICACIÓNYaruquíes 8.447e+05
## UBICACIÓNColta 6.420e+05
## UBICACIÓNGuano 6.784e+05
## UBICACIÓNGuamote 8.397e+05
## UBICACIÓNVeloz 5.065e+05
## UBICACIÓNMaldonado 8.631e+05
## UBICACIÓNCacha 7.648e+05
## UBICACIÓNLizarzaburu 8.691e+05
## UBICACIÓNSan Luis 6.693e+05
## UBICACIÓNPenipe 5.926e+05
## UBICACIÓNSan Juan 6.988e+05
## UBICACIÓNPallatanga 7.136e+05
## UBICACIÓNPunín 6.916e+05
## UBICACIÓNPungalá 5.968e+05
## UBICACIÓNChambo 7.511e+05
## UBICACIÓNAlausí 7.220e+05
## UBICACIÓNLicán 7.520e+05
## UBICACIÓNCumanda 7.883e+05
## UBICACIÓNLicto 6.557e+05
## UBICACIÓNCalpi 6.414e+05
## UBICACIÓNChunchi 6.446e+05
## UBICACIÓNQuimiag 6.336e+05
## UBICACIÓNCubijíes 5.904e+05
## UBICACIÓNFlores 6.166e+05
## `¿En qué cantón de la provincia de Chimborazo vive?`Chambo NA
## `¿En qué cantón de la provincia de Chimborazo vive?`Chunchi NA
## `¿En qué cantón de la provincia de Chimborazo vive?`Colta NA
## `¿En qué cantón de la provincia de Chimborazo vive?`Cumanda NA
## `¿En qué cantón de la provincia de Chimborazo vive?`Guamote NA
## `¿En qué cantón de la provincia de Chimborazo vive?`Guano NA
## `¿En qué cantón de la provincia de Chimborazo vive?`Pallatanga NA
## `¿En qué cantón de la provincia de Chimborazo vive?`Penipe NA
## `¿En qué cantón de la provincia de Chimborazo vive?`Riobamba NA
## `Edad del encuestado` 2.768e+04
## SexoFemenino 6.335e+04
## `Estado civil`Divorciado 1.281e+05
## `Estado civil`Otro 1.316e+05
## `Estado civil`Soltero 9.205e+04
## `Estado civil`Viudo 8.610e+04
## `Cantidad de personas que conforman el hogar` 2.662e+04
## `¿Qué parentesco tiene con el representante del hogar?`Otros parientes 6.794e+04
## `¿Qué parentesco tiene con el representante del hogar?`Yo soy el representante del Hogar 3.921e+05
## `¿Qué parentesco tiene con el representante del hogar?`Hijo o hija 1.032e+05
## `¿Qué parentesco tiene con el representante del hogar?`Padres o suegros 4.259e+05
## `¿Qué parentesco tiene con el representante del hogar?`Otros no parientes 1.060e+05
## `¿Cuál es el tipo de vivienda?`Departamento en casa o edificio 3.348e+05
## `¿Cuál es el tipo de vivienda?`Suite de lujo 1.476e+05
## `¿Cuál es el tipo de vivienda?`Cuarto(s) en casa de inquilinato 1.321e+05
## `¿Cuál es el tipo de vivienda?`Choza/Covacha/Otro 1.056e+05
## `¿Cuál es el tipo de vivienda?`Rancho 2.235e+05
## `¿Cuál es el tipo de vivienda?` Mediagua 7.106e+05
## `El material predominante de las paredes exteriores de la vivienda es de:`Caña no revestida/ Otros materiales 4.441e+05
## `El material predominante de las paredes exteriores de la vivienda es de:`Ladrillo o bloque 3.022e+05
## `El material predominante de las paredes exteriores de la vivienda es de:`Hormigón 3.575e+05
## `El material predominante de las paredes exteriores de la vivienda es de:` Caña revestida o bahareque/Madera 3.671e+05
## `EL material predominante del piso de la vivienda es de:` Cerámica, baldosa, vinil o marmetón 7.132e+04
## `EL material predominante del piso de la vivienda es de:`Tabla sin tratar 7.114e+04
## `EL material predominante del piso de la vivienda es de:`Tierra/ Caña/ Otros materiales 3.473e+05
## `EL material predominante del piso de la vivienda es de:`Ladrillo o cemento 1.727e+05
## `¿Cuántos cuartos de baño con ducha de uso exclusivo tiene este hogar?`Tiene 2 cuartos de baño exclusivos con ducha 2.961e+05
## `¿Cuántos cuartos de baño con ducha de uso exclusivo tiene este hogar?`Tiene 1 cuarto de baño exclusivo con ducha 2.724e+05
## `¿Cuántos cuartos de baño con ducha de uso exclusivo tiene este hogar?`Tiene 3 o más cuartos de baño exclusivos con ducha 3.274e+05
## `Tipo de servicio higiénico con que cuenta este hogar`Conectado a pozo séptico 3.177e+05
## `Tipo de servicio higiénico con que cuenta este hogar`Tiene 1 cuarto de baño exclusivo con ducha Letrina 2.386e+05
## `Tipo de servicio higiénico con que cuenta este hogar`Conectado a pozo ciego 2.207e+05
## `Tipo de servicio higiénico con que cuenta este hogar`Conectado a red pública de alcantarillado 2.399e+05
## `Tipo de servicio higiénico con que cuenta este hogar`Sin Marca 3.466e+05
## `¿Cuenta con servicio de internet?`Si 1.678e+05
## `¿Cuenta con computadora de escritorio?`Si 4.014e+04
## `¿Cuenta con computadora portátil?`Si 4.705e+04
## `¿Cuántos celulares activados tiene en este hogar?`Tiene 3 celulares 2.401e+05
## `¿Cuántos celulares activados tiene en este hogar?`Tiene 1 celular 1.641e+05
## `¿Cuántos celulares activados tiene en este hogar?`Tiene 4 ó más celulares 1.685e+05
## `¿Cuántos celulares activados tiene en este hogar?`Tiene 2 celulares 1.988e+05
## `¿Tiene su hogar servicio de teléfono convencional?`Si 6.474e+04
## `¿Tiene en su hogar cocina con horno?`Si 1.010e+05
## `¿Tiene refrigeradora?`Si 9.428e+04
## `¿Tiene lavadora?`Si 8.380e+04
## `¿Tiene equipo de sonido?`Si 5.683e+04
## `¿Cuántas TV a color tienen en este hogar?`Tiene 3 ó más TV a color 8.535e+04
## `¿Cuántas TV a color tienen en este hogar?`Tiene 1 TV a color 1.134e+05
## `¿Cuántas TV a color tienen en este hogar?`Tiene 2 TV a color 1.266e+05
## `¿Cuántos vehículos de uso exclusivo tiene este hogar?`Tiene 3 ó más vehículo exclusivo 1.071e+05
## `¿Cuántos vehículos de uso exclusivo tiene este hogar?`Tiene 2 vehículo exclusivo 8.389e+04
## `¿Cuántos vehículos de uso exclusivo tiene este hogar?`Tiene un vehículo exclusivo 6.501e+04
## `¿Alguien en el hogar compra vestimenta en centros comerciales?`Si 1.132e+05
## `¿En el hogar alguien ha usado internet en los últimos 6 meses?`Si 1.889e+05
## `¿En el hogar alguien utiliza correo electrónico que no es del trabajo?`Si 6.740e+04
## `¿En el hogar alguien está registrado en una red social?`Si 1.096e+05
## `Exceptuando los libros de texto o manuales de estudio y lecturas de trabajo¿Alguien del hogar ha leído algún libro completo en los últimos 3 meses?`Si 5.813e+04
## `¿Cuál es el nivel de instrucción del Jefe del hogar?`Hasta 3 años de educación superior (incompleta) 1.159e+05
## `¿Cuál es el nivel de instrucción del Jefe del hogar?`Postgrado 1.886e+05
## `¿Cuál es el nivel de instrucción del Jefe del hogar?`Primaria completa 1.706e+05
## `¿Cuál es el nivel de instrucción del Jefe del hogar?`Primaria incompleta 2.012e+05
## `¿Cuál es el nivel de instrucción del Jefe del hogar?`Secundaria completa 1.605e+05
## `¿Cuál es el nivel de instrucción del Jefe del hogar?`Secundaria incompleta 1.938e+05
## `¿Cuál es el nivel de instrucción del Jefe del hogar?`Sin estudios 2.974e+05
## `¿Alguien en el hogar está afiliado o cubierto por el seguro del IESS (general, voluntario o campesino) y/o seguro del ISSFA o ISSPOL?`Si 9.771e+04
## `¿Alguien en el hogar tiene seguro de salud privada con hospitalización, seguro de salud privada sin hospitalización, seguro internacional,seguros municipales y de Consejos Provinciales y/o seguro de vida?`No 2.940e+04
## `¿Cuál es la ocupación del Representante del hogar?` Empleados de oficina 1.441e+05
## `¿Cuál es la ocupación del Representante del hogar?`Fuerzas Armadas 1.946e+05
## `¿Cuál es la ocupación del Representante del hogar?`Inactivos 3.185e+05
## `¿Cuál es la ocupación del Representante del hogar?` Oficiales operarios y artesanos 9.517e+04
## `¿Cuál es la ocupación del Representante del hogar?` Operadores de instalaciones y máquinas 1.366e+05
## `¿Cuál es la ocupación del Representante del hogar?`Personal directivo de la Administración Pública y de empresas 1.279e+05
## `¿Cuál es la ocupación del Representante del hogar?` Profesionales científicos e intelectuales 1.828e+05
## `¿Cuál es la ocupación del Representante del hogar?`Técnicos y profesionales de nivel medio 1.501e+05
## `¿Cuál es la ocupación del Representante del hogar?`Trabajador calificados agropecuarios y pesqueros 8.851e+04
## `¿Cuál es la ocupación del Representante del hogar?`Trabajador de los servicios y comerciantes 1.095e+05
## `¿Cuál es la ocupación del Representante del hogar?`Trabajadores no calificados 6.899e+04
## NSE_2023_CUALIBajo 4.249e+05
## NSE_2023_CUALIMedio Alto 1.452e+05
## NSE_2023_CUALIMedio Bajo 4.509e+05
## NSE_2023_CUALIMedio Típico 2.910e+05
## `¿Para qué sector trabaja?`Sector privado 1.012e+05
## `¿Para qué sector trabaja?`Ninguno 8.923e+04
## `¿Consume usted arroz?`Si 2.029e+05
## `¿Consume usted queso industrial (CON MARCA)`Si 7.722e+04
## `¿Consume frutas?`Si 8.590e+04
## z value
## (Intercept) 0.001
## UBICACIÓNYaruquíes -0.001
## UBICACIÓNColta -0.001
## UBICACIÓNGuano 0.000
## UBICACIÓNGuamote 0.000
## UBICACIÓNVeloz -0.001
## UBICACIÓNMaldonado 0.000
## UBICACIÓNCacha 0.000
## UBICACIÓNLizarzaburu 0.000
## UBICACIÓNSan Luis 0.000
## UBICACIÓNPenipe -0.001
## UBICACIÓNSan Juan 0.000
## UBICACIÓNPallatanga 0.000
## UBICACIÓNPunín 0.000
## UBICACIÓNPungalá 0.000
## UBICACIÓNChambo 0.000
## UBICACIÓNAlausí 0.000
## UBICACIÓNLicán 0.000
## UBICACIÓNCumanda 0.000
## UBICACIÓNLicto -0.001
## UBICACIÓNCalpi 0.000
## UBICACIÓNChunchi 0.000
## UBICACIÓNQuimiag 0.000
## UBICACIÓNCubijíes 0.000
## UBICACIÓNFlores 0.000
## `¿En qué cantón de la provincia de Chimborazo vive?`Chambo NA
## `¿En qué cantón de la provincia de Chimborazo vive?`Chunchi NA
## `¿En qué cantón de la provincia de Chimborazo vive?`Colta NA
## `¿En qué cantón de la provincia de Chimborazo vive?`Cumanda NA
## `¿En qué cantón de la provincia de Chimborazo vive?`Guamote NA
## `¿En qué cantón de la provincia de Chimborazo vive?`Guano NA
## `¿En qué cantón de la provincia de Chimborazo vive?`Pallatanga NA
## `¿En qué cantón de la provincia de Chimborazo vive?`Penipe NA
## `¿En qué cantón de la provincia de Chimborazo vive?`Riobamba NA
## `Edad del encuestado` 0.002
## SexoFemenino 0.001
## `Estado civil`Divorciado -0.001
## `Estado civil`Otro 0.000
## `Estado civil`Soltero 0.000
## `Estado civil`Viudo 0.001
## `Cantidad de personas que conforman el hogar` -0.002
## `¿Qué parentesco tiene con el representante del hogar?`Otros parientes 0.001
## `¿Qué parentesco tiene con el representante del hogar?`Yo soy el representante del Hogar 0.000
## `¿Qué parentesco tiene con el representante del hogar?`Hijo o hija 0.000
## `¿Qué parentesco tiene con el representante del hogar?`Padres o suegros -0.001
## `¿Qué parentesco tiene con el representante del hogar?`Otros no parientes 0.000
## `¿Cuál es el tipo de vivienda?`Departamento en casa o edificio 0.000
## `¿Cuál es el tipo de vivienda?`Suite de lujo 0.002
## `¿Cuál es el tipo de vivienda?`Cuarto(s) en casa de inquilinato 0.000
## `¿Cuál es el tipo de vivienda?`Choza/Covacha/Otro 0.000
## `¿Cuál es el tipo de vivienda?`Rancho 0.001
## `¿Cuál es el tipo de vivienda?` Mediagua -0.001
## `El material predominante de las paredes exteriores de la vivienda es de:`Caña no revestida/ Otros materiales -0.002
## `El material predominante de las paredes exteriores de la vivienda es de:`Ladrillo o bloque -0.001
## `El material predominante de las paredes exteriores de la vivienda es de:`Hormigón -0.001
## `El material predominante de las paredes exteriores de la vivienda es de:` Caña revestida o bahareque/Madera -0.001
## `EL material predominante del piso de la vivienda es de:` Cerámica, baldosa, vinil o marmetón 0.000
## `EL material predominante del piso de la vivienda es de:`Tabla sin tratar 0.000
## `EL material predominante del piso de la vivienda es de:`Tierra/ Caña/ Otros materiales 0.000
## `EL material predominante del piso de la vivienda es de:`Ladrillo o cemento 0.000
## `¿Cuántos cuartos de baño con ducha de uso exclusivo tiene este hogar?`Tiene 2 cuartos de baño exclusivos con ducha -0.001
## `¿Cuántos cuartos de baño con ducha de uso exclusivo tiene este hogar?`Tiene 1 cuarto de baño exclusivo con ducha -0.001
## `¿Cuántos cuartos de baño con ducha de uso exclusivo tiene este hogar?`Tiene 3 o más cuartos de baño exclusivos con ducha -0.001
## `Tipo de servicio higiénico con que cuenta este hogar`Conectado a pozo séptico 0.001
## `Tipo de servicio higiénico con que cuenta este hogar`Tiene 1 cuarto de baño exclusivo con ducha Letrina -0.001
## `Tipo de servicio higiénico con que cuenta este hogar`Conectado a pozo ciego -0.001
## `Tipo de servicio higiénico con que cuenta este hogar`Conectado a red pública de alcantarillado -0.002
## `Tipo de servicio higiénico con que cuenta este hogar`Sin Marca 0.000
## `¿Cuenta con servicio de internet?`Si -0.002
## `¿Cuenta con computadora de escritorio?`Si -0.002
## `¿Cuenta con computadora portátil?`Si 0.001
## `¿Cuántos celulares activados tiene en este hogar?`Tiene 3 celulares 0.000
## `¿Cuántos celulares activados tiene en este hogar?`Tiene 1 celular 0.000
## `¿Cuántos celulares activados tiene en este hogar?`Tiene 4 ó más celulares 0.000
## `¿Cuántos celulares activados tiene en este hogar?`Tiene 2 celulares 0.000
## `¿Tiene su hogar servicio de teléfono convencional?`Si 0.000
## `¿Tiene en su hogar cocina con horno?`Si 0.001
## `¿Tiene refrigeradora?`Si 0.000
## `¿Tiene lavadora?`Si 0.002
## `¿Tiene equipo de sonido?`Si 0.000
## `¿Cuántas TV a color tienen en este hogar?`Tiene 3 ó más TV a color 0.002
## `¿Cuántas TV a color tienen en este hogar?`Tiene 1 TV a color 0.002
## `¿Cuántas TV a color tienen en este hogar?`Tiene 2 TV a color 0.001
## `¿Cuántos vehículos de uso exclusivo tiene este hogar?`Tiene 3 ó más vehículo exclusivo 0.000
## `¿Cuántos vehículos de uso exclusivo tiene este hogar?`Tiene 2 vehículo exclusivo -0.001
## `¿Cuántos vehículos de uso exclusivo tiene este hogar?`Tiene un vehículo exclusivo -0.001
## `¿Alguien en el hogar compra vestimenta en centros comerciales?`Si 0.001
## `¿En el hogar alguien ha usado internet en los últimos 6 meses?`Si 0.001
## `¿En el hogar alguien utiliza correo electrónico que no es del trabajo?`Si 0.002
## `¿En el hogar alguien está registrado en una red social?`Si -0.001
## `Exceptuando los libros de texto o manuales de estudio y lecturas de trabajo¿Alguien del hogar ha leído algún libro completo en los últimos 3 meses?`Si -0.001
## `¿Cuál es el nivel de instrucción del Jefe del hogar?`Hasta 3 años de educación superior (incompleta) 0.000
## `¿Cuál es el nivel de instrucción del Jefe del hogar?`Postgrado 0.000
## `¿Cuál es el nivel de instrucción del Jefe del hogar?`Primaria completa 0.000
## `¿Cuál es el nivel de instrucción del Jefe del hogar?`Primaria incompleta -0.001
## `¿Cuál es el nivel de instrucción del Jefe del hogar?`Secundaria completa 0.000
## `¿Cuál es el nivel de instrucción del Jefe del hogar?`Secundaria incompleta 0.001
## `¿Cuál es el nivel de instrucción del Jefe del hogar?`Sin estudios 0.001
## `¿Alguien en el hogar está afiliado o cubierto por el seguro del IESS (general, voluntario o campesino) y/o seguro del ISSFA o ISSPOL?`Si 0.001
## `¿Alguien en el hogar tiene seguro de salud privada con hospitalización, seguro de salud privada sin hospitalización, seguro internacional,seguros municipales y de Consejos Provinciales y/o seguro de vida?`No -0.003
## `¿Cuál es la ocupación del Representante del hogar?` Empleados de oficina -0.002
## `¿Cuál es la ocupación del Representante del hogar?`Fuerzas Armadas -0.001
## `¿Cuál es la ocupación del Representante del hogar?`Inactivos -0.001
## `¿Cuál es la ocupación del Representante del hogar?` Oficiales operarios y artesanos 0.002
## `¿Cuál es la ocupación del Representante del hogar?` Operadores de instalaciones y máquinas 0.000
## `¿Cuál es la ocupación del Representante del hogar?`Personal directivo de la Administración Pública y de empresas -0.001
## `¿Cuál es la ocupación del Representante del hogar?` Profesionales científicos e intelectuales -0.001
## `¿Cuál es la ocupación del Representante del hogar?`Técnicos y profesionales de nivel medio -0.002
## `¿Cuál es la ocupación del Representante del hogar?`Trabajador calificados agropecuarios y pesqueros -0.002
## `¿Cuál es la ocupación del Representante del hogar?`Trabajador de los servicios y comerciantes -0.001
## `¿Cuál es la ocupación del Representante del hogar?`Trabajadores no calificados -0.003
## NSE_2023_CUALIBajo 0.000
## NSE_2023_CUALIMedio Alto 0.000
## NSE_2023_CUALIMedio Bajo 0.000
## NSE_2023_CUALIMedio Típico 0.000
## `¿Para qué sector trabaja?`Sector privado 0.002
## `¿Para qué sector trabaja?`Ninguno 0.001
## `¿Consume usted arroz?`Si -0.001
## `¿Consume usted queso industrial (CON MARCA)`Si -0.001
## `¿Consume frutas?`Si 0.002
## Pr(>|z|)
## (Intercept) 0.999
## UBICACIÓNYaruquíes 1.000
## UBICACIÓNColta 1.000
## UBICACIÓNGuano 1.000
## UBICACIÓNGuamote 1.000
## UBICACIÓNVeloz 0.999
## UBICACIÓNMaldonado 1.000
## UBICACIÓNCacha 1.000
## UBICACIÓNLizarzaburu 1.000
## UBICACIÓNSan Luis 1.000
## UBICACIÓNPenipe 0.999
## UBICACIÓNSan Juan 1.000
## UBICACIÓNPallatanga 1.000
## UBICACIÓNPunín 1.000
## UBICACIÓNPungalá 1.000
## UBICACIÓNChambo 1.000
## UBICACIÓNAlausí 1.000
## UBICACIÓNLicán 1.000
## UBICACIÓNCumanda 1.000
## UBICACIÓNLicto 0.999
## UBICACIÓNCalpi 1.000
## UBICACIÓNChunchi 1.000
## UBICACIÓNQuimiag 1.000
## UBICACIÓNCubijíes 1.000
## UBICACIÓNFlores 1.000
## `¿En qué cantón de la provincia de Chimborazo vive?`Chambo NA
## `¿En qué cantón de la provincia de Chimborazo vive?`Chunchi NA
## `¿En qué cantón de la provincia de Chimborazo vive?`Colta NA
## `¿En qué cantón de la provincia de Chimborazo vive?`Cumanda NA
## `¿En qué cantón de la provincia de Chimborazo vive?`Guamote NA
## `¿En qué cantón de la provincia de Chimborazo vive?`Guano NA
## `¿En qué cantón de la provincia de Chimborazo vive?`Pallatanga NA
## `¿En qué cantón de la provincia de Chimborazo vive?`Penipe NA
## `¿En qué cantón de la provincia de Chimborazo vive?`Riobamba NA
## `Edad del encuestado` 0.998
## SexoFemenino 1.000
## `Estado civil`Divorciado 0.999
## `Estado civil`Otro 1.000
## `Estado civil`Soltero 1.000
## `Estado civil`Viudo 0.999
## `Cantidad de personas que conforman el hogar` 0.999
## `¿Qué parentesco tiene con el representante del hogar?`Otros parientes 0.999
## `¿Qué parentesco tiene con el representante del hogar?`Yo soy el representante del Hogar 1.000
## `¿Qué parentesco tiene con el representante del hogar?`Hijo o hija 1.000
## `¿Qué parentesco tiene con el representante del hogar?`Padres o suegros 0.999
## `¿Qué parentesco tiene con el representante del hogar?`Otros no parientes 1.000
## `¿Cuál es el tipo de vivienda?`Departamento en casa o edificio 1.000
## `¿Cuál es el tipo de vivienda?`Suite de lujo 0.999
## `¿Cuál es el tipo de vivienda?`Cuarto(s) en casa de inquilinato 1.000
## `¿Cuál es el tipo de vivienda?`Choza/Covacha/Otro 1.000
## `¿Cuál es el tipo de vivienda?`Rancho 0.999
## `¿Cuál es el tipo de vivienda?` Mediagua 0.999
## `El material predominante de las paredes exteriores de la vivienda es de:`Caña no revestida/ Otros materiales 0.998
## `El material predominante de las paredes exteriores de la vivienda es de:`Ladrillo o bloque 0.999
## `El material predominante de las paredes exteriores de la vivienda es de:`Hormigón 0.999
## `El material predominante de las paredes exteriores de la vivienda es de:` Caña revestida o bahareque/Madera 0.999
## `EL material predominante del piso de la vivienda es de:` Cerámica, baldosa, vinil o marmetón 1.000
## `EL material predominante del piso de la vivienda es de:`Tabla sin tratar 1.000
## `EL material predominante del piso de la vivienda es de:`Tierra/ Caña/ Otros materiales 1.000
## `EL material predominante del piso de la vivienda es de:`Ladrillo o cemento 1.000
## `¿Cuántos cuartos de baño con ducha de uso exclusivo tiene este hogar?`Tiene 2 cuartos de baño exclusivos con ducha 1.000
## `¿Cuántos cuartos de baño con ducha de uso exclusivo tiene este hogar?`Tiene 1 cuarto de baño exclusivo con ducha 0.999
## `¿Cuántos cuartos de baño con ducha de uso exclusivo tiene este hogar?`Tiene 3 o más cuartos de baño exclusivos con ducha 1.000
## `Tipo de servicio higiénico con que cuenta este hogar`Conectado a pozo séptico 0.999
## `Tipo de servicio higiénico con que cuenta este hogar`Tiene 1 cuarto de baño exclusivo con ducha Letrina 1.000
## `Tipo de servicio higiénico con que cuenta este hogar`Conectado a pozo ciego 0.999
## `Tipo de servicio higiénico con que cuenta este hogar`Conectado a red pública de alcantarillado 0.998
## `Tipo de servicio higiénico con que cuenta este hogar`Sin Marca 1.000
## `¿Cuenta con servicio de internet?`Si 0.999
## `¿Cuenta con computadora de escritorio?`Si 0.999
## `¿Cuenta con computadora portátil?`Si 0.999
## `¿Cuántos celulares activados tiene en este hogar?`Tiene 3 celulares 1.000
## `¿Cuántos celulares activados tiene en este hogar?`Tiene 1 celular 1.000
## `¿Cuántos celulares activados tiene en este hogar?`Tiene 4 ó más celulares 1.000
## `¿Cuántos celulares activados tiene en este hogar?`Tiene 2 celulares 1.000
## `¿Tiene su hogar servicio de teléfono convencional?`Si 1.000
## `¿Tiene en su hogar cocina con horno?`Si 0.999
## `¿Tiene refrigeradora?`Si 1.000
## `¿Tiene lavadora?`Si 0.999
## `¿Tiene equipo de sonido?`Si 1.000
## `¿Cuántas TV a color tienen en este hogar?`Tiene 3 ó más TV a color 0.998
## `¿Cuántas TV a color tienen en este hogar?`Tiene 1 TV a color 0.999
## `¿Cuántas TV a color tienen en este hogar?`Tiene 2 TV a color 0.999
## `¿Cuántos vehículos de uso exclusivo tiene este hogar?`Tiene 3 ó más vehículo exclusivo 1.000
## `¿Cuántos vehículos de uso exclusivo tiene este hogar?`Tiene 2 vehículo exclusivo 0.999
## `¿Cuántos vehículos de uso exclusivo tiene este hogar?`Tiene un vehículo exclusivo 0.999
## `¿Alguien en el hogar compra vestimenta en centros comerciales?`Si 0.999
## `¿En el hogar alguien ha usado internet en los últimos 6 meses?`Si 0.999
## `¿En el hogar alguien utiliza correo electrónico que no es del trabajo?`Si 0.999
## `¿En el hogar alguien está registrado en una red social?`Si 0.999
## `Exceptuando los libros de texto o manuales de estudio y lecturas de trabajo¿Alguien del hogar ha leído algún libro completo en los últimos 3 meses?`Si 1.000
## `¿Cuál es el nivel de instrucción del Jefe del hogar?`Hasta 3 años de educación superior (incompleta) 1.000
## `¿Cuál es el nivel de instrucción del Jefe del hogar?`Postgrado 1.000
## `¿Cuál es el nivel de instrucción del Jefe del hogar?`Primaria completa 1.000
## `¿Cuál es el nivel de instrucción del Jefe del hogar?`Primaria incompleta 0.999
## `¿Cuál es el nivel de instrucción del Jefe del hogar?`Secundaria completa 1.000
## `¿Cuál es el nivel de instrucción del Jefe del hogar?`Secundaria incompleta 0.999
## `¿Cuál es el nivel de instrucción del Jefe del hogar?`Sin estudios 0.999
## `¿Alguien en el hogar está afiliado o cubierto por el seguro del IESS (general, voluntario o campesino) y/o seguro del ISSFA o ISSPOL?`Si 0.999
## `¿Alguien en el hogar tiene seguro de salud privada con hospitalización, seguro de salud privada sin hospitalización, seguro internacional,seguros municipales y de Consejos Provinciales y/o seguro de vida?`No 0.998
## `¿Cuál es la ocupación del Representante del hogar?` Empleados de oficina 0.998
## `¿Cuál es la ocupación del Representante del hogar?`Fuerzas Armadas 0.999
## `¿Cuál es la ocupación del Representante del hogar?`Inactivos 1.000
## `¿Cuál es la ocupación del Representante del hogar?` Oficiales operarios y artesanos 0.998
## `¿Cuál es la ocupación del Representante del hogar?` Operadores de instalaciones y máquinas 1.000
## `¿Cuál es la ocupación del Representante del hogar?`Personal directivo de la Administración Pública y de empresas 0.999
## `¿Cuál es la ocupación del Representante del hogar?` Profesionales científicos e intelectuales 0.999
## `¿Cuál es la ocupación del Representante del hogar?`Técnicos y profesionales de nivel medio 0.998
## `¿Cuál es la ocupación del Representante del hogar?`Trabajador calificados agropecuarios y pesqueros 0.999
## `¿Cuál es la ocupación del Representante del hogar?`Trabajador de los servicios y comerciantes 0.999
## `¿Cuál es la ocupación del Representante del hogar?`Trabajadores no calificados 0.997
## NSE_2023_CUALIBajo 1.000
## NSE_2023_CUALIMedio Alto 1.000
## NSE_2023_CUALIMedio Bajo 1.000
## NSE_2023_CUALIMedio Típico 1.000
## `¿Para qué sector trabaja?`Sector privado 0.999
## `¿Para qué sector trabaja?`Ninguno 0.999
## `¿Consume usted arroz?`Si 0.999
## `¿Consume usted queso industrial (CON MARCA)`Si 1.000
## `¿Consume frutas?`Si 0.998
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 2.9677e+02 on 382 degrees of freedom
## Residual deviance: 2.0657e-07 on 272 degrees of freedom
## AIC: 222
##
## Number of Fisher Scoring iterations: 25TABLA ANOVA
Consumo de huevos
gfit02 <- glm(`¿Consume usted huevos?`~1, data = Hogar_clean, family = binomial)
anova(gfit02, gfit3, test = "Chisq")## Analysis of Deviance Table
##
## Model 1: `¿Consume usted huevos?` ~ 1
## Model 2: `¿Consume usted huevos?` ~ UBICACIÓN + `¿En qué cantón de la provincia de Chimborazo vive?` +
## `Edad del encuestado` + Sexo + `Estado civil` + `Cantidad de personas que conforman el hogar` +
## `¿Qué parentesco tiene con el representante del hogar?` +
## `¿Cuál es el tipo de vivienda?` + `El material predominante de las paredes exteriores de la vivienda es de:` +
## `EL material predominante del piso de la vivienda es de:` +
## `¿Cuántos cuartos de baño con ducha de uso exclusivo tiene este hogar?` +
## `Tipo de servicio higiénico con que cuenta este hogar` +
## `¿Cuenta con servicio de internet?` + `¿Cuenta con computadora de escritorio?` +
## `¿Cuenta con computadora portátil?` + `¿Cuántos celulares activados tiene en este hogar?` +
## `¿Tiene su hogar servicio de teléfono convencional?` +
## `¿Tiene en su hogar cocina con horno?` + `¿Tiene refrigeradora?` +
## `¿Tiene lavadora?` + `¿Tiene equipo de sonido?` + `¿Cuántas TV a color tienen en este hogar?` +
## `¿Cuántos vehículos de uso exclusivo tiene este hogar?` +
## `¿Alguien en el hogar compra vestimenta en centros comerciales?` +
## `¿En el hogar alguien ha usado internet en los últimos 6 meses?` +
## `¿En el hogar alguien utiliza correo electrónico que no es del trabajo?` +
## `¿En el hogar alguien está registrado en una red social?` +
## `Exceptuando los libros de texto o manuales de estudio y lecturas de trabajo¿Alguien del hogar ha leído algún libro completo en los últimos 3 meses?` +
## `¿Cuál es el nivel de instrucción del Jefe del hogar?` +
## `¿Alguien en el hogar está afiliado o cubierto por el seguro del IESS (general, voluntario o campesino) y/o seguro del ISSFA o ISSPOL?` +
## `¿Alguien en el hogar tiene seguro de salud privada con hospitalización, seguro de salud privada sin hospitalización, seguro internacional,seguros municipales y de Consejos Provinciales y/o seguro de vida?` +
## `¿Cuál es la ocupación del Representante del hogar?` +
## NSE_2023_CUALI + `¿Para qué sector trabaja?` + `¿Consume usted arroz?` +
## `¿Consume usted queso industrial (CON MARCA)` + `¿Consume frutas?`
## Resid. Df Resid. Dev Df Deviance Pr(>Chi)
## 1 382 296.77
## 2 272 0.00 110 296.77 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1Interpretación
Modelos: Se comparan dos modelos de regresión logística.
Modelo 1: Predice “¿Consume usted huevos?” utilizando un solo predictor
Modelo 2: Predice “¿Consume usted huevos?” utilizando un conjunto de predictores, incluyendo variables relacionadas con la ubicación, demografía, vivienda, acceso a tecnología, hábitos de consumo y factores socioeconómicos.
Resid. Df (Grados de Libertad Residual): Los grados de libertad residual para cada modelo.
Resid. Dev (Desviación Residual): Una medida de qué tan bien el modelo se ajusta a los datos. Valores más bajos indican un mejor ajuste.
Df Deviance: El cambio en los grados de libertad de un modelo a otro.
Pr(>Chi): El valor p asociado con el cambio en la desviación. Se utiliza para evaluar si los predictores adicionales en el Modelo 2 mejoran significativamente el ajuste en comparación con el Modelo 1.
Para el Modelo 1 (un solo predictor), hay 382 grados de libertad y una desviación residual de 296.77.
Para el Modelo 2 (múltiples predictores), hay 272 grados de libertad, una desviación residual de 0.00 y una mejora significativa en el ajuste en comparación con el Modelo 1 (valor p < 2.2e-16).
‘’: Altamente significativo (valor p < 0.001) ’’: Significativo (valor p < 0.01) ’’: Marginalmente significativo (valor p < 0.05) ‘.’: Significativo en el límite (valor p < 0.1) ’ ’: No significativo (valor p >= 0.1)
En conclusión, el conjunto de predictores en el Modelo 2 mejora significativamente la capacidad de predecir el consumo de huevos en comparación con un modelo con un solo predictor.
SEGUNDO ANOVA
Consumo de huevo
anova(gfit3, test = "Chisq")## Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
## Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
## Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
## Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
## Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
## Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
## Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
## Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
## Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
## Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
## Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
## Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
## Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
## Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
## Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
## Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
## Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
## Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
## Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
## Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
## Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
## Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
## Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
## Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
## Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
## Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
## Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred## Warning: glm.fit: algorithm did not converge## Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
## Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred## Warning: glm.fit: algorithm did not converge## Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred## Warning: glm.fit: algorithm did not converge## Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred
## Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred## Analysis of Deviance Table
##
## Model: binomial, link: logit
##
## Response: ¿Consume usted huevos?
##
## Terms added sequentially (first to last)
##
##
## Df
## NULL
## UBICACIÓN 24
## `¿En qué cantón de la provincia de Chimborazo vive?` 0
## `Edad del encuestado` 1
## Sexo 1
## `Estado civil` 4
## `Cantidad de personas que conforman el hogar` 1
## `¿Qué parentesco tiene con el representante del hogar?` 5
## `¿Cuál es el tipo de vivienda?` 6
## `El material predominante de las paredes exteriores de la vivienda es de:` 4
## `EL material predominante del piso de la vivienda es de:` 4
## `¿Cuántos cuartos de baño con ducha de uso exclusivo tiene este hogar?` 3
## `Tipo de servicio higiénico con que cuenta este hogar` 5
## `¿Cuenta con servicio de internet?` 1
## `¿Cuenta con computadora de escritorio?` 1
## `¿Cuenta con computadora portátil?` 1
## `¿Cuántos celulares activados tiene en este hogar?` 4
## `¿Tiene su hogar servicio de teléfono convencional?` 1
## `¿Tiene en su hogar cocina con horno?` 1
## `¿Tiene refrigeradora?` 1
## `¿Tiene lavadora?` 1
## `¿Tiene equipo de sonido?` 1
## `¿Cuántas TV a color tienen en este hogar?` 3
## `¿Cuántos vehículos de uso exclusivo tiene este hogar?` 3
## `¿Alguien en el hogar compra vestimenta en centros comerciales?` 1
## `¿En el hogar alguien ha usado internet en los últimos 6 meses?` 1
## `¿En el hogar alguien utiliza correo electrónico que no es del trabajo?` 1
## `¿En el hogar alguien está registrado en una red social?` 1
## `Exceptuando los libros de texto o manuales de estudio y lecturas de trabajo¿Alguien del hogar ha leído algún libro completo en los últimos 3 meses?` 1
## `¿Cuál es el nivel de instrucción del Jefe del hogar?` 7
## `¿Alguien en el hogar está afiliado o cubierto por el seguro del IESS (general, voluntario o campesino) y/o seguro del ISSFA o ISSPOL?` 1
## `¿Alguien en el hogar tiene seguro de salud privada con hospitalización, seguro de salud privada sin hospitalización, seguro internacional,seguros municipales y de Consejos Provinciales y/o seguro de vida?` 1
## `¿Cuál es la ocupación del Representante del hogar?` 11
## NSE_2023_CUALI 4
## `¿Para qué sector trabaja?` 2
## `¿Consume usted arroz?` 1
## `¿Consume usted queso industrial (CON MARCA)` 1
## `¿Consume frutas?` 1
## Deviance
## NULL
## UBICACIÓN 71.08
## `¿En qué cantón de la provincia de Chimborazo vive?` 0.00
## `Edad del encuestado` 7.23
## Sexo 2.59
## `Estado civil` 2.94
## `Cantidad de personas que conforman el hogar` 0.87
## `¿Qué parentesco tiene con el representante del hogar?` 11.98
## `¿Cuál es el tipo de vivienda?` 6.56
## `El material predominante de las paredes exteriores de la vivienda es de:` 7.61
## `EL material predominante del piso de la vivienda es de:` 0.37
## `¿Cuántos cuartos de baño con ducha de uso exclusivo tiene este hogar?` 5.12
## `Tipo de servicio higiénico con que cuenta este hogar` 9.33
## `¿Cuenta con servicio de internet?` 11.98
## `¿Cuenta con computadora de escritorio?` 2.62
## `¿Cuenta con computadora portátil?` 0.20
## `¿Cuántos celulares activados tiene en este hogar?` 7.25
## `¿Tiene su hogar servicio de teléfono convencional?` 0.74
## `¿Tiene en su hogar cocina con horno?` 1.33
## `¿Tiene refrigeradora?` 0.07
## `¿Tiene lavadora?` 0.99
## `¿Tiene equipo de sonido?` 0.01
## `¿Cuántas TV a color tienen en este hogar?` 2.25
## `¿Cuántos vehículos de uso exclusivo tiene este hogar?` 3.04
## `¿Alguien en el hogar compra vestimenta en centros comerciales?` 0.22
## `¿En el hogar alguien ha usado internet en los últimos 6 meses?` 7.93
## `¿En el hogar alguien utiliza correo electrónico que no es del trabajo?` 3.53
## `¿En el hogar alguien está registrado en una red social?` 0.90
## `Exceptuando los libros de texto o manuales de estudio y lecturas de trabajo¿Alguien del hogar ha leído algún libro completo en los últimos 3 meses?` 0.80
## `¿Cuál es el nivel de instrucción del Jefe del hogar?` 10.27
## `¿Alguien en el hogar está afiliado o cubierto por el seguro del IESS (general, voluntario o campesino) y/o seguro del ISSFA o ISSPOL?` 6.98
## `¿Alguien en el hogar tiene seguro de salud privada con hospitalización, seguro de salud privada sin hospitalización, seguro internacional,seguros municipales y de Consejos Provinciales y/o seguro de vida?` 1.38
## `¿Cuál es la ocupación del Representante del hogar?` 0.00
## NSE_2023_CUALI 72.09
## `¿Para qué sector trabaja?` 1946.36
## `¿Consume usted arroz?` 0.00
## `¿Consume usted queso industrial (CON MARCA)` 0.00
## `¿Consume frutas?` 1946.36
## Resid. Df
## NULL 382
## UBICACIÓN 358
## `¿En qué cantón de la provincia de Chimborazo vive?` 358
## `Edad del encuestado` 357
## Sexo 356
## `Estado civil` 352
## `Cantidad de personas que conforman el hogar` 351
## `¿Qué parentesco tiene con el representante del hogar?` 346
## `¿Cuál es el tipo de vivienda?` 340
## `El material predominante de las paredes exteriores de la vivienda es de:` 336
## `EL material predominante del piso de la vivienda es de:` 332
## `¿Cuántos cuartos de baño con ducha de uso exclusivo tiene este hogar?` 329
## `Tipo de servicio higiénico con que cuenta este hogar` 324
## `¿Cuenta con servicio de internet?` 323
## `¿Cuenta con computadora de escritorio?` 322
## `¿Cuenta con computadora portátil?` 321
## `¿Cuántos celulares activados tiene en este hogar?` 317
## `¿Tiene su hogar servicio de teléfono convencional?` 316
## `¿Tiene en su hogar cocina con horno?` 315
## `¿Tiene refrigeradora?` 314
## `¿Tiene lavadora?` 313
## `¿Tiene equipo de sonido?` 312
## `¿Cuántas TV a color tienen en este hogar?` 309
## `¿Cuántos vehículos de uso exclusivo tiene este hogar?` 306
## `¿Alguien en el hogar compra vestimenta en centros comerciales?` 305
## `¿En el hogar alguien ha usado internet en los últimos 6 meses?` 304
## `¿En el hogar alguien utiliza correo electrónico que no es del trabajo?` 303
## `¿En el hogar alguien está registrado en una red social?` 302
## `Exceptuando los libros de texto o manuales de estudio y lecturas de trabajo¿Alguien del hogar ha leído algún libro completo en los últimos 3 meses?` 301
## `¿Cuál es el nivel de instrucción del Jefe del hogar?` 294
## `¿Alguien en el hogar está afiliado o cubierto por el seguro del IESS (general, voluntario o campesino) y/o seguro del ISSFA o ISSPOL?` 293
## `¿Alguien en el hogar tiene seguro de salud privada con hospitalización, seguro de salud privada sin hospitalización, seguro internacional,seguros municipales y de Consejos Provinciales y/o seguro de vida?` 292
## `¿Cuál es la ocupación del Representante del hogar?` 281
## NSE_2023_CUALI 277
## `¿Para qué sector trabaja?` 275
## `¿Consume usted arroz?` 274
## `¿Consume usted queso industrial (CON MARCA)` 273
## `¿Consume frutas?` 272
## Resid. Dev
## NULL 296.77
## UBICACIÓN 225.69
## `¿En qué cantón de la provincia de Chimborazo vive?` 225.69
## `Edad del encuestado` 218.46
## Sexo 215.86
## `Estado civil` 212.93
## `Cantidad de personas que conforman el hogar` 212.06
## `¿Qué parentesco tiene con el representante del hogar?` 200.08
## `¿Cuál es el tipo de vivienda?` 193.51
## `El material predominante de las paredes exteriores de la vivienda es de:` 185.90
## `EL material predominante del piso de la vivienda es de:` 185.53
## `¿Cuántos cuartos de baño con ducha de uso exclusivo tiene este hogar?` 180.42
## `Tipo de servicio higiénico con que cuenta este hogar` 171.09
## `¿Cuenta con servicio de internet?` 159.10
## `¿Cuenta con computadora de escritorio?` 156.48
## `¿Cuenta con computadora portátil?` 156.28
## `¿Cuántos celulares activados tiene en este hogar?` 149.04
## `¿Tiene su hogar servicio de teléfono convencional?` 148.29
## `¿Tiene en su hogar cocina con horno?` 146.97
## `¿Tiene refrigeradora?` 146.90
## `¿Tiene lavadora?` 145.91
## `¿Tiene equipo de sonido?` 145.90
## `¿Cuántas TV a color tienen en este hogar?` 143.65
## `¿Cuántos vehículos de uso exclusivo tiene este hogar?` 140.61
## `¿Alguien en el hogar compra vestimenta en centros comerciales?` 140.39
## `¿En el hogar alguien ha usado internet en los últimos 6 meses?` 132.46
## `¿En el hogar alguien utiliza correo electrónico que no es del trabajo?` 128.93
## `¿En el hogar alguien está registrado en una red social?` 128.04
## `Exceptuando los libros de texto o manuales de estudio y lecturas de trabajo¿Alguien del hogar ha leído algún libro completo en los últimos 3 meses?` 127.23
## `¿Cuál es el nivel de instrucción del Jefe del hogar?` 116.97
## `¿Alguien en el hogar está afiliado o cubierto por el seguro del IESS (general, voluntario o campesino) y/o seguro del ISSFA o ISSPOL?` 109.98
## `¿Alguien en el hogar tiene seguro de salud privada con hospitalización, seguro de salud privada sin hospitalización, seguro internacional,seguros municipales y de Consejos Provinciales y/o seguro de vida?` 108.61
## `¿Cuál es la ocupación del Representante del hogar?` 2018.44
## NSE_2023_CUALI 1946.36
## `¿Para qué sector trabaja?` 0.00
## `¿Consume usted arroz?` 1658.01
## `¿Consume usted queso industrial (CON MARCA)` 1946.36
## `¿Consume frutas?` 0.00
## Pr(>Chi)
## NULL
## UBICACIÓN 1.497e-06
## `¿En qué cantón de la provincia de Chimborazo vive?`
## `Edad del encuestado` 0.0071792
## Sexo 0.1072274
## `Estado civil` 0.5679467
## `Cantidad de personas que conforman el hogar` 0.3522739
## `¿Qué parentesco tiene con el representante del hogar?` 0.0349988
## `¿Cuál es el tipo de vivienda?` 0.3632100
## `El material predominante de las paredes exteriores de la vivienda es de:` 0.1070125
## `EL material predominante del piso de la vivienda es de:` 0.9848739
## `¿Cuántos cuartos de baño con ducha de uso exclusivo tiene este hogar?` 0.1633846
## `Tipo de servicio higiénico con que cuenta este hogar` 0.0966240
## `¿Cuenta con servicio de internet?` 0.0005364
## `¿Cuenta con computadora de escritorio?` 0.1053394
## `¿Cuenta con computadora portátil?` 0.6568938
## `¿Cuántos celulares activados tiene en este hogar?` 0.1234739
## `¿Tiene su hogar servicio de teléfono convencional?` 0.3888840
## `¿Tiene en su hogar cocina con horno?` 0.2494557
## `¿Tiene refrigeradora?` 0.7967669
## `¿Tiene lavadora?` 0.3197317
## `¿Tiene equipo de sonido?` 0.9066424
## `¿Cuántas TV a color tienen en este hogar?` 0.5229684
## `¿Cuántos vehículos de uso exclusivo tiene este hogar?` 0.3858273
## `¿Alguien en el hogar compra vestimenta en centros comerciales?` 0.6396242
## `¿En el hogar alguien ha usado internet en los últimos 6 meses?` 0.0048517
## `¿En el hogar alguien utiliza correo electrónico que no es del trabajo?` 0.0603784
## `¿En el hogar alguien está registrado en una red social?` 0.3431070
## `Exceptuando los libros de texto o manuales de estudio y lecturas de trabajo¿Alguien del hogar ha leído algún libro completo en los últimos 3 meses?` 0.3705163
## `¿Cuál es el nivel de instrucción del Jefe del hogar?` 0.1738813
## `¿Alguien en el hogar está afiliado o cubierto por el seguro del IESS (general, voluntario o campesino) y/o seguro del ISSFA o ISSPOL?` 0.0082310
## `¿Alguien en el hogar tiene seguro de salud privada con hospitalización, seguro de salud privada sin hospitalización, seguro internacional,seguros municipales y de Consejos Provinciales y/o seguro de vida?` 0.2405060
## `¿Cuál es la ocupación del Representante del hogar?` 1.0000000
## NSE_2023_CUALI 8.225e-15
## `¿Para qué sector trabaja?` < 2.2e-16
## `¿Consume usted arroz?` 1.0000000
## `¿Consume usted queso industrial (CON MARCA)` 1.0000000
## `¿Consume frutas?` < 2.2e-16
##
## NULL
## UBICACIÓN ***
## `¿En qué cantón de la provincia de Chimborazo vive?`
## `Edad del encuestado` **
## Sexo
## `Estado civil`
## `Cantidad de personas que conforman el hogar`
## `¿Qué parentesco tiene con el representante del hogar?` *
## `¿Cuál es el tipo de vivienda?`
## `El material predominante de las paredes exteriores de la vivienda es de:`
## `EL material predominante del piso de la vivienda es de:`
## `¿Cuántos cuartos de baño con ducha de uso exclusivo tiene este hogar?`
## `Tipo de servicio higiénico con que cuenta este hogar` .
## `¿Cuenta con servicio de internet?` ***
## `¿Cuenta con computadora de escritorio?`
## `¿Cuenta con computadora portátil?`
## `¿Cuántos celulares activados tiene en este hogar?`
## `¿Tiene su hogar servicio de teléfono convencional?`
## `¿Tiene en su hogar cocina con horno?`
## `¿Tiene refrigeradora?`
## `¿Tiene lavadora?`
## `¿Tiene equipo de sonido?`
## `¿Cuántas TV a color tienen en este hogar?`
## `¿Cuántos vehículos de uso exclusivo tiene este hogar?`
## `¿Alguien en el hogar compra vestimenta en centros comerciales?`
## `¿En el hogar alguien ha usado internet en los últimos 6 meses?` **
## `¿En el hogar alguien utiliza correo electrónico que no es del trabajo?` .
## `¿En el hogar alguien está registrado en una red social?`
## `Exceptuando los libros de texto o manuales de estudio y lecturas de trabajo¿Alguien del hogar ha leído algún libro completo en los últimos 3 meses?`
## `¿Cuál es el nivel de instrucción del Jefe del hogar?`
## `¿Alguien en el hogar está afiliado o cubierto por el seguro del IESS (general, voluntario o campesino) y/o seguro del ISSFA o ISSPOL?` **
## `¿Alguien en el hogar tiene seguro de salud privada con hospitalización, seguro de salud privada sin hospitalización, seguro internacional,seguros municipales y de Consejos Provinciales y/o seguro de vida?`
## `¿Cuál es la ocupación del Representante del hogar?`
## NSE_2023_CUALI ***
## `¿Para qué sector trabaja?` ***
## `¿Consume usted arroz?`
## `¿Consume usted queso industrial (CON MARCA)`
## `¿Consume frutas?` ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1PREDICCIÓN DEL MODELO
Consumo de huevos
PREDICCIÓN 1
gfit13 <- glm(`¿Consume usted huevos?`~., data =Train.hogar, family = binomial)## Warning: glm.fit: algorithm did not converge## Warning: glm.fit: fitted probabilities numerically 0 or 1 occurredcbind(gfit3$coefficients, gfit13$coefficients)## Warning in cbind(gfit3$coefficients, gfit13$coefficients): number of rows of
## result is not a multiple of vector length (arg 2)## [,1]
## (Intercept) 1003.5357482
## UBICACIÓNYaruquíes -452.1716362
## UBICACIÓNColta -388.9915422
## UBICACIÓNGuano -214.9765273
## UBICACIÓNGuamote -295.2299533
## UBICACIÓNVeloz -493.1298575
## UBICACIÓNMaldonado -395.2810564
## UBICACIÓNCacha -303.5550843
## UBICACIÓNLizarzaburu -321.9854120
## UBICACIÓNSan Luis -73.4498109
## UBICACIÓNPenipe -388.7004568
## UBICACIÓNSan Juan -39.4309344
## UBICACIÓNPallatanga 329.5591920
## UBICACIÓNPunín -136.3148594
## UBICACIÓNPungalá -202.5801712
## UBICACIÓNChambo 210.0926437
## UBICACIÓNAlausí -305.3353137
## UBICACIÓNLicán -111.0802608
## UBICACIÓNCumanda -290.9002096
## UBICACIÓNLicto -538.3475940
## UBICACIÓNCalpi -197.8608295
## UBICACIÓNChunchi 86.2127240
## UBICACIÓNQuimiag -273.5903094
## UBICACIÓNCubijíes -115.1969102
## UBICACIÓNFlores -16.4481511
## `¿En qué cantón de la provincia de Chimborazo vive?`Chambo NA
## `¿En qué cantón de la provincia de Chimborazo vive?`Chunchi NA
## `¿En qué cantón de la provincia de Chimborazo vive?`Colta NA
## `¿En qué cantón de la provincia de Chimborazo vive?`Cumanda NA
## `¿En qué cantón de la provincia de Chimborazo vive?`Guamote NA
## `¿En qué cantón de la provincia de Chimborazo vive?`Guano NA
## `¿En qué cantón de la provincia de Chimborazo vive?`Pallatanga NA
## `¿En qué cantón de la provincia de Chimborazo vive?`Penipe NA
## `¿En qué cantón de la provincia de Chimborazo vive?`Riobamba NA
## `Edad del encuestado` 62.3668434
## SexoFemenino 38.4398411
## `Estado civil`Divorciado -149.7468552
## `Estado civil`Otro -43.6481668
## `Estado civil`Soltero -42.5966240
## `Estado civil`Viudo 83.9708491
## `Cantidad de personas que conforman el hogar` -41.3158088
## `¿Qué parentesco tiene con el representante del hogar?`Otros parientes 43.7695147
## `¿Qué parentesco tiene con el representante del hogar?`Yo soy el representante del Hogar -58.0361455
## `¿Qué parentesco tiene con el representante del hogar?`Hijo o hija -19.6942381
## `¿Qué parentesco tiene con el representante del hogar?`Padres o suegros -527.6199136
## `¿Qué parentesco tiene con el representante del hogar?`Otros no parientes -7.7589230
## `¿Cuál es el tipo de vivienda?`Departamento en casa o edificio 3.8594464
## `¿Cuál es el tipo de vivienda?`Suite de lujo 260.4234902
## `¿Cuál es el tipo de vivienda?`Cuarto(s) en casa de inquilinato 17.3889919
## `¿Cuál es el tipo de vivienda?`Choza/Covacha/Otro 52.0628951
## `¿Cuál es el tipo de vivienda?`Rancho 264.4383419
## `¿Cuál es el tipo de vivienda?` Mediagua -1060.0353609
## `El material predominante de las paredes exteriores de la vivienda es de:`Caña no revestida/ Otros materiales -930.3073496
## `El material predominante de las paredes exteriores de la vivienda es de:`Ladrillo o bloque -416.0878243
## `El material predominante de las paredes exteriores de la vivienda es de:`Hormigón -396.6483835
## `El material predominante de las paredes exteriores de la vivienda es de:` Caña revestida o bahareque/Madera -371.5733717
## `EL material predominante del piso de la vivienda es de:` Cerámica, baldosa, vinil o marmetón 16.2709700
## `EL material predominante del piso de la vivienda es de:`Tabla sin tratar -22.3094482
## `EL material predominante del piso de la vivienda es de:`Tierra/ Caña/ Otros materiales -69.8291260
## `EL material predominante del piso de la vivienda es de:`Ladrillo o cemento 7.9340496
## `¿Cuántos cuartos de baño con ducha de uso exclusivo tiene este hogar?`Tiene 2 cuartos de baño exclusivos con ducha -165.1241378
## `¿Cuántos cuartos de baño con ducha de uso exclusivo tiene este hogar?`Tiene 1 cuarto de baño exclusivo con ducha -208.7125584
## `¿Cuántos cuartos de baño con ducha de uso exclusivo tiene este hogar?`Tiene 3 o más cuartos de baño exclusivos con ducha -205.0595029
## `Tipo de servicio higiénico con que cuenta este hogar`Conectado a pozo séptico 297.9756556
## `Tipo de servicio higiénico con que cuenta este hogar`Tiene 1 cuarto de baño exclusivo con ducha Letrina -130.7663725
## `Tipo de servicio higiénico con que cuenta este hogar`Conectado a pozo ciego -258.5343872
## `Tipo de servicio higiénico con que cuenta este hogar`Conectado a red pública de alcantarillado -535.3769743
## `Tipo de servicio higiénico con que cuenta este hogar`Sin Marca 3.3214819
## `¿Cuenta con servicio de internet?`Si -278.0482007
## `¿Cuenta con computadora de escritorio?`Si -73.9614669
## `¿Cuenta con computadora portátil?`Si 70.0687147
## `¿Cuántos celulares activados tiene en este hogar?`Tiene 3 celulares -35.3517296
## `¿Cuántos celulares activados tiene en este hogar?`Tiene 1 celular -23.3406684
## `¿Cuántos celulares activados tiene en este hogar?`Tiene 4 ó más celulares -1.9175757
## `¿Cuántos celulares activados tiene en este hogar?`Tiene 2 celulares -30.0204661
## `¿Tiene su hogar servicio de teléfono convencional?`Si -13.6601890
## `¿Tiene en su hogar cocina con horno?`Si 134.2464224
## `¿Tiene refrigeradora?`Si -41.1344946
## `¿Tiene lavadora?`Si 135.7025477
## `¿Tiene equipo de sonido?`Si 4.9209250
## `¿Cuántas TV a color tienen en este hogar?`Tiene 3 ó más TV a color 196.0391249
## `¿Cuántas TV a color tienen en este hogar?`Tiene 1 TV a color 211.8977903
## `¿Cuántas TV a color tienen en este hogar?`Tiene 2 TV a color 86.4226243
## `¿Cuántos vehículos de uso exclusivo tiene este hogar?`Tiene 3 ó más vehículo exclusivo 17.2919895
## `¿Cuántos vehículos de uso exclusivo tiene este hogar?`Tiene 2 vehículo exclusivo -121.5428443
## `¿Cuántos vehículos de uso exclusivo tiene este hogar?`Tiene un vehículo exclusivo -67.7655569
## `¿Alguien en el hogar compra vestimenta en centros comerciales?`Si 86.0913725
## `¿En el hogar alguien ha usado internet en los últimos 6 meses?`Si 183.2000218
## `¿En el hogar alguien utiliza correo electrónico que no es del trabajo?`Si 112.8879524
## `¿En el hogar alguien está registrado en una red social?`Si -77.0013565
## `Exceptuando los libros de texto o manuales de estudio y lecturas de trabajo¿Alguien del hogar ha leído algún libro completo en los últimos 3 meses?`Si -35.8415959
## `¿Cuál es el nivel de instrucción del Jefe del hogar?`Hasta 3 años de educación superior (incompleta) 31.9738600
## `¿Cuál es el nivel de instrucción del Jefe del hogar?`Postgrado -81.8970497
## `¿Cuál es el nivel de instrucción del Jefe del hogar?`Primaria completa -0.6146997
## `¿Cuál es el nivel de instrucción del Jefe del hogar?`Primaria incompleta -168.1956249
## `¿Cuál es el nivel de instrucción del Jefe del hogar?`Secundaria completa -14.8852161
## `¿Cuál es el nivel de instrucción del Jefe del hogar?`Secundaria incompleta 154.4027208
## `¿Cuál es el nivel de instrucción del Jefe del hogar?`Sin estudios 266.9028297
## `¿Alguien en el hogar está afiliado o cubierto por el seguro del IESS (general, voluntario o campesino) y/o seguro del ISSFA o ISSPOL?`Si 108.6163736
## `¿Alguien en el hogar tiene seguro de salud privada con hospitalización, seguro de salud privada sin hospitalización, seguro internacional,seguros municipales y de Consejos Provinciales y/o seguro de vida?`No -88.5286532
## `¿Cuál es la ocupación del Representante del hogar?` Empleados de oficina -307.2312164
## `¿Cuál es la ocupación del Representante del hogar?`Fuerzas Armadas -276.2040019
## `¿Cuál es la ocupación del Representante del hogar?`Inactivos -162.6072618
## `¿Cuál es la ocupación del Representante del hogar?` Oficiales operarios y artesanos 197.1054854
## `¿Cuál es la ocupación del Representante del hogar?` Operadores de instalaciones y máquinas 24.6911913
## `¿Cuál es la ocupación del Representante del hogar?`Personal directivo de la Administración Pública y de empresas -177.6500865
## `¿Cuál es la ocupación del Representante del hogar?` Profesionales científicos e intelectuales -260.2877387
## `¿Cuál es la ocupación del Representante del hogar?`Técnicos y profesionales de nivel medio -331.0137539
## `¿Cuál es la ocupación del Representante del hogar?`Trabajador calificados agropecuarios y pesqueros -149.8220426
## `¿Cuál es la ocupación del Representante del hogar?`Trabajador de los servicios y comerciantes -148.9473831
## `¿Cuál es la ocupación del Representante del hogar?`Trabajadores no calificados -241.4548760
## NSE_2023_CUALIBajo 145.0582791
## NSE_2023_CUALIMedio Alto -4.0648739
## NSE_2023_CUALIMedio Bajo 172.4233715
## NSE_2023_CUALIMedio Típico 26.2509135
## `¿Para qué sector trabaja?`Sector privado 180.0823779
## `¿Para qué sector trabaja?`Ninguno 83.2836853
## `¿Consume usted arroz?`Si -216.6355587
## `¿Consume usted queso industrial (CON MARCA)`Si -44.7173126
## `¿Consume frutas?`Si 211.2713576
## [,2]
## (Intercept) 87.9529638
## UBICACIÓNYaruquíes 47.5325175
## UBICACIÓNColta -90.8749936
## UBICACIÓNGuano 12.6378352
## UBICACIÓNGuamote 20.8629408
## UBICACIÓNVeloz -45.3992403
## UBICACIÓNMaldonado -79.0348111
## UBICACIÓNCacha 9.9943066
## UBICACIÓNLizarzaburu 70.0942804
## UBICACIÓNSan Luis -57.7237363
## UBICACIÓNPenipe 20.3053525
## UBICACIÓNSan Juan 22.7421091
## UBICACIÓNPallatanga 28.7819830
## UBICACIÓNPunín 14.1223994
## UBICACIÓNPungalá 122.2111018
## UBICACIÓNChambo 52.2506736
## UBICACIÓNAlausí -17.6966578
## UBICACIÓNLicán -80.8288624
## UBICACIÓNCumanda -18.4770927
## UBICACIÓNLicto 86.5879646
## UBICACIÓNCalpi 16.3267852
## UBICACIÓNChunchi 13.9667748
## UBICACIÓNQuimiag -54.6976504
## UBICACIÓNCubijíes NA
## UBICACIÓNFlores NA
## `¿En qué cantón de la provincia de Chimborazo vive?`Chambo NA
## `¿En qué cantón de la provincia de Chimborazo vive?`Chunchi NA
## `¿En qué cantón de la provincia de Chimborazo vive?`Colta NA
## `¿En qué cantón de la provincia de Chimborazo vive?`Cumanda NA
## `¿En qué cantón de la provincia de Chimborazo vive?`Guamote NA
## `¿En qué cantón de la provincia de Chimborazo vive?`Guano NA
## `¿En qué cantón de la provincia de Chimborazo vive?`Pallatanga NA
## `¿En qué cantón de la provincia de Chimborazo vive?`Penipe 29.7709046
## `¿En qué cantón de la provincia de Chimborazo vive?`Riobamba -19.3824229
## `Edad del encuestado` -18.6952475
## SexoFemenino 86.1936757
## `Estado civil`Divorciado -4.6766223
## `Estado civil`Otro 48.8163235
## `Estado civil`Soltero -15.5632579
## `Estado civil`Viudo 20.8180237
## `Cantidad de personas que conforman el hogar` -18.9561345
## `¿Qué parentesco tiene con el representante del hogar?`Otros parientes -74.6090024
## `¿Qué parentesco tiene con el representante del hogar?`Yo soy el representante del Hogar 5.3378418
## `¿Qué parentesco tiene con el representante del hogar?`Hijo o hija -52.1897870
## `¿Qué parentesco tiene con el representante del hogar?`Padres o suegros 14.6972956
## `¿Qué parentesco tiene con el representante del hogar?`Otros no parientes 8.6787128
## `¿Cuál es el tipo de vivienda?`Departamento en casa o edificio 36.7392789
## `¿Cuál es el tipo de vivienda?`Suite de lujo 168.3033982
## `¿Cuál es el tipo de vivienda?`Cuarto(s) en casa de inquilinato -72.6263269
## `¿Cuál es el tipo de vivienda?`Choza/Covacha/Otro -42.3039750
## `¿Cuál es el tipo de vivienda?`Rancho 10.0389602
## `¿Cuál es el tipo de vivienda?` Mediagua 11.2876817
## `El material predominante de las paredes exteriores de la vivienda es de:`Caña no revestida/ Otros materiales 17.2961411
## `El material predominante de las paredes exteriores de la vivienda es de:`Ladrillo o bloque 4.0989692
## `El material predominante de las paredes exteriores de la vivienda es de:`Hormigón 5.2574198
## `El material predominante de las paredes exteriores de la vivienda es de:` Caña revestida o bahareque/Madera 60.7310301
## `EL material predominante del piso de la vivienda es de:` Cerámica, baldosa, vinil o marmetón -1.0260604
## `EL material predominante del piso de la vivienda es de:`Tabla sin tratar 36.5426707
## `EL material predominante del piso de la vivienda es de:`Tierra/ Caña/ Otros materiales 29.0491131
## `EL material predominante del piso de la vivienda es de:`Ladrillo o cemento 17.2945560
## `¿Cuántos cuartos de baño con ducha de uso exclusivo tiene este hogar?`Tiene 2 cuartos de baño exclusivos con ducha 98.2523908
## `¿Cuántos cuartos de baño con ducha de uso exclusivo tiene este hogar?`Tiene 1 cuarto de baño exclusivo con ducha -81.4329966
## `¿Cuántos cuartos de baño con ducha de uso exclusivo tiene este hogar?`Tiene 3 o más cuartos de baño exclusivos con ducha -71.7919952
## `Tipo de servicio higiénico con que cuenta este hogar`Conectado a pozo séptico -7.0158421
## `Tipo de servicio higiénico con que cuenta este hogar`Tiene 1 cuarto de baño exclusivo con ducha Letrina -96.1818342
## `Tipo de servicio higiénico con que cuenta este hogar`Conectado a pozo ciego -48.5434568
## `Tipo de servicio higiénico con que cuenta este hogar`Conectado a red pública de alcantarillado -64.9408343
## `Tipo de servicio higiénico con que cuenta este hogar`Sin Marca -8.9686753
## `¿Cuenta con servicio de internet?`Si 0.6216823
## `¿Cuenta con computadora de escritorio?`Si 32.2532849
## `¿Cuenta con computadora portátil?`Si 31.4722458
## `¿Cuántos celulares activados tiene en este hogar?`Tiene 3 celulares 15.9008256
## `¿Cuántos celulares activados tiene en este hogar?`Tiene 1 celular 5.8474507
## `¿Cuántos celulares activados tiene en este hogar?`Tiene 4 ó más celulares 20.5188681
## `¿Cuántos celulares activados tiene en este hogar?`Tiene 2 celulares -30.0427874
## `¿Tiene su hogar servicio de teléfono convencional?`Si 17.0187849
## `¿Tiene en su hogar cocina con horno?`Si -3.3350471
## `¿Tiene refrigeradora?`Si 22.2318611
## `¿Tiene lavadora?`Si 52.3376123
## `¿Tiene equipo de sonido?`Si -8.6726963
## `¿Cuántas TV a color tienen en este hogar?`Tiene 3 ó más TV a color 10.8261879
## `¿Cuántas TV a color tienen en este hogar?`Tiene 1 TV a color -20.3053157
## `¿Cuántas TV a color tienen en este hogar?`Tiene 2 TV a color -12.0263015
## `¿Cuántos vehículos de uso exclusivo tiene este hogar?`Tiene 3 ó más vehículo exclusivo 23.5550526
## `¿Cuántos vehículos de uso exclusivo tiene este hogar?`Tiene 2 vehículo exclusivo 20.0733232
## `¿Cuántos vehículos de uso exclusivo tiene este hogar?`Tiene un vehículo exclusivo 9.4285126
## `¿Alguien en el hogar compra vestimenta en centros comerciales?`Si -22.9843714
## `¿En el hogar alguien ha usado internet en los últimos 6 meses?`Si -22.0078898
## `¿En el hogar alguien utiliza correo electrónico que no es del trabajo?`Si 20.0904580
## `¿En el hogar alguien está registrado en una red social?`Si -31.2954453
## `Exceptuando los libros de texto o manuales de estudio y lecturas de trabajo¿Alguien del hogar ha leído algún libro completo en los últimos 3 meses?`Si 82.5005040
## `¿Cuál es el nivel de instrucción del Jefe del hogar?`Hasta 3 años de educación superior (incompleta) 14.8526339
## `¿Cuál es el nivel de instrucción del Jefe del hogar?`Postgrado 60.5508124
## `¿Cuál es el nivel de instrucción del Jefe del hogar?`Primaria completa 89.1188339
## `¿Cuál es el nivel de instrucción del Jefe del hogar?`Primaria incompleta 78.1020001
## `¿Cuál es el nivel de instrucción del Jefe del hogar?`Secundaria completa 4.6293272
## `¿Cuál es el nivel de instrucción del Jefe del hogar?`Secundaria incompleta 16.1409698
## `¿Cuál es el nivel de instrucción del Jefe del hogar?`Sin estudios -54.0857286
## `¿Alguien en el hogar está afiliado o cubierto por el seguro del IESS (general, voluntario o campesino) y/o seguro del ISSFA o ISSPOL?`Si -64.0413203
## `¿Alguien en el hogar tiene seguro de salud privada con hospitalización, seguro de salud privada sin hospitalización, seguro internacional,seguros municipales y de Consejos Provinciales y/o seguro de vida?`No -85.2407042
## `¿Cuál es la ocupación del Representante del hogar?` Empleados de oficina 17.6759560
## `¿Cuál es la ocupación del Representante del hogar?`Fuerzas Armadas -38.4250796
## `¿Cuál es la ocupación del Representante del hogar?`Inactivos -47.7689925
## `¿Cuál es la ocupación del Representante del hogar?` Oficiales operarios y artesanos -105.2009344
## `¿Cuál es la ocupación del Representante del hogar?` Operadores de instalaciones y máquinas -56.9378520
## `¿Cuál es la ocupación del Representante del hogar?`Personal directivo de la Administración Pública y de empresas -33.7263695
## `¿Cuál es la ocupación del Representante del hogar?` Profesionales científicos e intelectuales 7.1654629
## `¿Cuál es la ocupación del Representante del hogar?`Técnicos y profesionales de nivel medio -24.7943704
## `¿Cuál es la ocupación del Representante del hogar?`Trabajador calificados agropecuarios y pesqueros -204.5213795
## `¿Cuál es la ocupación del Representante del hogar?`Trabajador de los servicios y comerciantes -37.8824637
## `¿Cuál es la ocupación del Representante del hogar?`Trabajadores no calificados -121.0767340
## NSE_2023_CUALIBajo -74.1843112
## NSE_2023_CUALIMedio Alto 20.7076548
## NSE_2023_CUALIMedio Bajo -7.0310411
## NSE_2023_CUALIMedio Típico -36.7297557
## `¿Para qué sector trabaja?`Sector privado -7.9846763
## `¿Para qué sector trabaja?`Ninguno 99.3975603
## `¿Consume usted arroz?`Si 87.9529638
## `¿Consume usted queso industrial (CON MARCA)`Si 47.5325175
## `¿Consume frutas?`Si -90.8749936CURVA ROC
library(pROC)## Type 'citation("pROC")' for a citation.##
## Attaching package: 'pROC'## The following objects are masked from 'package:stats':
##
## cov, smooth, vartest_prob3 <- predict(gfit3, newdata = val.hogar, type = "response")## Warning in predict.lm(object, newdata, se.fit, scale = 1, type = if (type == :
## prediction from rank-deficient fit; attr(*, "non-estim") has doubtful casestest_roc3 <- roc(val.hogar$`¿Consume usted huevos?` ~ test_prob3, plot = TRUE, print.auc = TRUE)## Setting levels: control = No, case = Si## Setting direction: controls < cases
Imterpretación
Se utiliza cuando se tiene datos binarios en el eje x se tiene la especificación mide los valores negativos en el eje y se mide la sensibilidad como el comportamiento de los valores positivos y el AUC es de 1 que significa que es bueno y el modelo de predición logística de predicción es bueno. Esta siempre entre 0 y 1 se acerca a 1 se estabiliza de forma horizontal.
MAQUINA DE VECTOR SOPORTE (SVM)
Las máquinas de soporte vectorial (Support Vector Machines o SVM, por sus siglas en inglés) son un tipo de algoritmo de aprendizaje supervisado utilizado para clasificación y regresión. Están fundamentadas en la teoría del aprendizaje estadístico y la optimización convexa. La idea principal detrás de las SVM es encontrar el hiperplano que mejor separa las clases en un espacio de características de alta dimensión.
Otro modelo de clasificación binaria es el conocido como Support Vector Machine (Máquina de vestor de soporte). Los paquetes que se ajustan los modelos de SVM es e1071 y Liblinear. Ajustamos al modelo para nuestros datos de Kernel radial que es el que se usa por defecto.
# install.packages("e1071")
library(e1071)
fitsvm3 <- svm(`¿Consume usted huevos?`~., data = Train.hogar)## Warning in svm.default(x, y, scale = scale, ..., na.action = na.action):
## Variable(s) '`¿Qué parentesco tiene con el representante del hogar?`Yo soy el
## representante del Hogar' and '`¿Qué parentesco tiene con el representante del
## hogar?`Yerno o nuera' constant. Cannot scale data.summary(fitsvm3)##
## Call:
## svm(formula = `¿Consume usted huevos?` ~ ., data = Train.hogar)
##
##
## Parameters:
## SVM-Type: C-classification
## SVM-Kernel: radial
## cost: 1
##
## Number of Support Vectors: 102
##
## ( 65 37 )
##
##
## Number of Classes: 2
##
## Levels:
## No Sifitsvm3##
## Call:
## svm(formula = `¿Consume usted huevos?` ~ ., data = Train.hogar)
##
##
## Parameters:
## SVM-Type: C-classification
## SVM-Kernel: radial
## cost: 1
##
## Number of Support Vectors: 102PREDECCIÓN DE VECTOR DE SOPORTE (SVM) RADIAL
A continuación predecimos los valores de la respuesta y calculamos:
# Ajustar un modelo SVM con kernel Radial
library(caret)## Loading required package: latticepredictedSVM2 <- predict(fitsvm3, val.hogar)
matrizSVM12 <- confusionMatrix(val.hogar$`¿Consume usted huevos?`, predictedSVM2)
matrizSVM12## Confusion Matrix and Statistics
##
## Reference
## Prediction No Si
## No 0 13
## Si 0 102
##
## Accuracy : 0.887
## 95% CI : (0.8145, 0.9384)
## No Information Rate : 1
## P-Value [Acc > NIR] : 1.0000000
##
## Kappa : 0
##
## Mcnemar's Test P-Value : 0.0008741
##
## Sensitivity : NA
## Specificity : 0.887
## Pos Pred Value : NA
## Neg Pred Value : NA
## Prevalence : 0.000
## Detection Rate : 0.000
## Detection Prevalence : 0.113
## Balanced Accuracy : NA
##
## 'Positive' Class : No
## Interpretación
En resumen, el modelo parece tener un problema significativo ya que no ha realizado ninguna predicción para la clase “No” y se puede decir que el modelo es bastante malo
SEGUNDA PREDICCIÓN KERNEL POLINÓMICO
# Ajustar un modelo SVM con kernel polinómico
fitsvm23 <- svm(`¿Consume usted huevos?` ~ ., data = Train.hogar, kernel = "polynomial")## Warning in svm.default(x, y, scale = scale, ..., na.action = na.action):
## Variable(s) '`¿Qué parentesco tiene con el representante del hogar?`Yo soy el
## representante del Hogar' and '`¿Qué parentesco tiene con el representante del
## hogar?`Yerno o nuera' constant. Cannot scale data.summary(fitsvm23)##
## Call:
## svm(formula = `¿Consume usted huevos?` ~ ., data = Train.hogar,
## kernel = "polynomial")
##
##
## Parameters:
## SVM-Type: C-classification
## SVM-Kernel: polynomial
## cost: 1
## degree: 3
## coef.0: 0
##
## Number of Support Vectors: 137
##
## ( 100 37 )
##
##
## Number of Classes: 2
##
## Levels:
## No SipredictedSVM23 <- predict(fitsvm23, newdata = val.hogar)
matrizSVM23 <- confusionMatrix(val.hogar$`¿Consume usted huevos?`, predictedSVM23)
print(matrizSVM23)## Confusion Matrix and Statistics
##
## Reference
## Prediction No Si
## No 0 13
## Si 0 102
##
## Accuracy : 0.887
## 95% CI : (0.8145, 0.9384)
## No Information Rate : 1
## P-Value [Acc > NIR] : 1.0000000
##
## Kappa : 0
##
## Mcnemar's Test P-Value : 0.0008741
##
## Sensitivity : NA
## Specificity : 0.887
## Pos Pred Value : NA
## Neg Pred Value : NA
## Prevalence : 0.000
## Detection Rate : 0.000
## Detection Prevalence : 0.113
## Balanced Accuracy : NA
##
## 'Positive' Class : No
## Interpretación
En resumen, el modelo parece tener un problema significativo ya que no ha realizado ninguna predicción para la clase “No” y se puede decir que el modelo es bastante malo
TERCERA PREDICCIÓN KERNEL LINEAL
library(e1071)
# Ajustar un modelo SVM con kernel lineal
fitsvm33 <- svm(`¿Consume usted huevos?` ~ ., data = Train.hogar, kernel = "linear")## Warning in svm.default(x, y, scale = scale, ..., na.action = na.action):
## Variable(s) '`¿Qué parentesco tiene con el representante del hogar?`Yo soy el
## representante del Hogar' and '`¿Qué parentesco tiene con el representante del
## hogar?`Yerno o nuera' constant. Cannot scale data.summary(fitsvm33)##
## Call:
## svm(formula = `¿Consume usted huevos?` ~ ., data = Train.hogar,
## kernel = "linear")
##
##
## Parameters:
## SVM-Type: C-classification
## SVM-Kernel: linear
## cost: 1
##
## Number of Support Vectors: 87
##
## ( 58 29 )
##
##
## Number of Classes: 2
##
## Levels:
## No Silibrary(caret)
predictedSVM33 <- predict(fitsvm33, newdata = val.hogar)
matrizSVM33 <- confusionMatrix(val.hogar$`¿Consume usted huevos?`, predictedSVM33)
print(matrizSVM33)## Confusion Matrix and Statistics
##
## Reference
## Prediction No Si
## No 4 9
## Si 15 87
##
## Accuracy : 0.7913
## 95% CI : (0.7056, 0.8615)
## No Information Rate : 0.8348
## P-Value [Acc > NIR] : 0.9131
##
## Kappa : 0.1337
##
## Mcnemar's Test P-Value : 0.3074
##
## Sensitivity : 0.21053
## Specificity : 0.90625
## Pos Pred Value : 0.30769
## Neg Pred Value : 0.85294
## Prevalence : 0.16522
## Detection Rate : 0.03478
## Detection Prevalence : 0.11304
## Balanced Accuracy : 0.55839
##
## 'Positive' Class : No
## Interpretación
En resumen, el modelo tiene una exactitud razonable, pero la sensibilidad para la clase “No” es baja, lo que indica que el modelo tiene dificultades para identificar correctamente los casos negativos de la clase “No”.
CUARTA PREDICCIÓN KERNEL SIGMOIDAL
# Ajustar un modelo SVM con kernel sigmoidal
fitsvm43 <- svm(`¿Consume usted huevos?` ~ ., data = Train.hogar, kernel = "sigmoid")## Warning in svm.default(x, y, scale = scale, ..., na.action = na.action):
## Variable(s) '`¿Qué parentesco tiene con el representante del hogar?`Yo soy el
## representante del Hogar' and '`¿Qué parentesco tiene con el representante del
## hogar?`Yerno o nuera' constant. Cannot scale data.summary(fitsvm43)##
## Call:
## svm(formula = `¿Consume usted huevos?` ~ ., data = Train.hogar,
## kernel = "sigmoid")
##
##
## Parameters:
## SVM-Type: C-classification
## SVM-Kernel: sigmoid
## cost: 1
## coef.0: 0
##
## Number of Support Vectors: 96
##
## ( 59 37 )
##
##
## Number of Classes: 2
##
## Levels:
## No SipredictedSVM43 <- predict(fitsvm43, newdata = val.hogar)
matrizSVM43 <- confusionMatrix(val.hogar$`¿Consume usted huevos?`, predictedSVM43)
print(matrizSVM43)## Confusion Matrix and Statistics
##
## Reference
## Prediction No Si
## No 0 13
## Si 0 102
##
## Accuracy : 0.887
## 95% CI : (0.8145, 0.9384)
## No Information Rate : 1
## P-Value [Acc > NIR] : 1.0000000
##
## Kappa : 0
##
## Mcnemar's Test P-Value : 0.0008741
##
## Sensitivity : NA
## Specificity : 0.887
## Pos Pred Value : NA
## Neg Pred Value : NA
## Prevalence : 0.000
## Detection Rate : 0.000
## Detection Prevalence : 0.113
## Balanced Accuracy : NA
##
## 'Positive' Class : No
## Interpretación
En resumen, el modelo tiene una alta exactitud debido a que predice correctamente la clase “No”, pero no predice la clase “Si” en absoluto el modelo es bastante malo para predecir.
COMPARACIÓN
CONSUMO DE HUEVOS
Accuracy <- c(matrizSVM12$overall[1], matrizSVM23$overall[1], matrizSVM43$overall[1], matrizSVM33$overall[1])
names(Accuracy)= c("SVM-RADIAL", "SVM-POLYNOMIAL", "SVM-SIGMODAL", "SVM-LINEAL")
Accuracy## SVM-RADIAL SVM-POLYNOMIAL SVM-SIGMODAL SVM-LINEAL
## 0.8869565 0.8869565 0.8869565 0.7913043Interpretación
SVM-RADIAL: La precisión para SVM con núcleo radial es del 88.70%.
SVM-POLYNOMIAL: La precisión para SVM con núcleo polinómico es del 88.70%.
SVM-SIGMOIDAL: La precisión para SVM con núcleo sigmoide es del 88.70%.
SVM-LINEAL: La precisión para SVM con núcleo lineal es del 79.13%.
Estos valores indican la proporción de predicciones correctas para cada tipo de núcleo. En este caso, parece haber una consistencia en la precisión entre los diferentes tipos de núcleos, con una disminución en la precisión para SVM con núcleo lineal.
ÁRBOL DE DECISION
CONSUMO DE HUEVOS
library(tree)
X1 <- Train.hogar$`¿Consume usted huevos?`
X2 <- Train.hogar$UBICACIÓN
X3 <- Train.hogar$`¿En qué cantón de la provincia de Chimborazo vive?`
X4 <- Train.hogar$`Edad del encuestado`
X5 <- Train.hogar$Sexo
X6 <- Train.hogar$`Estado civil`
X7 <- Train.hogar$`Cantidad de personas que conforman el hogar`
X8 <- Train.hogar$`¿Qué parentesco tiene con el representante del hogar?`
X9 <- Train.hogar$`¿Cuál es el tipo de vivienda?`
X10 <- Train.hogar$`El material predominante de las paredes exteriores de la vivienda es de:`
X11 <- Train.hogar$`EL material predominante del piso de la vivienda es de:`
X12 <- Train.hogar$`¿Cuántos cuartos de baño con ducha de uso exclusivo tiene este hogar?`
X13 <- Train.hogar$`Tipo de servicio higiénico con que cuenta este hogar`
X14 <- Train.hogar$`¿Cuenta con servicio de internet?`
X15 <- Train.hogar$`¿Cuenta con computadora de escritorio?`
X16 <- Train.hogar$`¿Cuenta con computadora portátil?`
X17 <- Train.hogar$`¿Cuántos celulares activados tiene en este hogar?`
X18 <- Train.hogar$`¿Tiene su hogar servicio de teléfono convencional?`
X19 <- Train.hogar$`¿Tiene en su hogar cocina con horno?`
X20 <- Train.hogar$`¿Tiene refrigeradora?`
X21 <- Train.hogar$`¿Tiene lavadora?`
X22 <- Train.hogar$`¿Tiene equipo de sonido?`
X23 <- Train.hogar$`¿Cuántas TV a color tienen en este hogar?`
X24 <- Train.hogar$`¿Cuántos vehículos de uso exclusivo tiene este hogar?`
X25 <- Train.hogar$`¿Alguien en el hogar compra vestimenta en centros comerciales?`
X26 <- Train.hogar$`¿En el hogar alguien ha usado internet en los últimos 6 meses?`
X27 <- Train.hogar$`¿En el hogar alguien utiliza correo electrónico que no es del trabajo?`
X28 <- Train.hogar$`¿En el hogar alguien está registrado en una red social?`
X29 <- Train.hogar$`Exceptuando los libros de texto o manuales de estudio y lecturas de trabajo¿Alguien del hogar ha leído algún libro completo en los últimos 3 meses?`
X30 <- Train.hogar$`¿Cuál es el nivel de instrucción del Jefe del hogar?`
X31 <- Train.hogar$`¿Alguien en el hogar está afiliado o cubierto por el seguro del IESS (general, voluntario o campesino) y/o seguro del ISSFA o ISSPOL?`
X32 <- Train.hogar$`¿Alguien en el hogar tiene seguro de salud privada con hospitalización, seguro de salud privada sin hospitalización, seguro internacional,seguros municipales y de Consejos Provinciales y/o seguro de vida?`
X33 <- Train.hogar$`¿Cuál es la ocupación del Representante del hogar?`
X34 <- Train.hogar$NSE_2023_CUALI
X35 <- Train.hogar$`¿Para qué sector trabaja?`
X36 <- Train.hogar$`¿Consume usted arroz?`
X37 <- Train.hogar$`¿Consume usted queso industrial (CON MARCA)`
X38 <- Train.hogar$`¿Consume frutas?`
# Crear el árbol de decisión
treeA1 <- tree(X1 ~ X2 + X3 + X4 + X5 + X6 + X7 + X8 + X9 + X10 + X11 + X12 + X13 + X14 + X15 + X16 + X17 + X18 + X19 + X20 + X21 + X22 + X23 + X24 + X25 + X26 + X27 + X28 + X29 + X30 + X31 + X32 + X33 + X34 + X35 + X36 + X37 + X38, data = Train.hogar)
# Imprimir resumen del árbol
summary(treeA1)##
## Classification tree:
## tree(formula = X1 ~ X2 + X3 + X4 + X5 + X6 + X7 + X8 + X9 + X10 +
## X11 + X12 + X13 + X14 + X15 + X16 + X17 + X18 + X19 + X20 +
## X21 + X22 + X23 + X24 + X25 + X26 + X27 + X28 + X29 + X30 +
## X31 + X32 + X33 + X34 + X35 + X36 + X37 + X38, data = Train.hogar)
## Variables actually used in tree construction:
## [1] "X2" "X4" "X31" "X11" "X23" "X26" "X33" "X18" "X5" "X30" "X35"
## Number of terminal nodes: 17
## Residual mean deviance: 0.1928 = 48.4 / 251
## Misclassification error rate: 0.04104 = 11 / 268El modelo de árbol de clasificación ajustado que incluye múltiples variables independientes (X2, X4, X31, X11, X23, X26, X33, X18, X5, X30, X35). Lo que sugiere que está teniendo en cuenta múltiples características en la toma de decisiones. La desviación media residual es 0.1928, lo cual es una medida de cuán bien el modelo se ajusta a los datos. El modelo tiene un rendimiento sólido con una tasa de error de clasificación baja y un número moderado de nodos terminales.
Gráfica
plot(treeA1)
text(treeA1, pretty =1)## Warning in FUN(X[[i]], ...): abreviatura utilizada con caracteres no ASCII
## Warning in FUN(X[[i]], ...): abreviatura utilizada con caracteres no ASCII
## Warning in FUN(X[[i]], ...): abreviatura utilizada con caracteres no ASCII
## Warning in FUN(X[[i]], ...): abreviatura utilizada con caracteres no ASCII
## Warning in FUN(X[[i]], ...): abreviatura utilizada con caracteres no ASCII
## Warning in FUN(X[[i]], ...): abreviatura utilizada con caracteres no ASCII
## Warning in FUN(X[[i]], ...): abreviatura utilizada con caracteres no ASCII
## Warning in FUN(X[[i]], ...): abreviatura utilizada con caracteres no ASCII
## Warning in FUN(X[[i]], ...): abreviatura utilizada con caracteres no ASCII
## Warning in FUN(X[[i]], ...): abreviatura utilizada con caracteres no ASCII
## Warning in FUN(X[[i]], ...): abreviatura utilizada con caracteres no ASCII
## Warning in FUN(X[[i]], ...): abreviatura utilizada con caracteres no ASCII
## Warning in FUN(X[[i]], ...): abreviatura utilizada con caracteres no ASCII
## Warning in FUN(X[[i]], ...): abreviatura utilizada con caracteres no ASCII
## Warning in FUN(X[[i]], ...): abreviatura utilizada con caracteres no ASCII
## Warning in FUN(X[[i]], ...): abreviatura utilizada con caracteres no ASCII
## Warning in FUN(X[[i]], ...): abreviatura utilizada con caracteres no ASCII
## Warning in FUN(X[[i]], ...): abreviatura utilizada con caracteres no ASCII
## Warning in FUN(X[[i]], ...): abreviatura utilizada con caracteres no ASCII
## Warning in FUN(X[[i]], ...): abreviatura utilizada con caracteres no ASCII
## Warning in FUN(X[[i]], ...): abreviatura utilizada con caracteres no ASCII
Predicción
library(caret)
X1_test <- Train.hogar$`¿Consume usted huevos?`
X2_test <- Train.hogar$UBICACIÓN
X3_test <- Train.hogar$`¿En qué cantón de la provincia de Chimborazo vive?`
X4_test <- Train.hogar$`Edad del encuestado`
X5_test <- Train.hogar$Sexo
X6_test <- Train.hogar$`Estado civil`
X7_test <- Train.hogar$`Cantidad de personas que conforman el hogar`
X8_test <- Train.hogar$`¿Qué parentesco tiene con el representante del hogar?`
X9_test <- Train.hogar$`¿Cuál es el tipo de vivienda?`
X10_test <- Train.hogar$`El material predominante de las paredes exteriores de la vivienda es de:`
X11_test <- Train.hogar$`EL material predominante del piso de la vivienda es de:`
X12_test <- Train.hogar$`¿Cuántos cuartos de baño con ducha de uso exclusivo tiene este hogar?`
X13_test <- Train.hogar$`Tipo de servicio higiénico con que cuenta este hogar`
X14_test <- Train.hogar$`¿Cuenta con servicio de internet?`
X15_test <- Train.hogar$`¿Cuenta con computadora de escritorio?`
X16_test <- Train.hogar$`¿Cuenta con computadora portátil?`
X17_test <- Train.hogar$`¿Cuántos celulares activados tiene en este hogar?`
X18_test <- Train.hogar$`¿Tiene su hogar servicio de teléfono convencional?`
X19_test <- Train.hogar$`¿Tiene en su hogar cocina con horno?`
X20_test <- Train.hogar$`¿Tiene refrigeradora?`
X21_test <- Train.hogar$`¿Tiene lavadora?`
X22_test <- Train.hogar$`¿Tiene equipo de sonido?`
X23_test <- Train.hogar$`¿Cuántas TV a color tienen en este hogar?`
X24_test <- Train.hogar$`¿Cuántos vehículos de uso exclusivo tiene este hogar?`
X25_test <- Train.hogar$`¿Alguien en el hogar compra vestimenta en centros comerciales?`
X26_test <- Train.hogar$`¿En el hogar alguien ha usado internet en los últimos 6 meses?`
X27_test <- Train.hogar$`¿En el hogar alguien utiliza correo electrónico que no es del trabajo?`
X28_test <- Train.hogar$`¿En el hogar alguien está registrado en una red social?`
X29_test <- Train.hogar$`Exceptuando los libros de texto o manuales de estudio y lecturas de trabajo¿Alguien del hogar ha leído algún libro completo en los últimos 3 meses?`
X30_test <- Train.hogar$`¿Cuál es el nivel de instrucción del Jefe del hogar?`
X31_test <- Train.hogar$`¿Alguien en el hogar está afiliado o cubierto por el seguro del IESS (general, voluntario o campesino) y/o seguro del ISSFA o ISSPOL?`
X32_test <- Train.hogar$`¿Alguien en el hogar tiene seguro de salud privada con hospitalización, seguro de salud privada sin hospitalización, seguro internacional,seguros municipales y de Consejos Provinciales y/o seguro de vida?`
X33_test <- Train.hogar$`¿Cuál es la ocupación del Representante del hogar?`
X34_test <- Train.hogar$NSE_2023_CUALI
X35_test <- Train.hogar$`¿Para qué sector trabaja?`
X36_test <- Train.hogar$`¿Consume usted arroz?`
X37_test <- Train.hogar$`¿Consume usted queso industrial (CON MARCA)`
X38_test <- Train.hogar$`¿Consume frutas?`
predictions <- predict(treeA1, newdata = data.frame(X2 = X2_test, X3 = X3_test, X4 = X4_test, X5 = X5_test, X6 = X6_test, X7 = X7_test, X8 = X8_test, X9 = X9_test, X10 = X10_test, X11 = X11_test, X12 = X12_test, X13 = X13_test, X14 = X14_test, X15 = X15_test, X16 = X16_test, X17 = X17_test, X18 = X18_test, X19 = X19_test, X20 = X20_test, X21 = X21_test, X22 = X22_test, X23 = X23_test, X24 = X24_test, X25 = X25_test, X26 = X26_test, X27 = X27_test, X28 = X28_test, X29 = X29_test, X30 = X30_test, X31 = X31_test, X32 = X32_test, X33 = X33_test, X34 = X34_test, X35 = X35_test, X36 = X36_test, X37 = X37_test, X38 = X38_test), type = "class")
confusionMatrix(predictions, X1_test)## Confusion Matrix and Statistics
##
## Reference
## Prediction No Si
## No 33 7
## Si 4 224
##
## Accuracy : 0.959
## 95% CI : (0.9277, 0.9793)
## No Information Rate : 0.8619
## P-Value [Acc > NIR] : 1.276e-07
##
## Kappa : 0.8332
##
## Mcnemar's Test P-Value : 0.5465
##
## Sensitivity : 0.8919
## Specificity : 0.9697
## Pos Pred Value : 0.8250
## Neg Pred Value : 0.9825
## Prevalence : 0.1381
## Detection Rate : 0.1231
## Detection Prevalence : 0.1493
## Balanced Accuracy : 0.9308
##
## 'Positive' Class : No
## Interpretación
El modelo muestra un rendimiento muy sólido con una precisión del 95.9%, una sensibilidad del 89.19% y una especificidad del 96.97%. El valor de Kappa sugiere una concordancia sustancial entre las predicciones del modelo y las clases reales. El modelo parece ser capaz de identificar de manera efectiva ambas clases (“Si” y “No”). En general, el modelo es bastante bueno para predecir, y la tasa de error de clasificación es baja.
CONCLUSIÓN
En este caso en la comparación de los dos modelos, tanto las máquinas de vector de soporte y árboles de desición de las variables de CONSUMO DE HUEVOS, indican un modelo adecuado en los árboles de desición ya que es más precisos e indican una mejor predicción, mientras que el vector de soporte sus predicciónson es bastante malo y se puede decir que no son adecuados para predecir.