gibbsmat              package:mixstock              R Documentation

_P_l_o_t _G_i_b_b_s _s_a_m_p_l_e_r _o_u_t_p_u_t

_D_e_s_c_r_i_p_t_i_o_n:

     plot or report on output from Gibbs sampler

_U_s_a_g_e:

     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,...)

_A_r_g_u_m_e_n_t_s:

       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

_V_a_l_u_e:

     sets up an array of plots on the current output device

