mia icon indicating copy to clipboard operation
mia copied to clipboard

readQZA fails

Open antagomir opened this issue 1 year ago • 4 comments

I tried to load example files "distance.qza" and "ordination.qza" from QIIME/DEICODE tutorial.

After downloading them I did in R:

distance.qiime <- mia::readQZA("distance.qza")
ordination.qiime <- mia::readQZA("ordination.qza")

Both throw the error:

Error: Only files in format of 'BIOMV210DirFmt', 'TSVTaxonomyDirectoryFormat', NewickDirectoryFormat' and 'DNASequencesDirectoryFormat' are supported.

-> We might need to add support for more formats? I am not sure how critical feature this is.

antagomir avatar Feb 13 '24 21:02 antagomir

Interestingly, the qiime2R package seems to work. Something to learn from there?

remotes::install_github("jbisanz/qiime2R")
library(qiime2R)
d <- read_qza("distance.qza")
o <- read_qza("ordination.qza")

antagomir avatar Feb 13 '24 22:02 antagomir

We are now adding TSE support for qiime2R; perhaps this solves the issue and we could switch to using qiime2R functions instead of maintaining parallel mia implementation?

This is currently limited by the fact that qiime2R is not available via CRAN or Bioc. There is an open issue on that, however, perhaps the problem will be solved: https://github.com/jbisanz/qiime2R/issues/3

antagomir avatar Apr 21 '24 08:04 antagomir

I'll keep an eye on this issue. Might be good to add an example in the book once it's easier to load in qza files and transform them into a tse object.

Insaynoah avatar Apr 24 '24 19:04 Insaynoah

Support for TreeSE has been added now in qiime2R PR 64.

However that pkg is not yet in Bioc or CRAN

antagomir avatar Jul 24 '24 08:07 antagomir