netboost 2.8.0
The netboost package, implements a three-step dimension reduction technique. First, a boosting-based filter is combined with the topological overlap measure to identify the essential edges of the network. Second, sparse hierarchical clustering is applied on the selected edges to identify modules and finally module information is aggregated by the first principal components. The primary analysis is then carried out on these summary measures instead of the original data.
The package comes with an example dataset included. We import the acute myeloid leukemia patient data from The Cancer Genome Atlas public domain database. The dataset consists of one thousand DNA methylation sites and gene expression levels on chromosome 18 for 80 patients.
require("netboost")
## Loading required package: netboost
## Warning: replacing previous import 'utils::findMatches' by
## 'S4Vectors::findMatches' when loading 'AnnotationDbi'
##
## netboost 2.8.0 loadedDefault CPU cores: 1
data("tcga_aml_meth_rna_chr18", package = "netboost")
dim(tcga_aml_meth_rna_chr18)
## [1] 80 500
The netboost()
function integrates all major analysis steps and generates
multiple plots. In this step we also set analysis parameters:
stepno
defines the number of boosting steps taken
soft_power
(if null, automatically chosen) the exponent in the transformation
of the correlation
min_cluster_size
the minimal size of clusters, n_pc
the number of maximally
computed principal components
scale
if data should be scaled and centered prior to analysis
ME_diss_thres
defines the merging threshold for identified clusters.
For details on the options please see ?netboost
and the corresponding paper
Schlosser et al. 2020.
results <- netboost(datan = tcga_aml_meth_rna_chr18, stepno = 20L,
soft_power = 3L, min_cluster_size = 10L, n_pc = 2, scale = TRUE, ME_diss_thres = 0.25)
## idx: 1 (0.2%) - Tue Apr 25 17:55:03 2023
##
## Netboost extracted 10 modules (including background) with an average size of 17.5555555555556 (excluding background) from Tree 1.
##
## Netboost detected 9 modules and 1 background modules in 1 trees resulting in 15 aggreagate measures.
## Average size of the modules was 17.5555555555556.
## 342 of 500 features (68.4%) were not assigned to modules.
For each detected independent tree in the dataset (here one) the first graph shows a dendrogram of initial modules and at which level they are merged, the second graph a module dendrogram after merging and the third the dendrogram of features including the module-color-code.
results
contains the dendrograms (dendros), feature identifier (names) matched
to module assignment (colors), the aggregated dataset (MEs), the rotation matrix
to compute the aggregated dataset (rotation) and the proportion of variance
explained by the aggregate measures (var_explained).
Dependent on the minimum proportion of variance explained set in the
netboost()
call (default 0.5) up to n_pc
principal components are exported.
names(results)
## [1] "dendros" "names" "colors" "MEs"
## [5] "rotation" "var_explained" "filter"
colnames(results$MEs)
## [1] "ME0_1_pc1" "ME0_1_pc2" "ME7_pc1" "ME1_pc1" "ME1_pc2" "ME2_pc1"
## [7] "ME2_pc2" "ME8_pc1" "ME5_pc1" "ME3_pc1" "ME3_pc2" "ME4_pc1"
## [13] "ME4_pc2" "ME9_pc1" "ME6_pc1"
As you see for most modules the first principal component already explained more than 50% of the variance in the original features of this module. ME0_X_pcY denotes the background module (unclustered features) of the independent tree X.
Explained variance is reported by a matrix for the first n_pc
principal
components. Here we list the first 5 modules:
results$var_explained[,1:5]
## ME0_1 ME7 ME1 ME2 ME8
## PC1 0.06700469 0.61004480 0.4403646 0.49237958 0.59426699
## PC2 0.05502278 0.07484705 0.1174992 0.07779346 0.08341562
results$colors
use a numeric coding for the modules which matches their module
name. To list features of module ME8 we can extract them by:
results$names[results$colors==8]
## [1] "cg00027037" "cg00034852" "cg00220661" "cg00228017" "cg00366917"
## [6] "cg00430895" "cg00474194" "cg00481457" "cg00511081" "cg00539368"
## [11] "cg00576121" "cg00615915" "cg00736530" "cg00917154" "cg00940278"
## [16] "cg00955482"
The final dendrogram including all trees can be plotted including labels (results$names
) for
individual features. colorsrandom
controls if module-color matching should be
randomized to get a clearly differentiable pattern of the potentially many
modules. Labels are only suitable in applications with few features or with a
appropriately large pdf device.
set.seed(123)
nb_plot_dendro(nb_summary = results, labels = FALSE, colorsrandom = TRUE)
Next the primary analysis on the aggregated dataset (results$MEs
) can be
computed.
We also implemented a convenience function to transfer a clustering to a new
dataset. Here, we transfer the clustering to the same dataset resulting in
identical aggregate measures.
ME_transfer <- nb_transfer(nb_summary = results,
new_data = tcga_aml_meth_rna_chr18, scale = TRUE)
all(round(results$MEs, 10) == round(ME_transfer, 10))
## [1] TRUE
Netboost now also has a fully non-parametric implementation. Code is not run here to showcase the multicore option (Bioconductor vignette builder does not allow for multicore execution). Adjust cores
to your machine:
#results <- netboost(datan = tcga_aml_meth_rna_chr18,cores=10L,
#soft_power = 3L, min_cluster_size = 10L, n_pc = 2, qc_plot = FALSE,
#filter_method = "spearman", robust_PCs = TRUE, method = "spearman")
sessionInfo()
## R version 4.3.0 RC (2023-04-13 r84269)
## Platform: x86_64-pc-linux-gnu (64-bit)
## Running under: Ubuntu 22.04.2 LTS
##
## Matrix products: default
## BLAS: /home/biocbuild/bbs-3.17-bioc/R/lib/libRblas.so
## LAPACK: /usr/lib/x86_64-linux-gnu/lapack/liblapack.so.3.10.0
##
## locale:
## [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C
## [3] LC_TIME=en_GB LC_COLLATE=C
## [5] LC_MONETARY=en_US.UTF-8 LC_MESSAGES=en_US.UTF-8
## [7] LC_PAPER=en_US.UTF-8 LC_NAME=C
## [9] LC_ADDRESS=C LC_TELEPHONE=C
## [11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C
##
## time zone: America/New_York
## tzcode source: system (glibc)
##
## attached base packages:
## [1] stats graphics grDevices utils datasets methods base
##
## other attached packages:
## [1] netboost_2.8.0 BiocStyle_2.28.0
##
## loaded via a namespace (and not attached):
## [1] DBI_1.1.3 bitops_1.0-7 gridExtra_2.3
## [4] rlang_1.1.0 magrittr_2.0.3 matrixStats_0.63.0
## [7] compiler_4.3.0 RSQLite_2.3.1 png_0.1-8
## [10] vctrs_0.6.2 stringr_1.5.0 pkgconfig_2.0.3
## [13] crayon_1.5.2 fastmap_1.1.1 backports_1.4.1
## [16] magick_2.7.4 XVector_0.40.0 utf8_1.2.3
## [19] rmarkdown_2.21 preprocessCore_1.62.0 bit_4.0.5
## [22] xfun_0.39 zlibbioc_1.46.0 cachem_1.0.7
## [25] GenomeInfoDb_1.36.0 jsonlite_1.8.4 blob_1.2.4
## [28] highr_0.10 parallel_4.3.0 cluster_2.1.4
## [31] R6_2.5.1 bslib_0.4.2 stringi_1.7.12
## [34] rpart_4.1.19 jquerylib_0.1.4 Rcpp_1.0.10
## [37] bookdown_0.33 iterators_1.0.14 knitr_1.42
## [40] WGCNA_1.72-1 base64enc_0.1-3 R.utils_2.12.2
## [43] IRanges_2.34.0 Matrix_1.5-4 splines_4.3.0
## [46] nnet_7.3-18 tidyselect_1.2.0 rstudioapi_0.14
## [49] yaml_2.3.7 doParallel_1.0.17 codetools_0.2-19
## [52] lattice_0.21-8 tibble_3.2.1 Biobase_2.60.0
## [55] KEGGREST_1.40.0 evaluate_0.20 foreign_0.8-84
## [58] survival_3.5-5 RcppParallel_5.1.7 Biostrings_2.68.0
## [61] pillar_1.9.0 BiocManager_1.30.20 checkmate_2.1.0
## [64] foreach_1.5.2 stats4_4.3.0 generics_0.1.3
## [67] RCurl_1.98-1.12 S4Vectors_0.38.0 ggplot2_3.4.2
## [70] munsell_0.5.0 scales_1.2.1 glue_1.6.2
## [73] Hmisc_5.0-1 tools_4.3.0 data.table_1.14.8
## [76] fastcluster_1.2.3 grid_4.3.0 impute_1.74.0
## [79] AnnotationDbi_1.62.0 colorspace_2.1-0 GenomeInfoDbData_1.2.10
## [82] htmlTable_2.4.1 Formula_1.2-5 cli_3.6.1
## [85] fansi_1.0.4 dplyr_1.1.2 gtable_0.3.3
## [88] R.methodsS3_1.8.2 dynamicTreeCut_1.63-1 sass_0.4.5
## [91] digest_0.6.31 BiocGenerics_0.46.0 htmlwidgets_1.6.2
## [94] memoise_2.0.1 htmltools_0.5.5 R.oo_1.25.0
## [97] lifecycle_1.0.3 httr_1.4.5 GO.db_3.17.0
## [100] bit64_4.0.5
warnings()