Contents

In the mosbi package, similarities between biclusters are computed using different possible similarity metrics.

This vignette gives an overview about the implemented metrics.

library(mosbi)

The following similarity metrics are currently implemented:

# Bray-Curtis similarity
bray_curtis <- function(s1, s2, overlap) {
    return(((2 * overlap) / (s1 + s2)))
}

# Jaccard index
jaccard <- function(s1, s2, overlap) {
    return(((overlap) / (s1 + s2 - overlap)))
}

# overlap coefficient
overlap <- function(s1, s2, overlap) {
    return((overlap / min(s1, s2)))
}

# Fowlkes–Mallows index
folkes_mallows <- function(s1, s2, overlap) {
    tp <- choose(overlap, 2)
    fp <- choose(s1 - overlap, 2)
    fn <- choose(s2 - overlap, 2)

    return(sqrt((tp / (tp + fp)) * (tp / (tp + fn))))
}

The behavior of the similarity metrics will be evaluated for two scenarios:

# Scenario 1 - two biclusters of the same size
size1_1 <- rep(1000, 1000)
size2_1 <- rep(1000, 1000)
overlap_1 <- seq(1, 1000)

# Scenario 2 - two biclusters one of size 500, the other of size 1000
size1_2 <- rep(1000, 500)
size2_2 <- rep(500, 500)
overlap_2 <- seq(1, 500)

Two biclusters of the same size:

plot(overlap_1, bray_curtis(size1_1, size2_1, overlap_1),
    col = "red", type = "l", xlab = "Overlap", ylab = "Similarity", 
    ylim = c(0, 1)
)
lines(overlap_1, jaccard(size1_1, size2_1, overlap_1), col = "blue")
lines(overlap_1, overlap(size1_1, size2_1, overlap_1), col = "green", lty = 2)
lines(overlap_1, folkes_mallows(size1_1, size2_1, overlap_1), col = "orange")
legend(
    x = .8, legend = c("Bray-Curtis", "Jaccard", "Overlap", "Fowlkes–Mallows"),
    col = c("red", "blue", "green", "orange"),
    lty = 1, cex = 0.8, title = "Similarity metrics"
)

Two biclusters of different sizes:

plot(overlap_2, bray_curtis(size1_2, size2_2, overlap_2),
    col = "red", type = "l", xlab = "Overlap", ylab = "Similarity", 
    ylim = c(0, 1)
)
lines(overlap_2, jaccard(size1_2, size2_2, overlap_2), col = "blue")
lines(overlap_2, overlap(size1_2, size2_2, overlap_2), col = "green")
lines(overlap_2, folkes_mallows(size1_2, size2_2, overlap_2), col = "orange")
legend(
    x = .8, legend = c("Bray-Curtis", "Jaccard", "Overlap", "Fowlkes–Mallows"),
    col = c("red", "blue", "green", "orange"),
    lty = 1, cex = 0.8, title = "Similarity metrics"
)

1 Session Info

sessionInfo()
#> R version 4.2.0 RC (2022-04-21 r82226)
#> Platform: x86_64-pc-linux-gnu (64-bit)
#> Running under: Ubuntu 20.04.4 LTS
#> 
#> Matrix products: default
#> BLAS:   /home/biocbuild/bbs-3.16-bioc/R/lib/libRblas.so
#> LAPACK: /home/biocbuild/bbs-3.16-bioc/R/lib/libRlapack.so
#> 
#> 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       
#> 
#> attached base packages:
#> [1] stats     graphics  grDevices utils     datasets  methods   base     
#> 
#> other attached packages:
#> [1] mosbi_1.3.0      BiocStyle_2.25.0
#> 
#> loaded via a namespace (and not attached):
#>  [1] tidyselect_1.1.2      modeltools_0.2-23     xfun_0.30            
#>  [4] bslib_0.3.1           purrr_0.3.4           lattice_0.20-45      
#>  [7] flexclust_1.4-1       generics_0.1.2        colorspace_2.0-3     
#> [10] vctrs_0.4.1           htmltools_0.5.2       stats4_4.2.0         
#> [13] yaml_2.3.5            utf8_1.2.2            rlang_1.0.2          
#> [16] jquerylib_0.1.4       pillar_1.7.0          DBI_1.1.2            
#> [19] glue_1.6.2            BiocGenerics_0.43.0   RColorBrewer_1.1-3   
#> [22] QUBIC_1.25.0          lifecycle_1.0.1       stringr_1.4.0        
#> [25] munsell_0.5.0         gtable_0.3.0          evaluate_0.15        
#> [28] Biobase_2.57.0        knitr_1.38            fastmap_1.1.0        
#> [31] parallel_4.2.0        class_7.3-20          fansi_1.0.3          
#> [34] biclust_2.0.3         highr_0.9             Rcpp_1.0.8.3         
#> [37] scales_1.2.0          BiocManager_1.30.17   additivityTests_1.1-4
#> [40] RcppParallel_5.1.5    magick_2.7.3          jsonlite_1.8.0       
#> [43] fabia_2.43.0          ggplot2_3.3.5         digest_0.6.29        
#> [46] stringi_1.7.6         dplyr_1.0.8           bookdown_0.26        
#> [49] grid_4.2.0            BH_1.78.0-0           cli_3.3.0            
#> [52] tools_4.2.0           magrittr_2.0.3        sass_0.4.1           
#> [55] tibble_3.1.6          tidyr_1.2.0           pkgconfig_2.0.3      
#> [58] crayon_1.5.1          isa2_0.3.5            MASS_7.3-57          
#> [61] ellipsis_0.3.2        assertthat_0.2.1      rmarkdown_2.14       
#> [64] R6_2.5.1              igraph_1.3.1          compiler_4.2.0