Nothing
context("lmer and blmer")
test_that("blmer matches lmer exactly", {
source(system.file("common", "lmmData.R", package = "blme"))
control <- lmerControl(optimizer = "bobyqa")
lmerFit <- lmer(y ~ x.1 + x.2 + (1 + x.1 | g.1) + (1 + x.1 + x.2 | g.2), testData, control = control)
blmerFit <- blmer(y ~ x.1 + x.2 + (1 + x.1 | g.1) + (1 + x.1 + x.2 | g.2), testData, control = control,
cov.prior = NULL, fixef.prior = NULL, resid.prior = NULL)
expect_equal(lmerFit@pp$theta, blmerFit@pp$theta)
expect_equal(lmerFit@pp$u(1.0), blmerFit@pp$u(1.0))
expect_equal(lmerFit@pp$beta(1.0), blmerFit@pp$beta(1.0))
})
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.