tests/testthat/test-mvnpdf.R

In borishejblum/mypkgr: Example Of A Package For Teaching HeavyR

test_that("correct result for univariate gaussian", {
  expect_equal(mvnpdf(x=matrix(1.96), mean=0, varcovM=diag(1), Log=FALSE)$y,
               dnorm(1.96))
  expect_equal(mvnpdf(x=matrix(c(1.96, -0.5), ncol = 2), mean=0,
                      varcovM=diag(1), Log=FALSE)$y, dnorm(c(1.96, -0.5)))
})


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(mvnpdfsmartC(x=matrix(rep(1.96,2), nrow=2, ncol=1),
                            mean=c(0, 0), varcovM=diag(2), Log=FALSE)[1,1],
               mvtnorm::dmvnorm(rep(1.96, 2)))
})