Модель логистической регрессии здесь
glm.fit <- glm(activity~orbsize, data=data, family = "binomial")
summary(glm.fit)
Я пробую новые данные, и все работает нормально
prob <- predict(glm.fit, newdata, type = "response")
pred <- ifelse(prob > 0.5, "Nocturnal", "Diurnal")
Но как мне это построить?
xline <- seq(min(newdata$orbsize), max(data$orbsize), 5)
yline <- predict(glm.fit, newdata, type = "response")
lines(xline, yline, col = "red")
Вы должны использовать значения
xline
внутриnewdata
data.frame, чтобы делать прогнозыy
для этих значенийx
.