I am curious how to calculate probability from stan parameter estimate. In one dimensional it is easy - which(theta>=0.19&theta<0.21) /length(theta) gives probability that theta is within some interval. Any advice how to do it effectively when theta has multiple dimensions (around 40). My ultimate goal is to calculate probability that f(theta) is within hyper-rectangle.
Thanks for the advice.