tests/testthat/test.fit_mvnorm.R

In ashiklom/mvtraits: Fit multivariate models to partially missing trait data

context("Fit multivariate model")

dat <- as.matrix(iris[, -5])
niter <- 10000
nchains <- 8
keep_samples <- Inf
max_attempts <- 20
autofit <- TRUE

test_that("Simple multivariate fit, sequentially", {
  result <- fit_mvnorm(
    dat, niter = niter, nchains = nchains,
    autofit = autofit, keep_samples = keep_samples, max_attempts = max_attempts
  )
})

test_that("Multivariate fit, in parallel", {

  future::plan("multiprocess")
  result <- fit_mvnorm(
    dat, niter = niter, nchains = nchains,
    autofit = autofit, keep_samples = keep_samples, max_attempts = max_attempts
  )

})