dataset.jobs <- subset(dataset, !is.na(JOBSATIS))
dataset.jobs$is.satisfied <- dataset.jobs$JOBSATIS == 1 | dataset.jobs$JOBSATIS == 2
model2 <- glm(is.satisfied ~ DGRDG + NDGMEMG + HD03Y5 + NOCPRMG + EMSEC + WAPRSM + GENDER + MINRTY + MGRNAT + MGROTH + WASCSM + SATADV + SATBEN + SATCHAL + SATIND + SATLOC + SATRESP + SATSAL + SATSEC + SATSOC, data = dataset.jobs, family = "binomial")
library(MASS)
model2.aic <- stepAIC(model2)
summary(model2.aic)
library(InformationValue)
model.probs <- model2.aic$fitted.values
plotROC(dataset.jobs$is.satisfied == TRUE, model.probs)