The bioBakery help forum

Merged abundance table to biom format

Hi community!!!
I want to import metaphlan output of 60 samples to phyloseq. For that, I want to import the biom format of merged abundance table output from merge_metaphlan_tables.py. Can anyone please tell me how can I convert this merged table into biom file so that I can import it into phyloseq?

Thanks,
DC7

There’s no need to convert to biom, you can load the table with the merged profiles as an otu_table object. There’s a full and exhaustive tutorial on the phyloseq github page

1 Like

Thanks a lot @fbeghini. I am following the tutorial and successfully imported it in phyloseq.
But I am facing one important “ERROR” there which says the abundance proportion is not accepted their. They accepts read count. Have you ever tried phyloseq and encountered such error? Do you have any other way to do the analysis within phyloseq bypassing this problem? Here is the error message:

My code:

alpha_meas = c(“Observed”, “Chao1”, “ACE”, “Shannon”, “Simpson”, “InvSimpson”)
p <- plot_richness(physeq1, “type”, measures=alpha_meas, color=“time”)

ERROR returned:

Error in h(simpleError(msg, call)) : 
  error in evaluating the argument 'x' in selecting a method for function 't': function accepts only integers (counts)
In addition: Warning message:
In estimate_richness(physeq, split = TRUE, measures = measures) :
  The data you have provided does not have
any singletons. This is highly suspicious. Results of richness
estimates (for example) are probably unreliable, or wrong, if you have already
trimmed low-abundance taxa from the data.

We recommended that you find the un-trimmed data and retry.

Thanks and regards,
DC7

It may be possible that the column type is not a numeric and a string/factor instead.
You should cast the numeric values in the otu table to numeric with as.numeric.

This is a support page for the BioBakery tools, for issues with other tools like phyloseq, I suggest you to post in their specific support page

1 Like