View source: R/all_postprodn_fns.R
summary.angmcmc | R Documentation |
Summary statistics for parameters from an angmcmc object
## S3 method for class 'angmcmc'
summary(object, par.name, comp.label, chain.no, ...)
object |
angular MCMC object. |
par.name |
vector of names of parameters for which point estimates are to be computed. If |
comp.label |
vector of component labels (positive integers, e.g., |
chain.no |
vector of chain numbers whose samples are to be be used. in the estimation. By default all chains are used. |
... |
additional arguments affecting the summary produced. |
Computes (after thinning and discarding burn-in) point estimates with 95% posterior credible sets for all components and all parameters, together with the sample averages of log likelihood and log posterior density.
Returns a list with elements estimate, lower, upper, llik
and lpd
. estimate
is itself a list with three elements: mean
, median
and mode
providing the
sample mean, sample median and (sample) MAP estimates.
Note that summary.angmcmc
has its own print method, providing a table the estimated mean and 95% credible intervals for each parameter
# illustration only - more iterations needed for convergence
fit.vmsin.20 <- fit_vmsinmix(tim8, ncomp = 3, n.iter = 20,
n.chains = 1)
summary(fit.vmsin.20)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.