mfbox_fastica_run.m

toolbox for spm 5 for data, model free analysis

function [params,A,W,S]=mfbox_fastica_run(X,mask,params,runflag)
% Copyright by Ingo R. Keck
% Signal Processing & Information Theory group
% Institute of Biophysics, University of Regensburg, Germany
% Homepage: http://ingokeck.de/
%           http://www-aglang.uni-regensburg.de
%
% This file is free software, subject to the 
% GNU GENERAL PUBLIC LICENSE, see gpl.txt

error(nargchk(1,4,nargin));
error(nargchk(1,4,nargout));

A = [];
W = [];
S = [];

s = size(X);
dim = s(1:(end-1));
if (length(dim)==1), dim = [dim,1]; end
n = s(end);

if (nargin<2), mask = logical(ones(dim)); end
if (nargin<3), params = []; end
if (nargin<4), runflag = 1; end

params = mfbox_checkparam(params,'mf','fastica', ...
    struct('approach','symm','numComponents',10,'numComponents_min',1, ...
    'numComponents_max',200,'numComponents_type','int','g','tanh',...
    'stabilization','on','firstEig',1,'lastEig',200,'verbose','off'));

%params

if (abs(runflag-0.5)<1)
    if (exist('OCTAVE_HOME')~=5)% matlab
        params = mfbox_fasticag(X,mask,params,runflag);
    elseif (exist('OCTAVE_HOME')==5) %octave
        [params,runflag] = mfbox_getparam(params,runflag);
    end
end

if (runflag>0 && isstruct(params))
    [S,A,W] = mfbox_fastica(X',params);
    S=S';
end