This model aims to seek posterior of beta0 and sigma. In this model, f1 and f2 are known (deterministic) functions.
My question is about z. The “z=pyro.sample()” sentence will assume z has a prior which is a multi Gaussian, and this model will seek posterior of z, too. But in my case, z_loc and z_cov are deterministic given beta0, and I don’t want any inference on z; it should be a forward-pass instead. How can I deal with it?

if z is low dimensional, e.g. 1-4 dimension, you can use gauss-hermite quadrature

if it’s higher dimensional it becomes much harder and you probably just need to explicitly do inference over z, which effectively does the integral via sampling

Unfortunately, the dimension of z is rather high (e.g. over 100). I am trying to by pass the inference over z because the computation time is much longer than I desire. May be I should firstly try to reduce the dimension of z.