Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

mztab export fails for particular dataset from user. #427

Open
ypriverol opened this issue Oct 9, 2024 · 0 comments
Open

mztab export fails for particular dataset from user. #427

ypriverol opened this issue Oct 9, 2024 · 0 comments
Assignees
Labels
bug Something isn't working

Comments

@ypriverol
Copy link
Member

Description of the bug

Reported here: https://nfcore.slack.com/archives/C02Q3FL29PD/p1728424167422049

Command error:
  2024-10-08 21:43:49,054 [mztab_prh] - Calculating protein coverage (bottleneck)...
  2024-10-08 21:44:04,074 [mztab_prh] - Getting ambiguity members...
  2024-10-08 21:44:04,087 [mztab_prh] - Matching PRH to best search engine score...
  2024-10-08 21:44:04,246 [mztab_prh] - Matching PRH to modifications...
  2024-10-08 21:44:04,263 [mztab_prh] - Matching PRH to protein quantification...
  2024-10-08 21:44:04,323 [mztab_peh] - Constructing PEH sub-table...
  2024-10-08 21:44:04,323 [mztab_peh] - report.shape: (23503, 23),  pr.shape: (7377, 15), len(precursor_list): 8833, index_ref.shape: (5, 6)
  2024-10-08 21:44:04,326 [mztab_peh] - Finding modifications...
  2024-10-08 21:44:04,376 [mztab_peh] - Extracting sequence...
  2024-10-08 21:44:04,425 [mztab_peh] - Checking accession uniqueness...
  2024-10-08 21:44:04,464 [mztab_peh] - Matching precursor IDs...
  2024-10-08 21:44:04,469 [mztab_peh] - Getting scores per run
  2024-10-08 21:44:04,493 [mztab_peh] - Getting peptide abundances per study variable
  2024-10-08 21:44:04,524 [mztab_peh] - Getting peptide properties...
  2024-10-08 21:44:04,536 [mztab_peh] - Re-ordering columns...
  2024-10-08 21:44:04,569 [mztab_psh] - Constructing PSH sub-table
  /usr/local/lib/python3.10/site-packages/quantmsutils/diann/diann2mztab.py:1089: FutureWarning: In a future version of pandas, a length 1 tuple will be returned when iterating over a groupby with a grouper equal to a list of length 1. Don't supply a list with a single grouper to avoid this warning.
    for n, group in report.groupby(["Run"]):
  Traceback (most recent call last):
    File "/usr/local/bin/quantmsutilsc", line 10, in <module>
      sys.exit(main())
    File "/usr/local/lib/python3.10/site-packages/quantmsutils/quantmsutilsc.py", line 38, in main
      cli()
    File "/usr/local/lib/python3.10/site-packages/click/core.py", line 1157, in __call__
      return self.main(*args, **kwargs)
    File "/usr/local/lib/python3.10/site-packages/click/core.py", line 1078, in main
      rv = self.invoke(ctx)
    File "/usr/local/lib/python3.10/site-packages/click/core.py", line 1688, in invoke
      return _process_result(sub_ctx.command.invoke(sub_ctx))
    File "/usr/local/lib/python3.10/site-packages/click/core.py", line 1434, in invoke
      return ctx.invoke(self.callback, **ctx.params)
    File "/usr/local/lib/python3.10/site-packages/click/core.py", line 783, in invoke
      return __callback(*args, **kwargs)
    File "/usr/local/lib/python3.10/site-packages/click/decorators.py", line 33, in new_func
      return f(get_current_context(), *args, **kwargs)
    File "/usr/local/lib/python3.10/site-packages/quantmsutils/diann/diann2mztab.py", line 181, in diann2mztab
      diann_directory.convert_to_mztab(
    File "/usr/local/lib/python3.10/site-packages/quantmsutils/diann/diann2mztab.py", line 431, in convert_to_mztab
      psh = mztab_psh(report, str(self.base_path), database)
    File "/usr/local/lib/python3.10/site-packages/quantmsutils/diann/diann2mztab.py", line 1116, in mztab_psh
      pd.merge_asof(group, target, on="RT.Start", direction="nearest"),
    File "/usr/local/lib/python3.10/site-packages/pandas/core/reshape/merge.py", line 600, in merge_asof
      return op.get_result()
    File "/usr/local/lib/python3.10/site-packages/pandas/core/reshape/merge.py", line 1742, in get_result
      join_index, left_indexer, right_indexer = self._get_join_info()
    File "/usr/local/lib/python3.10/site-packages/pandas/core/reshape/merge.py", line 1026, in _get_join_info
      (left_indexer, right_indexer) = self._get_join_indexers()
    File "/usr/local/lib/python3.10/site-packages/pandas/core/reshape/merge.py", line 2042, in _get_join_indexers
      raise ValueError(f"Merge keys contain null values on {side} side")
  ValueError: Merge keys contain null values on right side

Command used and terminal output

No response

Relevant files

No response

System information

No response

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
bug Something isn't working
Projects
None yet
Development

No branches or pull requests

2 participants