tests/testthat/test-InputPrimaryUnCor.R

In CPBayes: Bayesian Meta Analysis for Studying Cross-Phenotype Genetic Associations

context("Primary Inputs of function for Uncorrelated version")

test_that("Absent of BetaHat and SE will throw error", {
  expect_error(cpbayes_uncor(1:10), "BetaHat or SE vector is missing!")
  expect_error(cpbayes_uncor(Burnin=10), "BetaHat or SE vector is missing!")
})

test_that("Throws error if BetaHat is not a numeric vector of length more than 1 and no missing value", {
  expect_error(cpbayes_uncor(1, 1:10), "Number of elements in the BetaHat vector must be more than 1!")
  expect_error(cpbayes_uncor(NA, 1:10), "BetaHat must be a numeric vector.")
  expect_error(cpbayes_uncor(c(NA,1), 1:10), "BetaHat for one or more phenotypes are missing!")
  expect_error(cpbayes_uncor("AB", 1:10), "BetaHat must be a numeric vector.")
  expect_error(cpbayes_uncor(c("A", 1), 1:10), "BetaHat must be a numeric vector.")
  expect_error(cpbayes_uncor(matrix(0,2,2), 1:10), "BetaHat must be a vector.")
  expect_error(cpbayes_uncor(data.frame(rep(1:10)), 1:10), "BetaHat must be a vector.")
  expect_warning(cpbayes_uncor(matrix(1:10, 10, 1), 1:10), "BetaHat is a matrix!")
})

test_that("Throws error if SE is not a positive numeric vector of length more than 1 and no missing value", {
  expect_error(cpbayes_uncor(1:10, 1), "Number of elements in the SE vector must be more than 1!")
  expect_error(cpbayes_uncor(1:10, NA), "SE must be a numeric vector.")
  expect_error(cpbayes_uncor(1:10, c(NA,1)), "SE for one or more phenotypes are missing!")
  expect_error(cpbayes_uncor(1:10, "AB"), "SE must be a numeric vector.")
  expect_error(cpbayes_uncor(1:10, c("A", 1)), "SE must be a numeric vector.")
  expect_error(cpbayes_uncor(1:10, matrix(0,2,2)), "SE must be a vector.")
  expect_error(cpbayes_uncor(1:10, data.frame(rep(1:10))), "SE must be a vector.")
  expect_error(cpbayes_uncor(1:10, c(-1, 2,3)), "One or more elements in the SE vector are not positive!")
  expect_warning(cpbayes_uncor(1:10, matrix(1:10, 10, 1)), "SE is a matrix!")
})

test_that("Throws error if Beta and SE are not of same length", {
  expect_error(cpbayes_uncor(1:10, 1:15), "BetaHat and SE vectors must have the same number of elements!")
  expect_error(cpbayes_uncor(1:100, 1:15), "BetaHat and SE vectors must have the same number of elements!")
})