Translated more of UpdateGamma() (#11)

ocbe-uio · Sep 16, 2024 · 58bf5d3 · 58bf5d3
1 parent 8df46e8
commit 58bf5d3
Showing 1 changed file with 3 additions and 1 deletion.
diff --git a/src/UpdateGamma_cpp.cpp b/src/UpdateGamma_cpp.cpp
@@ -87,8 +87,10 @@ Rcpp::List UpdateGamma_cpp(
     );
     return out;
   } else {
-    // post.gamma <- rep(list(rep(0, p)), S)
     Rcpp::List post_gamma = Rcpp::List::create(S);
+    for (arma::uword g = 0; g < S; g++) {
+      post_gamma[g] = arma::zeros<arma::vec>(p);
+    }
 
     if (MRF_G) {
       for (arma::uword g = 0; g < S; g++) { // loop through subgroups