kordk / torch-ecpg

(GPU accelerated) eCpG mapper
BSD 3-Clause "New" or "Revised" License
2 stars 0 forks source link

mlr error: IndexError: index 1 is out of bounds for axis 0 with size 1 #19

Closed rituroy closed 1 year ago

rituroy commented 1 year ago

ritu@pnldev [11:35:45] /songlab/proj/cbi/torch_ecpg $

pwd /songlab/proj/cbi/torch_ecpg (py310)

ritu@pnldev [11:36:02] /songlab/proj/cbi/torch_ecpg $

tecpg run mlr -r 10 -p 0.1 2> mlr.err 1> mlr.out &

ritu@pnldev [11:37:22] /songlab/proj/cbi/torch_ecpg $

cat mlr.err Traceback (most recent call last): File "/home/ritu/anaconda3/envs/py310/bin/tecpg", line 33, in sys.exit(load_entry_point('tecpg', 'console_scripts', 'tecpg')()) File "/home/ritu/anaconda3/envs/py310/bin/tecpg", line 25, in importlib_load_entry_point return next(matches).load() File "/home/ritu/anaconda3/envs/py310/lib/python3.10/importlib/metadata/init.py", line 171, in load module = import_module(match.group('module')) File "/home/ritu/anaconda3/envs/py310/lib/python3.10/importlib/init.py", line 126, in import_module return _bootstrap._gcd_import(name[level:], package, level) File "", line 1050, in _gcd_import File "", line 1027, in _find_and_load File "", line 1006, in _find_and_load_unlocked File "", line 688, in _load_unlocked File "", line 883, in exec_module File "", line 241, in _call_with_frames_removed File "/data/songlab/proj/cbi/torch_ecpg/torch-ecpg/tecpg/main.py", line 9, in main() File "/data/songlab/proj/cbi/torch_ecpg/torch-ecpg/tecpg/main.py", line 6, in main start() File "/data/songlab/proj/cbi/torch_ecpg/torch-ecpg/tecpg/cli.py", line 495, in start cli(obj={}) File "/home/ritu/anaconda3/envs/py310/lib/python3.10/site-packages/click/core.py", line 1130, in call return self.main(args, kwargs) File "/home/ritu/anaconda3/envs/py310/lib/python3.10/site-packages/click/core.py", line 1055, in main rv = self.invoke(ctx) File "/home/ritu/anaconda3/envs/py310/lib/python3.10/site-packages/click/core.py", line 1657, in invoke return _process_result(sub_ctx.command.invoke(sub_ctx)) File "/home/ritu/anaconda3/envs/py310/lib/python3.10/site-packages/click/core.py", line 1657, in invoke return _process_result(sub_ctx.command.invoke(sub_ctx)) File "/home/ritu/anaconda3/envs/py310/lib/python3.10/site-packages/click/core.py", line 1404, in invoke return ctx.invoke(self.callback, ctx.params) File "/home/ritu/anaconda3/envs/py310/lib/python3.10/site-packages/click/core.py", line 760, in invoke return __callback(args, kwargs) File "/home/ritu/anaconda3/envs/py310/lib/python3.10/site-packages/click/decorators.py", line 26, in new_func return f(get_current_context(), *args, *kwargs) File "/data/songlab/proj/cbi/torch_ecpg/torch-ecpg/tecpg/cli.py", line 336, in mlr output = regression_single(args, logger) File "/data/songlab/proj/cbi/torch_ecpg/torch-ecpg/tecpg/regression_single.py", line 216, in regression_single if not filter_p or p_value_np[1] >= p_thresh: IndexError: index 1 is out of bounds for axis 0 with size 1 (py310)

ritu@pnldev [11:38:29] /songlab/proj/cbi/torch_ecpg $

ritu@pnldev [11:44:27] /songlab/proj/cbi/torch_ecpg $

cat mlr.out [INFO] CUDA GPU detected. This device supports CUDA. [INFO] Reading 3 dataframes... [INFOTIMER] Reading 1/3: C.csv [INFO] Reading csv file /data/songlab/proj/cbi/torch_ecpg/data/C.csv with separator , [INFOTIMER] Read 1/3 in 0.003 seconds [INFOTIMER] Reading 2/3: M.csv [INFO] Reading csv file /data/songlab/proj/cbi/torch_ecpg/data/M.csv with separator , [INFOTIMER] Read 2/3 in 0.0045 seconds [INFOTIMER] Reading 3/3: G.csv [INFO] Reading csv file /data/songlab/proj/cbi/torch_ecpg/data/G.csv with separator , [INFOTIMER] Read 3/3 in 0.0042 seconds [INFOTIMER] Finished reading 3 dataframes in 0.0119 seconds. [INFO] Use CPU not supplied. Checking if CUDA is available. [INFO] Using CUDA [INFO] Clearing output directory [INFO] Removing directory /data/songlab/proj/cbi/torch_ecpg/output... [INFO] Creating directory /data/songlab/proj/cbi/torch_ecpg/output... [INFO] Running full regression in expression only mode with "all" region filtration. [INFOTIMER] Converted C to tensor in 1.3383 seconds [INFOTIMER] Created root oneX tensor in 0.0005 seconds [INFOTIMER] Set up output dataframe [INFO] Running with 96 degrees of freedom [INFO] Calculating chunks...

kordk commented 1 year ago

Does this error occur only when using the -p flag? Does it also occur using a higher p-value threshold (e.g., - p 0.8)?

liamgd commented 1 year ago

Fixed in 83539af.

The argument order in cli.py swapped the update and expression_only arguments, so expression_only was always one. See the commit differences for what changed.