tests/testthat/test-api_germline_gwas_results_queries.R

In CRI-iAtlas/iatlas-graphql-client: iatlas API Client

test_that("query_germline_gwas_results",{
  expected_columns <- c(
    "dataset_name",
    "dataset_display",
    "feature_name",
    "feature_display",
    "feature_germline_module",
    "feature_germline_category",
    "snp_name",
    "snp_rsid",
    "snp_chr",
    "snp_bp",
    "p_value",
    "maf"
  )
  result1 <- query_germline_gwas_results(
    datasets = "TCGA",
    feature = "Attractors_G_SIGLEC9",
    snps = "2:206655924:C:T",
    query_dir = query_dir
  )

  expect_named(result1, expected_columns)
  expect_true(nrow(result1) > 0)


  result2 <- query_germline_gwas_results(
    datasets = "TCGA",
    min_p_value = 4.24e-25,
    max_p_value = 5.25e-25,
    query_dir = query_dir
  )
  expect_named(result2, expected_columns)
  expect_true(nrow(result2) > 0)
  expect_true(all(result2$p_value >= 4.24e-25))
  expect_true(all(result2$p_value <= 5.25e-25))

  result3 <- query_germline_gwas_results(
    datasets = "TCGA",
    snps = "none",
    query_dir = query_dir
  )
  expect_named(result3, expected_columns)
  expect_equal(nrow(result3), 0)
})