ENH: Check input to binary models

I don't think that's what we want. It still works for continuous in interval [0,1].
My guess is that the right check is `not (endog < 0 & endog > 1).any()`

but I haven't gone through the code. S-shaped continuous variables, fractions seem to work well.