Nothing
context("blmer numerical results with residual variance prior")
test_that("blmer fits the eight schools example correctly", {
# eight schools
y <- c(28, 8, -3, 7, -1, 1, 18, 12)
sigma <- c(15, 10, 16, 11, 9, 11, 10, 18)
y.z <- (y - mean(y)) / sigma
g <- 1:8
control <- lmerControl(check.conv.singular = "ignore",
check.nobs.vs.nRE = "ignore",
check.nobs.vs.nlev = "ignore")
eightSchools <- blmer(y.z ~ 1 + (1 | g), REML = FALSE,
control = control,
resid.prior = point,
cov.prior = NULL, fixef.prior = NULL)
expect_equal(eightSchools@theta, 0, tolerance = 1.0e-7)
})
test_that("blmer fits test data with inv.gamma prior, matching previous version", {
source(system.file("common", "lmmData.R", package = "blme"))
control <- lmerControl(optimizer = "bobyqa")
fit <- blmer(y ~ x.1 + x.2 + (1 + x.1 | g.1) + (1 + x.1 + x.2 | g.2), testData, control = control,
cov.prior = NULL, resid.prior = invgamma(2, 1.0))
expect_equal(fit@pp$theta, c(0.725321928185923, -0.251272308917427, 1.5828609906233, 0.946932542474828, 0.467970716580088, -0.183212783510381, 1.07158442297183, 0.122067368879505, 0.223238050522642),
tolerance = 1.0e-5)
})
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.