Import_biom with tree problem
0
0
Entering edit mode
@andrebolerbarros-16788
Last seen 5 months ago
Portugal

Hey!

I'm trying to import my mothur data in phyloseq, by first of all import the list, group and tree files:

   phyl<-import_mothur(mothur_list_file = "total.final.pick.opti_mcc.list", mothur_group_file="total.final.pick.groups", mothur_tree_file = "total.final.pick.tre")

However, I get an error message:

`Error in 3:as.integer(rawlines[2]) :

NA/NaN argument In addition:

Warning messages:

1: In asMethod(object) : NAs introduced by coercion

2: In import_mothur_otulist(mothur_list_file, cutoff) : NAs introduced by coercion`

Apparently, there's a problem with the group file, since I can't really import it whatsoever (gives that error message). So, I've tried with import_biom:

   phyl<-import_biom(BIOMfilename = "total.final.pick.opti_mcc.0.03.biom", treefilename = "total.final.pick.tre")

Which, in turn, gives the error output:

`Error in validObject(.Object):

invalid class “phyloseq” object: Component taxa/OTU names do not match.

Taxa indices are critical to analysis. Try taxa_names()`

Is there any suggestion about importing mothur files "directly from the pipeline?

Note: I've seen post suggestion using different tree-producing softwares, but those also require several modifications of the trees. I would like to know if there's a more straightforward solution

phyloseq mothur 16s • 1.8k views
ADD COMMENT
0
Entering edit mode

Solved in phyloseq github; Issue #408:
Removing the label-row from list file solved the problem

ADD REPLY

Login before adding your answer.

Traffic: 700 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6