I am having a distribution with rather small standard deviation, if I use pm.summary() the sd’s(hpd) will be rounded down to zero. How do i cope with this?
You can use
az.summary(idata, round_to="none") to display the raw numbers (defaults to 2 decimals to round results) – see here.
Hope this helps