Gneiss: metadata for filtered samples still in ols regression results

I’ve filtered samples from my dataset so only some of the metadata categories remain (eg., filtered Category 4 samples so only Categories 1-3 remain). After running OLS regression in gneiss, the resulting OLS visualization artifact still contains all of the categories instead of just the filtered ones. Moreover, the “Regression Coefficients Summary” interactive table shows p-values (and corrected p-values) for those “filtered” categories as 0 for all balances, and the coefficients are very large (eg., 1e18).

It appears that gneiss is using all factor levels in the full (not-filtered) metadata table.

It seems that I’m going to have to filter my metadata table to match the filtered feature table prior to running qiime gneiss ols-regression. Is there any good way of doing this from the command line (or Jupyter) other than:

  1. export filtered feature table to biom
  2. convert biom to OTU table (or summarize samples)
  3. load metadata table + list of samples into R or python
  4. filter the metadata table based on the list of samples
  5. write out the filtered metadata table

This is a lot of steps just to filter the metadata table.

Hi @nick-youngblut. Could you post your qza files so that we could reproduce this error?

OK. Attached is an artifact with only one taxon in the denominator. The denominator.csv file produced by qiime tools export doesn't include the Feature ID field in the header.

y42.qzv (40.7 KB)

Hi @nick-youngblut. Could you post the full metadata / table? Just want to reproduce the results, and sanity check to make sure that the metadata is really filtered.

This topic was automatically closed 31 days after the last reply. New replies are no longer allowed.