Skip to content

Commit

Permalink
up
Browse files Browse the repository at this point in the history
  • Loading branch information
TuomasBorman committed Nov 27, 2024
1 parent d9b92e0 commit f9fe87f
Show file tree
Hide file tree
Showing 15 changed files with 350 additions and 679 deletions.
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
2 changes: 1 addition & 1 deletion docs/index.html
Original file line number Diff line number Diff line change
Expand Up @@ -107,7 +107,7 @@ <h3 class="anchored" data-anchor-id="bioconductor">Bioconductor</h3>
<ul>
<li><a href="https://microbiome.github.io/outreach/intro.html" target="_blank">Intro</a></li>
<li><a href="https://microbiome.github.io/outreach/bioconductor.html" target="_blank">Bioconductor project</a></li>
<li><a href="https://microbiome.github.io/outreach/oma.html" target="_blank">Orchestrating microbiome Analysis with Bioconductor</a></li>
<li><a href="https://microbiome.github.io/outreach/oma.html" target="_blank">Orchestrating Microbiome Analysis with Bioconductor</a></li>
</ul>
</section>
<section id="quarto" class="level3">
Expand Down
191 changes: 22 additions & 169 deletions docs/oma.html

Large diffs are not rendered by default.

Binary file modified docs/oma_files/figure-revealjs/bioc_packages-1.png
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file modified docs/oma_files/figure-revealjs/data_and_data_container-1.png
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file modified docs/oma_files/figure-revealjs/data_container-1.png
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file modified docs/oma_files/figure-revealjs/pages_and_book-1.png
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file modified docs/oma_files/figure-revealjs/se_field-1.png
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file removed docs/oma_files/figure-revealjs/show_alpha-1.png
Binary file not shown.
Binary file removed docs/oma_files/figure-revealjs/show_prevalence-1.png
Binary file not shown.
Binary file removed docs/oma_files/figure-revealjs/sow_pcoa-1.png
Binary file not shown.
Binary file modified docs/pcoa_files/figure-revealjs/pca-load-1.png
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
2 changes: 1 addition & 1 deletion docs/quarto.html
Original file line number Diff line number Diff line change
Expand Up @@ -516,7 +516,7 @@ <h2>Example 2.3: YAML Parameters</h2>
<p>editor: visual</p>
<p>smaller: true</p>
<p>author: Escherichia coli</p>
<p>date: 2024-11-24</p>
<p>date: 2024-11-27</p>
<p><code>---</code></p>
</section>
<section id="exercise-2" class="slide level2">
Expand Down
692 changes: 318 additions & 374 deletions docs/search.json

Large diffs are not rendered by default.

142 changes: 8 additions & 134 deletions quarto/oma.qmd
Original file line number Diff line number Diff line change
Expand Up @@ -30,7 +30,6 @@ format:

- ~2,300 R packages
- Review, testing, documentation
- Genomics, transcriptomics, microbiomics, ...

```{r}
#| label: bioc_packages
Expand Down Expand Up @@ -115,7 +114,7 @@ p1 <- ggplot(pkgs_date, aes(x = Date, y = N, fill = Field)) +
p1
```

## Data containers form the foundation {.smaller}
## Data containers form the core {.smaller}

```{r}
#| label: data_container
Expand All @@ -130,7 +129,7 @@ ellipse_data <- data.frame(
y = c(2, 1, 0), # Centers of ellipses
a = c(4, 3, 2), # Widths of ellipses
b = c(3, 2, 1), # Heights of ellipses
label = c("COMMUNITY", "PACKAGES", "DATA"), # Labels for each ellipse
label = c("COMMUNITY", "METHODS", "DATA CONTAINER"), # Labels for each ellipse
label_y = c(4, 1.75, 0) # Adjusted vertical positions for labels
)
Expand Down Expand Up @@ -306,17 +305,17 @@ _Reduce overlapping efforts, improve interoperability, ensure sustainability._

- Extension to SummarizedExperiment
- Optimal for microbiome data
- Links microbiome field to larger SE family
- Links microbiome field to larger SummarizedExperiment family

##

![](images/SE.png){fig-alt="SummarizedExperiment class" fig-align="center" width=10%}

## {transition="fade" transition-speed="slow"}
##

![](images/paste-14DB8F76.png){fig-alt="TreeSummarizedExperiment class" fig-align="center" width=10%}

## MIcrobiome Analysis (mia) {transition="none"}
## MIcrobiome Analysis (mia)

```{r}
#| label: mia_stats
Expand Down Expand Up @@ -375,133 +374,8 @@ perc <- paste0(round(which(rownames(df) == "mia") / nrow(df), 3)*100, "%")
- Scalable & optimized for large datasets
- Comprehensive documentation

_Allows us to develop modular and efficient workflows_
_Allows us to develop efficient microbiome data science workflows_

## {auto-animate="true"}

```r
# Load package
library(mia)
# Load example dataset
data("peerj13075")
tse <- peerj13075
```

```{r}
#| label: show_treese
# Load package
library(mia)
# Load example dataset
data("peerj13075")
tse <- peerj13075
tse
```


## {auto-animate="true"}

```r
# Agglomerate to genus level
tse <- agglomerateByRank(tse, rank = "genus")
```

## {auto-animate="true"}

```r
# Agglomerate to genus level
tse <- agglomerateByRank(tse, rank = "genus")

# Add relative abundances
tse <- transformAssay(tse, method = "relabundance")
```

```{r}
#| label: show_transform
# Agglomerate to genus level
tse <- agglomerateByRank(tse, rank = "genus")
# Add relative abundances
tse <- transformAssay(tse, method = "relabundance")
```

## {auto-animate="true"}

```r
# Load visualization package
library(miaViz)
# Summarize abundance of top taxa
plotAbundanceDensity(tse, assay.type = "relabundance")
```

```{r}
#| label: show_prevalence
# Load visualization package
library(miaViz)
# Summarize abundance of top taxa
plotAbundanceDensity(tse, assay.type = "relabundance")
```

## {auto-animate="true"}

```r
# Calculate alpha diversity indices
tse <- addAlpha(tse, index = "shannon")
```

## {auto-animate="true"}

```r
# Calculate alpha diversity indices
tse <- addAlpha(tse, index = "shannon")

# Load single-cell analysis package that has useful, complementary tools
library(scater)
# Plot alpha diversity
plotColData(tse, x = "Geographical_location", y = "shannon")
```

```{r}
#| label: show_alpha
#| fig-height: 4
# Calculate alpha diversity indices
tse <- addAlpha(tse, index = "shannon")
# Load single-cell analysis package that has useful, complementary tools
library(scater)
# Plot alpha diversity
plotColData(tse, x = "Geographical_location", y = "shannon")
```

## {auto-animate="true"}

```r
# Perform PCoA
tse <- runMDS(tse, assay.type = "relabundance", FUN = getDissimilarity, method = "bray")
```

## {auto-animate="true"}

```r
# Perform PCoA
tse <- runMDS(tse, assay.type = "relabundance", FUN = getDissimilarity, method = "bray")
# Plot PCoA
plotReducedDim(tse, dimred = "MDS", colour_by = "Geographical_location")
```

```{r}
#| label: show_pcoa
# Perform PCoA
tse <- runMDS(tse, assay.type = "relabundance", FUN = getDissimilarity, method = "bray")
# Plot PCoA
plotReducedDim(tse, dimred = "MDS", colour_by = "Geographical_location")
```

## Orchestrating Microbiome Analysis with Bioconductor

Expand All @@ -515,7 +389,7 @@ plotReducedDim(tse, dimred = "MDS", colour_by = "Geographical_location")
[microbiome.github.io/OMA](https://microbiome.github.io/OMA/docs/devel/){preview-link="true"}
:::

## Summary
## Poem of the day

::: columns
::: {.column width="60%"}
Expand All @@ -525,7 +399,7 @@ _mia framework’s waiting here._


_OMA’s guide will light the way,_
_Helping you every step of the day._
_helping you every step of the day._

:::

Expand Down

0 comments on commit f9fe87f

Please sign in to comment.