Fixed variance of some parameters

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?