stan2mcmc: Convert a Stanfit object to MCMC object

Description Usage Arguments Details Value Author(s) Examples

Description

Converts a large stanfit object into a MCMC object for easier analysis, it extracts the relevant MCMC samples of the model from stanFit including the baseline mean epg, dispersion, pre- and post-treatment zero-inflation parameters and the calculated reduction.

Usage

1
stan2mcmc(stanFit)

Arguments

stanFit

A stanfit object from the output of either fecr_stan() or fec_stan()

Details

The output can be analyzed as a typical MCMC object with the functions from the coda package. NOTE: The resulting MCMC object does not contain warm-up samples and is already thinned.

Value

A MCMC object with a list of relevant parameters depending on the chosen model.

Author(s)

Craig Wang craig.wang@uzh.ch

Examples

1
2
3
4
5
6
7
8
9
## Not run: 
data(epgs)

## apply zero-infation model for the paired design to the data vector
model <- fecr_stan(epgs[,1], epgs[,2], rawCounts=FALSE, preCF=10, paired=TRUE, zeroInflation=TRUE)
samples <- stan2mcmc(model$stan.samples)
summary(samples)

## End(Not run)


Search within the eggCounts package
Search all R packages, documentation and source code

Questions? Problems? Suggestions? or email at ian@mutexlabs.com.

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.