tests/testthat/tests-increment_IGscale.R

In emittman/dpOMP: Parallelized Gibbs Sampler for DPGMM Model

context("Increment IG scale")

len <- 10

Gk <- sample(1:10, 1)

xTy <- rnorm(len)

xTx <- crossprod(matrix(runif(len*len),len,len))

beta <- rnorm(len)

incR <- 2*t(xTy)%*%beta - Gk * t(beta) %*% xTx %*% beta

incC <- increment_IGscale(xTy, xTx, beta, Gk)

test_that("Correct output", {
  expect_equal(drop(incR), incC)
})