title: Error code for tool subset_taxa; Error in dimnames(x) <- dn : length of 'dimnames' [1] not equal to array extent
body: <p>I load .biom data into R using library phyloseq.</p>
  <pre><code>Data2020 &lt;- import_biom(&quot;XXX&quot;)
  mapfile2020 &lt;- import_qiime_sample_data(&quot;XXX&quot;)
  tree &lt;- read_tree(&quot;XXX&quot;)
  
  mydata &lt;- merge_phyloseq(Data2020, mapfile2020, tree)
  colnames(tax_table(mydata)) &lt;- c(&quot;Kingdom&quot;, &quot;Phylum&quot;, &quot;Class&quot;, &quot;Order&quot;, &quot;Family&quot;, &quot;Genus&quot;, &quot;Species&quot;) 
  
  ps &lt;- mydata 
  ps1 &lt;- subset_taxa(ps, Phylum != &quot;NA&quot;)
  </code></pre>
  <p>Subset to the data I want.</p>
  <pre><code>WRB_data &lt;- subset_samples(ps1, Site==&quot;WRB&quot;)
  </code></pre>
  <p>Then Try to Subset by a taxa.</p>
  <pre><code>WRB_data2 &lt;- subset_taxa(WRB_data, Phylum==&quot;Acidobacteria&quot;)
  </code></pre>
  <p>Always results in this error. Please Help!?</p>
  <pre><code>Error in dimnames(x) &lt;- dn : length of 'dimnames' [1] not equal to array extent
  </code></pre>