andersonwinkler / PALM

PALM: Permutation Analysis of Linear Models
70 stars 28 forks source link

Error with -accel tail: «The coordinates of the input points must be finite values; Inf and NaN are not permitted» #52

Closed danielrmrtns closed 8 months ago

danielrmrtns commented 8 months ago

Dear PALM team,

I've been having some troubles with the usage of -accel tail. Sometimes, depending on the input data, the p-value maps are not calculated and the following error gets caught:

Computing p-values.
Error using interp1 (/usr/local/MATLAB/R2018a/toolbox/matlab/polyfun/interp1.p:151)
The coordinates of the input points must be finite values; Inf and NaN are not permitted.
Error in andersondarling (/mnt/data1/fconnect/software/MATLAB/palm-alpha119/palm_pareto.m:196->interp1)
Error in palm_pareto (/mnt/data1/fconnect/software/MATLAB/palm-alpha119/palm_pareto.m:128->andersondarling)
Error in palm_saveall (/mnt/data1/fconnect/software/MATLAB/palm-alpha119/palm_saveall.m:77->palm_pareto)
Error in palm_core (/mnt/data1/fconnect/software/MATLAB/palm-alpha119/palm_core.m:2469->palm_saveall)
Error in palm (/mnt/data1/fconnect/software/MATLAB/palm-alpha119/palm.m:81->palm_core)

Without -accel tail, this does not arise. Do you reckon what could be causing this?

Note: The design is of a standard 2-sample unpaired t-test, with n=30 binarized volumes in the 1st group against n=79 in the 2nd.

andersonwinkler commented 8 months ago

Hi Daniel,

This can happen with unusual data (like you seem to have indeed, binarized) and/or when too few permutations are done.

The Gamma approximation tends to be more robust in these cases ("-approx gamma"). I'd give it a try.

All the best,

Anderson