tests/testthat/test-calc_pattern_risk_ratio.R

In tbilab/meToolkit: Build and Customize Phewas Multimorbidity Explorer Apps.

context("test-calc_pattern_risk_ratio")

test_that("Gets expected RR interval", {
  set.seed(42)
  expected_results <- list(
    PE = 1.455625,
    lower = 0.6169498,
    upper = 3.434385
  )
  RR_results <- calc_pattern_risk_ratio(
    pattern_n=120, pattern_snp=5,
    other_n=2000, other_snp=56,
    CI_size = 0.95
  )
  expect_equivalent(RR_results, expected_results, tolerance = .002)
})