Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

dada2 ouput and importing sample data #1275

Closed
insectnate opened this issue Dec 2, 2019 · 4 comments
Closed

dada2 ouput and importing sample data #1275

insectnate opened this issue Dec 2, 2019 · 4 comments

Comments

@insectnate
Copy link

insectnate commented Dec 2, 2019

I am attempting to switch from using qiime2 and do things directly in R however following the dada2 pipeline I have a phyloseq obj that works but it lacks sample data. I have tried a number of solutions suggested on this and the dada2 forum to add the sample data but I get either an error that
Problem with OTU/taxa indices among those you provided.
or
Component sample names do not match.
I am using a metadata format (csv file) that has worked in phyloseq before but here there is some kind of problem with the merging of the sample data.
This is very confusing as the sample.id names match from the otu_table import and the metadata file. The phyloseq obj made from just the otu table and tax table works fine but without sample data.....

@spholmes
Copy link
Contributor

spholmes commented Dec 3, 2019 via email

@insectnate
Copy link
Author

insectnate commented Dec 3, 2019

Hi Susan,
I have looked for errant #, using unix vs mac line break invisibles, but found nothing that I tried had any effect so far on the behavior of the phyloseq merge where I try to add the sample data to the existing phyloseq object produced from the dada2 output (seqtab.nochim matrix and taxa matrix).

The weird part is I have used this sample data file successfully in phyloseq following the qiime2 to R methodology detailed by @ViridianaAvila previously.
#821

One big difference I notice is the use of the actual DNA sequence as the OTU ID.

@insectnate
Copy link
Author

Ahhh!!! I had to transpose the otu table from dada2 and then switch taxa_are_rows to TRUE. Then the phyloseq(otu_table(nochimTable, taxa_are_rows = TRUE), tax_table(taxa), sample_data(metadata)) seems to work to give a functional phyloseq object from the dada2 pipeline output.

I used import_qiime_sample_data function to import the metadata.
Cheers,
Nathan

@spholmes
Copy link
Contributor

spholmes commented Dec 3, 2019 via email

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants