tests/testthat/test-latent_linear.R

In inlabru: Bayesian Latent Gaussian Modelling using INLA and Extensions

test_that("bru: linear component", {
  skip_on_cran()
  local_bru_safe_inla()

  # Seed influences data as well as predict()!
  set.seed(123)

  input.df <- data.frame(x = cos(1:100))
  input.df <- within(input.df, y <- 5 + 2 * x + rnorm(100, mean = 0, sd = 0.1))

  fit <- bru(
    y ~ myLin(main = x, model = "linear") + Intercept(1),
    family = "gaussian",
    data = input.df
  )

  expect_equal(fit$summary.fixed["myLin", "mean"], 2.002273, tolerance = midtol)
  expect_equal(fit$summary.fixed["myLin", "sd"], 0.01323361, tolerance = hitol)

  skip_if_not_installed("sn")
  pr <- predict(
    fit,
    data.frame(x = c(1, 2)),
    ~ myLin + 2,
    n.samples = 5,
    seed = 1L
  )

  expect_equal(pr[, "mean"], c(4.005013, 6.010026), tolerance = midtol)
})