Revision 7

Email me for help: aso.toshihiko@gmail.com

Now the scripts use niimath if it exists. If not, you must specify the path to fslmaths/FSL5.
Updated to work without any previous subfunctions.

Example usage:

• Modify Line 56 to capture your fMRI directories and run
  Einsteining_v05( 'Sbj0001', '0.72', '1200', '9', '0.3', '0', '8', '0')
  or Einsteining_v05( 'Sbj0001', 0.72, 1200, 9, 0.3, 0, 8, 0)

function Einsteining_v05( Sdir, TR, Nvols, MaxLag, MinR, Fixed, Sm, onlyLag)
% Sdir: Subject directory name
% TR = in seconds
% Number of volumes per run
% MaxLag = 8 -> lag mapping up to -8TR / +8TR
% MinR = 0.2 -> Cross-correlogram peak < 0.2 will not be used (leaving such voxels as NaN)
% Fixed = 0 -> Recursive (flexible) LFO instead of Fixed-sLFO algorithm
% Spatial smoothing at 8mm FWHM
% onlyLag = 0 -> Lag mapping + deperfusioning (creating clean 4D dataset without perfusion lag structure)