NaNs produced

[ialmodovar]

You have addressed the underflow issue before and made the updated using version 0.2.2.1. I'm currently using pTFCE version 0.2.2.1 but I'm still getting NaNs. The files (at the bottom), which are of Zmaps and the mask, produce NaNs for 1470 voxels.

I did as follows,

library("oro.nifti")
library("pTFCE")
Z <- readNIfTI(fname="Pseudo-Zmap.nii.gz")
mask <- readNIfTI(fname="Pseudo-mask.nii.gz")
pp <- ptfce(-Z,mask)
* Estimating smoothness based on the data...
  |======================================================================| 100%
* Performing pTFCE...
  |======================================================================| 100%
There were 50 or more warnings (use warnings() to see the first 50)
> warnings()
1: In sqrt(d * (8 * s + d)) : NaNs produced
2: In sqrt(d * (8 * s + d)) : NaNs produced
3: In sqrt(d * (8 * s + d)) : NaNs produced
4: In sqrt(d * (8 * s + d)) : NaNs produced
5: In sqrt(d * (8 * s + d)) : NaNs produced
...

I was hoping you could help me with this?

Many thanks,

Israel

Pseudo-Zmap.nii.gz
Pseudo-mask.nii.gz

[maitra]

I can verify these warnings also. Strange.