DPMrandom {DPpackage} | R Documentation |
Extracts random effects from DPpackage objects: DPMlmm, DPMolmm, and DPMglmm.
DPMrandom(object,centered=FALSE,predictive=FALSE,ngrid=1000,gridl=NULL)
object |
DPM 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. |
ngrid |
number of grid points where the density estimate is evaluated. This is only used if dimension of the random effects is lower or equal than 2. The default value is 1000. |
gridl |
The limits of the interval or rectangle covered by the grid as c(xl,xu) or c(xl, xu, yl, yu), respectively. If not specified the grid is defined automatically. This is only used if dimension of the random effects is lower or equal than 2 and if predictive=TRUE. |
Alejandro Jara <ajarav@udec.cl>
## Not run: # School Girls Data Example data(schoolgirls) attach(schoolgirls) # Prior information prior<-list(alpha=1, tau1=0.01,tau2=0.01, nu0=4.01, tinv=diag(10,2), nub=4.01, tbinv=diag(10,2), mb=rep(0,2), Sb=diag(1000,2)) # Initial state state <- NULL # MCMC parameters nburn<-5000 nsave<-10000 nskip<-20 ndisplay<-1000 mcmc <- list(nburn=nburn,nsave=nsave,nskip=nskip,ndisplay=ndisplay) # Fitting the model fit1<-DPMlmm(fixed=height~1,random=~age|child,prior=prior,mcmc=mcmc, state=state,status=TRUE) fit1 # Extract random effects DPMrandom(fit1) DPMrandom(fit1,centered=TRUE) plot(DPMrandom(fit1)) plot(DPMrandom(fit1,centered=TRUE)) # Extract predictive information of random effects DPMrandom(fit1,predictive=TRUE) plot(DPMrandom(fit1,predictive=TRUE,gridl=c(75,89,3.8,7.5))) ## End(Not run)