Correct ELBO

statwangz · statwangz · commit 1aede41b34ec · 2025-04-18T16:37:06.000+08:00
diff --git a/R/mfairELBO.R b/R/mfairELBO.R
@@ -37,7 +37,7 @@ getELBO <- function(Y, object, obs_indices = NULL) {
           ) -
           projSparse(
             tcrossprod(as.matrix(mu_sq), as.matrix(nu_sq)), obs_indices
-          ))
+          )) / 2
     } else {
       elbo2 <- -tau * sum(
         (Y - tcrossprod(as.matrix(mu), as.matrix(nu)))^2 +
