I have polar angle data only (no eccentricity), and I'm trying to run @RetinoProc. My TR is 1.5 seconds, and this appears to be an issue for some reason.
The 3dRetinoPhase stage is where the problem occurs:
3dRetinoPhase -ccw ccw.rh.niml.dset -clw clw.rh.niml.dset -Tpol 24 -nwedges 1 -pre_stim 6 -phase_estimate FFT -prefix s1008-3.rh
** ERROR: Frequency indices 1116274688 and/or 68 outside max of 1
** ERROR: Failed to set bins.
Is your dataset's TR of 1.500000 sec valid?
++ Warning: 141898 bad floats computed in MAKER_4D_to_typed_fbuc
++ Output dataset ./s1008-3.rh.pol+.niml.dset
** ERROR: Frequency indices 1116274688 and/or 68 outside max of 1
** ERROR: Failed to set bins.
Is your dataset's TR of 1.500000 sec valid?
++ Warning: 141898 bad floats computed in MAKER_4D_to_typed_fbuc
++ Output dataset ./s1008-3.rh.pol-.niml.dset
++ Output dataset ./s1008-3.rh.pol.field.niml.dset
If I change the TR from 1.5 to 2 with 3drefit and adjust the -period_pol accordingly (from 24 to 32 seconds), it runs without complaining and gives me a usable Max.PwR subbrick, whereas MaxPwR can't be used to threshold because nothing survives a value greater than 0 when I keep the TR at 1.5.
Is there a reason why 3dRetinoPhase wouldn't like a 1.5 second TR, or is there a different problem that I'm missing?
Phil