Suponha que eu tenha os dados abaixo, aplique neles uma regressão logística e queira visualizar a função estimada que retorna a probabilidade.
#Simula os dados
nobs<-100
beta0=-10
beta1=0.006
saldo=runif(nobs,1300,2300)
p_x <- exp(beta0 + beta1 * saldo) / (1 + exp(beta0 + beta1 * saldo))
y <- rbinom(n=length(saldo), size=1, prob=p_x)
data <- data.frame(saldo, p_x, y)
#Regressão
default.glm = glm(y~saldo, data=data,family=binomial)
summary(default.glm)
O código abaixo não está plotando adequadamente a função y=probabilidade prevista e x= saldo.
plot(saldo,y)
lines(data$saldo, default.glm$fitted, type="l", col="red")