test_that("COLOC_run works", {
gwas.qtl_paths <- catalogueR::eQTLcatalogue_example_queries(
fnames = c("BST1__Alasoo_2018.macrophage_IFNg+Salmonella.tsv.gz",
"BST1__Alasoo_2018.macrophage_IFNg.tsv.gz")
)
for(x in names(gwas.qtl_paths)){
message(x)
dat <- data.table::fread(gwas.qtl_paths[x], nThread = 1)
testthat::expect_true(methods::is(dat,"data.table"))
testthat::expect_gte(nrow(dat),30000)
}
coloc_QTLs <- catalogueR::COLOC_run(gwas.qtl_paths = gwas.qtl_paths)
testthat::expect_true(methods::is(coloc_QTLs, "data.table"))
testthat::expect_gte(nrow(coloc_QTLs), 70000)
testthat::expect_gte(nrow(subset(coloc_QTLs, PP.H4 > .1)), 9300)
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.