context("create my_inla object")
describe("my_inla", {
data(cbpp, package = "lme4")
meta1 <- 1L
meta2 <- 2L
meta3 <- 3L
meta4 <- 4L
meta5 <- "Meta 5"
meta6 <- "Meta 6"
it("fills the object correctly", {
expect_that(
my_inla(
data = cbpp,
meta1 = meta1,
meta2 = meta2,
meta3 = meta3,
meta4 = meta4,
meta5 = meta5,
meta6 = meta6
)@Data,
is_identical_to(cbpp)
)
expect_that(
my_inla(
data = cbpp,
meta1 = meta1,
meta2 = meta2,
meta3 = meta3,
meta4 = meta4,
meta5 = meta5,
meta6 = meta6
)@Meta1,
is_identical_to(meta1)
)
expect_that(
my_inla(
data = cbpp,
meta1 = meta1,
meta2 = meta2,
meta3 = meta3,
meta4 = meta4,
meta5 = meta5,
meta6 = meta6
)@Status,
is_identical_to("new")
)
expect_that(
my_inla(
data = cbpp,
meta1 = meta1,
meta2 = meta2,
meta3 = meta3,
meta4 = meta4,
meta5 = meta5,
meta6 = meta6
)@Weight,
throws_error('no slot of name "Weight"')
)
})
object <- my_inla(
data = cbpp,
meta1 = meta1,
meta2 = meta2,
meta3 = meta3,
meta4 = meta4,
meta5 = meta5,
meta6 = meta6
)
require(INLA)
model <- inla(incidence ~ period + f(herd, model = "iid"), family = "poisson", data = cbpp)
it("creates object from the correct class", {
expect_that(
object,
is_a("my_inla")
)
expect_that(
my_inla(data = object, model.fit = model, status = "converged"),
is_a("my_inla")
)
})
it("fills the object correctly", {
expect_that(
my_inla(data = object, model.fit = model, status = "converged")@Model,
is_identical_to(model)
)
})
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.