gneiss ols-regression error "SVD did not converge"

Hello,

Can someone enlighten me on the meaning of this error, and what I might do to resolve it?

Thanks.

CLI command (qiime2-2019.7):
qiime gneiss ols-regression --verbose --i-table input/M_balances.qza --i-tree input/M_hierarchy.qza --m-metadata-file input/GCMP_EMP_map_r28.txt --p-formula “temperature+host_scientific_name” --o-visualization output/M_regression_summary.qzv

Red text:
Plugin error from gneiss:
SVD did not converge

Verbose output:
Traceback (most recent call last):
File “/home/dylan/miniconda3/envs/qiime2-2019.7/lib/python3.6/site-packages/q2cli/commands.py”, line 327, in call
results = action(**arguments)
File “</home/dylan/miniconda3/envs/qiime2-2019.7/lib/python3.6/site-packages/decorator.py:decorator-gen-278>”, line 2, in ols_regression
File “/home/dylan/miniconda3/envs/qiime2-2019.7/lib/python3.6/site-packages/qiime2/sdk/action.py”, line 240, in bound_callable
output_types, provenance)
File “/home/dylan/miniconda3/envs/qiime2-2019.7/lib/python3.6/site-packages/qiime2/sdk/action.py”, line 445, in callable_executor
ret_val = self._callable(output_dir=temp_dir, **view_args)
File “/home/dylan/miniconda3/envs/qiime2-2019.7/lib/python3.6/site-packages/q2_gneiss/regression/_regression.py”, line 34, in ols_regression
ols_summary(output_dir, res, tree)
File “/home/dylan/miniconda3/envs/qiime2-2019.7/lib/python3.6/site-packages/gneiss/plot/_regression_plot.py”, line 288, in ols_summary
_k, _l = model.kfold(), model.lovo()
File “/home/dylan/miniconda3/envs/qiime2-2019.7/lib/python3.6/site-packages/gneiss/regression/_ols.py”, line 472, in lovo
res_i.fit(**kwargs)
File “/home/dylan/miniconda3/envs/qiime2-2019.7/lib/python3.6/site-packages/gneiss/regression/_ols.py”, line 174, in fit
inv = np.linalg.pinv(np.dot(X.T, X))
File “<array_function internals>”, line 6, in pinv
File “/home/dylan/miniconda3/envs/qiime2-2019.7/lib/python3.6/site-packages/numpy/linalg/linalg.py”, line 1970, in pinv
u, s, vt = svd(a, full_matrices=False, hermitian=hermitian)
File “<array_function internals>”, line 6, in svd
File “/home/dylan/miniconda3/envs/qiime2-2019.7/lib/python3.6/site-packages/numpy/linalg/linalg.py”, line 1636, in svd
u, s, vh = gufunc(a, signature=signature, extobj=extobj)
File “/home/dylan/miniconda3/envs/qiime2-2019.7/lib/python3.6/site-packages/numpy/linalg/linalg.py”, line 106, in _raise_linalgerror_svd_nonconvergence
raise LinAlgError(“SVD did not converge”)
numpy.linalg.LinAlgError: SVD did not converge

2 Likes

Hi @dylan, when this error shows up, it implies that there is something weird going on with the metadata.

How many samples are there? How many hosts are there? Are there any categorical values that snuck into your temperature column? Are there missing entries in your metadata?

There is also a relevant post here that may help with debugging: Gneiss Singular matrix error

3 Likes

Thanks for your help!

Temperature was not in fact all numeric, so I removed the offenders and it fixed the issue.

5 Likes