tests/testthat/test-det_AtA.R

In jianan/qtlpvl: Test of Pleiotropy vs. Close Linkage of QTL

context("tests of det in Rcpp")

test_that("tests of det",{
  set.seed(74695464)
  n <- 500
  p <- 50

  X <- matrix(rnorm(n*p), n, p)
  
  det.r <- determinant(crossprod(X))$modulus
  det.c <- det_AtA(X)
  expect_equal(det.r, det.c,  check.attributes=FALSE)
})