Bioconductor has a rich ecosystem of metadata around packages, usage, and build status. This package is a simple collection of functions to access that metadata from R in a tidy data format. The goal is to expose metadata for data mining and value-added functionality such as package searching, text mining, and analytics on packages.
BiocPkgTools 1.10.2
Bioconductor has a rich ecosystem of metadata around packages, usage, and build status. This package is a simple collection of functions to access that metadata from R in a tidy data format. The goal is to expose metadata for data mining and value-added functionality such as package searching, text mining, and analytics on packages.
Functionality includes access to :
The Bioconductor build reports are available online as HTML pages.
However, they are not very computable.
The biocBuildReport function does some heroic parsing of the HTML
to produce a tidy data.frame for further processing in R.
library(BiocPkgTools)
head(biocBuildReport())## # A tibble: 6 x 9
## pkg author version git_last_commit git_last_commit_date node stage result
## <chr> <chr> <chr> <chr> <dttm> <chr> <chr> <chr>
## 1 a4 Laure C~ 1.40.0 f5e7837 2021-05-19 00:00:00 nebb~ inst~ OK
## 2 a4 Laure C~ 1.40.0 f5e7837 2021-05-19 00:00:00 nebb~ buil~ OK
## 3 a4 Laure C~ 1.40.0 f5e7837 2021-05-19 00:00:00 nebb~ chec~ OK
## 4 a4 Laure C~ 1.40.0 f5e7837 2021-05-19 00:00:00 toka~ inst~ OK
## 5 a4 Laure C~ 1.40.0 f5e7837 2021-05-19 00:00:00 toka~ buil~ OK
## 6 a4 Laure C~ 1.40.0 f5e7837 2021-05-19 00:00:00 toka~ chec~ OK
## # ... with 1 more variable: bioc_version <chr>Because developers may be interested in a quick view of their own
packages, there is a simple function, problemPage, to produce an HTML report of
the build status of packages matching a given author regex supplied to the authorPattern argument.
The default is to report only “problem” build statuses (ERROR, WARNING).
problemPage(authorPattern = "V.*Carey")In similar fashion, maintainers of packages that have many downstream packages that
depend on them may wish to check that a change they introduced hasn’t suddenly broken
a large number of these. You can use the dependsOn argument to produce the summary report
of those packages that “depend on” the given package.
problemPage(dependsOn = "limma")When run in an interactive environment, the problemPage function
will open a browser window for user interaction. Note that if you want
to include all your package results, not just the broken ones, simply
specify includeOK = TRUE.
Bioconductor supplies download stats for all packages. The biocDownloadStats
function grabs all available download stats for all packages in all
Experiment Data, Annotation Data, and Software packages. The results
are returned as a tidy data.frame for further analysis.
head(biocDownloadStats())## # A tibble: 6 x 7
## Package Year Month Nb_of_distinct_IPs Nb_of_downloads repo Date
## <chr> <int> <chr> <int> <int> <chr> <date>
## 1 ABarray 2021 Jan 60 127 Software 2021-01-01
## 2 ABarray 2021 Feb 51 139 Software 2021-02-01
## 3 ABarray 2021 Mar 77 146 Software 2021-03-01
## 4 ABarray 2021 Apr 75 145 Software 2021-04-01
## 5 ABarray 2021 May 73 122 Software 2021-05-01
## 6 ABarray 2021 Jun 56 101 Software 2021-06-01The download statistics reported are for all available versions of a package. There are no separate, publicly available statistics broken down by version.
The majority of Bioconductor Software packages are also available through other channels
such as Anaconda, who also provided download statistics for packages installed from
their repositories. Access to these counts is provided by the anacondaDownloadStats function:
head(anacondaDownloadStats())## # A tibble: 6 x 7
## Package Year Month Nb_of_distinct_IPs Nb_of_downloads repo Date
## <chr> <chr> <chr> <int> <dbl> <chr> <date>
## 1 ABAData 2018 Apr NA 8 Anaconda 2018-04-01
## 2 ABAData 2018 Aug NA 5 Anaconda 2018-08-01
## 3 ABAData 2018 Dec NA 133 Anaconda 2018-12-01
## 4 ABAData 2018 Jul NA 6 Anaconda 2018-07-01
## 5 ABAData 2018 Jun NA 18 Anaconda 2018-06-01
## 6 ABAData 2018 Mar NA 13 Anaconda 2018-03-01Note that Anaconda do not provide counts for distinct IP addresses, but this column is included for compatibility with the Bioconductor count tables.
The R DESCRIPTION file contains a plethora of information regarding package
authors, dependencies, versions, etc. In a repository such as Bioconductor, these
details are available in bulk for all inclucded packages. The biocPkgList returns
a data.frame with a row for each package. Tons of information are avaiable, as
evidenced by the column names of the results.
bpi = biocPkgList()
colnames(bpi)## [1] "Package" "Version" "Depends"
## [4] "Suggests" "License" "MD5sum"
## [7] "NeedsCompilation" "Title" "Description"
## [10] "biocViews" "Author" "Maintainer"
## [13] "git_url" "git_branch" "git_last_commit"
## [16] "git_last_commit_date" "Date/Publication" "source.ver"
## [19] "win.binary.ver" "mac.binary.ver" "vignettes"
## [22] "vignetteTitles" "hasREADME" "hasNEWS"
## [25] "hasINSTALL" "hasLICENSE" "Rfiles"
## [28] "dependencyCount" "Imports" "Enhances"
## [31] "dependsOnMe" "VignetteBuilder" "suggestsMe"
## [34] "LinkingTo" "Archs" "URL"
## [37] "SystemRequirements" "BugReports" "importsMe"
## [40] "Video" "linksToMe" "License_restricts_use"
## [43] "PackageStatus" "organism" "OS_type"
## [46] "License_is_FOSS"Some of the variables are parsed to produce list columns.
head(bpi)## # A tibble: 6 x 46
## Package Version Depends Suggests License MD5sum NeedsCompilation Title
## <chr> <chr> <list> <list> <chr> <chr> <chr> <chr>
## 1 a4 1.40.0 <chr [5]> <chr [6]> GPL-3 e2ac1~ no Autom~
## 2 a4Base 1.40.0 <chr [2]> <chr [4]> GPL-3 aa474~ no Autom~
## 3 a4Classif 1.40.0 <chr [2]> <chr [4]> GPL-3 18794~ no Autom~
## 4 a4Core 1.40.0 <chr [1]> <chr [2]> GPL-3 8a703~ no Autom~
## 5 a4Preproc 1.40.0 <chr [1]> <chr [4]> GPL-3 4e2db~ no Autom~
## 6 a4Reporting 1.40.0 <chr [1]> <chr [2]> GPL-3 db956~ no Autom~
## # ... with 38 more variables: Description <chr>, biocViews <list>,
## # Author <list>, Maintainer <list>, git_url <chr>, git_branch <chr>,
## # git_last_commit <chr>, git_last_commit_date <chr>, Date/Publication <chr>,
## # source.ver <chr>, win.binary.ver <chr>, mac.binary.ver <chr>,
## # vignettes <list>, vignetteTitles <list>, hasREADME <chr>, hasNEWS <chr>,
## # hasINSTALL <chr>, hasLICENSE <chr>, Rfiles <list>, dependencyCount <chr>,
## # Imports <list>, Enhances <list>, dependsOnMe <list>, ...As a simple example of how these columns can be used, extracting
the importsMe column to find the packages that import the
GEOquery package.
require(dplyr)
bpi = biocPkgList()
bpi %>%
filter(Package=="GEOquery") %>%
pull(importsMe) %>%
unlist()## [1] "bigmelon" "ChIPXpress" "coexnet"
## [4] "conclus" "crossmeta" "DExMA"
## [7] "EGAD" "GAPGOM" "MACPET"
## [10] "minfi" "MoonlightR" "phantasus"
## [13] "recount" "SRAdb" "BeadArrayUseCases"
## [16] "GSE13015" "geneExpressionFromGEO" "MetaIntegrator"For the end user of Bioconductor, an analysis often starts with finding a
package or set of packages that perform required tasks or are tailored
to a specific operation or data type. The biocExplore() function
implements an interactive bubble visualization with filtering based on
biocViews terms. Bubbles are sized based on download statistics. Tooltip
and detail-on-click capabilities are included. To start a local session:
biocExplore()The Bioconductor ecosystem is built around the concept of interoperability
and dependencies. These interdependencies are available as part of the
biocPkgList() output. The BiocPkgTools provides some convenience
functions to convert package dependencies to R graphs. A modular approach leads
to the following workflow.
data.frame of dependencies using buildPkgDependencyDataFrame.igraph object from the dependency data frame using buildPkgDependencyIgraphigraph functionality to perform arbitrary network operations.
Convenience functions, inducedSubgraphByPkgs and subgraphByDegree are available.A dependency graph for all of Bioconductor is a starting place.
library(BiocPkgTools)
dep_df = buildPkgDependencyDataFrame()
g = buildPkgDependencyIgraph(dep_df)
g## IGRAPH 40f807f DN-- 3672 34805 --
## + attr: name (v/c), edgetype (e/c)
## + edges from 40f807f (vertex names):
## [1] a4 ->a4Base a4 ->a4Preproc
## [3] a4 ->a4Classif a4 ->a4Core
## [5] a4 ->a4Reporting a4Base ->a4Preproc
## [7] a4Base ->a4Core a4Classif ->a4Core
## [9] a4Classif ->a4Preproc ABAEnrichment->R
## [11] abseqR ->R ABSSeq ->R
## [13] ABSSeq ->methods acde ->R
## [15] acde ->boot ACE ->R
## + ... omitted several edgeslibrary(igraph)
head(V(g))## + 6/3672 vertices, named, from 40f807f:
## [1] a4 a4Base a4Classif ABAEnrichment abseqR
## [6] ABSSeqhead(E(g))## + 6/34805 edges from 40f807f (vertex names):
## [1] a4 ->a4Base a4 ->a4Preproc a4 ->a4Classif
## [4] a4 ->a4Core a4 ->a4Reporting a4Base->a4PreprocSee inducedSubgraphByPkgs and subgraphByDegree to produce
subgraphs based on a subset of packages.
See the igraph documentation for more detail on graph analytics, setting vertex and edge attributes, and advanced subsetting.
The visNetwork package is a nice interactive visualization tool that implements graph plotting in a browser. It can be integrated into shiny applications. Interactive graphs can also be included in Rmarkdown documents (see vignette)
igraph_network = buildPkgDependencyIgraph(buildPkgDependencyDataFrame())The full dependency graph is really not that informative to look at, though doing so is possible. A common use case is to visualize the graph of dependencies “centered” on a package of interest. In this case, I will focus on the GEOquery package.
igraph_geoquery_network = subgraphByDegree(igraph_network, "GEOquery")The subgraphByDegree() function returns all nodes and connections within
degree of the named package; the default degree is 1.
The visNework package can plot igraph objects directly, but more flexibility
is offered by first converting the graph to visNetwork form.
library(visNetwork)
data <- toVisNetworkData(igraph_geoquery_network)The next few code chunks highlight just a few examples of the visNetwork capabilities, starting with a basic plot.
visNetwork(nodes = data$nodes, edges = data$edges, height = "500px")For fun, we can watch the graph stabilize during drawing, best viewed interactively.
visNetwork(nodes = data$nodes, edges = data$edges, height = "500px") %>%
visPhysics(stabilization=FALSE)Add arrows and colors to better capture dependencies.
data$edges$color='lightblue'
data$edges[data$edges$edgetype=='Imports','color']= 'red'
data$edges[data$edges$edgetype=='Depends','color']= 'green'
visNetwork(nodes = data$nodes, edges = data$edges, height = "500px") %>%
visEdges(arrows='from') Add a legend.
ledges <- data.frame(color = c("green", "lightblue", "red"),
label = c("Depends", "Suggests", "Imports"), arrows =c("from", "from", "from"))
visNetwork(nodes = data$nodes, edges = data$edges, height = "500px") %>%
visEdges(arrows='from') %>%
visLegend(addEdges=ledges)[Work in progress]
The biocViews package is a small ontology of terms describing Bioconductor packages. This is a work-in-progress section, but here is a small example of plotting the biocViews graph.
library(biocViews)
data(biocViewsVocab)
biocViewsVocab## A graphNEL graph with directed edges
## Number of Nodes = 490
## Number of Edges = 489library(igraph)
g = igraph.from.graphNEL(biocViewsVocab)
library(visNetwork)
gv = toVisNetworkData(g)
visNetwork(gv$nodes, gv$edges, width="100%") %>%
visIgraphLayout(layout = "layout_as_tree", circular=TRUE) %>%
visNodes(size=20) %>%
visPhysics(stabilization=FALSE)The dependency burden of a package, namely the amount of functionality that
a given package is importing, is an important parameter to take into account
during package development. A package may break because one or more of its
dependencies have changed the part of the API our package is importing or
this part has even broken. For this reason, it may be useful for package
developers to quantify the dependency burden of a given package. To do that
we should first gather all dependency information using the function
buildPkgDependencyDataFrame() but setting the arguments to work with
packages in Bioconductor and CRAN and dependencies categorised as Depends
or Imports, which are the ones installed by default for a given package.
library(BiocPkgTools)
depdf <- buildPkgDependencyDataFrame(repo=c("BioCsoft", "CRAN"),
dependencies=c("Depends", "Imports"))
depdf## # A tibble: 124,260 x 3
## Package dependency edgetype
## <chr> <chr> <chr>
## 1 a4 a4Base Depends
## 2 a4 a4Preproc Depends
## 3 a4 a4Classif Depends
## 4 a4 a4Core Depends
## 5 a4 a4Reporting Depends
## 6 a4Base a4Preproc Depends
## 7 a4Base a4Core Depends
## 8 a4Classif a4Core Depends
## 9 a4Classif a4Preproc Depends
## 10 ABAEnrichment R Depends
## # ... with 124,250 more rowsFinally, we call the function pkgDepMetrics() to obtain different metrics
on the dependency burden of a package we want to analyze, in the case below,
the package BiocPkgTools itself:
pkgDepMetrics("BiocPkgTools", depdf)## ImportedAndUsed Exported Usage DepOverlap DepGainIfExcluded
## utils 1 247 0.40 0.01 0
## rlang 4 476 0.84 0.02 0
## tools 1 118 0.85 0.01 0
## graph 1 116 0.86 0.08 0
## igraph 9 782 1.15 0.11 4
## RBGL 1 77 1.30 0.09 0
## htmltools 1 75 1.33 0.08 0
## xml2 1 66 1.52 0.01 0
## tidyr 1 62 1.61 0.24 1
## magrittr 1 41 2.44 0.01 0
## DT 1 39 2.56 0.23 5
## dplyr 8 285 2.81 0.23 0
## httr 3 91 3.30 0.10 0
## tidyselect 1 25 4.00 0.09 0
## tibble 2 44 4.55 0.17 0
## jsonlite 1 17 5.88 0.01 0
## htmlwidgets 1 14 7.14 0.11 0
## rvest 3 40 7.50 0.34 2
## gh 1 10 10.00 0.17 3
## stringr 5 49 10.20 0.08 0
## BiocManager 1 5 20.00 0.02 0
## BiocFileCache NA 29 NA 0.51 10
## biocViews NA 31 NA 0.18 6
## readr NA 114 NA 0.32 5
## rex NA 16 NA 0.02 1In this resulting table, rows correspond to dependencies and columns provide the following information:
ImportedAndUsed: number of functionality calls imported and used in
the package.Exported: number of functionality calls exported by the dependency.Usage: (ImportedAndUsedx 100) / Exported. This value provides an
estimate of what fraction of the functionality of the dependency is
actually used in the given package.DepOverlap: Similarity between the dependency graph structure of the
given package and the one of the dependency in the corresponding row,
estimated as the Jaccard index
between the two sets of vertices of the corresponding graphs. Its values
goes between 0 and 1, where 0 indicates that no dependency is shared, while
1 indicates that the given package and the corresponding dependency depend
on an identical subset of packages.DepGainIfExcluded: The ‘dependency gain’ (decrease in the total number
of dependencies) that would be obtained if this package was excluded
from the list of direct dependencies.The reported information is ordered by the Usage column to facilitate the
identification of dependencies for which the analyzed package is using a small
fraction of their functionality and therefore, it could be easier remove them.
To aid in that decision, the column DepOverlap reports the overlap of the
dependency graph of each dependency with the one of the analyzed package. Here
a value above, e.g., 0.5, could, albeit not necessarily, imply that removing
that dependency could substantially lighten the dependency burden of the analyzed
package.
An NA value in the ImportedAndUsed column indicates that the function
pkgDepMetrics() could not identify what functionality calls in the analyzed
package are made to the dependency. This may happen because pkgDepMetrics()
has failed to identify the corresponding calls, as it happens with imported
built-in constants such as DNA_BASES from Biostrings, or that although the
given package is importing that dependency, none of its functionality is actually
being used. In such a case, this dependency could be safely removed without
any further change in the analyzed package.
We can find out what actually functionality calls are we importing as follows:
imp <- pkgDepImports("BiocPkgTools")
imp %>% filter(pkg == "DT")## # A tibble: 1 x 2
## pkg fun
## <chr> <chr>
## 1 DT datatablesessionInfo()## R version 4.1.1 (2021-08-10)
## Platform: x86_64-w64-mingw32/x64 (64-bit)
## Running under: Windows Server 2012 R2 x64 (build 9600)
##
## Matrix products: default
##
## locale:
## [1] LC_COLLATE=C
## [2] LC_CTYPE=English_United States.1252
## [3] LC_MONETARY=English_United States.1252
## [4] LC_NUMERIC=C
## [5] LC_TIME=English_United States.1252
##
## attached base packages:
## [1] stats graphics grDevices utils datasets methods base
##
## other attached packages:
## [1] biocViews_1.60.0 visNetwork_2.0.9 igraph_1.2.6
## [4] dplyr_1.0.7 BiocPkgTools_1.10.2 htmlwidgets_1.5.4
## [7] knitr_1.34 BiocStyle_2.20.2
##
## loaded via a namespace (and not attached):
## [1] Rcpp_1.0.7 tidyr_1.1.3 assertthat_0.2.1
## [4] digest_0.6.27 utf8_1.2.2 BiocFileCache_2.0.0
## [7] R6_2.5.1 stats4_4.1.1 RSQLite_2.2.8
## [10] evaluate_0.14 httr_1.4.2 pillar_1.6.2
## [13] rlang_0.4.11 curl_4.3.2 jquerylib_0.1.4
## [16] blob_1.2.2 DT_0.19 RUnit_0.4.32
## [19] rmarkdown_2.11 readr_2.0.1 stringr_1.4.0
## [22] selectr_0.4-2 RCurl_1.98-1.5 bit_4.0.4
## [25] compiler_4.1.1 xfun_0.26 pkgconfig_2.0.3
## [28] BiocGenerics_0.38.0 htmltools_0.5.2 tidyselect_1.1.1
## [31] tibble_3.1.4 bookdown_0.24 XML_3.99-0.8
## [34] fansi_0.5.0 crayon_1.4.1 tzdb_0.1.2
## [37] dbplyr_2.1.1 rappdirs_0.3.3 bitops_1.0-7
## [40] RBGL_1.68.0 jsonlite_1.7.2 lifecycle_1.0.0
## [43] DBI_1.1.1 magrittr_2.0.1 graph_1.70.0
## [46] cli_3.0.1 stringi_1.7.4 cachem_1.0.6
## [49] rex_1.2.0 xml2_1.3.2 bslib_0.3.0
## [52] ellipsis_0.3.2 filelock_1.0.2 generics_0.1.0
## [55] vctrs_0.3.8 gh_1.3.0 tools_4.1.1
## [58] bit64_4.0.5 Biobase_2.52.0 glue_1.4.2
## [61] purrr_0.3.4 hms_1.1.0 parallel_4.1.1
## [64] fastmap_1.1.0 yaml_2.2.1 BiocManager_1.30.16
## [67] rvest_1.0.1 memoise_2.0.0 sass_0.4.0