Description Usage Arguments Author(s) Examples
Extracts random effects from PTglmm objects: PTlmm, PTolmm, and PTglmm.
1 |
object |
PT fitted model object from which random effects estimates can be extracted. |
centered |
logical variable indicating whether the random effects
should be extracted centered, |
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 and if predictive=TRUE. 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 <atjara@uc.cl>
Tim Hanson <hansont@stat.sc.edu>
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 | ## Not run:
# School Girls Data Example
data(schoolgirls)
attach(schoolgirls)
# Prior information
prior<-list(alpha=1,
M=4,
tau1=0.01,tau2=0.01,
nu0=4.01,
tinv=diag(10,2),
mub=rep(0,2),
Sb=diag(1000,2))
# Initial state
state <- NULL
# MCMC parameters
nburn<-5000
nsave<-5000
nskip<-0
ndisplay<-100
mcmc <- list(nburn=nburn,nsave=nsave,nskip=nskip,ndisplay=ndisplay,
tune1=1.5,tune2=1.1)
# Fitting the model
fit1<-PTlmm(fixed=height~1,random=~age|child,prior=prior,mcmc=mcmc,
state=state,status=TRUE)
fit1
# Extract random effects
PTrandom(fit1)
plot(PTrandom(fit1))
# Extract predictive information of random effects
PTrandom(fit1,predictive=TRUE)
plot(PTrandom(fit1,predictive=TRUE,gridl=c(75,89,3.8,7.5)))
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.