gibbsmat {mixstock}R Documentation

Plot Gibbs sampler output

Description

plot or report on output from Gibbs sampler

Usage

gibbsmat(x, burnin = 500, R = 2, H = 2, trans = TRUE)
gibbsrpt(x,burnin=500,R=2,H=2,trans=TRUE,plot=TRUE)
plotvar(vec,best,name,...)

Arguments

x results of an MCMC fit
burnin length of burn-in period for chain
R number of rookeries
H number of haplotypes
trans transform chain values?
plot plot posterior distributions for each variable?
vec numeric vector (posterior distribution chain)
best best-fit/point estimate value
name name of variable, for x label
... other arguments to plot

Value

sets up an array of plots on the current output device


[Package mixstock version 0.9 Index]