Thanks Ricardo, that makes sense.
To modify my code to get the log-likelihood of each observations, I’m guessing I need to change my_loglike function to not sum the likelihoods, and change the LogLike otype to pt.dvector, to get a vector of log likelihood values?