hap.py
hap.py copied to clipboard
som.py throws error `cannot call `vectorize` on size 0 inputs unless `otypes` is set`
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?