From 268aac5b4b651a1eab6da9562505c43976eb492b Mon Sep 17 00:00:00 2001 From: Dr Nathan Green Date: Sat, 24 Aug 2024 19:21:45 +0100 Subject: [PATCH] fixed bug where Id mixed up rep and each for number of groups and 2 * still need to double check this is actually correct --- R/default_prior_cure.R | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/R/default_prior_cure.R b/R/default_prior_cure.R index 5b4030d..4041ca4 100644 --- a/R/default_prior_cure.R +++ b/R/default_prior_cure.R @@ -13,6 +13,7 @@ default_prior_cure <- function(formula_cure, prior_cure = list(), bg_model = 2) { + browser() nTx <- formula_cure$fe_nlevels[1] n_groups <- formula_cure$n_groups nvars <- formula_cure$nvars @@ -63,7 +64,7 @@ default_prior_cure <- function(formula_cure, rep(formula_cure$n_groups) names(params_cf) <- - paste0(names(params_cf), "_", rep(1:2, each = formula_cure$n_groups)) + paste0(names(params_cf), "_", rep(1:formula_cure$n_groups, each = 2)) } else if (is_pooled_cf(formula_cure)) { params_cf <- list(a_cf = 3,