gibbsmat {mixstock} | R Documentation |
plot or report on output from Gibbs sampler
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,...)
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 |
sets up an array of plots on the current output device