mfbox_infomap.m

toolbox for spm 5 for data, model free analysis

function [grid,newdata]=mfbox_infomap(timeline,data,oldgrid,newsd,newtl,infometh,maxiter,nbhfunc,nbhdeg)
% usage: [grid,newdata]=mfbox_infomap(timeline,data,oldgrid,newsd,newtl,infometh,maxiter,nbhfunc,nbhdeg)
%
% Copyright by Peter Gruber
% 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(4,9,nargin));

if (nargin<4), newtl = linspace(timeline(1),timeline(end),length(timeline)); end
if (nargin<5), infometh = 'var'; end
if (nargin<6), maxiter = 1000; end
if (nargin<7), nbhfunc = 'exp'; end
if (nargin<8), nbhdeg = 2; end

sd = size(data);
data = reshape(data,[],sd(end));

if (length(timeline)~=size(data,2)), error('Wrong data dimension'); end

grid = diag(1./(newsd-1))*(mfbox_mkgrid(newsd)'-1);
minog = min(oldgrid,[],2);
maxog = max(oldgrid,[],2);
grid = diag(maxog-minog)*grid+repmat(minog,1,size(grid,2));

switch infometh
    case 'var'
        info = var(data,0,2);
        grid = mfbox_som(grid',newsd,oldgrid',info,'maxiter',maxiter)';
    case 'kurt'
        info = kurtosis(data');
        grid = mfbox_som(grid',newsd,oldgrid',info,'maxiter',maxiter)';
end

newdata = zeros(prod(newsd),length(newtl),'single');
for t=1:length(newtl)
    m = max(find(timeline<=newtl(t)));
    if (timeline(m)==newtl(t))
        newdata(:,t) = mfbox_transgrid(oldgrid,data(:,m),grid);
    else
        r = timeline(m:(m+1));
        og = [repmat(oldgrid,1,2);reshape(repmat([1,2],size(oldgrid,2),1),1,[])];
        ng = [grid;repmat((newtl(t)-r(1))/(r*[-1;1]),1,size(grid,2))];
        newdata(:,t) = mfbox_transgrid(og,reshape(data(:,m:(m+1)),[],1),ng);
    end
end

newdata = reshape(newdata,[newsd,length(newtl)]);

return