inst/tests/test-dataset.R

In poissonconsulting/jaggernaut: Bayesian Analyses using 'JAGS'

context("dataset")

test_that("test multi par", {
  
  model <- jags_model("
                      model {
                      alpha ~ dunif(-20, 20)
                      beta1 ~ dunif(-10, 10)
                      beta2 ~ dunif(-10, 10)
                      beta3 ~ dunif(-10, 10)
                      
                      for (i in 1:length(Year)) {
                      log(eC[i]) <- alpha + beta1 * Year[i] 
                      + beta2 * Year[i]^2 + beta3 * Year[i]^3
                      C[i] ~ dpois(eC[i])
                      }
                      }",
 derived_code = "data{
 for (i in 1:length(Year)) {
 log(prediction[i]) <- alpha + beta1 * Year[i]
 + beta2 * Year[i]^2 + beta3 * Year[i]^3
 }
 }",
select_data = c("C","Year*")
  )

data <- jaggernaut::falcon

data$C <- data$Pairs

analysis <- jags_analysis (model, data, mode = "test")

expect_equal(dataset(analysis), data)

analysis <- jags_analysis (combine(model,model), data, mode = "test")

expect_equal(dataset(analysis), data)
x <- dataset(analysis, converted = TRUE)
expect_is(x, "list")
expect_equal(names(x), c("Model1", "Model2"))
expect_equal(x[["Model1"]], x[["Model2"]])
expect_equal(names(x[["Model1"]]), c("C", "Year"))
expect_equal(names(x[["Model2"]]), c("C", "Year"))
})