Are you saying that
# n is the number of params
# point_est_of_params is some n-length array
# var_of_params is some n-length array
x = pm.Normal("x", mu=point_est_of_params, sigma=tt.sqrt(var_of_params), len=n)
won’t do what you want?
Are you saying that
# n is the number of params
# point_est_of_params is some n-length array
# var_of_params is some n-length array
x = pm.Normal("x", mu=point_est_of_params, sigma=tt.sqrt(var_of_params), len=n)
won’t do what you want?