@junpenglao Is there currently any support for dynamic shapes (i.e. using a tensorvariable as input to the shape argument)? I am trying to estimate the parameters of a Compound Poisson process with a Binomial distribution as the underlying distribution. The number of Binomial RVs follow a Poisson distribution, so I need a way to either dynamically set the shape argument or generate those RVs in a loop.