DPpackage (version 1.1-6)

DPcaterpillar: Caterpillar Plots for Random Effects

Description

This generic function produces Caterpillar Plots for Random Effects from DPrandom objects.

Usage

DPcaterpillar(object, midpoint="mean", hpd=TRUE , ask=TRUE,  nfigr=1, nfigc=1, ...)

Arguments

object

DPrandom object from which random effects estimates can be extracted.

midpoint

variable indicating whether the mean or median of the posterior distribution of random effects should be considered as "midpoint" in the caterpillar plot.

hpd

logical variable indicating whether the hpd (TRUE) or pd (FALSE) of random effects should be considered in the caterpillar plot.

ask

logical variable indicating whether the caterpillar plots should be display gradually (TRUE) or not (FALSE).

nfigr

integer variable indicating the number of caterpillar plots by row.

nfigc

integer variable indicating the number of caterpillar plots by column.

...

further arguments passed to or from other methods.

Examples

## Not run: 
#     # School Girls Data Example
# 
#       data(schoolgirls)
#       attach(schoolgirls)
# 
#     # Prior information
#     # Prior information
# 
#       tinv<-diag(10,2)
#       prior<-list(alpha=1,nu0=4.01,tau1=0.001,tau2=0.001,
#       tinv=tinv,mub=rep(0,2),Sb=diag(1000,2))
# 
#     # Initial state
#       state <- NULL
# 
#     # MCMC parameters
# 
#       nburn<-5000
#       nsave<-25000
#       nskip<-20
#       ndisplay<-1000
#       mcmc <- list(nburn=nburn,nsave=nsave,nskip=nskip,
#                    ndisplay=ndisplay)
# 
#     # Fit the model
#     
#       fit1<-DPlmm(fixed=height~1,random=~age|child,prior=prior,
#                   mcmc=mcmc,state=state,status=TRUE)
#       fit1
# 
# 
#     # Extract random effects
#     
#       DPcaterpillar(DPrandom(fit1))
# ## End(Not run)