Hi guys, recently i’ve learned about the powerful packages for bayesian inference in julia

i would like some tips on how to attack this problem, given this probabiliy distribution over the angles \theta:

where W and Z are defined as

finally the \Sigma_i are defined as

I have some observations of these sets of angles \vec \theta_1,\vec \theta_2 and i want to infer the posterior distribution of the matrix J and the vectors h_{1/2}. but i have no clue on how to leverage the marvelous tools in the julia ecosystem like Turing.jl etc… since the integral in Z will unavoidably appear in the `logpdf`

.

if this question is considered off-topic i will delete it. Thank you in advance