2 Regresión Logística
2 Regresión Logística
2 Regresión Logística
Regresin logstica
2. REGRESIN LOGSTICA 2.1. Introduccin 2.2. Modelo logit 2.3. Construccin de modelos 2.4. Interpretacin de Odd ratios 2.5. Ejercicio en R.
2.4. Interpretacin
Bajo el modelo de regresin logstica:
1 g ( x 1) g ( x)
La pendiente del coeficiente representa el cambio en el logit correspondiente al cambio en una unidad en x. Por lo tanto la interpretacin no es directa, la forma de interpretar los coeficientes de la regresin dependen del tipo de datos.
Odd ratios
Para interpretar los resultados se utilizan odd ratios o razn What??? de probabilidad:
Para los parmetros de la regresin logstica:
e
( EE ) e
1 =exp 1
(en Excel)
=exp EE1
(en Excel)
Interpretacin: Si el intervalo de confianza incluye el valor 1, entonces el efecto de la variable es irrelevante, si el intervalo no incluye a 1, entonces la variable tiene significado biolgico.
Ejemplo:
Ejemplo: si Y corresponde a la presencia o ausencia de parasitismo de nidos en un sitio, 2 sugiere que la probabilidad entonces de parasitismo es 2 veces mayor que la probabilidad de ausencia de parasitismo en el sitio.
determinado Y=1 con respecto al grupo de comparacin. Ejemplo: si y representa la presencia o ausencia de parasitismo en tres diferentes tipos de hbitat, 0.5 sugiere que la entonces un valor de probabilidad de presencia de parasitismo es menor en un tipo de hbitat con respecto al tipo de hbitat de comparacin.
Ejemplo con R
mydata<-read.csv("c://curso//ejerciciosr//binary.csv") attach(mydata) names(mydata) summary(gre) sd(gre) summary(gpa) sd(gpa) table(rank) table(admit) table(rank,admit) mylogit<- glm(admit~gre+gpa+as.factor(rank), family=binomial(link="logit"), na.action=na.pass) summary(mylogit) confint(mylogit)
exp(mylogit$coefficients) exp(confint(mylogit)) rank <- c(1,2,3,4) gre <- c(mean(mydata$gre)) gpa <- c(mean(mydata$gpa)) newdata1 <- data.frame(gre,gpa,rank) newdata1 newdata1$rankP <-predict(mylogit,newdata=newdata1,type="response") newdata1