Introducción a la Regresión Logística
Introducción a la Regresión Logística
MARCO TEÓRICO
La regresión logística es una técnica de aprendizaje automático que proviene del campo
de la estadística, que requiere datos binarios y supone una distribución de Bernoulli.
Este es un método para problemas de clasificación, en los que se obtienen un valor binario
entre 0 y 1. Por ejemplo, un problema de clasificación es identificar si una operación dada
es exitosa o no. Asociándole una etiqueta “éxito” a unos registros y “fracaso” a otros.
Simplificando mucho es identificar si al realizar una afirmación sobre registro esta es cierta
o no (Rodríguez, 2018).
Esta regresión permite medir la relación entre las variables dependientes (lo que se desea
predecir) con una o más variables independientes o explicativas con las cuales se pretende
dar cuenta del comportamiento de la variable respuesta.
Este modelo se hace visible a través de funciones matemáticas, que expresan la relación
entre las variables dependientes e independientes, que en su forma mas simple incluye una
sola variable explicativa X1, a saber:
𝟏
𝑷(𝒁 = 𝟏|𝑿) =
𝟏 + 𝒆𝒙𝒑[−(𝜷𝟎 + 𝜷𝟏 𝑿𝟏 )]
𝟏
𝑷(𝒁 = 𝟏|𝑿) =
𝟏 + 𝒆𝒙𝒑[−(𝜷𝟎 + 𝜷𝟏 𝑿𝟏 + 𝜷𝟐 𝑿𝟐 + 𝜷𝟑 𝑿𝟑 + ⋯ + 𝜷𝒑 𝑿𝒑 )]
Regresión Logística
Elkin morales
Ricardo Ortiz
Es importante anotar que las variables explicativas del modelo (independientes), pueden
ser del tipo nominal, ordinal o continuo.
P(Z|X)
Este modelo es el mas usado dada su sencillez y la eficacia que tiene a la hora de interpretar
los resultados ademas de que constituye la representación más parsimoniosa del problema
que se quiere resolver.
Este modelo de regresión se usa para modelar la relación entre una variable dependiente
dicotómica (éxito o fracaso) y variables independientes. Las variables independientes
pueden ser numéricas o categóricas y su dicotomía se puede expresar como estructuralmente
insegura (Z=1), o segura (Z=0), aunque esta codificación es arbitraria e irrelevante.
Una vez estimados los parámetros del modelo, es necesario revisar la importancia de cada
variable para el modelo, determinando cual de ellas es la mas influyente y adquiere mayor
relevancia en el mismo. Estas determinaciones deben realizarse a través de revisión de p-
valores para un nivel de significancia prestablecido, y del análisis de los errores
encontrados.
Regresión Logística
Elkin morales
Ricardo Ortiz
Para lo anterior, existen procedimientos tales como el stepwise, con el cual pueden
constituirse modelos de varias combinaciones de las variables independientes e ir
identificando y excluyendo aquellas que nos son influyentes o relevantes dentro del modelo.
Aunque esta regresión permite trabajar con mayor número de categorías, el modelo
logístico multivariado resulta menos preciso a la hora de predecir un resultado, dado que
ya no habría solo dos opciones (caso regresión binaria) sino, por ejemplo 5 opciones como
en el caso de la clasificación de niveles arriba mencionada.
El p-valor es el nivel de significación más pequeño para el que, con los datos de la muestra,
habría que rechazar Ho. Valores superiores al nivel de significancia especificado llevarían
a la aceptación de Ho, mientras que valores inferiores al nivel de significancia implicarían
el rechazo de Ho.
La odds ratio (oportunidad relativa) es la probabilidad que un evento ocurra dividido por la
probabilidad de que el evento no ocurra. Es muy utilizada estadística para comparar la
frecuencia de exposición de diferentes factores de riesgo.
Regresión Logística
Elkin morales
Ricardo Ortiz
EJEMPLOS DE APLICACIÓN
EJEMPLO 1:
(Logistic Regression Model for Evaluating Soil Liquefaction Probability Using CPT Data)
Es una modelación logística basada en 396 casos históricos, teniendo en cuenta que el
caso de licuefacción en suelos granulares saturados sometidos a cargas sísmicas es uno
de los tópicos más importantes en cuanto a la estabilidad de las estructuras que se
construyen en zonas de alto potencial de licuefacción.
Varios autores (Seed et al, Robertson et al) han desarrollado procesos simplificados para
la determinación de potencial de licuación basados en pruebas SPT, pero dadas las
características de este último ensayo, se hace necesario el análisis de las muestras
obtenidas en laboratorio (por ejemplo granulometría) para poder implementar los
procedimientos de cálculo.
Con el ensayo de Cono CPT se tiene la ventaja de obtener parámetros únicamente sitio y de
forma continua.
2)normalmente distribuidas
Otros autores sugieren que, aunque las variables explicativas no sean linealmente
dependientes del modelo y no estén normalmente distribuidas, la transformación logística
de PL (x) se sigue cumpliendo 1) 2) y 3).
Se tomaron 396 datos de campo históricos de varios terremotos ocurridos en China, Taiwan
y EEUU , cuyo procesamiento incluyó:
a) Manipulación de la información
Regresión Logística
Elkin morales
Ricardo Ortiz
b) Revisión de probabilidades mediante diferentes transformaciones de la variable
explicativa.
c) Procedimiento de clasificación del suelo usando datos del ensayo CPT para arenas,
mezclas de arenas y mezclas de limos.
Las variables explicativas pueden ser dividida en dos: la primera corresponde la RELACION
DE TENSION CICLICA AL INICIO DE LA LICUEFACCIÓN (CRR) evaluada por medio del ensayo
de cono CPT (qn) y la segunda corresponde a perturbación sísmica inducida por los
terremotos que es representada como un equivalente de la relación de carga cíclica CSR
Donde CSR7.5 corresponde a la relación de carga cíclica para un terremoto de magnitud 7.5
Regresión Logística
Elkin morales
Ricardo Ortiz
√𝑞 se observa mas normalmente distribuida que q y CSR es mas log- normalmente
distribuido que normalmente distribuido.
Finalmente √𝑞 y ln(CSR) son las variables explicativas más apropiadas del modelo, donde:
Luego:
Regresión Logística
Elkin morales
Ricardo Ortiz
Finalmente, se dicotomiza la variable dependiente PL = 0.5 para dividir la condición de suelo
licuable o no licuable usando el procedimiento simplificado de Seed and Idriss (1971) cuyo
criterio de licuación se de para un factor de seguridad FS >1.
El análisis de regresión logística es hecho entre PL y FS, el cual muestra una forma sigmoide:
Regresión Logística
Elkin morales
Ricardo Ortiz
EJEMPLO 2:
Para este trabajo inicialmente se desarrolló un sistema de información geográfica SIG con
información específica para el salvador, denominado SIGSAL, que permitió la recopilación
e integración de diferentes bases de datos, mapas e información de la zona de estudio.
-Modelo digital del Terreno (MDT): Realizado a través del proceso de digitalización y
estructuración de bases de datos, georreferenciación, edición y validación.
-LITOLOGIA
-ELEVACIÓN
-PENDIENTE
-ORIENTACIÓN
-RUGOSIDAD
Regresión Logística
Elkin morales
Ricardo Ortiz
-USOS DE SUELO
Mediante el uso de la regresión logística para un nivel de significancia del 99%, se encontró
que la mayoría de los factores tienen un p-valor <0.01, lo cual indica que tienen un alto nivel
de significación, a excepción del factor USO DE SUELO, que obtuvo un p-valor>0.01.
Con este ultimo modelo se mejoró el nivel de significancia del modelo ajustado y se
redujeron de forma notable las variables de la ecuación, resultando tan solo dos factores
influyentes: RUGOSIDAD(R) Y LITOLOGÍA (S), ajustándose a la siguiente ecuación:
Donde S1, S3 y S4 son categorías del suelo (roca dura, suelo consolidado y suelo no
consolidado respectivamente).
Regresión Logística
Elkin morales
Ricardo Ortiz
EJEMPLO 3:
Regresión Logística
Elkin morales
Ricardo Ortiz
Modelo de regresión ordinal
Para validar cada uno de los modelos individuales, se calculan los estadísticos de Wald y el
de líneas paralelas. La prueba de Wald permite evaluar la significación estadística de una
variable explicativa sobre la variable dependiente, bajo la hipótesis:
H0: βj = 0, vs Ha: βj ≠ 0
Para j = 1, 2, 3, 4, 5, 6, 7.
Regresión Logística
Elkin morales
Ricardo Ortiz
DESARROLLO DE UN MODELO DE REGRESIÓN LOGÍSTICA
Para el presente trabajo se desarrolló un modelo de regresión logística para el cálculo La
resistencia no-drenada es la resistencia del suelo cuando se carga hasta la falla en
condiciones no-drenadas o sea cuando las cargas que producen la falla se aplican sobre la
masa de suelo a una velocidad superior a la del drenaje del suelo. El caso más común de
resistencia no-drenada, se presenta en los depósitos naturales de arcilla saturada cuando
éstos son cargados o descargados en forma relativamente rápida, comparada con la rata
en la cual puede ocurrir drenaje y/o consolidación.
Cuando se presenta esta condición se asume que hay un fenómeno de resistencia no-
drenada; el contenido de agua y el volumen de la arcilla permanecen constantes durante la
carga no-drenada y se generan presiones de poros en exceso. El comportamiento no-
drenado de arcillas saturadas se analiza en términos de esfuerzos totales y la evaluación
de las presiones de poros es innecesaria. Bajo esta situación se asume un método de
análisis Ø = 0 y la resistencia no-drenada Cu es igual al valor de cohesión en la envolvente
de Mohr-Coulomb para esfuerzos totales.
Modelo de análisis
(Cu)/σ=0.11+0.0037IP
Donde:
• z= profundidad (m)
• σ = ɣd*z
Los datos observados con los cuales se quiere validar y ajustar el modelo de Matteo et al
(2009) corresponden a 59 ensayos realizados para arcillas de Bogotá por la Universidad
Nacional, los cuales fueron suministrados por el ingeniero Rubén Tovar y se resumen en el
archivo llamado [Link]
Para los suelos con resistencia al corte no drenado alto Y=1 (buena condición) y para suelos
con resistencia al corte no drenado bajo Y = 0 (mala condición).
0-15 bajo
<15 medio.
Regresión Logística
Elkin morales
Ricardo Ortiz
Datos de campo
Regresión Logística
Elkin morales
Ricardo Ortiz
Datos categorizados
Regresión Logística
Elkin morales
Ricardo Ortiz
MODELACION EN R
Respecto al modelo, el logaritmo de odds de que Cu esta positivamente relacionado con los
valores de los esfuerzos verticales (SIGMA) (coeficiente parcial = 0.25760), siendo
significativa esta relación (p-value = 0.00739). También existe una relación significativa
positiva entre el logaritmo de odds de Cu mayor a 150kpa y los rangos de IP(bajo y medio)
(p-value = 0.13707).
2. Intervalos de confianza
2.5 % 97.5 %
(Intercept) - 19.0269291 -2.9745124
SIGMA 0.1207142 0.5149528
IPmedio -0.4653233 9.9804677
Regresión Logística
Elkin morales
Ricardo Ortiz
intervalos basados en el error estándar
2.5 % 97.5 %
(Intercept) -16.1187035 - 1.1572319
SIGMA 0.0691228 0.4460754
IPmedio -1.1660534 8.4937908
Regresión Logística
Elkin morales
Ricardo Ortiz
4. Caculo de diferencia de residuos,
grados de libertad y p-value
[Link] de confusión
Predicciones
Observaciones 0 1
0 12 2
1 2 43
Regresión Logística
Elkin morales
Ricardo Ortiz
Predicciones
Observaciones 0 1
0 VN FP
1 FN VP
MÉTRICAS
Regresión Logística
Elkin morales
Ricardo Ortiz
CONCLUSIONES
El modelo logístico creado para predecir la probabilidad de que un suelo tenga una
resistencia no drenada superior a 150kpa a partir del índice plasticidad y de los esfuerzos
verticales es en conjunto significativo acorde al (p-value = 2.2596719942e-11). El p-value
de una de las variables es significativo (SIGMA= 0.00739), mientras que de la otra variable
no es tan bueno IP = 0.1307). La ratio de error obtenido empleando las observaciones con
las que se ha entrenado el modelo muestra un porcentaje de verdaderos positivos altos, lo
cual nos indica que el modelo predice de forma acertada nuevos valores.
BIBIOGRAFÍA
- Sheng-Yao Lai; Wen-Jong Chang; and Ping-Sien Lin. Logistic Regression Model for
Evaluating Soil Liquefaction Probability Using CPT Data.. J. Geotech. Geoenviron.
Eng., 2006, 132(6): 694-704.
- María José García Rodríguez, Universidad de Alcalá, 2008. EVALUACIÓN DE LA
SUSCEPTIBILIDAD Y PELIGROSIDAD DE DESLIZAMIENTOS DE LADERAS PARA EL
SALVADOR.
- ZÚÑIGA Gutiérrez Martín, CUEVAS Sandoval Alfredo, SÁNCHEZ Calvo Mateo,
BARRAGÁN Trinidad Raziel. Análisis de la Amenaza de Riesgo por Deslizamiento de
Laderas en la Ciudad de Chilpancingo, Gro. Articulo “Foro de estudios sobre
guerrero” Mayo 2015-Abril 2016 Vol. 2 No.3 224-243
- LADY JOHANNA SAUZA RODRÍGUEZ, DETERMINACIÓN DEL ESTADO DE REDES
LOCALES DE ALCANTARILLADO Y SU NECESIDAD O NO DE SER SOMETIDAS A
RENOVACIÓN/REHABILITACIÓN TENIENDO EN CUENTA CCTV EN UN NÚMERO
LIMITADO DE TUBERÍAS. CASO DE ESTUDIO ZONA 1 DE LA EAAB – ESP.
UNIVERSIDAD DE LOS ANDES, BOGOTÁ D.C., ENERO DE 2020.
- Das, Braja M. Fundamentos de ingeniería geotécnica. Cuarta Edición. ISBN:978-607-
519-372-4
Regresión Logística
Elkin morales
Ricardo Ortiz