biocore / mmvec

Neural networks for microbe-metabolite interaction analysis
BSD 3-Clause "New" or "Revised" License
118 stars 51 forks source link

ranks with full taxonomic hierarqui rather than featureid #158

Closed iratxezarra closed 2 years ago

iratxezarra commented 3 years ago

Hi

Is there a way to include the full hierarchical taxonomic names in the rank.qza file? ( rather than just the feature_id itself?) if i convert this rank file into qzv, then I visualize it in qiime view, and save them as .tsv I do not have a way to know which taxa they belong to...

thanks in advance

mortonjt commented 3 years ago

They are both metadata files, so it should be doable. You can definitely combine the taxonomy file with the metadata, but I don’t know of a qiime2 command that do this.

It’ll likely require something as follows in python

import qiime2 import pandas as pd taxa = qiime2.Artifact.load(“taxonomy.qza”).view(pd.DataFrame) ranks = qiime2.Artifact.load(“ranks.qza”).view(pd.DataFrame) combined = pd.concat((ranks.T, taxa), axis=1) combined.to_csv(“combined.txt”, sep=‘\t’)

In this way, you can view in excel, but it’ll be hard to sort the molecules ...

On Thu, Apr 22, 2021 at 7:33 AM iratxezarra @.***> wrote:

Hi

Is there a way to include the full hierarchical taxonomic names in the rank.qza file? ( rather than just the feature_id itself?) if i convert this rank file into qzv, then I visualize it in qiime view, and save them as .tsv I do not have a way to know which taxa they belong to...

thanks in advance

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub https://github.com/biocore/mmvec/issues/158, or unsubscribe https://github.com/notifications/unsubscribe-auth/AA75VXPXGG2KDZEANXDZQTTTKAXS7ANCNFSM43MTAIZQ .

iratxezarra commented 3 years ago

Great, Thanks so much. I will give it a try! :-) Iratxe Zarraonaindia Martinez Maria Goyri, Lab. 1.08. (Applied genomics and Bioinformatics) University of the Basque Country (UPV/EHU) Bº Sarriena s/n 48940, Leioa; Bilbao; Spain Tfno: +34-946015510 E-mail: @., @. http://www.genomic-resources.eu/

El jue, 22 abr 2021 a las 16:44, Jamie Morton @.***>) escribió:

They are both metadata files, so it should be doable. You can definitely combine the taxonomy file with the metadata, but I don’t know of a qiime2 command that do this.

It’ll likely require something as follows in python

import qiime2 import pandas as pd taxa = qiime2.Artifact.load(“taxonomy.qza”).view(pd.DataFrame) ranks = qiime2.Artifact.load(“ranks.qza”).view(pd.DataFrame) combined = pd.concat((ranks.T, taxa), axis=1) combined.to_csv(“combined.txt”, sep=‘\t’)

In this way, you can view in excel, but it’ll be hard to sort the molecules ...

On Thu, Apr 22, 2021 at 7:33 AM iratxezarra @.***> wrote:

Hi

Is there a way to include the full hierarchical taxonomic names in the rank.qza file? ( rather than just the feature_id itself?) if i convert this rank file into qzv, then I visualize it in qiime view, and save them as .tsv I do not have a way to know which taxa they belong to...

thanks in advance

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub https://github.com/biocore/mmvec/issues/158, or unsubscribe < https://github.com/notifications/unsubscribe-auth/AA75VXPXGG2KDZEANXDZQTTTKAXS7ANCNFSM43MTAIZQ

.

— You are receiving this because you authored the thread. Reply to this email directly, view it on GitHub https://github.com/biocore/mmvec/issues/158#issuecomment-824905314, or unsubscribe https://github.com/notifications/unsubscribe-auth/AKUMSD6WFEHBEDAF5VXFRSTTKAY6NANCNFSM43MTAIZQ .