lindenmp / rs-fmri Goto Github PK

Linden,

Apologies if this is a non-issue, but I'm getting values of "NaN" as output for QCFC_DistDep and QCFC_DistDep_Pval.

Not sure if its what it returns when there isn't anything of substance or the way my data is organized so I was wondering if you wouldn't mind taking a look when and if you have a minute.

Script is enclosed with the appropriate .mat files so it should run with minimal to no edits.

Troubleshoot_Nan_Script.zip

I think there may have been a small mix-up with the oversampling and maximum allowable frequency factor given in func/JP14_getTransform.m

if nargin < 5
		hifac = 8;
	end

	if nargin < 4
		ofac = 1;
	end

Shouldn't the hifac value be set to 1 here? Otherwise the sampling frequencies that are computed in the line below

f = (1/(T*ofac):1/(T*ofac):hifac*N/(2*T)).';

Would exceed the nyquist limit of the sampler. In Power's original code these are swapped (ofac=8, hfac=1)