Dear all,
I’m preparing case-control study with a binary outcome in which I’d like to model the numerical covariates with splines; I’ll have 3-5 controls per case (we didn’t decide finally yet). The model should account for the matched nature of the data, so normally I’d go for stan_clogit
and use its strata
argument. However, because I’d like to use splines, I was planning on using stan_gamm4
which doesn’t seem to support conditioning on stratum out-of-the-box.
Thus, I was wondering if others have done this before and might be able to give some advise on the best approach. I suppose it would be possible to simply pre-compute the data matrix for the smoothers and use them in stan_clogit
but it would be neat to be able to use the machinery of stan_gamm4
, e.g. plot_nonlinear
. Including stratum as a random-effects intercept might be another approach, but it would change the nature of the estimates.
Might it be possible to actually condition on strata with stan_gamm4
(perhaps in a future release)?
Cheers,
Ben
- Operating System: Mac OS X Catalina
- rstanarm Version: 2.19.2