mfbox_pos_informative_run.m

toolbox for spm 5 for data, model free analysis

function [A,W,S,params,reliability]=mfbox_pos_informative_run(A,W,S,X,params,runflag)
% select and label components
%
% Usage:
%  [A,W,S,params,reliability]=mfbox_pos_informative_run(A,W,S,X,params,runflag)
%
%  X           - cell with
%                (NxT) data
%                (Tx3) grid
%                timeline
%                mask
%                design
%                reference
%  A           - (TxD) mixing matrix
%  W           - (DxT) demixing matrix
%  S           - struct with
%                dat              - (MxD) components (M<=N unique datapoints)
%                mask             - mask selecting unique subset of grid
%                part             - cell array containing the different parts 
%                                   of the analysis
%                map              - (TxC) maps (timestep,component) to element in
%                                   1...D (D>=C)
%                spatialnoisevar  - (Mx1) spatial variance of the noise
%                temporalnoisevar - (Tx1) temporal variance of the noise
%                names            - (Cx1) cell with names
%  params      - struct
%  reliability - struct
%  runflag     - -1 get default parameter
%                 0 interactive ask parameters
%                 1 interactive ask parameters and run
%                 2 run
%
% Copyright by Peter Gruber and Fabian J. Theis
% Signal Processing & Information Theory group
% Institute of Biophysics, University of Regensburg, Germany
% Homepage: http://research.fabian.theis.name
%           http://www-aglang.uni-regensburg.de
%
% This file is free software, subject to the 
% GNU GENERAL PUBLIC LICENSE, see gpl.txt

error(nargchk(3,6,nargin));
error(nargchk(1,5,nargout));

if (nargin<4), X = []; end
if (nargin<5), params = []; end
if (nargin<6), runflag = 1; end

params = mfbox_checkparam(params,'pos','informative',struct());
reliability = struct();

if (runflag>0 && isstruct(params))
    timeline = [];
    if (length(X)>3)
        timeline = X{4};
    end
    design = [];
    if (length(X)>4)
        design = X{5};
    end
    nds = [];
    if (~isempty(design)) 
        sd = size(design);
        nds = zeros(length(timeline),sd(2));
        for i=1:sd(2)
            nds(:,i) = interp1(linspace(timeline(1),timeline(end), ...
                sd(1)),design(:,i),timeline)';
        end
        nds = nds-repmat(mean(nds),size(nds,1),1);
        nds = nds./repmat(std(nds),size(nds,1),1);
    end
    [reliability.mdllhdesign,reliability.mdllhtransient,reliability.mdllhnoise] = ...
        mfbox_getinformative(double(A),nds);
end