Robust Bayesian Longitudinal Regularized Semiparametric Mixed Models
Our recently developed fully robust Bayesian semiparametric mixed-effect model for high-dimensional longitudinal studies with heterogeneous observations can be implemented through this package. This model can distinguish between time-varying interactions and constant-effect-only cases to avoid model misspecifications. Facilitated by spike-and-slab priors, this model leads to superior performance in estimation, identification and statistical inference. In particular, robust Bayesian inferences in terms of valid Bayesian credible intervals on both parametric and nonparametric effects can be validated on finite samples. The Markov chain Monte Carlo algorithms of the proposed and alternative models are efficiently implemented in 'C++'.
- To install from github, run these two lines of code in R
install.packages("devtools")
devtools::install_github("kunfa/Blend")
- Released versions of Blend are available on CRAN (link), and can be installed within R via
install.packages("Blend")
library(Blend)
data(dat)
fit = Blend(y,x,t,J,kn,degree)
fit$coefficient
Coverage(fit)
plot_Blend(fit,sparse=TRUE)
fit = Blend(y,x,t,J,kn,degree,structural=FALSE)
fit = Blend(y,x,t,J,kn,degree, robust=FALSE)
fit = Blend(y,x,t,J,kn,degree, robust=FALSE, structural=FALSE)
- Updated README file.
- Fixed minor bugs.
This package provides implementation for methods proposed in
-Fan, K., Ren, J., Ma, Shuangge and Wu, C. (2025+). robust Bayesian Regularized Semiparametric Mixed Models in Longitudinal Studies. (submitted).