future.apply

DESCRIPTION

@@ -46,7 +46,8 @@ Imports:
     rstan (>= 2.26.0),
     rstantools (>= 2.1.1),
     R6,
-    assertthat
+    assertthat,
+    future.apply
 LinkingTo: 
     BH (>= 1.66.0),
     Rcpp (>= 0.12.0),

R/analyse.R

@@ -152,9 +152,9 @@ analyse <- function(imputations, fun = ancova, delta = NULL, ...) {
         )
     }
 
-    results <- lapply(
-        imputations$imputations,
-        function(x, ...) {
+    results <- future.apply::future_lapply(
+        X = imputations$imputations,
+        FUN = function(x, ...) {
             dat2 <- extract_imputed_df(x, imputations$data, delta)
             fun(dat2, ...)
         },

vignettes/quickstart.Rmd

@@ -411,3 +411,59 @@ poolObj <- pool(
     alternative = "two.sided"
 )
 ```
+
+
+```{r, eval = FALSE}
+library(future)
+library(lubridate)
+some_fun <- function(data) {
+    Sys.sleep(.5)
+    as_datetime(2)
+    NULL
+}
+
+my_fun <- function(
+    data,
+    vars,
+    visits = NULL,
+    weights = c("counterfactual", "equal", "proportional_em", "proportional")
+) {
+    some_fun(data)    
+    ancova(data, vars, visits, weights)
+}
+plan("multisession", workers = 4)
+system.time({
+anaObj_par <- analyse(
+    imputeObj,
+    my_fun,
+    vars = set_vars(
+        subjid = "PATIENT",
+        outcome = "CHANGE",
+        visit = "VISIT",
+        group = "THERAPY",
+        covariates = c("BASVAL")
+    )
+)
+anaObj_par
+})
+
+plan("sequential")
+system.time({
+anaObj_seq <- analyse(
+    imputeObj,
+    my_fun,
+    vars = set_vars(
+        subjid = "PATIENT",
+        outcome = "CHANGE",
+        visit = "VISIT",
+        group = "THERAPY",
+        covariates = c("BASVAL")
+    )
+)
+anaObj_seq
+})
+
+identical(anaObj_par, anaObj_seq)
+
+```
+