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 <ajarav@udec.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.0-7 Index]