Tag: phyloseq

Cannot install Phyloseq and dada2

Hello all, I have been having issues installing packages that I really need to use. Basically, I cannot download either Phyloseq or dada2 and I believe it’s because I don’t have GenomeInfoDbData. But at the same time, I cannot install GenomeInfoDbData because I can’t seem to update the dependencies (“fansi”…

Continue Reading Cannot install Phyloseq and dada2

combine OTU and tax table and replace actual sequences with OTU ids (Phyloseq/dada2)

This is to the part of the question “replace actual sequences with OTU ids (Phyloseq/dada2)?” I contacted the phyloseq/dada2 developers and based on Susan Holmes’ reply (github.com/joey711/phyloseq/issues/1030) I came up with this piece of code to replace the amplicon sequences with a numbered OTU header. Further discussion can be found…

Continue Reading combine OTU and tax table and replace actual sequences with OTU ids (Phyloseq/dada2)

Change size of label annotations in a ggplot

I am trying to change text label sizes inside my plot (not the axes, rather the label annotations) I am working with a phyloseq object but I don’t think that matters. Here is the code and the output. Any suggestions? plot_ordination(prokaryote_ra, ordBC, color = “Stage”, label=”SampleID”) + ggtitle(“PCoA: Bray-Curtis”) graph…

Continue Reading Change size of label annotations in a ggplot

Zooplankton diversity monitoring strategy for the urban coastal region using metabarcoding analysis

1. Eyun, S. Phylogenomic analysis of Copepoda (Arthropoda, Crustacea) reveals unexpected similarities with earlier proposed morphological phylogenies. BMC Evol. Biol. 17, 23 (2017). PubMed  PubMed Central  Google Scholar  2. Eyun, S. et al. Evolutionary history of chemosensory-related gene families across the Arthropoda. Mol. Biol. Evol. 34, 1838–1862 (2017). CAS  PubMed …

Continue Reading Zooplankton diversity monitoring strategy for the urban coastal region using metabarcoding analysis

t-SNE for microbiome analysis

t-SNE for microbiome analysis 1 Hi I wondered whether there is a good tutorial on how to make t-SNE plot on microbiome samples. I have a phyloseq object for which I have metadata age and gender. Ideally I want something like this: So different shapes for gender, and in my…

Continue Reading t-SNE for microbiome analysis

Phyloseq Objects for different time points

Phyloseq Objects for different time points 0 Hi I am very new to this microbiome analysis so this might be a very simple question… I have a question related to the microbiome analysis I’m doing. I have microbiome data from 3 different timepoints (T1, T7, and T13), for I made…

Continue Reading Phyloseq Objects for different time points

How to get normalized count table from DESeq?

How to get normalized count table from DESeq? 1 Hi, I’m using Deseq compare differential abundance. Here is my code: ds.all <- phyloseq_to_deseq2(ps0.infant.pbs, ~ sample_type) geoMeans <- apply(counts(ds.all),1,gm_mean) ds.all <- estimateSizeFactors(ds.all,geoMeans = geoMeans) dds.all <- DESeq(ds.all,fitType = “local”) Then as the results I got 8 ASVs that showed significantly different….

Continue Reading How to get normalized count table from DESeq?

Network plot using NetCoMi and igraph

Hi, I am trying to plot the Network plot as suggested here (github.com/stefpeschel/NetCoMi#single-association-network-on-genus-level) by using igraph and NetCoMi. But I am not getting the network plot as expected- I just want to label the hub genera and phyla. and want the network plot in spherical layout. # Agglomerate to genus…

Continue Reading Network plot using NetCoMi and igraph

Bioconductor – DESeq2

DOI: 10.18129/B9.bioc.DESeq2     This package is for version 3.10 of Bioconductor; for the stable, up-to-date release version, see DESeq2. Differential gene expression analysis based on the negative binomial distribution Bioconductor version: 3.10 Estimate variance-mean dependence in count data from high-throughput sequencing assays and test for differential expression based on…

Continue Reading Bioconductor – DESeq2

export ANCOM-BC normalisation matrix

Hello! I have similar doubt. Did you solve it? I run ancombc with my Phyloseq opbject and extract the data, but we are not certain about the normalization: out = ancombc(phyloseq = phyloseq_obj, formula = “condition_1”, p_adj_method = “BH”, zero_cut = 0.90, lib_cut = 1000, group = NULL, struc_zero =…

Continue Reading export ANCOM-BC normalisation matrix