DPpackage (version 1.1-0)

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

# 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))