replace for loop by vectorized operation

Author: NightlordTW

R/RcppExports.R

@@ -317,7 +317,7 @@ run_simulations_par <- function(nsim, n, muT, muR, SigmaT, SigmaR, lequi_tol, ue
 #' @param typey Integer vector indicating the classification of each endpoint, where `1` corresponds to a primary endpoint and `2` corresponds to a secondary endpoint.
 #' @param adseq Logical. If `TRUE`, applies sequential (hierarchical) testing.
 #' @param k Integer. Minimum number of endpoints required for equivalence.
-#' @param arm_seed arma::ivec. Random seed vector (one per simulation).
+#' @param arm_seed Integer vector. Random seed for each simulation.
 #'
 #' @details
 #' This function performs equivalence testing using either the Difference of Means (DOM) or Ratio of Means (ROM) approach.

man/run_simulations_2x2.Rd

@@ -123,14 +123,15 @@ arma::mat check_equivalence(const arma::uvec& typey,
 
   // If no primary endpoints exist, assume sequential test passes automatically
   bool sumpe = (num_primary == 0);
-  int sumtypey = 0;
+  int sumtypey = arma::accu(tbioq.cols(arma::find(typey == 1)));
 
+  //int sumtypey = 0;
   // Count number of primary endpoints that meet equivalence
-  for (size_t i = 0; i < typey.n_elem; i++) {
-    if (typey(i) == 1) {
-      sumtypey += tbioq(0, i);
-    }
-  }
+  //for (size_t i = 0; i < typey.n_elem; i++) {
+  //  if (typey(i) == 1) {
+  //    sumtypey += tbioq(0, i);
+  //  }
+  //}
 
   // Ensure all primary endpoints pass if sequential testing is enabled
   if (num_primary > 0) {
@@ -832,7 +833,7 @@ arma::mat run_simulations_par(const int nsim,
 //' @param typey Integer vector indicating the classification of each endpoint, where `1` corresponds to a primary endpoint and `2` corresponds to a secondary endpoint.
 //' @param adseq Logical. If `TRUE`, applies sequential (hierarchical) testing.
 //' @param k Integer. Minimum number of endpoints required for equivalence.
-//' @param arm_seed arma::ivec. Random seed vector (one per simulation).
+//' @param arm_seed Integer vector. Random seed for each simulation.
 //'
 //' @details
 //' This function performs equivalence testing using either the Difference of Means (DOM) or Ratio of Means (ROM) approach.
@@ -903,6 +904,7 @@ arma::mat run_simulations_2x2(const int nsim,
   return results.t(); // Transpose before returning
 }
 
+
 RCPP_MODULE(test)
 {
    function("rcpp_test_2x2_dom", &test_2x2_dom);

src/module.cpp

@@ -46,7 +46,7 @@ lequi_upper <- c(AUC = log(1.25), Cmax = log(1.25))
                  list_lequi.tol = list("T_vs_R" = lequi_lower),
                  list_uequi.tol = list("T_vs_R" = lequi_upper),
                  dtype = "2x2", ctype = "DOM", lognorm = FALSE, 
-                 adjust = "no", ncores = 1, nsim = 10000, seed = 1234))
+                 adjust = "no", ncores = 1, nsim = 5000, seed = 1234))
 ```
 The total sample size required is `r ss$response$n_total` subjects, which corresponds to the estimate obtained using the PASS software (n=37).