tests/testthat/test-sample_omega.R

In josbop/BayesBiSurv:

test_that("simulating omega", {
  omega <- c(1.4, 2.2, 1.7)
  t <- c(0.3, 1.9, 2.1)
  y <- c(1, 0, 3)
  x <- list(1, 0, 1)
  theta <- 0.45
  gamma <- 0.32
  p1 <- c(0, 1, 2, 3)
  lambda <- c(0.74, 0.64, 0.8)

  expect_type(
    sample_omega(omega, t, y, x, theta, gamma, p1, lambda),
    "double"
  )
  expect_length(
    sample_omega(omega, t, y, x, theta, gamma, p1, lambda),
    length(omega)
    )
})