eggCounts (version 2.4)

stan2mcmc: Convert a Stanfit object to a MCMC object

Description

Converts a stanfit object into a mcmc object for easier analysis.

Usage

stan2mcmc(stanFit)

Value

A MCMC object with a list of relevant parameters.

Arguments

stanFit

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

Author

Craig Wang

Details

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

Examples

data(epgs)

## apply zero-infation model for the paired design 
model <- fecr_stan(epgs$before, epgs$after, rawCounts = FALSE, 
            indEfficacy = FALSE, preCF = 10, 
            paired = TRUE, zeroInflation = TRUE)
samples <- stan2mcmc(model$stan.samples)
summary(samples)
plot(samples)