# NOT RUN {
data('fake_data')
data = fake_data
data$y = 1/(1+exp(-data$y))
data$y = ifelse(data$y > .5, 1, 0)
table(data$y)
fit5 = ConsReg(y~x1+x2+x3+x4, data = data,
family = 'binomial', penalty = 10000,
LOWER = -.5, UPPER = .2,
optimizer = 'gosolnp')
pr = predict(fit5, newdata = data[1:3,], type = 'probability')
pr
# }
Run the code above in your browser using DataLab