Hi all,
I am using RetroTS.py to turn respiratory and cardiac tsv.gz files into regressors for 1st level analysis.
I kept getting this error and any suggestion is greatly appreciated! 
[bash input: target physio file, json file [has p(frequency) = 250], v (TR) = 2, n (slice) = 72, output prefix]
RetroTS.py
-phys_file sub-664_ses-001_task-reward_recording-cardiac_physio.tsv.gz
-phys_json sub-664_ses-001_task-reward_recording-cardiac_physio.json
-v 2
-n 72
-prefix physio_card_regressors.1D
#output
Traceback (most recent call last):
File “/Users/xiy/abin/RetroTS.py”, line 687, in
retro_ts(
File “/Users/xiy/abin/RetroTS.py”, line 279, in retro_ts
cardiac_peak, error = peak_finder(cardiac_info,v=phys_dat[k])
File “/Users/xiy/abin/lib_RetroTS/PeakFinder.py”, line 517, in peak_finder
r[“p_trace_r”] = interp1d(
File “/Users/xiy/opt/anaconda3/lib/python3.8/site-packages/scipy/interpolate/interpolate.py”, line 467, in init
self._y = self._reshape_yi(self.y)
File "/Users/xiy/opt/anaconda3/lib/python3.8/site-packages/scipy/interpolate/polyint.py", line 110, in _reshape_yi
** return yi.reshape((yi.shape[0], -1))**
ValueError: cannot reshape array of size 0 into shape (0,newaxis)
- I am following the suggestions from this post: So what are people using to process their BIDS physio files in to regressors for Fitlins?
- Same error with respiratory files. Specifying p instead of json file yielded the same error.
Upon inspection, physio files look normal (two columns: timepoints, amplitude).
Please let me know if a physio file needs to be uploaded to replicate this error.
Thank you all!