Skip to content

issue with wilcoxon signed rank test #24

@lauramason326

Description

@lauramason326

Hi
I am trying to run lefse on an OTU dataset, but it keeps failing at the wilcoxon signed rank step.
Do you have any advice?

Thanks
Laura

lefse_run.py OTU_Sample2.in samples_16S_otu.res


Number of significantly discriminative features: 5307 ( 84201 ) before internal wilcoxon
R[write to console]: Error in .Primitive("[")(c(72519.717356904, 53305.4470592899, 43080.6615065647,  :
  subscript out of bounds
Traceback (most recent call last):
  File "/fs/ess/PAS1212/bioinformatic_tools/miniconda3/envs/lefse-1.1.2/bin/lefse_run.py", line 10, in <module>
    sys.exit(lefse_run())
  File "/fs/ess/PAS1212/bioinformatic_tools/miniconda3/envs/lefse-1.1.2/lib/python3.7/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 "/fs/ess/PAS1212/bioinformatic_tools/miniconda3/envs/lefse-1.1.2/lib/python3.7/site-packages/lefse/lefse.py", line 225, in test_lda_r
    res = dict([(pp,[float(ff) for ff in rres.rx(pp,True)] if pp in rowns else [0.0]*lenc ) for pp in [p[0],p[1]]])
  File "/fs/ess/PAS1212/bioinformatic_tools/miniconda3/envs/lefse-1.1.2/lib/python3.7/site-packages/lefse/lefse.py", line 225, in <listcomp>
    res = dict([(pp,[float(ff) for ff in rres.rx(pp,True)] if pp in rowns else [0.0]*lenc ) for pp in [p[0],p[1]]])
  File "/fs/ess/PAS1212/bioinformatic_tools/miniconda3/envs/lefse-1.1.2/lib/python3.7/site-packages/rpy2/robjects/vectors.py", line 79, in __call__
    res = fun(*conv_args, **kwargs)
  File "/fs/ess/PAS1212/bioinformatic_tools/miniconda3/envs/lefse-1.1.2/lib/python3.7/site-packages/rpy2/rinterface_lib/conversion.py", line 45, in _
    cdata = function(*args, **kwargs)
  File "/fs/ess/PAS1212/bioinformatic_tools/miniconda3/envs/lefse-1.1.2/lib/python3.7/site-packages/rpy2/rinterface.py", line 680, in __call__
    raise embedded.RRuntimeError(_rinterface._geterrmessage())
rpy2.rinterface_lib.embedded.RRuntimeError: Error in .Primitive("[")(c(72519.717356904, 53305.4470592899, 43080.6615065647,  :
  subscript out of bounds

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions