Extract samples of each parameter from the BKMR fit
Arguments
- fit
An object containing the results returned by a the
kmbayesfunction- sel
logical expression indicating samples to keep; defaults to keeping the second half of all samples
Value
a list where each component contains the posterior samples of one of the parameters (or vector of parameters) being estimated
Examples
## First generate dataset
set.seed(111)
dat <- SimData(n = 50, M = 4)
y <- dat$y
Z <- dat$Z
X <- dat$X
## Fit model with component-wise variable selection
## Using only 100 iterations to make example run quickly
## Typically should use a large number of iterations for inference
set.seed(111)
fitkm <- kmbayes(y = y, Z = Z, X = X, iter = 100, verbose = FALSE, varsel = TRUE)
#> Iteration: 10 (10% completed; 0.00486 secs elapsed)
#> Iteration: 20 (20% completed; 0.01027 secs elapsed)
#> Iteration: 30 (30% completed; 0.01583 secs elapsed)
#> Iteration: 40 (40% completed; 0.02201 secs elapsed)
#> Iteration: 50 (50% completed; 0.02795 secs elapsed)
#> Iteration: 60 (60% completed; 0.04045 secs elapsed)
#> Iteration: 70 (70% completed; 0.04589 secs elapsed)
#> Iteration: 80 (80% completed; 0.05089 secs elapsed)
#> Iteration: 90 (90% completed; 0.05624 secs elapsed)
#> Iteration: 100 (100% completed; 0.06148 secs elapsed)
samps <- ExtractSamps(fitkm)