tests/testthat/test-bivariate.R

In robingenuer/mypkgr: Example of a package for teaching advanced R

context("Bivariate gaussian test")

test_that("correct results for bivariate gaussian", {
  expect_equal(mvnpdf(x=matrix(rep(1.96,2), nrow=2, ncol=1),
                      mean=c(0, 0), varcovM=diag(2), Log=FALSE)$y,
               mvtnorm::dmvnorm(rep(1.96, 2)))
  expect_equal(mvnpdfsmart(x=matrix(rep(1.96,2), nrow=2, ncol=1),
                           mean=c(0, 0), varcovM=diag(2), Log=FALSE)$y,
               mvtnorm::dmvnorm(rep(1.96, 2)))
  expect_equal(mvnpdfoptim(x=matrix(rep(1.96,2), nrow=2, ncol=1),
                           mean=c(0, 0), varcovM=diag(2), Log=FALSE)$y,
               mvtnorm::dmvnorm(rep(1.96, 2)))
  expect_equal(mvnpdf_invC(x=matrix(rep(1.96,2), nrow=2, ncol=1),
                           mean=c(0, 0), varcovM=diag(2), Log=FALSE)$y,
               mvtnorm::dmvnorm(rep(1.96, 2)))
  expect_equal(mvnpdfC(x=matrix(rep(1.96,2), nrow=2, ncol=1),
                      mean=c(0, 0), varcovM=diag(2), Log=FALSE),
               mvtnorm::dmvnorm(rep(1.96, 2)))
})