MODELO ANDES

datos1 <- ANDES


head(datos1)

library(tidyverse)
library(ISLR)
default1 <-ANDES$Xanthomonas
balance1<- ANDES$Northing + ANDES$Easting

modelo_logistico1 <- glm(default1 ~ balance1, data = datos1, family = "binomial")

ggplot(data = datos1, aes(x = balance1, y = default1)) +
  geom_point(aes(color = as.factor(default1)), shape = 1) + 
  stat_function(fun = function(x){predict(modelo_logistico1,
                                          newdata = data.frame(balance1 = x),
                                          type = "response")}) +
  theme_bw() +
  labs(title = "Regresión logística ANDES",
       y = "Probabilidad default") +
  theme(legend.position = "none")