SCM Repository
[matrix] Diff of /pkg/man/mcmcsamp.Rd
Diff of /pkg/man/mcmcsamp.Rd
Parent Directory
|
Revision Log
|
Patch
48 |
\examples{ |
\examples{ |
49 |
require("lattice", quietly = TRUE, character = TRUE) |
require("lattice", quietly = TRUE, character = TRUE) |
50 |
(fm1 <- lmer(Reaction ~ Days + (Days|Subject), sleepstudy)) |
(fm1 <- lmer(Reaction ~ Days + (Days|Subject), sleepstudy)) |
51 |
|
set.seed(101) |
52 |
samp1 <- mcmcsamp(fm1, n = 1000) |
samp1 <- mcmcsamp(fm1, n = 1000) |
53 |
frm <- |
frm <- |
54 |
data.frame(vals = c(samp1), iter = rep(1:nrow(samp1), ncol(samp1)), |
data.frame(vals = c(samp1), iter = rep(1:nrow(samp1), ncol(samp1)), |
65 |
print(summary(samp1)) |
print(summary(samp1)) |
66 |
print(autocorr.diag(samp1)) |
print(autocorr.diag(samp1)) |
67 |
} |
} |
68 |
|
(eDF <- mean(samp1[,"deviance"]) - deviance(fm1)) # potentially useful approximate D.F. |
69 |
|
\dontshow{ |
70 |
|
stopifnot(abs(eDF - 15) < 1) |
71 |
|
} |
72 |
} |
} |
73 |
\keyword{methods} |
\keyword{methods} |
74 |
\keyword{datagen} |
\keyword{datagen} |
|
Legend:
Removed from v.1298 |
|
changed lines |
|
Added in v.1299 |
|
|