Closed sambrilleman closed 6 years ago
This line where the log likelihood for the random effects distribution is missing brackets around the term following the -0.5 multiplier. It should be...
ll_b <- -0.5 (c(determinant(Sigma, logarithm = TRUE)$modulus) + (b %% chol2inv(chol(Sigma)) %% b)[1] + length(b) log(2 * pi))
Need to check this updated on stan-dev branch too...
This line where the log likelihood for the random effects distribution is missing brackets around the term following the -0.5 multiplier. It should be...
ll_b <- -0.5 (c(determinant(Sigma, logarithm = TRUE)$modulus) + (b %% chol2inv(chol(Sigma)) %% b)[1] + length(b) log(2 * pi))
Need to check this updated on stan-dev branch too...