explicit args

DESCRIPTION

@@ -53,13 +53,17 @@ Depends:
   R (>= 4.0.0)
 Imports:
   brms,
+  coda,
   dplyr,
   emmeans,
   MASS,
+  posterior,
+  purrr,
   rlang,
   stats,
   tibble,
   tidyr,
+  tidyselect,
   trialr,
   utils
 Suggests:

NAMESPACE

@@ -12,15 +12,27 @@ importFrom(brms,brmsformula)
 importFrom(brms,get_prior)
 importFrom(brms,prior)
 importFrom(brms,unstr)
+importFrom(coda,as.mcmc)
 importFrom(dplyr,left_join)
 importFrom(emmeans,emmeans)
+importFrom(posterior,as_draws_df)
+importFrom(posterior,mcse_mean)
+importFrom(posterior,mcse_median)
+importFrom(posterior,mcse_quantile)
+importFrom(posterior,mcse_sd)
+importFrom(purrr,map2_df)
+importFrom(purrr,map_dbl)
+importFrom(purrr,map_df)
 importFrom(rlang,warn)
 importFrom(stats,as.formula)
+importFrom(stats,median)
 importFrom(stats,model.matrix)
 importFrom(stats,rnorm)
 importFrom(stats,runif)
+importFrom(stats,sd)
 importFrom(tibble,tibble)
 importFrom(tidyr,expand_grid)
+importFrom(tidyselect,any_of)
 importFrom(trialr,rlkjcorr)
 importFrom(utils,capture.output)
 importFrom(utils,globalVariables)

R/brm_package.R

@@ -20,13 +20,18 @@
 #'     doi:10.1177/009286150804200402
 #' @family help
 #' @importFrom brms brm brmsformula get_prior prior unstr
+#' @importFrom coda as.mcmc
 #' @importFrom dplyr left_join
 #' @importFrom emmeans emmeans
 #' @importFrom MASS mvrnorm
+#' @importFrom posterior as_draws_df mcse_mean mcse_median mcse_quantile
+#'   mcse_sd
+#' @importFrom purrr map_dbl map_df map2_df
 #' @importFrom rlang warn
-#' @importFrom stats as.formula model.matrix rnorm runif
+#' @importFrom stats as.formula median model.matrix rnorm runif sd
 #' @importFrom tibble tibble
 #' @importFrom tidyr expand_grid
+#' @importFrom tidyselect any_of
 #' @importFrom trialr rlkjcorr
 #' @importFrom utils capture.output globalVariables
 NULL

tests/testthat/test-brm_marginal_draws.R

@@ -1,6 +1,10 @@
 test_that("brm_marginal_draws() on response", {
   set.seed(0L)
-  sim <- brm_simulate()
+  sim <- brm_simulate(
+    n_group = 2L,
+    n_patient = 100L,
+    n_time = 4L
+  )
   data <- sim$data
   data$group <- paste("treatment", data$group)
   data$time <- paste("visit", data$time)
@@ -89,7 +93,11 @@ test_that("brm_marginal_draws() on response", {
 
 test_that("brm_marginal_draws() on change", {
   set.seed(0L)
-  sim <- brm_simulate()
+  sim <- brm_simulate(
+    n_group = 2L,
+    n_patient = 100L,
+    n_time = 4L
+  )
   data <- sim$data
   data$group <- paste("treatment", data$group)
   data$time <- paste("visit", data$time)

tests/testthat/test-brm_marginal_probabilities.R

@@ -1,6 +1,10 @@
 test_that("brm_marginal_probabilities() on response", {
   set.seed(0L)
-  sim <- brm_simulate()
+  sim <- brm_simulate(
+    n_group = 2L,
+    n_patient = 100L,
+    n_time = 4L
+  )
   data <- sim$data
   data$group <- paste("treatment", data$group)
   data$time <- paste("visit", data$time)
@@ -63,7 +67,11 @@ test_that("brm_marginal_probabilities() on response", {
 
 test_that("brm_marginal_probabilities() on change and multiple probs", {
   set.seed(0L)
-  sim <- brm_simulate()
+  sim <- brm_simulate(
+    n_group = 2L,
+    n_patient = 100L,
+    n_time = 4L
+  )
   data <- sim$data
   data$group <- paste("treatment", data$group)
   data$time <- paste("visit", data$time)

tests/testthat/test-brm_marginal_summaries.R

@@ -1,6 +1,10 @@
 test_that("brm_marginal_summaries() on response", {
   set.seed(0L)
-  sim <- brm_simulate()
+  sim <- brm_simulate(
+    n_group = 2L,
+    n_patient = 100L,
+    n_time = 4L
+  )
   data <- sim$data
   data$group <- paste("treatment", data$group)
   data$time <- paste("visit", data$time)
@@ -121,7 +125,11 @@ test_that("brm_marginal_summaries() on response", {
 
 test_that("brm_marginal_summaries() on change", {
   set.seed(0L)
-  sim <- brm_simulate()
+  sim <- brm_simulate(
+    n_group = 2L,
+    n_patient = 100L,
+    n_time = 4L
+  )
   data <- sim$data
   data$group <- paste("treatment", data$group)
   data$time <- paste("visit", data$time)

tests/testthat/test-brm_model.R

@@ -1,6 +1,10 @@
 test_that("brm_model() runs", {
   set.seed(0L)
-  data <- brm_simulate()$data
+  data <- brm_simulate(
+    n_group = 2L,
+    n_patient = 100L,
+    n_time = 4L
+  )$data
   formula <- brm_formula(
     response = "response",
     group = "group",