hap.py icon indicating copy to clipboard operation
hap.py copied to clipboard

som.py throws error `cannot call `vectorize` on size 0 inputs unless `otypes` is set`

Open marrip opened this issue 10 months ago • 0 comments

I am trying to run som.py on two vcfs using -R/-Toption to focus only on certain regions. I got:

2024-12-02 13:37:32,097 WARNING  No reference file found at default locations. You can set the environment variable 'HGREF' or 'HG19' to point to a suitable Fasta file.
2024-12-02 13:37:52,964 ERROR    cannot call `vectorize` on size 0 inputs unless `otypes` is set
2024-12-02 13:37:52,964 ERROR    Traceback (most recent call last):
2024-12-02 13:37:52,964 ERROR      File "/usr/local/bin/som.py", line 974, in <module>
2024-12-02 13:37:52,965 ERROR        main()
2024-12-02 13:37:52,965 ERROR      File "/usr/local/bin/som.py", line 865, in main
2024-12-02 13:37:52,965 ERROR        recall = binomialCI(res["tp"], res["tp"]+res["fn"], ci_alpha)
2024-12-02 13:37:52,965 ERROR      File "/usr/local/lib/python2.7/site-packages/numpy/lib/function_base.py", line 2091, in __call__
2024-12-02 13:37:52,966 ERROR        return self._vectorize_call(func=func, args=vargs)
2024-12-02 13:37:52,966 ERROR      File "/usr/local/lib/python2.7/site-packages/numpy/lib/function_base.py", line 2161, in _vectorize_call
2024-12-02 13:37:52,966 ERROR        ufunc, otypes = self._get_ufunc_and_otypes(func=func, args=args)
2024-12-02 13:37:52,966 ERROR      File "/usr/local/lib/python2.7/site-packages/numpy/lib/function_base.py", line 2117, in _get_ufunc_and_otypes
2024-12-02 13:37:52,966 ERROR        raise ValueError('cannot call `vectorize` on size 0 inputs '
2024-12-02 13:37:52,966 ERROR    ValueError: cannot call `vectorize` on size 0 inputs unless `otypes` is set

I did the same comparison with hap.py before and that worked fine. A I using the wrong option to supply the bedfile or what is the cause for this and how can I circumvent it?

marrip avatar Dec 02 '24 13:12 marrip