Skip to content

Commit

Permalink
Update iSEEtree part (#2)
Browse files Browse the repository at this point in the history
  • Loading branch information
RiboRings authored Sep 16, 2024
1 parent f9be1ec commit 1185b59
Showing 1 changed file with 20 additions and 9 deletions.
29 changes: 20 additions & 9 deletions vignettes/eccb_workflow.Rmd
Original file line number Diff line number Diff line change
Expand Up @@ -826,25 +826,36 @@ This task could be performed with the NetCoMi function `netCompare`, as describe

## Interactive microbiome data exploration with iSEEtree

```{r agg_iseetree}
tse_phylum <- agglomerateByRank(tse, rank = "phylum")
tse_phylum <- transformAssay(tse_phylum, method = "relabundance")
tse_phylum <- transformAssay(tse_phylum, method = "clr", pseudocount = 1)
```

[iSEEtree](https://microbiome.github.io/iSEEtree/) is a shiny app that provides
a Graphical User Interface (GUI) to visualize and interactively explore `TreeSE`
objects. We will demonstrate its capability by reproducing some of the plots
that were previously created in this workshop, this time without the need of
coding.
coding (small lie).

As we did before for `tse`, we add alpha and beta diversity metrics also to the agglomerated experiment `tse_family`, so we will be able to visualize these
metrics in the app.

```{r prep_iSEEtree}
# Add information on alpha diversity
tse_family <- addAlpha(tse_family,
index = c("coverage_diversity", "shannon_diversity"))
# Add information on beta diversity
tse_family <- runMDS(tse_family,
FUN = getDissimilarity,
method = "bray",
assay.type = "relabundance",
ncomp = 10,
name = "MDS_bray")
```

Launching the app is easy. All you need to do is run the command `iSEE` on the
`TreeSE` object that you want to explore.

```{r iseetree1}
# Launch app
if (interactive()) {
iSEE(tse_phylum)
iSEE(tse_family)
}
```

Expand All @@ -855,7 +866,7 @@ panel layout with the `initial` argument to generate the alpha diversity plot.
```{r iseetree2}
# Launch app with custom layout
if (interactive()) {
iSEE(tse_phylum, initial = c(ColumnDataTable(), ColumnDataPlot()))
iSEE(tse_family, initial = c(ColumnDataTable(), ColumnDataPlot()))
}
```

Expand Down

0 comments on commit 1185b59

Please sign in to comment.