Skip to content

Commit

Permalink
brm_marginal_summaries
Browse files Browse the repository at this point in the history
  • Loading branch information
wlandau-lilly committed Jun 8, 2023
1 parent 74c0627 commit 59ef917
Show file tree
Hide file tree
Showing 2 changed files with 31 additions and 39 deletions.
22 changes: 11 additions & 11 deletions R/brm_marginal_summaries.R
Original file line number Diff line number Diff line change
Expand Up @@ -31,15 +31,18 @@
#' for the credible intervals.
#' @examples
#' set.seed(0L)
#' sim <- brm_simulate()
#' data <- sim$data
#' data <- brm_data(
#' data = tibble::as_tibble(brm_simulate()$data),
#' outcome = "response",
#' role = "response",
#' group = "group",
#' time = "time",
#' patient = "patient"
#' )
#' data$group <- paste("treatment", data$group)
#' data$time <- paste("visit", data$time)
#' formula <- brm_formula(
#' response = "response",
#' group = "group",
#' time = "time",
#' patient = "patient",
#' data = data,
#' effect_base = FALSE,
#' interaction_base = FALSE
#' )
Expand All @@ -58,12 +61,9 @@
#' )
#' draws <- brm_marginal_draws(
#' model = model,
#' group = "group",
#' time = "time",
#' patient = "patient",
#' data = data,
#' control = "treatment 1",
#' baseline = "visit 1",
#' outcome = "response"
#' baseline = "visit 1"
#' )
#' suppressWarnings(brm_marginal_summaries(draws))
brm_marginal_summaries <- function(
Expand Down
48 changes: 20 additions & 28 deletions tests/testthat/test-brm_marginal_summaries.R
Original file line number Diff line number Diff line change
@@ -1,18 +1,17 @@
test_that("brm_marginal_summaries() on response", {
set.seed(0L)
sim <- brm_simulate(
n_group = 2L,
n_patient = 100L,
n_time = 4L
data <- brm_data(
data = tibble::as_tibble(brm_simulate()$data),
outcome = "response",
role = "response",
group = "group",
time = "time",
patient = "patient"
)
data <- sim$data
data$group <- paste("treatment", data$group)
data$time <- paste("visit", data$time)
formula <- brm_formula(
response = "response",
group = "group",
time = "time",
patient = "patient",
data = data,
effect_base = FALSE,
interaction_base = FALSE
)
Expand All @@ -31,12 +30,9 @@ test_that("brm_marginal_summaries() on response", {
)
draws <- brm_marginal_draws(
model = model,
group = "group",
time = "time",
patient = "patient",
data = data,
control = "treatment 1",
baseline = "visit 1",
outcome = "response"
baseline = "visit 1"
)
suppressWarnings(
x <- brm_marginal_summaries(
Expand Down Expand Up @@ -125,19 +121,18 @@ test_that("brm_marginal_summaries() on response", {

test_that("brm_marginal_summaries() on change", {
set.seed(0L)
sim <- brm_simulate(
n_group = 2L,
n_patient = 100L,
n_time = 4L
data <- brm_data(
data = tibble::as_tibble(brm_simulate()$data),
outcome = "response",
role = "change",
group = "group",
time = "time",
patient = "patient"
)
data <- sim$data
data$group <- paste("treatment", data$group)
data$time <- paste("visit", data$time)
formula <- brm_formula(
response = "response",
group = "group",
time = "time",
patient = "patient",
data = data,
effect_base = FALSE,
interaction_base = FALSE
)
Expand All @@ -156,12 +151,9 @@ test_that("brm_marginal_summaries() on change", {
)
draws <- brm_marginal_draws(
model = model,
group = "group",
time = "time",
patient = "patient",
data = data,
control = "treatment 1",
baseline = "visit 1",
outcome = "change"
baseline = "visit 1"
)
suppressWarnings(
x <- brm_marginal_summaries(
Expand Down

0 comments on commit 59ef917

Please sign in to comment.