# Introduction to iSEEde

Kevin Rue-Albrecht1*

1University of Oxford

* kevin.rue-albrecht@imm.ox.ac.uk

#### 24 October 2023

#### Package

iSEEde 1.0.0

# 1 Basics

## 1.1 Install `iSEEde`

is an open-source statistical environment which can be easily modified to enhance its functionality via packages. iSEEde is an
package available via the Bioconductor repository for packages.
can be installed on any operating system from CRAN after which you can install iSEEde by using the following commands in your
session:

    if (!requireNamespace("BiocManager", quietly = TRUE)) {
        install.packages("BiocManager")
    }

    BiocManager::install("iSEEde")

    ## Check that you have a valid Bioconductor installation
    BiocManager::valid()

## 1.2 Required knowledge

iSEEde is based on many other packages and in particular on those that have implemented the infrastructure needed for dealing with omics data and interactive visualisation. That is, packages like SummarizedExperiment, SingleCellExperiment, iSEE and shiny.

If you are asking yourself the question “Where do I start using Bioconductor?” you might be interested in this blog post.

## 1.3 Asking for help

As package developers, we try to explain clearly how to use our packages and in which order to use the functions. But
and `Bioconductor` have a steep learning curve so it is critical to learn where to ask for help. The blog post quoted above mentions some but we would like to highlight the Bioconductor support site as the main resource for getting help: remember to use the `iSEEde` tag and check the older posts. Other alternatives are available such as creating GitHub issues and tweeting. However, please note that if you want to receive help you should adhere to the posting guidelines. It is particularly critical that you provide a small reproducible example and your session information so package developers can track down the source of the error.

## 1.4 Citing `iSEEde`

We hope that iSEEde will be useful for your research. Please use the following information to cite the package and the overall approach. Thank you!

    ## Citation info
    citation("iSEEde")
    #> To cite package 'iSEEde' in publications use:
    #>
    #>   Rue-Albrecht K (2023). _iSEEde: iSEE extension for panels related to differential
    #>   expression analysis_. doi:10.18129/B9.bioc.iSEEde
    #>   <https://doi.org/10.18129/B9.bioc.iSEEde>, R package version 1.0.0,
    #>   <https://bioconductor.org/packages/iSEEde>.
    #>
    #> A BibTeX entry for LaTeX users is
    #>
    #>   @Manual{,
    #>     title = {iSEEde: iSEE extension for panels related to differential expression analysis},
    #>     author = {Kevin Rue-Albrecht},
    #>     year = {2023},
    #>     note = {R package version 1.0.0},
    #>     url = {https://bioconductor.org/packages/iSEEde},
    #>     doi = {10.18129/B9.bioc.iSEEde},
    #>   }

# 2 Quick start to using to `iSEEde`

The example below demonstrates how to use iSEEde functionality in a simple differential expression workflow.

Specifically, the DESeq2 package is used to perform a simple differential expression analysis. Then, the results of one contrast – extracted using the function `DESeq2::results()` function – are embedded into the SummarizedExperiment object `airway` using the function `iSEEde::embedContrastResults()`. This ensures that the differential expression results for that contrast are accessible to specialised interactive panels in iSEE applications.

    library("iSEEde")
    library("airway")
    library("DESeq2")
    library("iSEE")

    # Example data ----

    data("airway")
    airway$dex <- relevel(airway$dex, "untrt")

    dds <- DESeqDataSet(airway, ~ 0 + dex + cell)

    dds <- DESeq(dds)
    res_deseq2 <- results(dds, contrast = list("dextrt", "dexuntrt"))
    head(res_deseq2)
    #> log2 fold change (MLE): dextrt vs dexuntrt
    #> Wald test p-value: dextrt vs dexuntrt
    #> DataFrame with 6 rows and 6 columns
    #>                   baseMean log2FoldChange     lfcSE      stat      pvalue       padj
    #>                  <numeric>      <numeric> <numeric> <numeric>   <numeric>  <numeric>
    #> ENSG00000000003 708.602170     -0.3812539  0.100654 -3.787752 0.000152016 0.00128292
    #> ENSG00000000005   0.000000             NA        NA        NA          NA         NA
    #> ENSG00000000419 520.297901      0.2068127  0.112219  1.842944 0.065337213 0.19646961
    #> ENSG00000000457 237.163037      0.0379205  0.143445  0.264356 0.791505314 0.91141884
    #> ENSG00000000460  57.932633     -0.0881679  0.287142 -0.307054 0.758802543 0.89500551
    #> ENSG00000000938   0.318098     -1.3782416  3.499906 -0.393794 0.693733216         NA

    # iSEE / iSEEde ---

    airway <- embedContrastResults(res_deseq2, airway, name = "dex: trt vs untrt")
    contrastResults(airway)
    #> DataFrame with 63677 rows and 1 column
    #>                   dex: trt vs untrt
    #>                 <iSEEDESeq2Results>
    #> ENSG00000000003 <iSEEDESeq2Results>
    #> ENSG00000000005 <iSEEDESeq2Results>
    #> ENSG00000000419 <iSEEDESeq2Results>
    #> ENSG00000000457 <iSEEDESeq2Results>
    #> ENSG00000000460 <iSEEDESeq2Results>
    #> ...                             ...
    #> ENSG00000273489 <iSEEDESeq2Results>
    #> ENSG00000273490 <iSEEDESeq2Results>
    #> ENSG00000273491 <iSEEDESeq2Results>
    #> ENSG00000273492 <iSEEDESeq2Results>
    #> ENSG00000273493 <iSEEDESeq2Results>

    app <- iSEE(airway, initial = list(
      DETable(ContrastName="dex: trt vs untrt", HiddenColumns = c("baseMean",
        "lfcSE", "stat"), PanelWidth = 4L),
      VolcanoPlot(ContrastName="dex: trt vs untrt", PanelWidth = 4L),
      MAPlot(ContrastName="dex: trt vs untrt", PanelWidth = 4L)
    ))

    if (interactive()) {
        shiny::runApp(app)
    }

# 3 Reproducibility

The iSEEde package (Rue-Albrecht, 2023) was made possible thanks to:

- R (R Core Team, 2023)
- BiocStyle (Oleś, 2023)
- knitr (Xie, 2023)
- RefManageR (McLean, 2017)
- rmarkdown (Allaire, Xie, Dervieux et al., 2023)
- sessioninfo (Wickham, Chang, Flight et al., 2021)
- testthat (Wickham, 2011)

This package was developed using biocthis.

Code for creating the vignette

    ## Create the vignette
    library("rmarkdown")
    system.time(render("iSEEde.Rmd", "BiocStyle::html_document"))

    ## Extract the R code
    library("knitr")
    knit("iSEEde.Rmd", tangle = TRUE)

Date the vignette was generated.

    #> [1] "2023-10-24 17:41:02 EDT"

Wallclock time spent generating the vignette.

    #> Time difference of 19.404 secs

session information.

    #> ─ Session info ───────────────────────────────────────────────────────────────────────────────────────────────────────
    #>  setting  value
    #>  version  R version 4.3.1 (2023-06-16)
    #>  os       Ubuntu 22.04.3 LTS
    #>  system   x86_64, linux-gnu
    #>  ui       X11
    #>  language (EN)
    #>  collate  C
    #>  ctype    en_US.UTF-8
    #>  tz       America/New_York
    #>  date     2023-10-24
    #>  pandoc   2.7.3 @ /usr/bin/ (via rmarkdown)
    #>
    #> ─ Packages ───────────────────────────────────────────────────────────────────────────────────────────────────────────
    #>  package              * version   date (UTC) lib source
    #>  abind                  1.4-5     2016-07-21 [2] CRAN (R 4.3.1)
    #>  airway               * 1.21.0    2023-10-24 [2] Bioconductor
    #>  AnnotationDbi        * 1.64.0    2023-10-24 [2] Bioconductor
    #>  backports              1.4.1     2021-12-13 [2] CRAN (R 4.3.1)
    #>  beachmat               2.18.0    2023-10-24 [2] Bioconductor
    #>  bibtex                 0.5.1     2023-01-26 [2] CRAN (R 4.3.1)
    #>  Biobase              * 2.62.0    2023-10-24 [2] Bioconductor
    #>  BiocGenerics         * 0.48.0    2023-10-24 [2] Bioconductor
    #>  BiocManager            1.30.22   2023-08-08 [2] CRAN (R 4.3.1)
    #>  BiocParallel           1.36.0    2023-10-24 [2] Bioconductor
    #>  BiocStyle            * 2.30.0    2023-10-24 [2] Bioconductor
    #>  Biostrings             2.70.0    2023-10-24 [2] Bioconductor
    #>  bit                    4.0.5     2022-11-15 [2] CRAN (R 4.3.1)
    #>  bit64                  4.0.5     2020-08-30 [2] CRAN (R 4.3.1)
    #>  bitops                 1.0-7     2021-04-24 [2] CRAN (R 4.3.1)
    #>  blob                   1.2.4     2023-03-17 [2] CRAN (R 4.3.1)
    #>  bookdown               0.36      2023-10-16 [2] CRAN (R 4.3.1)
    #>  bslib                  0.5.1     2023-08-11 [2] CRAN (R 4.3.1)
    #>  cachem                 1.0.8     2023-05-01 [2] CRAN (R 4.3.1)
    #>  circlize               0.4.15    2022-05-10 [2] CRAN (R 4.3.1)
    #>  cli                    3.6.1     2023-03-23 [2] CRAN (R 4.3.1)
    #>  clue                   0.3-65    2023-09-23 [2] CRAN (R 4.3.1)
    #>  cluster                2.1.4     2022-08-22 [3] CRAN (R 4.3.1)
    #>  codetools              0.2-19    2023-02-01 [3] CRAN (R 4.3.1)
    #>  colorspace             2.1-0     2023-01-23 [2] CRAN (R 4.3.1)
    #>  colourpicker           1.3.0     2023-08-21 [2] CRAN (R 4.3.1)
    #>  ComplexHeatmap         2.18.0    2023-10-24 [2] Bioconductor
    #>  crayon                 1.5.2     2022-09-29 [2] CRAN (R 4.3.1)
    #>  DBI                    1.1.3     2022-06-18 [2] CRAN (R 4.3.1)
    #>  DelayedArray           0.28.0    2023-10-24 [2] Bioconductor
    #>  DelayedMatrixStats     1.24.0    2023-10-24 [2] Bioconductor
    #>  DESeq2               * 1.42.0    2023-10-24 [2] Bioconductor
    #>  digest                 0.6.33    2023-07-07 [2] CRAN (R 4.3.1)
    #>  doParallel             1.0.17    2022-02-07 [2] CRAN (R 4.3.1)
    #>  dplyr                  1.1.3     2023-09-03 [2] CRAN (R 4.3.1)
    #>  DT                     0.30      2023-10-05 [2] CRAN (R 4.3.1)
    #>  edgeR                * 4.0.0     2023-10-24 [2] Bioconductor
    #>  ellipsis               0.3.2     2021-04-29 [2] CRAN (R 4.3.1)
    #>  evaluate               0.22      2023-09-29 [2] CRAN (R 4.3.1)
    #>  fansi                  1.0.5     2023-10-08 [2] CRAN (R 4.3.1)
    #>  fastmap                1.1.1     2023-02-24 [2] CRAN (R 4.3.1)
    #>  fontawesome            0.5.2     2023-08-19 [2] CRAN (R 4.3.1)
    #>  foreach                1.5.2     2022-02-02 [2] CRAN (R 4.3.1)
    #>  generics               0.1.3     2022-07-05 [2] CRAN (R 4.3.1)
    #>  GenomeInfoDb         * 1.38.0    2023-10-24 [2] Bioconductor
    #>  GenomeInfoDbData       1.2.11    2023-10-17 [2] Bioconductor
    #>  GenomicRanges        * 1.54.0    2023-10-24 [2] Bioconductor
    #>  GetoptLong             1.0.5     2020-12-15 [2] CRAN (R 4.3.1)
    #>  ggplot2                3.4.4     2023-10-12 [2] CRAN (R 4.3.1)
    #>  ggrepel                0.9.4     2023-10-13 [2] CRAN (R 4.3.1)
    #>  GlobalOptions          0.1.2     2020-06-10 [2] CRAN (R 4.3.1)
    #>  glue                   1.6.2     2022-02-24 [2] CRAN (R 4.3.1)
    #>  gtable                 0.3.4     2023-08-21 [2] CRAN (R 4.3.1)
    #>  htmltools              0.5.6.1   2023-10-06 [2] CRAN (R 4.3.1)
    #>  htmlwidgets            1.6.2     2023-03-17 [2] CRAN (R 4.3.1)
    #>  httpuv                 1.6.12    2023-10-23 [2] CRAN (R 4.3.1)
    #>  httr                   1.4.7     2023-08-15 [2] CRAN (R 4.3.1)
    #>  igraph                 1.5.1     2023-08-10 [2] CRAN (R 4.3.1)
    #>  IRanges              * 2.36.0    2023-10-24 [2] Bioconductor
    #>  iSEE                 * 2.14.0    2023-10-24 [2] Bioconductor
    #>  iSEEde               * 1.0.0     2023-10-24 [1] Bioconductor
    #>  iterators              1.0.14    2022-02-05 [2] CRAN (R 4.3.1)
    #>  jquerylib              0.1.4     2021-04-26 [2] CRAN (R 4.3.1)
    #>  jsonlite               1.8.7     2023-06-29 [2] CRAN (R 4.3.1)
    #>  KEGGREST               1.42.0    2023-10-24 [2] Bioconductor
    #>  knitr                  1.44      2023-09-11 [2] CRAN (R 4.3.1)
    #>  later                  1.3.1     2023-05-02 [2] CRAN (R 4.3.1)
    #>  lattice                0.22-5    2023-10-24 [3] CRAN (R 4.3.1)
    #>  lifecycle              1.0.3     2022-10-07 [2] CRAN (R 4.3.1)
    #>  limma                * 3.58.0    2023-10-24 [2] Bioconductor
    #>  locfit                 1.5-9.8   2023-06-11 [2] CRAN (R 4.3.1)
    #>  lubridate              1.9.3     2023-09-27 [2] CRAN (R 4.3.1)
    #>  magrittr               2.0.3     2022-03-30 [2] CRAN (R 4.3.1)
    #>  Matrix                 1.6-1.1   2023-09-18 [3] CRAN (R 4.3.1)
    #>  MatrixGenerics       * 1.14.0    2023-10-24 [2] Bioconductor
    #>  matrixStats          * 1.0.0     2023-06-02 [2] CRAN (R 4.3.1)
    #>  memoise                2.0.1     2021-11-26 [2] CRAN (R 4.3.1)
    #>  mgcv                   1.9-0     2023-07-11 [3] CRAN (R 4.3.1)
    #>  mime                   0.12      2021-09-28 [2] CRAN (R 4.3.1)
    #>  miniUI                 0.1.1.1   2018-05-18 [2] CRAN (R 4.3.1)
    #>  munsell                0.5.0     2018-06-12 [2] CRAN (R 4.3.1)
    #>  nlme                   3.1-163   2023-08-09 [3] CRAN (R 4.3.1)
    #>  org.Hs.eg.db         * 3.18.0    2023-09-21 [2] Bioconductor
    #>  pillar                 1.9.0     2023-03-22 [2] CRAN (R 4.3.1)
    #>  pkgconfig              2.0.3     2019-09-22 [2] CRAN (R 4.3.1)
    #>  plyr                   1.8.9     2023-10-02 [2] CRAN (R 4.3.1)
    #>  png                    0.1-8     2022-11-29 [2] CRAN (R 4.3.1)
    #>  promises               1.2.1     2023-08-10 [2] CRAN (R 4.3.1)
    #>  R6                     2.5.1     2021-08-19 [2] CRAN (R 4.3.1)
    #>  RColorBrewer           1.1-3     2022-04-03 [2] CRAN (R 4.3.1)
    #>  Rcpp                   1.0.11    2023-07-06 [2] CRAN (R 4.3.1)
    #>  RCurl                  1.98-1.12 2023-03-27 [2] CRAN (R 4.3.1)
    #>  RefManageR           * 1.4.0     2022-09-30 [2] CRAN (R 4.3.1)
    #>  rintrojs               0.3.2     2022-08-09 [2] CRAN (R 4.3.1)
    #>  rjson                  0.2.21    2022-01-09 [2] CRAN (R 4.3.1)
    #>  rlang                  1.1.1     2023-04-28 [2] CRAN (R 4.3.1)
    #>  rmarkdown              2.25      2023-09-18 [2] CRAN (R 4.3.1)
    #>  RSQLite                2.3.1     2023-04-03 [2] CRAN (R 4.3.1)
    #>  S4Arrays               1.2.0     2023-10-24 [2] Bioconductor
    #>  S4Vectors            * 0.40.0    2023-10-24 [2] Bioconductor
    #>  sass                   0.4.7     2023-07-15 [2] CRAN (R 4.3.1)
    #>  scales                 1.2.1     2022-08-20 [2] CRAN (R 4.3.1)
    #>  scuttle              * 1.12.0    2023-10-24 [2] Bioconductor
    #>  sessioninfo          * 1.2.2     2021-12-06 [2] CRAN (R 4.3.1)
    #>  shape                  1.4.6     2021-05-19 [2] CRAN (R 4.3.1)
    #>  shiny                  1.7.5.1   2023-10-14 [2] CRAN (R 4.3.1)
    #>  shinyAce               0.4.2     2022-05-06 [2] CRAN (R 4.3.1)
    #>  shinydashboard         0.7.2     2021-09-30 [2] CRAN (R 4.3.1)
    #>  shinyjs                2.1.0     2021-12-23 [2] CRAN (R 4.3.1)
    #>  shinyWidgets           0.8.0     2023-08-30 [2] CRAN (R 4.3.1)
    #>  SingleCellExperiment * 1.24.0    2023-10-24 [2] Bioconductor
    #>  SparseArray            1.2.0     2023-10-24 [2] Bioconductor
    #>  sparseMatrixStats      1.14.0    2023-10-24 [2] Bioconductor
    #>  statmod                1.5.0     2023-01-06 [2] CRAN (R 4.3.1)
    #>  stringi                1.7.12    2023-01-11 [2] CRAN (R 4.3.1)
    #>  stringr                1.5.0     2022-12-02 [2] CRAN (R 4.3.1)
    #>  SummarizedExperiment * 1.32.0    2023-10-24 [2] Bioconductor
    #>  tibble                 3.2.1     2023-03-20 [2] CRAN (R 4.3.1)
    #>  tidyselect             1.2.0     2022-10-10 [2] CRAN (R 4.3.1)
    #>  timechange             0.2.0     2023-01-11 [2] CRAN (R 4.3.1)
    #>  utf8                   1.2.4     2023-10-22 [2] CRAN (R 4.3.1)
    #>  vctrs                  0.6.4     2023-10-12 [2] CRAN (R 4.3.1)
    #>  vipor                  0.4.5     2017-03-22 [2] CRAN (R 4.3.1)
    #>  viridisLite            0.4.2     2023-05-02 [2] CRAN (R 4.3.1)
    #>  xfun                   0.40      2023-08-09 [2] CRAN (R 4.3.1)
    #>  xml2                   1.3.5     2023-07-06 [2] CRAN (R 4.3.1)
    #>  xtable                 1.8-4     2019-04-21 [2] CRAN (R 4.3.1)
    #>  XVector                0.42.0    2023-10-24 [2] Bioconductor
    #>  yaml                   2.3.7     2023-01-23 [2] CRAN (R 4.3.1)
    #>  zlibbioc               1.48.0    2023-10-24 [2] Bioconductor
    #>
    #>  [1] /tmp/RtmppS7J1M/Rinst2b06eb1d35a3e
    #>  [2] /home/biocbuild/bbs-3.18-bioc/R/site-library
    #>  [3] /home/biocbuild/bbs-3.18-bioc/R/library
    #>
    #> ──────────────────────────────────────────────────────────────────────────────────────────────────────────────────────

# 4 Bibliography

This vignette was generated using BiocStyle (Oleś, 2023) with knitr (Xie, 2023) and rmarkdown (Allaire, Xie, Dervieux et al., 2023) running behind the scenes.

Citations made with RefManageR (McLean, 2017).

[1] J. Allaire, Y. Xie, C. Dervieux, et al. rmarkdown: Dynamic Documents for R. R package version 2.25. 2023. URL: https://github.com/rstudio/rmarkdown.

[2] M. W. McLean. “RefManageR: Import and Manage BibTeX and BibLaTeX References in R”. In: The Journal of Open Source Software (2017). DOI: 10.21105/joss.00338.

[3] A. Oleś. BiocStyle: Standard styles for vignettes and other Bioconductor documents. R package version 2.30.0. 2023. DOI: 10.18129/B9.bioc.BiocStyle. URL: https://bioconductor.org/packages/BiocStyle.

[4] R Core Team. R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing. Vienna, Austria, 2023. URL: https://www.R-project.org/.

[5] K. Rue-Albrecht. iSEEde: iSEE extension for panels related to differential expression analysis. R package version 1.0.0. 2023. DOI: 10.18129/B9.bioc.iSEEde. URL: https://bioconductor.org/packages/iSEEde.

[6] H. Wickham. “testthat: Get Started with Testing”. In: The R Journal 3 (2011), pp. 5–10. URL: https://journal.r-project.org/archive/2011-1/RJournal_2011-1_Wickham.pdf.

[7] H. Wickham, W. Chang, R. Flight, et al. sessioninfo: R Session Information. R package version 1.2.2. 2021. URL: https://CRAN.R-project.org/package=sessioninfo.

[8] Y. Xie. knitr: A General-Purpose Package for Dynamic Report Generation in R. R package version 1.44. 2023. URL: https://yihui.org/knitr/.