MODELOS SONSON

datos1 <- SONSON_1


head(datos1)

library(tidyverse)

package 㤼㸱tidyverse㤼㸲 was built under R version 3.6.3Registered S3 methods overwritten by 'dbplyr':
  method         from
  print.tbl_lazy     
  print.tbl_sql      
-- Attaching packages --------------------------------------- tidyverse 1.3.0 --
v ggplot2 3.3.2     v purrr   0.3.4
v tibble  3.0.3     v dplyr   1.0.2
v tidyr   1.1.2     v stringr 1.4.0
v readr   1.3.1     v forcats 0.5.0
package 㤼㸱ggplot2㤼㸲 was built under R version 3.6.3package 㤼㸱tibble㤼㸲 was built under R version 3.6.3package 㤼㸱tidyr㤼㸲 was built under R version 3.6.3package 㤼㸱readr㤼㸲 was built under R version 3.6.3package 㤼㸱purrr㤼㸲 was built under R version 3.6.3package 㤼㸱dplyr㤼㸲 was built under R version 3.6.3package 㤼㸱stringr㤼㸲 was built under R version 3.6.3package 㤼㸱forcats㤼㸲 was built under R version 3.6.3-- Conflicts ------------------------------------------ tidyverse_conflicts() --
x dplyr::filter() masks stats::filter()
x dplyr::lag()    masks stats::lag()

library(ISLR)

package 㤼㸱ISLR㤼㸲 was built under R version 3.6.3

default1 <- SONSON_1$Xanthomonas
balance1<- SONSON_1$Northing + SONSON_1$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 SONSON 1",
       y = "Probabilidad default") +
  theme(legend.position = "none")

datos2 <- SONSON_2


head(datos2)

default2 <- SONSON_2$Xanthomonas
balance2<- SONSON_2$Northing + SONSON_2$Easting

modelo_logistico2 <- glm(default2 ~ balance2, data = datos2, family = "binomial")

ggplot(data = datos2, aes(x = balance2, y = default2)) +
  geom_point(aes(color = as.factor(default2)), shape = 1) + 
  stat_function(fun = function(x){predict(modelo_logistico2,
                                          newdata = data.frame(balance2 = x),
                                          type = "response")}) +
  theme_bw() +
  labs(title = "Regresión logística SONSON 2",
       y = "Probabilidad default") +
  theme(legend.position = "none")

datos3 <- SONSON_3


head(datos3)

default3 <- SONSON_3$Xanthomonas
balance3<- SONSON_3$Northing + SONSON_3$Easting

modelo_logistico3 <- glm(default3 ~ balance3, data = datos3, family = "binomial")

ggplot(data = datos3, aes(x = balance3, y = default3)) +
  geom_point(aes(color = as.factor(default3)), shape = 1) + 
  stat_function(fun = function(x){predict(modelo_logistico3,
                                          newdata = data.frame(balance3 = x),
                                          type = "response")}) +
  theme_bw() +
  labs(title = "Regresión logística SONSON 3",
       y = "Probabilidad default") +
  theme(legend.position = "none")