Hi,
I’m not sure I understand exactly what you want to do, but:
- If you want to compute the standard deviation during sampling, you need to use theano operators, as the variables are tensors, not scalars.
- If you want to compute it after sampling, you can just use the trace and do numpy operations on it, as the trace contains scalars.
Hope this helps 
PS: Any reason why you’re not using the built-in Beta-Binomial distribution?