ENH: Small p-values are no longer truncated, plus masking on thresholded output
This is an update to 9d9c6025 fixing the following issues:
- A 4D input will be read as 3D, but a warning message will be output to the terminal
- Very small ( e.g. 1x10^-16 ) p-values will be preserved in the thresholded output
- The thresholded output will be masked by both the original data mask and the mask of significant voxels ( in either p or 1-p mode ). Voxels outside the mask will be set to zero
Note the change in the thresholding in line 98 will change the output slightly - but this should be in-line with the improved behaviour for extreme p-values.
Edited by Matthew Webster