paraplot: visual display of the posterior inferences of the parameters...

In zoib: Bayesian Inference for Beta Regression and Zero-or-One Inflated Beta Regression

visual display of the posterior inferences of the parameters from a zoib model

Description

paraplot plots the posterior mode, mean, or median with Bayesian credible intervals for the parameters from a zoib model. plot.para can accommodate up to 4 groups (models) of inferences.

Usage

paraplot(para1, para2=NULL, para3=NULL, para4=NULL,
          tickx=NULL, jitter=NULL, pch = 1:4, col=1:4, 
          legpos=NULL, legtext, annotate=FALSE)

Arguments

para1	a data frame or matrix with posterior inferences from approach (model) 1.
para2	a data frame or matrix with posterior inferences from approach (model) 4; specify para2 only needed.
para3	a data frame or matrix with posterior inferences from approach (model) 3; specify para3 only needed.
para4	a data frame or matrix with posterior inferences from approach (model) 2; specify para4 only needed.
tickx	customization of the tickx of the x-axis.
jitter	jittering distance between the plotting positions of different groups to enhance plot readability.
pch	point types for different groups.
col	colors for different types.
legpos	the positions of legend in the format of c(x,y).
legtext	contents of the legtext, such as the names of the approaches.
annotate	whether to label parameters.

Details

Each of the applicabale para objects contains 3 columns. The 1st column is the point estimate (posterior mean, median or mode), the 2nd and 3rd columns contains the lower and upper bounds of the Bayesian credible intervals. The rows represent the parameters from an approach (model). It is not necessary that all para objects have the same set of parameters.

Author(s)

Fang Liu ([e-mail])

Examples

  ## Not run: 
  set.seed(12) 
  x=rnorm(4); para1 = cbind(x, x-1,x+1); rownames(para1) = c("a","b","c","d")
  x=rnorm(4)+1; para2 = cbind(x, x-1,x+1); rownames(para2) = c("a","b","e","f")
  x=rnorm(3)+2; para3 = cbind(x, x-1,x+1); rownames(para3) = c("a","b","d")
  paraplot(para1, para2, para3, para4=NULL, legpos=c(-1.5,6),
            legtext=c("model 1","model 2","model 3"),annotate=TRUE)
  paraplot(para1, legpos=c(-2,3), legtext="m1", annotate=TRUE)
	
## End(Not run)

zoib documentation built on May 31, 2023, 7:49 p.m.