R: Caterpillar Plots for Random Effects

DPcaterpillar {DPpackage} | R Documentation

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

### Author(s)

Alejandro Jara <atjara@uc.cl>

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

[Package DPpackage version 1.1-7.4 Index]