could you please guide me why I am getting an error?
I am working with the latest version of lefse in conda environment
I already have worked with the tutorial command (below attached from your shared workflow example)
without any error.
Could you please advise if it is preferable to convert my data into relative abundance?
# (updatedR) r01mt19@kitcat:~/intestineFinal/intestineSmartData/Final Final: STEPS/Chall ds mr vs veg/tables$ lefse_format_input.py lefseConda.txt lefseConda.in -c 1 -s 2 -u 3 -o 1000000
# (updatedR) r01mt19@kitcat:~/intestineFinal/intestineSmartData/Final Final: STEPS/Chall ds mr vs veg/tables$ lefse_run.py lefseConda.in lefseConda.res
# Number of significantly discriminative features: 1 ( 1 ) before internal wilcoxon
# R[write to console]: Error in lda.default(x, grouping, ...) :
# variable 1 appears to be constant within groups
#
# R[write to console]: In addition:
# R[write to console]: Warning messages:
#
# R[write to console]: 1: package ‘survival’ was built under R version 4.1.1
#
# R[write to console]: 2: package ‘mvtnorm’ was built under R version 4.1.1
#
# R[write to console]: 3: package ‘coin’ was built under R version 4.1.1
#
# Traceback (most recent call last):
# File "/home/r01mt19/.conda/envs/updatedR/bin/lefse_run.py", line 10, in <module>
# sys.exit(lefse_run())
# File "/home/r01mt19/.conda/envs/updatedR/lib/python3.9/site-packages/lefse/lefse_run.py", line 90, in lefse_run
# 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 "/home/r01mt19/.conda/envs/updatedR/lib/python3.9/site-packages/lefse/lefse.py", line 206, in test_lda_r
# z = robjects.r('z <- suppressWarnings(lda(as.formula('+f+'),data=sub_d,tol='+str(tol_min)+'))')
# File "/home/r01mt19/.conda/envs/updatedR/lib/python3.9/site-packages/rpy2/robjects/__init__.py", line 438, in __call__
# res = self.eval(p)
# File "/home/r01mt19/.conda/envs/updatedR/lib/python3.9/site-packages/rpy2/robjects/functions.py", line 198, in __call__
# return (super(SignatureTranslatedFunction, self)
# File "/home/r01mt19/.conda/envs/updatedR/lib/python3.9/site-packages/rpy2/robjects/functions.py", line 125, in __call__
# res = super(Function, self).__call__(*new_args, **new_kwargs)
# File "/home/r01mt19/.conda/envs/updatedR/lib/python3.9/site-packages/rpy2/rinterface_lib/conversion.py", line 45, in _
# cdata = function(*args, **kwargs)
# File "/home/r01mt19/.conda/envs/updatedR/lib/python3.9/site-packages/rpy2/rinterface.py", line 680, in __call__
# raise embedded.RRuntimeError(_rinterface._geterrmessage())
# rpy2.rinterface_lib.embedded.RRuntimeError: Error in lda.default(x, grouping, ...) :
# variable 1 appears to be constant within groups
lefseConda.txt (223.2 KB)