Extract Posterior Samples
Extract posterior samples of specified parameters.
## S3 method for class 'brmsfit' posterior_samples( x, pars = NA, fixed = FALSE, add_chain = FALSE, subset = NULL, as.matrix = FALSE, as.array = FALSE, ... ) posterior_samples(x, pars = NA, ...) ## S3 method for class 'brmsfit' as.data.frame(x, row.names = NULL, optional = TRUE, ...) ## S3 method for class 'brmsfit' as.matrix(x, ...) ## S3 method for class 'brmsfit' as.array(x, ...)
x |
An |
pars |
Names of parameters for which posterior samples should be returned, as given by a character vector or regular expressions. By default, all posterior samples of all parameters are extracted. |
fixed |
Indicates whether parameter names
should be matched exactly ( |
add_chain |
A flag indicating if the returned |
subset |
A numeric vector indicating the rows
(i.e., posterior samples) to be returned.
If |
as.matrix |
Should the output be a |
as.array |
Should the output be an |
... |
For |
row.names, optional |
See |
Currently there are methods for brmsfit objects.
as.data.frame.brmsfit, as.matrix.brmsfit, and
as.array.brmsfit are basically aliases of
posterior_samples.brmsfit and differ from
each other only in type of the returned object.
A data.frame (matrix or array) containing the posterior samples, with one column per parameter. In case an array is returned, it contains one additional dimension for the chains.
## Not run:
fit <- brm(rating ~ treat + period + carry + (1|subject),
data = inhaler, family = "cumulative")
# extract posterior samples of population-level effects
samples1 <- posterior_samples(fit, "^b")
head(samples1)
# extract posterior samples of group-level standard deviations
samples2 <- posterior_samples(fit, "^sd_")
head(samples2)
## End(Not run)Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.