in heavy_likelihood function, if I use t(t()) two times, there is error. I will continue checking it.