Does Julia have a function to calculate the density points where p% of the distribution is included?
Something like the scipy.stats norm.ppf function mentioned in this answer
Example: 2-sided 95% confidence interval:
> norm.ppf(1-(1-0.95)/2)
1.96
> norm.ppf(1-(1+0.95)/2)
-1.96