There are a few points that come to my mind.
- You only have 16 samples (since a paired test is equivalent to a 1-sample test using A-B as the inputs). The program allows that, but it is really too few for good statistics in the permutation/randomization. That may be one cause of your lack of results.
- Was any blurring applied to the data during the processing that lead to the 3dttest++ run? That is, if you used afni_proc.py, did you specify blurring in the time series analysis? I ask because the amount of blurring you use is very low (0 and 2 mm) -- if this is the only blurring ever applied, it is too small to show FMRI results in most cases.
- I would use a wider range of pthr values, for example 'pthr=0.01/0.001/19' to specify 19 thresholds evenly spaced from 0.01 to 0.001, with step size (0.01-0.001)/(19-1)=0.0005. Using a range of pthr values is the main point of ETAC, and is what I was thinking of when I developed the idea while hiking in the Grand Canyon.
- Although using multiple hpow values seemed like a good idea to me at the time, it doesn't actually seem to make much difference. I now usually just use 'hpow=2'
I hope this helps, and apologize for the delay. If you have any followup questions, please email me. (The homework is to somehow find my email address.)