mcmcapply | R Documentation |

Conveniently calculates statistics for mcmc.list objects.

mcmcapply(x, FUN, ...) ## S3 method for class 'mcmc.list' stack(x, ...)

`x` |
Objects of class |

`FUN` |
A function to be used in the calculations, returning a single value. |

`...` |
Other arguments passed to |

`mcmcapply`

returns a certain statistics based on `FUN`

after coercing into a matrix. `FUN`

can be missing,
in this case `mcmcapply`

is equivalent
to calling `as.matrix`

on an 'mcmc.list' object.

`stack`

can be used to concatenates 'mcmc.list'
objects into a single vector
along with index variables indicating where each observation originated
from (e.g. iteration, variable, chain).

`mcmcapply`

returns statistic value for each variable
based on `FUN`

, using all values in all chains of the MCMC object.

`stack`

returns a data frame with columns:
iter, variable, chain, value.

Peter Solymos, solymos@ualberta.ca

data(regmod) mcmcapply(regmod, mean) mcmcapply(regmod, sd) x <- stack(regmod) head(x) summary(x) library(lattice) xyplot(value ~ iter | variable, data=x, type="l", scales = "free", groups=chain)

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.