Construct a phyloseq object directly from the dada2 outputs
Entering edit mode
Manijeh • 0
Last seen 13 months ago

I have a group of samples, 8 of them control, and the rest have a disease. we want to analyze their microbiome in these two groups. I can do the analysis in DADA2 and I got response from ASSIGN TAXONOMY. But for the phyloseq section, I cant define "ps" or "phyloseq object" for further analysis. I don't know what should I replace instead of these commands in the tutorial.

``samples.out <- rownames(seqtab.nochim) subject <- sapply(strsplit(samples.out, "D"),[, 1) gender <- substr(subject,1,1) subject <- substr(subject,2,999) day <- as.integer(sapply(strsplit(samples.out, "D"),[`, 2)) samdf <- data.frame(Subject=subject, Gender=gender, Day=day) samdf$When <- "Early" samdf$When[samdf$Day>100] <- "Late" rownames(samdf) <- samples.out


phyloseq dada2 • 302 views

Login before adding your answer.

Traffic: 262 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6