tests/testthat/test-derive.R

In audrey-b/simanalyse: Facilitate Model Comparisons and Simulation Studies

library(testthat)
library(simanalyse)

context("function-derive")

test_that("sma_derive works", {
  set.seed(10L)
  code <- "for(i in 1:10){x[i] ~ dnorm(0,1/variance)}"
  parameters = nlist(variance=4)
  dat <- sims::sims_simulate(code, parameters = parameters, nsims=2)
  #sma_derive(dat, "y=x+c", values=list(c=2))
  res <- sma_analyse(dat, code, code.add = "variance ~ dunif(0,10)", 
                              mode=sma_set_mode("quick"), monitor="variance",
                              deviance=FALSE)
  
  result <- sma_derive(res, "sd=sqrt(variance*c)", values=list(c=2))
  
  result1 <- readRDS("Derive_results/result1.rds")
  
  expect_equal(result, result1)
  
  result <- sma_derive(parameters, "sd=sqrt(variance)")
  
  result2 <- readRDS("Derive_results/result2.rds")
  
  expect_equal(result, result2)
  
})