Hi @junpenglao! Thanks for answering.
See my code below:
var_alpha = theano.shared(value = 1., borrow = False)
phi_e = pm.Uniform(‘phi_e’, lower = lb_phi_e, upper = ub_phi_e, shape = (ub_phi_e.size))
S1_phi_e = (phi_e/phi_e.sum()*np.log(phi_e/phi_e.sum())).sum()
pm.DensityDist(‘phi_e_potential’, logp = -var_alpha*S1_phi_e)
The error I receive with this approach is “TypeError: ‘TensorVariable’ object is not callable”