Hello,
I am trying to follow the example given here https://github.com/joey711/phyloseq/issues/683, using DeSeq2 to examine differential abundance at Genus rather than OTU level, however When I import my data into R studio doing the following I only have one taxonomic rank (Kingdom).
map<-import_qiime_sample_data(file.choose())
tree<-read_tree(file.choose())
biomfile<-import_biom(file.choose(),parseFunction=parse_taxonomy_greengenes)
physeq<-merge_phyloseq(biomfile,tree,map)
physeq
phyloseq-class experiment-level object
otu_table() OTU Table: [ 2206 taxa and 66 samples ]
sample_data() Sample Data: [ 66 samples by 11 sample variables ]
tax_table() Taxonomy Table: [ 2206 taxa by 1 taxonomic ranks ]
phy_tree() Phylogenetic Tree: [ 2206 tips and 2204 internal nodes ]
Previously I had imported my data using the following code and I see 7 taxonomic ranks, however when I try to use this data with the code supplied by joey711 (link above) I get the error below
otufile=("otu_table_m1000_minus_cyano_chloro_json_w_consensuslineage.txt")
mapfile=("Mapping_file_LS.txt")
trefile=("rep_set_tree.tre")
qiimedata = import_qiime(otufile, mapfile, trefile)
qiimedata
phyloseq-class experiment-level object
otu_table() OTU Table: [ 2206 taxa and 66 samples ]
sample_data() Sample Data: [ 66 samples by 11 sample variables ]
tax_table() Taxonomy Table: [ 2206 taxa by 7 taxonomic ranks ]
phy_tree() Phylogenetic Tree: [ 2206 tips and 2204 internal nodes ]
GenusFiltered <- subset_taxa(qiimedata, Genus != "NA" & Species == "NA")
Error in rownames<-(*tmp*, value = c("sp1", "sp0")) :
length of 'dimnames' [1] not equal to array extent
The package version of phyloseq that I am using in R studio is 1.19.1
Any help to fix this would be much appreciated
Secondly does anyone know if there is way to use the DeSeq2 function in Qiime to compare higher taxonomic ranks, eg Genus? I am aware that I would have to take the results with a pinch of salt.
If there isnt a way is there a way in which I can use Qiime to do permutation based tests, other than at OTU level?
Cheers,
Sarah
I am trying to follow the example given here https://github.com/joey711/phyloseq/issues/683, using DeSeq2 to examine differential abundance at Genus rather than OTU level, however When I import my data into R studio doing the following I only have one taxonomic rank (Kingdom).
map<-import_qiime_sample_data(file.choose())
tree<-read_tree(file.choose())
biomfile<-import_biom(file.choose(),parseFunction=parse_taxonomy_greengenes)
physeq<-merge_phyloseq(biomfile,tree,map)
physeq
phyloseq-class experiment-level object
otu_table() OTU Table: [ 2206 taxa and 66 samples ]
sample_data() Sample Data: [ 66 samples by 11 sample variables ]
tax_table() Taxonomy Table: [ 2206 taxa by 1 taxonomic ranks ]
phy_tree() Phylogenetic Tree: [ 2206 tips and 2204 internal nodes ]
Previously I had imported my data using the following code and I see 7 taxonomic ranks, however when I try to use this data with the code supplied by joey711 (link above) I get the error below
otufile=("otu_table_m1000_minus_cyano_chloro_json_w_consensuslineage.txt")
mapfile=("Mapping_file_LS.txt")
trefile=("rep_set_tree.tre")
qiimedata = import_qiime(otufile, mapfile, trefile)
qiimedata
phyloseq-class experiment-level object
otu_table() OTU Table: [ 2206 taxa and 66 samples ]
sample_data() Sample Data: [ 66 samples by 11 sample variables ]
tax_table() Taxonomy Table: [ 2206 taxa by 7 taxonomic ranks ]
phy_tree() Phylogenetic Tree: [ 2206 tips and 2204 internal nodes ]
GenusFiltered <- subset_taxa(qiimedata, Genus != "NA" & Species == "NA")
Error in rownames<-(*tmp*, value = c("sp1", "sp0")) :
length of 'dimnames' [1] not equal to array extent
The package version of phyloseq that I am using in R studio is 1.19.1
Any help to fix this would be much appreciated
Secondly does anyone know if there is way to use the DeSeq2 function in Qiime to compare higher taxonomic ranks, eg Genus? I am aware that I would have to take the results with a pinch of salt.
If there isnt a way is there a way in which I can use Qiime to do permutation based tests, other than at OTU level?
Cheers,
Sarah
Comment