## testing getME function
require(robustlmm)
formatNum <- function(x, ...)
format(round(x, 8), trim = TRUE, drop0trailing = TRUE, ...)
options(str = strOptions(formatNum = formatNum))
sleepstudy2 <- within(sleepstudy, Group <- letters[1:4])
rfm <- rlmer(Reaction ~ Days + (Days|Subject) + (1|Group), sleepstudy2,
rho.e = smoothPsi, rho.b=smoothPsi,
rho.sigma.e = psi2propII(smoothPsi, k=2.28),
rho.sigma.b = psi2propII(smoothPsi, k=2.28),
init = lmer(Reaction ~ Days + (Days|Subject) + (1|Group), sleepstudy2,
control=lmerControl(optimizer="bobyqa")),
doFit=FALSE)
(nmME <- eval(formals(robustlmm:::getME.rlmerMod)$name))
for (nm in nmME) {
cat("\nName:", nm, "\n")
value <- getME(rfm, name=nm)
if (nm == "theta") {
value <- value + 1
} else if (nm == "A") {
value@x <- value@x + 1
}
if (substr(nm, 1, 3) == "rho") {
print(value)
} else {
str(value)
}
}
g.all <- getME(rfm, "ALL")
g.all[grepl("^rho", names(g.all))] <- NULL
g.all[["theta"]] <- g.all[["theta"]] + 1
str(g.all, max.level = 2)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.