Contrasts

[fweber144]

At least for brms-fitted reference models, factor predictor columns in the $data element of the brmsfit get an attribute contrasts assigned. This makes sense and in the future, this should probably also be respected by projpred when fitting the submodels. For now, I think the warnings produced by some predict() methods (see below) should be suppressed and some documentation should be added, stating that any contrasts attributes from the reference model's dataset are ignored.

Reprex (a bit artificial since options(projpred.glm_fitter = "fit_glm_callback") is currently a "hidden" feature, but I encountered a similar issue while trying to solve #70):

options(mc.cores = parallel::detectCores(logical = FALSE))
dat <- data.frame(
  group = gl(n = 3, k = floor(41 / 3), length = 41,
             labels = paste0("gr", seq_len(3)))
)
set.seed(457211)
dat$y <- rnorm(nrow(dat), mean = -4.2, sd = 4)

bfit <- brms::brm(y ~ group,
                  data = dat,
                  backend = "cmdstanr",
                  refresh = 0,
                  seed = 1140350788)

library(projpred)
options(projpred.glm_fitter = "fit_glm_callback")
prj <- project(bfit,
               solution_terms = c("group"),
               nclusters = 3,
               seed = 46782345)
## --> Throws:
# Warning messages:
# 1: contrasts dropped from factor group
# 2: contrasts dropped from factor group
# 3: contrasts dropped from factor group
##

# Explanation:
refm <- get_refmodel(bfit)
str(refm$fetch_data())
## --> Note the `contrasts` attribute here.

[fweber144]

Removed the fixed in develop label because #284 was just a quick fix suppressing the warnings. As explained above, in the future, contrasts attributes should also be respected by projpred when fitting the submodels.