DPMrandom {DPpackage} | R Documentation |
Extracts Random Effects
Description
Extracts random effects from DPpackage objects: DPMlmm, DPMolmm, and DPMglmm.
Usage
DPMrandom(object,centered=FALSE,predictive=FALSE,ngrid=1000,gridl=NULL)
Arguments
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. |
Examples
## 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)
[Package
DPpackage version 1.0-5
Index]