Cannot fit bayesian model with rethinking::quap for small (1e-4) parameter values

Question

I\'m trying to calculate posterior for p using rethinking package version 2.13:

rethinking::quap(
  alist(
    n ~ dbinom(N, p),