Thank you for the reply. Well, I tried to use the adapt_diag but it always keeps on giving me mass matrix error around the degree of freedom and scale when it enters data fusion part.
So, is there a way I could get the posterior where my parameters are around the range of (Ex: 0.00xxx, xx.x*10-x)?