Skip to content

Commit

Permalink
Merge branch 'dev' of github.com:adrientaudiere/MiscMetabar into dev
Browse files Browse the repository at this point in the history
  • Loading branch information
adrientaudiere committed Dec 12, 2023
2 parents 254965f + 33249e3 commit b17e816
Show file tree
Hide file tree
Showing 5 changed files with 13,962 additions and 6 deletions.
4 changes: 2 additions & 2 deletions DESCRIPTION
Original file line number Diff line number Diff line change
Expand Up @@ -95,8 +95,8 @@ Remotes:
github::adrientaudiere/lulu,
github::david-barnett/microViz,
github::erocoar/gghalves,
bioc::release/phyloseq,
bioc::release/dada2,
bioc::phyloseq,
bioc::dada2,
github::brendanf/FUNGuildR
Config/testthat/edition: 3
Config/testthat/parallel: true
Expand Down
4 changes: 2 additions & 2 deletions paper/paper.bib
Original file line number Diff line number Diff line change
@@ -1,4 +1,4 @@
@article{altschul1990,
@article{altschul1990,
doi = {10.32388/rhq6vj},
title = {Basic local alignment search tool},
author = {Altschul, Stephen F and Gish, Warren and Miller, Webb and Myers, Eugene W and Lipman, David J},
Expand Down Expand Up @@ -300,4 +300,4 @@ @article{xu2023
number = {2},
year = {2023},
publisher = {Elsevier}
}
}
3 changes: 1 addition & 2 deletions paper/paper.md
Original file line number Diff line number Diff line change
Expand Up @@ -62,7 +62,6 @@ The filter functions `subset_samples_pq()` and `subset_taxa_pq()` complement `su

I also implement a function to filter taxa based on their blast to a custom database (`filter_asv_blast()`). This function uses the blastn software [@altschul1990] to compare ASV sequences to a database and filter out species that are below a given threshold of e-value and/or bit-score.


## Exploration

`MiscMetabar` provides a large number of facilities to explore the biological diversity in a phyloseq object. In most functions, a parameter enables the effect of the number of reads (sampling depth) to be controlled by rarefaction or other statistical methods, depending on the function. For example, the alpha diversity analysis (function `hill_pq()`) uses the HSD-Tuckey test on a linear model that includes the square roots of the number of reads as the first explanatory variable.
Expand All @@ -79,4 +78,4 @@ The targets R package [@landau2021] improves the efficiency and reproducibility

I thank Will Landau, Paul McMurdie, and Benjamin Callahan for their excellent R packages on which `MiscMetabar` rests. I also want to acknowledge Franck Richard, Lise Roy, Élisa Taschen and the [Mycea](https://mycea.fr/) team for the discussion and work around metabarcoding.

# References
# References
Loading

0 comments on commit b17e816

Please sign in to comment.