phyloseq is a set of classes, wrappers, and tools (in R) to make it easier to import, store, and analyze phylogenetic sequencing data; and to reproducibly share that data and analysis with others. See the phyloseq front page:
Thank you for creating phyloseq, your contributions are, dare I say invaluable and have helped me immensely!
I am trying to find bootstrap in my tree but I haven't been successful yet,
your advice is appreciated.
I created the tree based on this paper: https://f1000research.com/articles/5-1492
as following:
seqs <- getSequences(seqtab)
names(seqs) <- seqs # This propagates to the tip labels of the tree
alignment <- AlignSeqs(DNAStringSet(seqs), anchor=NA)
phang.align <- phyDat(as(alignment, "matrix"), type="DNA")
dm <- dist.ml(phang.align)
treeNJ <- NJ(dm) # Note, tip order != sequence order
fit = pml(treeNJ, data=phang.align)
a quick example of the tree I viewed using plot_tree
here the structure of the final tree created before merging it into the phyloseq object
dismiss the $data because it takes a space out of the 24 list shown below
I don't know which one of these is considered the bootstrap, if this is to be figured out, I think bootstrap would be easily integrated into the tree plot
Thank you for creating phyloseq, your contributions are, dare I say invaluable and have helped me immensely!
I am trying to find bootstrap in my tree but I haven't been successful yet, your advice is appreciated. I created the tree based on this paper: https://f1000research.com/articles/5-1492 as following:
a quick example of the tree I viewed using plot_tree
here the structure of the final tree created before merging it into the phyloseq object dismiss the $data because it takes a space out of the 24 list shown below I don't know which one of these is considered the bootstrap, if this is to be figured out, I think bootstrap would be easily integrated into the tree plot
Originally posted by @marwa38 in https://github.com/joey711/phyloseq/issues/1354#issuecomment-969286418