I’m going to write a power-law distribution ( y = x^-2) in which x comes from a uniform distribution in [0,100] interval. According to this topic I could write the following code:
import theano.tensor as tt
log_density = -alpha*tt.log(x)
(alpha is 2). but I can’t incorporate the limitation on x!. How can I do this?