Skip to content

issue with wilcoxon signed rank test #24

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

    Fields

    No fields configured for issues without a type.

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions