inst/tests/test-control-mode.R

In andrewzm/INLA: Functions which allow to perform full Bayesian analysis of latent Gaussian models using Integrated Nested Laplace Approximaxion

context("test 'control.mode'")

test_that("Case 1", {
    set.seed(123)
    n = 100
    x = arima.sim(n, model=list(ar=0.9))
    x = x/sd(x)
    y = 1 + x + rnorm(n, sd=0.01)
    
    idx = 1:n
    formula = y ~ 1 + f(idx,  model="ar1")
    
    inla.setOption(inla.call='inla.work')
    r = inla(formula, data = data.frame(y, idx),
            control.inla = list(int.strategy = "eb"))
    rr = inla(formula, data = data.frame(y, idx),
            control.mode = list(result = r,  fixed=TRUE))
    expect_true(all(abs(r$summary.random[[1]] - rr$summary.random[[1]]) < 1e-8))
    expect_true(all(abs(r$summary.fixed - rr$summary.fixed) < 1e-8))
})