mfbox_rel_paramsweepg.m

toolbox for spm 5 for data, model free analysis

function varargout=mfbox_rel_paramsweepg(varargin)
% MFBOX_REL_PARAMSWEEPG M-file for mfbox_rel_paramsweepg.fig
%
% See also: GUIDE, GUIDATA, GUIHANDLES
%
% Copyright by Peter Gruber, Fabian J. Theis and Ingo R. Keck
% 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
%

% Begin initialization code - DO NOT EDIT
gui_Singleton = 1;
gui_State = struct('gui_Name',mfilename, ...
    'gui_Singleton',gui_Singleton, ...
    'gui_OpeningFcn',@mfbox_rel_paramsweepg_OpeningFcn, ...
    'gui_OutputFcn',@mfbox_rel_paramsweepg_OutputFcn, ...
    'gui_LayoutFcn',[], ...
    'gui_Callback',[]);
if (nargin&&(ischar(varargin{1})))
    gui_State.gui_Callback = str2func(varargin{1});
end

if (nargout)
    [varargout{1:nargout}] = gui_mainfcn(gui_State,varargin{:});
else
    gui_mainfcn(gui_State,varargin{:});
end
% End initialization code - DO NOT EDIT


% --- Executes just before mfbox_rel_paramsweepg is made visible.
function mfbox_rel_paramsweepg_OpeningFcn(hObject,eventdata,handles,varargin)
% varargin{1} contains the mfbss struct
% varargin{2} (optional) contains the parameters
% varargin{3} (optional) contains runflag (0: param edit only, 1: param+run, 2: run only)

handles.output = hObject;
handles.mfbss = varargin{1};

tmethod = {'none','variance'};
if isfield(handles.mfbss,'design')
    if (~isempty(handles.mfbss.design)), tmethod{end+1} = 'design'; end
end
if isfield(handles.mfbss,'reference')
    if (~isempty(handles.mfbss.reference)), tmethod{end+1} = 'reference'; end
end
smethod = {'none','spatial','temporal','temporalspatial', ...
    'spatialvar','temporalvar','temporalspatialvar'};

if (nargin>4), handles.params = varargin{2};
else handles.params = [];
end
if (nargin>5), handles.runflag = varargin{3};
else handles.runflag = 1;
end

set(handles.steps,'String',sprintf('%1.2f',handles.params.steps));
set(handles.range_from1,'String',sprintf('%1.2f',handles.params.range_from1));
set(handles.range_to1,'String',sprintf('%1.2f',handles.params.range_to1));
set(handles.range_from2,'String',sprintf('%1.2f',handles.params.range_from2));
set(handles.range_to2,'String',sprintf('%1.2f',handles.params.range_to2));
set(handles.range_from3,'String',sprintf('%1.2f',handles.params.range_from3));
set(handles.range_to3,'String',sprintf('%1.2f',handles.params.range_to3));
t = find(strcmp(tmethod,handles.params.threshold_type));
if (isempty(t)), t = 1; end
set(handles.thresholdtype,'String',tmethod,'Value',t);
set(handles.threshold,'String',sprintf('%1.4f',handles.params.threshold));
t = find(strcmp(smethod,handles.params.sortby));
if (isempty(t)), t = 1; end
set(handles.sortby,'String',smethod,'Value',t);
switch handles.params.threshold_type
    case 'none'
        set(handles.threshold,'Enable','off');
    otherwise
        set(handles.threshold,'Enable','on');
end

handles.pp = handles.mfbss.params.(handles.mfbss.name);
tf = fieldnames(handles.pp);

m1 = {''};
mp1 = 1;
for i=1:length(tf)
    l = [tf{i},'_type'];
    if (isfield(handles.pp,l)&&any(strcmp(handles.pp.(l),{'int','float','boolean'})))
        m1{end+1} = tf{i};
    end
    if (strcmp(tf{i},handles.params.parameter1))
        mp1 = length(m1);
        tf{i} = '_';
    end
end
set(handles.parameter1,'String',m1,'Value',mp1);
m2 = {''};
mp2 = 1;
for i=1:length(tf)
    l = [tf{i},'_type'];
    if (isfield(handles.pp,l)&&any(strcmp(handles.pp.(l),{'int','float','boolean'})))
        m2{end+1} = tf{i};
    end
    if (strcmp(tf{i},handles.params.parameter2))
        mp2 = length(m2);
        tf{i} = '_';
    end
end
set(handles.parameter2,'String',m2,'Value',mp2);
m3 = {''};
mp3 = 1;
for i=1:length(tf)
    l = [tf{i},'_type'];
    if (isfield(handles.pp,l)&&any(strcmp(handles.pp.(l),{'int','float','boolean'})))
        m3{end+1} = tf{i};
    end
    if (strcmp(tf{i},handles.params.parameter3))
        mp3 = length(m3);
        tf{i} = '_';
    end
end
set(handles.parameter3,'String',m3,'Value',mp3);

handles.ppn1 = m1{mp1};
if (length(handles.ppn1)>0)
    handles.pptype1 = handles.pp.([handles.ppn1,'_type']);
    if (any(strcmp(handles.pptype1,{'float','int'})))
        handles.ppmin1 = handles.pp.([handles.ppn1,'_min']);
        handles.ppmax1 = handles.pp.([handles.ppn1,'_max']);
    else
        handles.ppmin1 = 0;
        handles.ppmax1 = 1;
    end
else
    handles.ppmin1 = 1;
    handles.ppmax1 = 1;
    handles.pptype1 = '';
end
handles.ppn2 = m2{mp2};
if (length(handles.ppn2)>0)
    handles.pptype2 = handles.pp.([handles.ppn2,'_type']);
    if (any(strcmp(handles.pptype2,{'float','int'})))
        handles.ppmin2 = handles.pp.([handles.ppn2,'_min']);
        handles.ppmax2 = handles.pp.([handles.ppn2,'_max']);
    else
        handles.ppmin2 = 0;
        handles.ppmax2 = 1;
    end
else
    handles.ppmin2 = 1;
    handles.ppmax2 = 1;
    handles.pptype2 = '';
end
handles.ppn3 = m3{mp3};
if (length(handles.ppn3)>0)
    handles.pptype3 = handles.pp.([handles.ppn3,'_type']);
    if (any(strcmp(handles.pptype3,{'float','int'})))
        handles.ppmin3 = handles.pp.([handles.ppn3,'_min']);
        handles.ppmax3 = handles.pp.([handles.ppn3,'_max']);
    else
        handles.ppmin3 = 0;
        handles.ppmax3 = 1;
    end
else
    handles.ppmin3 = 1;
    handles.ppmax3 = 1;
    handles.pptype3 = '';
end

guidata(hObject,handles);

uiwait(handles.figure1);


% --- Outputs from this function are returned to the command line.
function varargout=mfbox_rel_paramsweepg_OutputFcn(hObject,eventdata,handles) 
% varargout{1} contains the (modified) parameters

if (isstruct(handles))
    close(handles.figure1);
    varargout{1} = handles.params;
else
    varargout{1} = [];
end


% --- Executes on button press in run.
function run_Callback(hObject,eventdata,handles)

handles.params.steps = str2double(get(handles.steps,'String'));

contents = get(handles.parameter1,'String');
if (iscell(contents))
    handles.params.parameter1 = contents{get(handles.parameter1,'Value')};
    handles.params.range_from1 = str2double(get(handles.range_from1,'String'));
    handles.params.range_to1 = str2double(get(handles.range_to1,'String'));
else
    handles.params.parameter1 = '';
    handles.params.range_from1 = 0;
    handles.params.range_to1 = 0;
end
contents = get(handles.parameter2,'String');
if (iscell(contents))
    handles.params.parameter2 = contents{get(handles.parameter2,'Value')};
    handles.params.range_from2 = str2double(get(handles.range_from2,'String'));
    handles.params.range_to2 = str2double(get(handles.range_to2,'String'));
else
    handles.params.parameter2 = '';
    handles.params.range_from2 = 0;
    handles.params.range_to2 = 0;
end
contents = get(handles.parameter3,'String');
if (iscell(contents))
    handles.params.parameter3 = contents{get(handles.parameter3,'Value')};
    handles.params.range_from3 = str2double(get(handles.range_from3,'String'));
    handles.params.range_to3 = str2double(get(handles.range_to3,'String'));
else
    handles.params.parameter3 = '';
    handles.params.range_from3 = 0;
    handles.params.range_to3 = 0;
end
t = get(handles.sortby,'String');
handles.params.sortby = t{get(handles.sortby,'Value')};
handles.params.threshold = str2double(get(handles.threshold,'String'));
t = get(handles.thresholdtype,'String');
handles.params.threshold_type = t{get(handles.thresholdtype,'Value')};

guidata(hObject,handles);

uiresume(handles.figure1);


% --- Executes on range_from1.
function range_from1_Callback(hObject,eventdata,handles)

val = str2double(get(hObject,'String'));
max_val = str2double(get(handles.range_to1,'String'));
val = min(val,max_val);
val = min(handles.ppmax1,max(handles.ppmin1,val));
if (strcmp(handles.pptype1,'int')) val = floor(val);
elseif (strcmp(handles.pptype1,'float')) val = val;
else val = handles.ppmin1;
end
set(hObject,'String',mat2str(val));

guidata(hObject,handles);


% --- Executes on range_from1.
function range_to1_Callback(hObject, eventdata, handles)

val = str2double(get(hObject,'String'));
min_val = str2double(get(handles.range_from1,'String'));
val = max(val,min_val);
val = min(handles.ppmax1,max(handles.ppmin1,val));
if (strcmp(handles.pptype1,'int')) val = floor(val);
elseif (strcmp(handles.pptype1,'float')) val = val;
else val = handles.ppmax1;
end
set(hObject,'String',mat2str(val));

guidata(hObject,handles);