reportSAM {numOSL} | R Documentation |
Summarizing distributions of parameters simulated from statistical age models.
reportSAM(obj, burn = 10000, thin = 5, plot = TRUE, outfile = NULL, ...)
obj |
list(required): an object of S3 class |
burn |
integer(with default): number of iterations (i.e., the initial, non-stationary |
thin |
integer(with default): take every |
plot |
logical(with default): plot the MCMC output or not |
outfile |
character(optional): if specified, simulated parameters will be written to a CSV file named |
... |
do not use |
Function reportSAM summarizes the output of a Markov chain (the mean values, the standard deviations, and the 95 percent confidence intervals). The initial i (burn=i
) samples may have been affected by the inital state and has to be discarded ("burn-in"). Autocorrelation of simulated samples can be reduced by taking every j-th (thin=j
) iteration ("thining").
Return a list that contains the following elements:
pars |
parameters and relevant statistics |
quantile |
quantiles of simulated parameters |
maxlik |
calculated maximum logged likelihood value |
Lunn D, Jackson C, Best N, Thomas A, Spiegelhalter D, 2013. The BUGS book: a practical introduction to bayesian analysis. Chapman & Hall/CRC Press.
Gelman A, Carlin JB, Stern HS, Dunson DB, Vehtari A, Rubin DB, 2013. Bayesian data analysis. Chapman & Hall/CRC Press.