`varsel()` for GAMMs: `object not found` error

[fweber144]

On branch develop (commit 96d8e11), this reprex:

options(mc.cores = parallel::detectCores(logical = FALSE))
data("kidiq", package = "rstanarm")
kidiq_gr <- within(kidiq, {
  mom_age_gr <- cut(mom_age,
                    breaks = unique(quantile(mom_age, probs = seq(0, 1, 0.1))),
                    include.lowest = TRUE)
})
fit_gauss_gamm <- rstanarm::stan_gamm4(
  kid_score ~ s(mom_iq, by = mom_age_gr),
  random = ~ (mom_iq | mom_age_gr),
  data = kidiq_gr,
  iter = 500,
  seed = 734572
)

library(projpred)
### To avoid issue #144:
library(lme4)
###
vs_gauss_gamm <- varsel(fit_gauss_gamm, nclusters = 3, nclusters_pred = 5)

throws the following error:

Error in eval(formula[[2]], data, environment(formula)) : 
  object 'mom_age_gr' not found

[AlejandroCatalina]

This happens because we have not implemented parsing support for by syntax in smooth terms. It is desirable to have but we haven’t had the time yet. Best, Alejandro

…

On 10 Jun 2021, 4:09 PM +0300, Frank Weber ***@***.***>, wrote: On branch develop (commit 96d8e11), this reprex: options(mc.cores = parallel::detectCores(logical = FALSE)) data("kidiq", package = "rstanarm") kidiq_gr <- within(kidiq, { mom_age_gr <- cut(mom_age, breaks = unique(quantile(mom_age, probs = seq(0, 1, 0.1))), include.lowest = TRUE) }) fit_gauss_gamm <- rstanarm::stan_gamm4( kid_score ~ s(mom_iq, by = mom_age_gr), random = ~ (mom_iq | mom_age_gr), data = kidiq_gr, iter = 500, seed = 734572 ) library(projpred) ### To avoid issue #144: library(lme4) ### vs_gauss_gamm <- varsel(fit_gauss_gamm, nclusters = 3, nclusters_pred = 5) throws the following error: Error in eval(formula[[2]], data, environment(formula)) : object 'mom_age_gr' not found — You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub, or unsubscribe.