An example corresponding to bins=np.linspace(lo, hi,nbins+1) would be really helpful. Automated binning is great for exploratory visualization but you often want a consistent binning that is independent of data variance in a production environment.