Those are not the points, but the log likelihood at the points you specified.
You don’t need to use the sigmoid transformation theta, you can more simply define the Bernoulli with the logit values via the argument logit_p=mu.
If I had to guess the -inf probability comes from an overflow of the sigmoid to zero or one, together with an observation that’s the opposite value. Using the logit_p parametrization will avoid this.