DPrandom              package:DPpackage              R Documentation

_E_x_t_r_a_c_t_s _R_a_n_d_o_m _E_f_f_e_c_t_s

_D_e_s_c_r_i_p_t_i_o_n:

     This generic function extracts Random Effects' information from
     DPpackage model objects.

_U_s_a_g_e:

     DPrandom(object,centered=FALSE,predictive=FALSE)

_A_r_g_u_m_e_n_t_s:

  object: DP fitted model object from which random effects estimates
          can be extracted.

centered: logical variable indicating whether the random effects should
          be extracted centered, 'bi', or  uncentered 'thetai'.

predictive: logical variable indicating whether actual or predictive
          information of the random effects should be extracted.

_E_x_a_m_p_l_e_s:

     ## 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
         
           DPrandom(fit1)
           DPrandom(fit1,centered=TRUE)
           
           plot(DPrandom(fit1))
           plot(DPrandom(fit1,centered=TRUE))

         # Extract predictive information of random effects
         
           DPrandom(fit1,predictive=TRUE)
           DPrandom(fit1,centered=TRUE,predictive=TRUE)
           
           plot(DPrandom(fit1,predictive=TRUE))
           plot(DPrandom(fit1,centered=TRUE,predictive=TRUE))
     ## End(Not run)

