Returns the active status of genesets, which are specified by
their collection,name compound keys. This function is vectorized and
supports query of multiple gene sets at a time. If a requested
collection,name gene set doesn't exist, this throws an error.
is.active(x, i, j)logical indicating if geneset is active. throws an error if
any requested geneset does not exist in x.
dge.stats <- exampleDgeResult()
y <- exampleExpressionSet(do.voom = FALSE)
gdb <- conform(exampleGeneSetDb(), y, min.gs.size = 10)
#> Warning: Deactivating 5 gene sets because conformation of GeneSetDb to the target creates gene sets smaller than 10 or greater than Inf
# size 9 geneset:
geneSet(gdb, "c2", "BYSTRYKH_HEMATOPOIESIS_STEM_CELL_IL3RA")
#> collection name active N n feature_id
#> 1 c2 BYSTRYKH_HEMATOPOIESIS_STEM_CELL_IL3RA FALSE 9 6 10134
#> 2 c2 BYSTRYKH_HEMATOPOIESIS_STEM_CELL_IL3RA FALSE 9 6 29950
#> 3 c2 BYSTRYKH_HEMATOPOIESIS_STEM_CELL_IL3RA FALSE 9 6 4297
#> 4 c2 BYSTRYKH_HEMATOPOIESIS_STEM_CELL_IL3RA FALSE 9 6 5879
#> 5 c2 BYSTRYKH_HEMATOPOIESIS_STEM_CELL_IL3RA FALSE 9 6 7423
#> 6 c2 BYSTRYKH_HEMATOPOIESIS_STEM_CELL_IL3RA FALSE 9 6 896
is.active(gdb, "c2", "BYSTRYKH_HEMATOPOIESIS_STEM_CELL_IL3RA")
#> [1] FALSE
# geneset with >100 genes
is.active(gdb, "c7", "GSE3982_MAC_VS_NEUTROPHIL_LPS_STIM_DN")
#> [1] TRUE