stan2mcmc: Convert a Stanfit object to MCMC object

Description Usage Arguments Details Value Author(s) Examples

View source: R/stan2mcmc.R

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)

eggCounts documentation built on May 29, 2017, 5:59 p.m.

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