In adrientaudiere/MiscMetabar: Miscellaneous Functions for Metabarcoding Analysis

knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>",
  message = FALSE
)

library(MiscMetabar)
data(data_fungi)
data(data_fungi_sp_known)

Export phyloseq object

You can export a phyloseq object to csv (and txt for phylogenetic tree) files in a folder. It is possible to export each table into one file or to merge all slot (except phytree) in one file (args one_file = TRUE). Finally, if rdata is set to TRUE, a rdata file containing the phyloseq object is also writed.

write_pq(data_fungi, path = "fungi_phyloseq")
write_pq(data_fungi, path = "fungi_phyloseq", one_file = TRUE)
write_pq(data_fungi, path = "fungi_phyloseq", rdata = TRUE)

Finally, you can use the function save_pq() to write the phyloseq object in all 3 versions (one table for each slot, a file merging each slot and an Rdata file).

save_pq(data_fungi)

Import

To import a Rdata file, just use load() base function. In order to import phyloseq object from a folder create using write_pq() or save_pq(), please use the function read_pq().

d <- read_pq(path = "fungi_phyloseq")

Tracking sequences, clusters and samples

In bioinformatic pipeline, we often need to track the number of samples, sequences and clusters across step in the pipeline. MiscMetabar propose two utilities to achieve this goal : track_wkflow() and a version to compute value per samples : track_wkflow_samples(). The function track_wkflow() can deal with (i) fastq and fastg.gz files, dada-class object, derep-class object, matrix of samples x clusters (e.g. otu_table) and phyloseq-class object.

track_wkflow(list(data_fungi, data_fungi_sp_known))

Session information

sessionInfo()

adrientaudiere/MiscMetabar documentation built on Feb. 7, 2025, 5:47 p.m.