Error in Step B LDA Effect Size (LEfSe)

Downstream analysis and statistics LEfSe
1

Hi,

I’m trying to proceed to step B. After searching around, I did found someone else having similar issue. The following error was received:

/galaxy_venv/local/lib/python2.7/site-packages/rpy2/rinterface/init.py:185: RRuntimeWarning: Error in La.svd(x, nu, nv) : error code 1 from Lapack routine ‘dgesdd’

warnings.warn(x, RRuntimeWarning)
Traceback (most recent call last):
File “/shed_tools/testtoolshed.g2.bx.psu.edu/repos/george-weingart/lefse/a6284ef17bf3/lefse/run_lefse.py”, line 89, in
if params[‘rank_tec’] == ‘lda’: lda_res,lda_res_th = test_lda_r(cls,feats,class_sl,params[‘n_boots’],params[‘f_boots’],params[‘lda_abs_th’],0.0000000001,params[‘nlogs’])
File “/export/shed_tools/testtoolshed.g2.bx.psu.edu/repos/george-weingart/lefse/a6284ef17bf3/lefse/lefse.py”, line 189, in test_lda_r
z = robjects.r(‘z ← suppressWarnings(lda(as.formula(’+f+’),data=sub_d,tol=’+str(tol_min)+’))’)
File “/galaxy_venv/local/lib/python2.7/site-packages/rpy2/robjects/init.py”, line 359, in call
res = self.eval(p)
File “/galaxy_venv/local/lib/python2.7/site-packages/rpy2/robjects/functions.py”, line 178, in call
return super(SignatureTranslatedFunction, self).call(*args, **kwargs)
File “/galaxy_venv/local/lib/python2.7/site-packages/rpy2/robjects/functions.py”, line 106, in call
res = super(Function, self).call(*new_args, **new_kwargs)
rpy2.rinterface.RRuntimeError: Error in La.svd(x, nu, nv) : error code 1 from Lapack routine ‘dgesdd’

Attached is the data I use.
bf_ptmid_lefse_nosqrt campy.txt (762.5 KB)

Appreciate any help and input.

Thanks,
Xinran

2

Hello Xinran,
LEfSe will often have errors if special characters are used in labels, such as “+” or “-” like you have in your data. I recommend replacing those labels with alphanumeric ones (“positive” and “negative” for instance) and trying to re-run.
Best,
Meg

3

Hi Meg,

Thank you for helping. I tried as you suggested, but I still had the error. Then I figured that the format of cells containing relative abundance data in excel should be all scientific. Where in my data I have some cells in scientific format and others in regular format. Once I uniformed the style, the error seems to be solved.

Thanks,
Xinran

4

Xinran,
Good to hear, and thank you for sharing what fixed the problem!
Best,
Meg