Potential bug in Predictive.mvnormal()

[sergioventurini]

Hello,

I'm struggling with the code for the Predictive.mvnormal() function. In particular, I was wondering why in the last if() statement that computes the gamma_contrib scalar you used twice 'prod(vapply(seq_along(d) ...))', where seq_along(d) is always equal to 1. Can you please check it?

Many thanks in advance for your response.