Nothing
## for old versions of lme4, the refit for g/lmer doesn't move very
## far from the fit so we just suppress the test
lme4Version <- packageVersion("lme4")
if (lme4Version >= "1.1-6") {
context("refit generic for blmerMod and bglmerMod classes")
test_that("refit for blmer matches original, not lmer", {
source(system.file("common", "lmmData.R", package = "blme"))
control <- lmerControl(optimizer = "Nelder_Mead")
cov.prior <- "g.1 ~ wishart(scale = 2)"
fit <- blmer(y ~ x.1 + x.2 + (1 + x.1 | g.1), testData, cov.prior = cov.prior, control = control)
blmerRefit <- refit(fit)
lmerRefit <- getS3method("refit", "merMod")(fit)
expect_equal(fit@theta, blmerRefit@theta, tolerance = 1.0e-02)
expect_equal(fit@beta, blmerRefit@beta, tolerance = 1.0e-03)
expect_false(all(abs(fit@theta - lmerRefit@theta) <= 1.0e-02))
expect_false(all(abs(fit@beta - lmerRefit@beta) <= 1.0e-03))
})
test_that("refit for bglmer matches original, not glmer", {
source(system.file("common", "glmmData.R", package = "blme"))
control <- if (lme4Version >= "1.1-8")
glmerControl(optimizer = "Nelder_Mead", nAGQ0initStep = FALSE)
else
glmerControl(optimizer = "Nelder_Mead")
fit <- bglmer(y ~ x.1 + x.2 + (1 | g), testData, family = binomial(), control = control,
cov.prior = wishart)
bglmerRefit <- refit(fit)
glmerRefit <- getS3method("refit", "merMod")(fit)
expect_equal(fit@theta, bglmerRefit@theta, tolerance = 1.0e-3)
expect_equal(fit@beta, bglmerRefit@beta, tolerance = 1.0e-3)
expect_false(all(abs(fit@theta - glmerRefit@theta) <= 1.0e-3))
expect_false(all(abs(fit@beta - glmerRefit@beta) <= 1.0e-3))
})
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.