Dear all,
when i was working on my soil 16 amplicon sequencing data (12 samples taken for 2 treatements), by looking at the Pcoa plot of weighted UniFrac distances, i discovered one sample isolated from all the others and noticed this sample has very high levels of phosphorus compared to the others.
I then tried DA analysis (ALDEx2 in R studio ggpicrust2) on the q2 picrust2 functional annotated data to confirm my hypothesis but i think i cannot perform it becouse it is just one sample and i think i need at least 3 per group. Can anyone suggest me a way to discover differentally abundant OTUs or functional annotations for 1 sample vs all? or other ways to analyze the same data?
thanks in advance,
Michele
i will report my error here just in case:
computing center with all features
Error in apply(data[, grp2], 1, sd) : dim(X) must have a positive length