biplot_pq
,
circle_pq
, upset_pq
,
ggvenn_pq
)For an introduction to metabarcoding in R, Please visite the state
of the field vignettes. The import,
export and track vignette explains how import and export
phyloseq
object. Its also show how to summarize useful
information (number of sequences, samples and clusters) accross
bioinformatic pipelines.
If you are interested in ecological metrics, see the vignettes describing alpha-diversity and beta-diversity analysis. The vignette filter taxa and samples describes some data-filtering processes using MiscMetabar and the reclustering tutorial introduces the different way of clustering already-clustered OTU/ASV. The vignette tengeler explore the dataset from Tengeler et al. (2020) using some MiscMetabar functions.
For developers, I also wrote a vignette describing som rules of codes.
library("MiscMetabar")
library("phyloseq")
library("magrittr")
data("data_fungi")
summary_plot_pq(data_fungi)
Tengeler, A.C., Dam, S.A., Wiesmann, M. et al. Gut microbiota from persons with attention-deficit/hyperactivity disorder affects the brain in mice. Microbiome 8, 44 (2020). https://doi.org/10.1186/s40168-020-00816-x
sessionInfo()
#> R version 4.4.1 (2024-06-14)
#> Platform: x86_64-pc-linux-gnu
#> Running under: Ubuntu 24.04.1 LTS
#>
#> Matrix products: default
#> BLAS: /usr/lib/x86_64-linux-gnu/openblas-pthread/libblas.so.3
#> LAPACK: /usr/lib/x86_64-linux-gnu/openblas-pthread/libopenblasp-r0.3.26.so; LAPACK version 3.12.0
#>
#> locale:
#> [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C
#> [3] LC_TIME=en_US.UTF-8 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
#>
#> time zone: Etc/UTC
#> tzcode source: system (glibc)
#>
#> attached base packages:
#> [1] stats graphics grDevices utils datasets methods base
#>
#> other attached packages:
#> [1] magrittr_2.0.3 MiscMetabar_0.9.3 purrr_1.0.2 dplyr_1.1.4
#> [5] dada2_1.33.0 Rcpp_1.0.13 ggplot2_3.5.1 phyloseq_1.49.0
#> [9] rmarkdown_2.28
#>
#> loaded via a namespace (and not attached):
#> [1] bitops_1.0-8 deldir_2.0-4
#> [3] permute_0.9-7 rlang_1.1.4
#> [5] ade4_1.7-22 matrixStats_1.4.0
#> [7] compiler_4.4.1 mgcv_1.9-1
#> [9] png_0.1-8 vctrs_0.6.5
#> [11] reshape2_1.4.4 stringr_1.5.1
#> [13] pwalign_1.1.0 pkgconfig_2.0.3
#> [15] crayon_1.5.3 fastmap_1.2.0
#> [17] XVector_0.45.0 labeling_0.4.3
#> [19] utf8_1.2.4 Rsamtools_2.21.1
#> [21] UCSC.utils_1.1.0 xfun_0.47
#> [23] zlibbioc_1.51.1 cachem_1.1.0
#> [25] GenomeInfoDb_1.41.1 jsonlite_1.8.8
#> [27] biomformat_1.33.0 highr_0.11
#> [29] rhdf5filters_1.17.0 DelayedArray_0.31.11
#> [31] Rhdf5lib_1.27.0 BiocParallel_1.39.0
#> [33] jpeg_0.1-10 parallel_4.4.1
#> [35] cluster_2.1.6 R6_2.5.1
#> [37] bslib_0.8.0 stringi_1.8.4
#> [39] RColorBrewer_1.1-3 ComplexUpset_1.3.3
#> [41] GenomicRanges_1.57.1 jquerylib_0.1.4
#> [43] SummarizedExperiment_1.35.1 iterators_1.0.14
#> [45] knitr_1.48 IRanges_2.39.2
#> [47] Matrix_1.7-0 splines_4.4.1
#> [49] igraph_2.0.3 tidyselect_1.2.1
#> [51] abind_1.4-5 yaml_2.3.10
#> [53] vegan_2.6-8 codetools_0.2-20
#> [55] hwriter_1.3.2.1 lattice_0.22-6
#> [57] tibble_3.2.1 plyr_1.8.9
#> [59] Biobase_2.65.1 withr_3.0.1
#> [61] ShortRead_1.63.1 evaluate_0.24.0
#> [63] survival_3.7-0 RcppParallel_5.1.9
#> [65] Biostrings_2.73.1 pillar_1.9.0
#> [67] MatrixGenerics_1.17.0 DT_0.33
#> [69] foreach_1.5.2 stats4_4.4.1
#> [71] generics_0.1.3 S4Vectors_0.43.2
#> [73] munsell_0.5.1 scales_1.3.0
#> [75] glue_1.7.0 maketools_1.3.0
#> [77] tools_4.4.1 interp_1.1-6
#> [79] sys_3.4.2 data.table_1.16.0
#> [81] GenomicAlignments_1.41.0 buildtools_1.0.0
#> [83] rhdf5_2.49.0 grid_4.4.1
#> [85] tidyr_1.3.1 ape_5.8
#> [87] crosstalk_1.2.1 latticeExtra_0.6-30
#> [89] colorspace_2.1-1 patchwork_1.2.0
#> [91] networkD3_0.4 nlme_3.1-166
#> [93] GenomeInfoDbData_1.2.12 cli_3.6.3
#> [95] fansi_1.0.6 S4Arrays_1.5.7
#> [97] gtable_0.3.5 sass_0.4.9
#> [99] digest_0.6.37 BiocGenerics_0.51.1
#> [101] SparseArray_1.5.31 htmlwidgets_1.6.4
#> [103] farver_2.1.2 htmltools_0.5.8.1
#> [105] multtest_2.61.0 lifecycle_1.0.4
#> [107] httr_1.4.7 MASS_7.3-61