MAP

% multiThreshold = 'experimenter GUI for multiThreshold

% allows the experimenter to design experiments.

% The *running* of experiments is left to subjGUI.m

%

% There are three kinds of experiments known as 'earOptions':

% 1. Measurements using real listeners

%  'left', 'right',  'diotic', 'dichoticLeft', 'dichoticRight'

% 2. Measurements using the MAP model as the subject

%   'MAPmodelListen',  'MAPmodelMultiCh', 'MAPmodelSingleCh'

% 3. Monte Carlo simulations

%   'statsModelLogistic','statsModelRareEvent'

%

% There are many stimulus configurations relating to different measurement

% requirements. These configurations are defined in paradigm files located

% in the 'paradigms' folder with names of the form paradigm_<name>.m. These

% files specify values in the stimulusParameter and experiment structures.

% Some minor parameters are specified in the intialiseGUI function below.

% Each configuration is only a start up arrangement and the user can modify

% the parameters on the GUI itself.

%

% the 'RUN' button initiates the measurements and hands control over to the

% subjGUI program. When the measurements are complete control is handed

% back and the stack unwinds without any further action

function varargout = multiThreshold(varargin)

%MULTITHRESHOLD M-file for multiThreshold.fig

%      MULTITHRESHOLD, by itself, creates a new MULTITHRESHOLD or raises the existing

%      singleton*.

%

%      H = MULTITHRESHOLD returns the handle to a new MULTITHRESHOLD or the handle to

%      the existing singleton*.

%

%      MULTITHRESHOLD('Property','Value',...) creates a new MULTITHRESHOLD using the

%      given property value pairs. Unrecognized properties are passed via

%      varargin to multiThreshold_OpeningFcn.  This calling syntax produces a

%      warning when there is an existing singleton*.

%

%      MULTITHRESHOLD('CALLBACK') and MULTITHRESHOLD('CALLBACK',hObject,...) call the

%      local function named CALLBACK in MULTITHRESHOLD.M with the given input

%      arguments.

%

%      *See GUI Options on GUIDE's Tools menu.  Choose "GUI allows only one

%      instance to run (singleton)".

%

% See also: GUIDE, GUIDATA, GUIHANDLES

% Edit the above text to modify the response to help multiThreshold

% Last Modified by GUIDE v2.5 25-Oct-2011 07:51:14

% Begin initialization code - DO NOT EDIT

gui_Singleton = 1;

gui_State = struct('gui_Name',       mfilename, ...

    'gui_Singleton',  gui_Singleton, ...

    'gui_OpeningFcn', @multiThreshold_OpeningFcn, ...

    'gui_OutputFcn',  @multiThreshold_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

% --------------------------------------------------  multiThreshold_OpeningFcn

function multiThreshold_OpeningFcn(hObject, eventdata, handles, varargin)

% Choose default command line output for multiThreshold

handles.output = hObject;

cla(handles.axes1)

cla(handles.axes2)

cla(handles.axes4)

cla(handles.axes5)

% Update handles structure

guidata(hObject, handles);

function varargout = multiThreshold_OutputFcn(hObject, eventdata, handles)

% Get default command line output from handles structure

initializeGUI(handles)

varargout{1} = handles.output;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% -----------------------------------------------------initializeGUI

function initializeGUI(handles)

% Populate the edit boxes and popup menus on the GUI

% Then wait for user action

global stimulusParameters experiment betweenRuns

global targetTypes maskerTypes backgroundTypes

global variableNames paradigmNames threshEstNames  cueNames betweenRunsVariables

% Specify order of fields in main structures

% identify as empty values or empty strings only

orderGlobals

addpath ('paradigms')   %paradigm informations is stored here

addpath (['..' filesep 'testPrograms']) % model physiology tests

% specify all variables that  need to be set on the GUI

variableNames={'stimulusDelay','maskerDuration','maskerLevel',...

    'maskerRelativeFrequency', 'targetFrequency', 'gapDuration',...

    'targetDuration','targetLevel','rampDuration',...

    'cueTestDifference', 'WRVstartValues', 'WRVsteps', 'WRVlimits',...

    'OHIOnTones'};

% Variable variables

%  (names of variable that can changed between runs)

betweenRunsVariables={'stimulusDelay','maskerDuration','maskerLevel',...

    'maskerRelativeFrequency','targetFrequency', 'gapDuration',...

    'targetDuration','targetLevel','OHIOnTones'};

% populate the 'between runs variable' menus

set(handles.popupmenuVaryParameter1,'string',betweenRunsVariables)

set(handles.popupmenuVaryParameter2,'string',betweenRunsVariables)

% Trial presentation order - randomize at startup

presentationOrderNames={'randomize within blocks', 'fixed sequence', ...

    'randomize across blocks'};

set(handles.popupmenuRandomize,'string', presentationOrderNames)

% targetType- value must be set in paradigm

targetTypes={'tone','noise', 'pinkNoise','whiteNoise','24TalkerBabble',...

    'speech', 'Twister','digitStrings'};

set(handles.popupmenuTargetType, 'string', targetTypes);

% maskerType - value must be set in paradigm

maskerTypes={'tone','noise', 'pinkNoise','TEN','whiteNoise','24TalkerBabble', ...

    'speech'};

set(handles.popupmenuMaskerType, 'string', maskerTypes);

% background Type- value must be set in paradigm (default = 1, 'none')

backgroundTypes={'none','noise', 'pinkNoise', 'TEN','noiseDich',...

    'pinkNoiseDich','whiteNoise','24TalkerBabble',...

    '16TalkerBabble','8TalkerBabble','4TalkerBabble',...

    '4TalkerReversedBabble','2TalkerBabble','1TalkerBabble'};

set(handles.popupmenuBackgroundType, 'string', backgroundTypes);

set(handles.editBackgroundLevel,'string', '0')

% Establish available paradigms by scanning paradigms folder

paradigmNames= what('paradigms');

paradigmNames= paradigmNames.m; % select m files only

idx=strmatch('paradigm_', paradigmNames); % with 'paradigm_'

paradigmNames=paradigmNames(idx);

for i=1:length(paradigmNames)   % strip off file extension

    paradigmNames{i}=paradigmNames{i}(10:end-2);

end

set(handles.popupmenuParadigm,'string', paradigmNames)

% startup paradigm, 'training' (could be anywhere on the list)

startupParadigm='training';

idx= find(strcmp(paradigmNames,startupParadigm));

if ~isempty(idx)

    set(handles.popupmenuParadigm,'value', idx)

else

    % training paradigm must exist

    error(['expGUI_MT\initializeGUI: No ' startupParadigm...

        ' paradigm found in paradigms folder'])

end

earOptions={'left', 'right',  'diotic', 'dichoticLeft', 'dichoticRight',...

    'MAPmodelListen',  'MAPmodelMultiCh', 'MAPmodelSingleCh'...

    'statsModelLogistic','statsModelRareEvent'};

set(handles.popupmenuEar,'string', earOptions)

defaultOption=1;

experiment.ear=earOptions{defaultOption};

set(handles.popupmenuEar,'value', defaultOption)    % 'left' is deafult

set(handles.pushbuttonSingleShot, 'visible', 'off') % use only for MAP

% phase

phaseOptions={'sin','cos','alt','rand'};

set(handles.popupmenuPhase,'string', phaseOptions)

set(handles.popupmenuPhase,'value', 1)

% Cue

cueNames={'cued', 'noCue'};

set(handles.popupmenuCueNoCue, 'string', cueNames);

set(handles.popupmenuCueNoCue, 'value', 1);

% threshold assessment method - value must be set in paradigm

threshEstNames={'oneIntervalUpDown', 'MaxLikelihood', ...

    '2I2AFC++', '2I2AFC+++'};

set(handles.popupmenuThreshEst, 'string', threshEstNames);

experiment.stopCriteria2IFC=[75 3 5];

experiment.stopCriteriaSI=20;

% ** editBoxes that are only set by hand

% music (relative) level, 'tada!' (manual setting only)

% increase for the hard of hearing

set(handles.editMusicLevel,'string','0')

% Catch Trial Rate

set(handles.editCatchTrialRate,'string','0.2   0.1  ');

% calibration

stimulusParameters.restoreCalibration=7;

set(handles.editcalibrationdB,'string',...

    stimulusParameters.restoreCalibration)

% a MAPplot

experiment.MAPplot=0;	%default

set(handles.editMAPplot,'string',num2str(experiment.MAPplot))

% saveData

experiment.saveData=1;

set(handles.editSaveData,'string',num2str(experiment.saveData))

set(handles.pushbuttonSave,'visible','off')

% printTracks

experiment.printTracks=0;

set(handles.editPrintTracks,'string',num2str(experiment.printTracks))

% standard delay between button press and initiation of next stimulus

experiment.clickToStimulusPause=1;

% buttonBoxType: NB no problem if mouse is used instead

experiment.buttonBoxType='square';

% estimation of the mean

% {'logisticLS', 'logisticML', 'rareEvent'}

experiment.functionEstMethod='logisticLS';

% message box

set(handles.textMSG,'backgroundcolor', 'w', 'ForegroundColor', 'b', 'string', '')

set(handles.editMsgFont,'string','10')

set(handles.editSubjectFont,'string','14')

% default psychometric bin size and logistic slopes

experiment.psyBinWidth=1;  % dB

maxLogisticK=2; % steepest slope contemplated

experiment.maxLogisticK=maxLogisticK;

experiment.numPossLogisticK=100;

experiment.possLogSlopes= ...

    0.01: maxLogisticK/experiment.numPossLogisticK: maxLogisticK;

experiment.meanSearchStep=0.25; % dB

% indicate that this is the first run and the GUI should be located in

% default location.

experiment.justInitialized=1;

% set up GUI based on training paradigm

aParadigmSelection(handles);

earSetUp(handles)

aThresholdAssessmentMethod(handles)

aShowRelevantObjects(handles)

% Done. Now wait for action

% -------------------------------------------- popupmenuMaskerType_Callback

function popupmenuMaskerType_Callback(hObject, eventdata, handles)

% show or remove masker frequency box

aShowRelevantObjects(handles)

% -------------------------------------------- popupmenuTargetType_Callback

function popupmenuTargetType_Callback(hObject, eventdata, handles)

% show or remove target frequency box

aShowRelevantObjects(handles)

% -------------------------------------------- popupmenuParadigm_Callback

function popupmenuParadigm_Callback(hObject, eventdata, handles)

% Any change to the paradigm selection causes all boxes to be shown

% showParameters(handles)

aParadigmSelection(handles);

% -------------------------------------------- aParadigmSelection

function aParadigmSelection(handles)

global experiment stimulusParameters betweenRuns  paradigmNames

global variableNames

% identify paradigm selected

chosenOption=get(handles.popupmenuParadigm,'value');

paradigm=paradigmNames{chosenOption};

experiment.paradigm=paradigm;

%Paradigm: read in all relevant parameters

% a file must exist with this name 'paradigm_<paradigm>'

% 'handles' are only occasionally used

addpath ('paradigms')

cmd=['paradigm_' paradigm '(handles);'];

try

    eval(cmd)

catch

    error(['ExpGUI\aParadigmSelection:'...

        'paradigm file not found or error in file'])

end

rmpath ('paradigms')

if ~isfield(experiment,'maskerInUse')

    error('selected paradigm does not specify if masker is used')

end

if ~experiment.maskerInUse

    stimulusParameters.maskerType='tone';

    stimulusParameters.maskerPhase='sin';

    stimulusParameters.maskerDuration=0.0;

    stimulusParameters.maskerLevel= -50;

    stimulusParameters.maskerRelativeFrequency= 1 ;

end

% if a variable is subject to change, specify list of values here

% eg. stimulusParameters.targetFrequency=betweenRuns.variableList1;

cmd=['stimulusParameters.' ...

    betweenRuns.variableName1 '=betweenRuns.variableList1;'];

eval (cmd);

cmd=['stimulusParameters.' ...

    betweenRuns.variableName2 '=betweenRuns.variableList2;'];

eval (cmd);

% establish popup menus on the basis of the paradigm file

% set(handles.popupmenuRandomize,'value', betweenRuns.randomizeSequence)

sequenceOptions=get(handles.popupmenuRandomize,'string');

idx=find(strcmp(betweenRuns.randomizeSequence, sequenceOptions)==1);

set(handles.popupmenuRandomize,'value', idx)

phaseOptions=get(handles.popupmenuPhase,'string');

idx=find(strcmp(stimulusParameters.maskerPhase, phaseOptions)==1);

set(handles.popupmenuPhase,'value', idx)

if stimulusParameters.includeCue

    set(handles.popupmenuCueNoCue,'value', 1)

else

    set(handles.popupmenuCueNoCue,'value', 2)

end

set(handles.text34, 'string', stimulusParameters.WRVname)

% Put the new data values into the edit boxes on the GUI

for i=1:length(variableNames)

    cmd=(['set(handles.edit' variableNames{i} ...

        ',''string'', num2str(stimulusParameters.' ...

        variableNames{i} '));']);

    eval(cmd);

end

% backgroundLevel is not a variableName (?!)

set(handles.editBackgroundLevel,'string', num2str...

    (stimulusParameters.backgroundLevel))

% values related to assessment method

switch experiment.threshEstMethod

    case {'MaxLikelihood','oneIntervalUpDown'}

        set(handles.editstopCriteriaBox, 'string', ...

            num2str(experiment.singleIntervalMaxTrials))

    case {'2I2AFC++','2I2AFC+++'}

        set(handles.editstopCriteriaBox, 'string', ...

            num2str(experiment.stopCriteria2IFC))

    otherwise

        error([' aResetPopupMenus:  threshEstMethod not recognised -> ' ...

            experiment.threshEstMethod])

end

% assessment method popup may be changed between paradigms

%  e.g. SRT must be one interval

x=get(handles.popupmenuThreshEst, 'string');

set(handles.popupmenuThreshEst, 'value', ...

    strmatch(experiment.threshEstMethod, x));

% on RUN the sample rate will be picked from the text box

% However, MAP overrules the sample rate and sets its own

aSetSampleRate(stimulusParameters.subjectSampleRate, handles);

% used for plotting functions (NB affected by paradigm settings)

experiment.predictionLevels=stimulusParameters.WRVlimits(1):...

    experiment.meanSearchStep:stimulusParameters.WRVlimits(2);

experiment.possLogSlopes=abs(experiment.possLogSlopes)*...

    sign(experiment.psyFunSlope);

aResetPopupMenus(handles)

% ------------------------------------------------------ aResetPopupMenus

function aResetPopupMenus(handles)

global   stimulusParameters betweenRuns variableNames

global targetTypes maskerTypes experiment backgroundTypes betweenRunsVariables

switch experiment.threshEstMethod

    case {'MaxLikelihood','oneIntervalUpDown'}

        set(handles.editstopCriteriaBox, 'string', ...

            num2str(experiment.singleIntervalMaxTrials))

    case {'2I2AFC++','2I2AFC+++'}

        set(handles.editstopCriteriaBox, 'string', ...

            num2str(experiment.stopCriteria2IFC))

    otherwise

        error([' aResetPopupMenus:  threshEstMethod not recognised -> ' ...

            experiment.threshEstMethod])

end

% forced noCue

switch experiment.paradigm

    case 'discomfort'

        set(handles.popupmenuCueNoCue,'value', 2)

end

%set variables popupmenus as specified in betweenRuns

variableParameter1ID=0; variableParameter2ID=0;

for i=1:length(betweenRunsVariables) %  variableNames

    if strcmp(betweenRunsVariables{i},betweenRuns.variableName1)

        variableParameter1ID=i;

    end

    if strcmp(betweenRunsVariables{i},betweenRuns.variableName2)

        variableParameter2ID=i;

    end

end

if variableParameter1ID==0 || variableParameter2ID==0;

    Error('a ResetPopMenu: variableParameter not identified')

end

% display popupmenus

set(handles.popupmenuVaryParameter1, 'value',round(variableParameter1ID))

set(handles.popupmenuVaryParameter2, 'value',round(variableParameter2ID))

% targetType

idx= find(strcmp(stimulusParameters.targetType, targetTypes));

set(handles.popupmenuTargetType,'value', idx)

% paradigm selection may alter the maskerType

idx= find(strcmp(stimulusParameters.maskerType, maskerTypes));

set(handles.popupmenuMaskerType,'value', idx)

aShowRelevantObjects(handles)

% % backgroundType popup

idx= find(strcmp(stimulusParameters.backgroundType, backgroundTypes));

set(handles.popupmenuBackgroundType,'value', idx)

set(handles.editBackgroundLevel,'string', ...

    num2str(stimulusParameters.backgroundLevel))

% ---------------------------------------------- aShowRelevantObjects

function aShowRelevantObjects(handles)

global experiment stimulusParameters

% called from aShowRelevantObjects

% always on

set(handles.edittargetLevel, 'visible', 'on')

set(handles.edittargetDuration, 'visible', 'on')

set(handles.edittargetFrequency, 'visible', 'on')

switch experiment.paradigm(1:3)

    case 'OHI'

        set(handles.editOHIOnTones, 'visible', 'on')

        set(handles.textOHIOnTones, 'visible', 'on')

    otherwise

        set(handles.editOHIOnTones, 'visible', 'off')

        set(handles.textOHIOnTones, 'visible', 'off')

end

switch experiment.ear

    case {'statsModelLogistic', 'statsModelRareEvent'}

        set(handles.editStatsModel, 'visible', 'on')

        set(handles.textStatsModel, 'visible', 'on')

        set(handles.pushbuttonStop, 'visible', 'on')

        showModelPushButtons(handles, 0)

        set(handles.editCatchTrialRate, 'visible', 'off')

        set(handles.textCatchTrials, 'visible', 'off')

        set(handles.editcalibrationdB, 'visible', 'off')

        set(handles.textcalibration, 'visible', 'off')

        set(handles.popupmenuCueNoCue, 'visible', 'off')

        set(handles.textCue, 'visible', 'off')

        set(handles.editMusicLevel,'visible', 'off')

        set(handles.textMusicLevel,'visible', 'off')

    case {'MAPmodel',  'MAPmodelListen', 'MAPmodelMultiCh', 'MAPmodelSingleCh'}

        set(handles.popupmenuCueNoCue, 'visible', 'off')

        set(handles.editStatsModel, 'visible', 'off')

        set(handles.textStatsModel, 'visible', 'off')

        set(handles.pushbuttonStop, 'visible', 'on')

        showModelPushButtons(handles, 1)

        set(handles.editcalibrationdB, 'visible', 'off')

        set(handles.textcalibration, 'visible', 'off')

        set(handles.textCue, 'visible', 'off')

        set(handles.editMusicLevel,'visible', 'off')

        set(handles.textMusicLevel,'visible', 'off')

    otherwise

        % i.e. using real subjects (left, right, diotic, dichotic)

        set(handles.editStatsModel, 'visible', 'off')

        set(handles.textStatsModel, 'visible', 'off')

        set(handles.pushbuttonStop, 'visible', 'off')

        showModelPushButtons(handles, 0)

        set(handles.editCatchTrialRate, 'visible', 'on')

        set(handles.textCatchTrials, 'visible', 'on')

        set(handles.editcalibrationdB, 'visible', 'on')

        set(handles.textcalibration, 'visible', 'on')

        set(handles.popupmenuCueNoCue, 'visible', 'on')

        set(handles.textCue, 'visible', 'on')

        set(handles.editMusicLevel,'visible', 'on')

        set(handles.textMusicLevel,'visible', 'on')

end

switch experiment.threshEstMethod

    case {'MaxLikelihood','oneIntervalUpDown'}

        set(handles.popupmenuCueNoCue,'visible', 'on')

        set(handles.textCue,'visible', 'on')

        set(handles.editstopCriteriaBox, 'string', ...

            num2str(experiment.singleIntervalMaxTrials))

        if stimulusParameters.includeCue==0

            set(handles.editcueTestDifference,'visible', 'off')

            set(handles.textcueTestDifference,'visible', 'off')

        else

            set(handles.editcueTestDifference,'visible', 'on')

            set(handles.textcueTestDifference,'visible', 'on')

        end

    case {'2I2AFC++','2I2AFC+++'}

        set(handles.editCatchTrialRate, 'visible', 'off')

        set(handles.textCatchTrials, 'visible', 'off')

        set(handles.popupmenuCueNoCue,'visible', 'off')

        set(handles.textCue,'visible', 'off')

        set(handles.editcueTestDifference,'visible', 'off')

        set(handles.textcueTestDifference,'visible', 'off')

end

% show/ remove masker related boxes

if ~experiment.maskerInUse

    set(handles.editmaskerDuration,'visible', 'off')

    set(handles.editmaskerLevel,'visible', 'off')

    set(handles.editmaskerRelativeFrequency,'visible', 'off')

    set(handles.editgapDuration,'visible', 'off')

    set(handles.textmaskerDuration,'visible', 'off')

    set(handles.textmaskerLevel,'visible', 'off')

    set(handles.textmaskerRelativeFrequency,'visible', 'off')

    set(handles.textgapDuration,'visible', 'off')

    set(handles.popupmenuMaskerType,'visible', 'off')

    set(handles.textMaskerType,'visible', 'off')

    % paradigms with maskers

else

    set(handles.editmaskerDuration,'visible', 'on')

    set(handles.editmaskerLevel,'visible', 'on')

    set(handles.editmaskerRelativeFrequency,'visible', 'on')

    set(handles.editgapDuration,'visible', 'on')

    set(handles.popupmenuMaskerType,'visible', 'on')

    set(handles.textmaskerDuration,'visible', 'on')

    set(handles.textmaskerLevel,'visible', 'on')

    set(handles.textmaskerRelativeFrequency,'visible', 'on')

    set(handles.textgapDuration,'visible', 'on')

    set(handles.popupmenuMaskerType,'visible', 'on')

    set(handles.textMaskerType,'visible', 'on')

    % masker relative frequency /type

    chosenOption=get(handles.popupmenuMaskerType,'value');

    maskerTypes=get(handles.popupmenuMaskerType,'string');

    maskerType=maskerTypes{chosenOption};

    switch maskerType

        case 'tone'

            set(handles.editmaskerRelativeFrequency,'visible', 'on')

        otherwise

            set(handles.editmaskerRelativeFrequency,'visible', 'off')

    end

end

eval(['set(handles.edit' stimulusParameters.WRVname ...

    ',''visible'',''off'' )'])

% target type

chosenOption=get(handles.popupmenuTargetType,'value');

targetTypes=get(handles.popupmenuTargetType,'string');

targetType=targetTypes{chosenOption};

switch targetType

    case 'tone'

        set(handles.edittargetFrequency,'visible', 'on')

    otherwise

        set(handles.edittargetFrequency,'visible', 'off')

end

if strcmp(stimulusParameters.backgroundType,'none')

    set(handles.popupmenuBackgroundType, 'visible', 'on');

    set(handles.editBackgroundLevel,'visible','off')

    set(handles.textBGlevel,'visible','off')

else

    set(handles.popupmenuBackgroundType, 'visible', 'on');

    set(handles.editBackgroundLevel,'visible','on')

    set(handles.textBGlevel,'visible','on')

end

% ------------------------------------------------ showModelPushButtons

function showModelPushButtons(handles,on)

if on

    set(handles.pushbuttonOME, 'visible', 'on')

    set(handles.pushbuttonBM, 'visible', 'on')

    set(handles.pushbuttonRP, 'visible', 'on')

    set(handles.pushbuttonAN, 'visible', 'on')

    set(handles.pushbuttonPhLk, 'visible', 'on')

    set(handles.pushbuttonSYN, 'visible', 'on')

    set(handles.pushbuttonFM, 'visible', 'on')

    set(handles.pushbuttonParams, 'visible', 'on')

    set(handles.pushbuttonSingleShot, 'visible', 'on')

    set(handles.editMAPplot,'visible', 'on')

    set(handles.textMAPplot,'visible', 'on')

else

    set(handles.pushbuttonOME, 'visible', 'off')

    set(handles.pushbuttonBM, 'visible', 'off')

    set(handles.pushbuttonRP, 'visible', 'off')

    set(handles.pushbuttonAN, 'visible', 'off')

    set(handles.pushbuttonPhLk, 'visible', 'off')

    set(handles.pushbuttonSYN, 'visible', 'off')

    set(handles.pushbuttonFM, 'visible', 'off')

    set(handles.pushbuttonParams, 'visible', 'off')

    set(handles.pushbuttonSingleShot, 'visible', 'off')

    set(handles.editMAPplot,'visible', 'off')

    set(handles.textMAPplot,'visible', 'off')

end

% ------------------------------------------------ pushbuttonRun_Callback

function pushbuttonRun_Callback(hObject, eventdata, handles)

global checkForPreviousGUI % holds screen positioning across repeated calls

global experiment stimulusParameters betweenRuns paradigmNames errormsg

checkForPreviousGUI.GUIposition=get(handles.figure1,'position');

experiment.singleShot=0;

experiment.stop=0;

switch experiment.paradigm

    case 'profile'

        %% special option for  successive and linked measurements

        global resultsTable

        experiment.paradigm='threshold_16ms';

        set(handles.edittargetDuration,'string', num2str(0.25))

        set(handles.editstopCriteriaBox,'string','10') % nTrials

        run (handles)

        if experiment.stop

            optionNo=strmatch('profile',paradigmNames);

            set(handles.popupmenuParadigm,'value',optionNo);

            experiment.paradigm='profile';

            aParadigmSelection(handles)

            return

        end

        longTone=resultsTable(2:end,2:end);

        set(handles.edittargetDuration,'string', num2str(0.016))

        set(handles.editstopCriteriaBox,'string','30') % nTrials

        run (handles)

        if experiment.stop

            disp(errormsg)

            optionNo=strmatch('profile',paradigmNames);

            set(handles.popupmenuParadigm,'value',optionNo);

            experiment.paradigm='profile';

            experiment.stop=-0;

            aParadigmSelection(handles)

            return

        end

        shortTone=resultsTable(2:end,2:end);

        % use 16ms thresholds for TMC

        thresholds16ms=betweenRuns.thresholds;

        optionNo=strmatch('TMC',paradigmNames);

        set(handles.popupmenuParadigm,'value',optionNo);

        aParadigmSelection(handles)

        set(handles.edittargetLevel,'string', thresholds16ms+10);

        set(handles.editstopCriteriaBox,'string','30')  % nTrials

        pause(.1)

        run (handles)

        if experiment.stop

            disp(errormsg)

            optionNo=strmatch('profile',paradigmNames);

            set(handles.popupmenuParadigm,'value',optionNo);

            experiment.paradigm='profile';

            experiment.stop=-0;

            aParadigmSelection(handles)

            return

        end

        TMC=resultsTable(2:end,2:end);

        gaps=resultsTable(2:end,1);

        BFs=resultsTable(1, 2:end);

        % use 16ms threshold for IFMC

        optionNo=strmatch('IFMC',paradigmNames);

        set(handles.popupmenuParadigm,'value',optionNo);

        aParadigmSelection(handles)

        set(handles.edittargetLevel,'string', thresholds16ms+10);

        set(handles.editstopCriteriaBox,'string','30')  % nTrials

        pause(.1)

        run (handles)

        if experiment.stop

            disp(errormsg)

            optionNo=strmatch('profile',paradigmNames);

            set(handles.popupmenuParadigm,'value',optionNo);

            experiment.paradigm='profile';

            experiment.stop=-0;

            aParadigmSelection(handles)

            return

        end

        IFMCs=resultsTable(2:end,2:end);

        offBFs=resultsTable(2:end,1);

        % reset original paradigm

        optionNo=strmatch('profile',paradigmNames);

        set(handles.popupmenuParadigm,'value',optionNo);

        aParadigmSelection(handles)

%% save data and plot profile

%         save profile longTone shortTone gaps BFs TMC offBFs IFMCs

fileName=['MTprofile' Util_timeStamp];

profile2mFile(longTone, shortTone, gaps, BFs, TMC, offBFs, IFMCs,...

    fileName)

while ~fopen(fileName), end

plotProfile(fileName, 'profile_CMA_L')

%% xx

        if strcmp(errormsg,'manually stopped')

            disp(errormsg)

            optionNo=strmatch('profile',paradigmNames);

            set(handles.popupmenuParadigm,'value',optionNo);

            experiment.paradigm='profile';

            experiment.stop=-0;

            aParadigmSelection(handles)

            return

        end

    otherwise

        run (handles)

        experiment.stop=0;

end

function run (handles)

global experiment expGUIhandles stimulusParameters

tic

expGUIhandles=handles;

set(handles.pushbuttonSave,'visible','off')

set(handles.pushbuttonStop, 'backgroundColor', [.941 .941 .941])

% set(handles.editparamChanges,'visible','off')

% message box white (removes any previous error message)

set(handles.textMSG,...

    'backgroundcolor', 'w', 'ForegroundColor', 'b', 'string', '')

errorMsg=aReadAndCheckParameterBoxes(handles);

if ~isempty(errorMsg)

    append=0;

    warning=1;

    addToMsg(['error message:' errorMsg], append, warning)

    return

end

if isnan(stimulusParameters.targetLevel)

    addToMsg('Error: targetlevel not set', 1)

    error('Error: targetlevel not set')

end

set(handles.textMSG,'backgroundcolor', 'w')

set(handles.textMSG,'string', ' ')

% leave the program and start up multiThreshold

subjGUI_MT % leave the program and start up multiThreshold

experiment.justInitialized=0;% prevents moving subjectGUI

toc

% --- Executes on button press in pushbuttonSingleShot.

function pushbuttonSingleShot_Callback(hObject, eventdata, handles)

global experiment paradigmNames

experiment.singleShot=1;

% startup paradigm, 'training' (could be anywhere on the list)

startupParadigm='training';

idx= find(strcmp(paradigmNames,startupParadigm));

set(handles.popupmenuParadigm,'value', idx)

aParadigmSelection(handles);

% special test for spontaneous activity

x=get(handles.editWRVstartValues, 'string');

y=get(handles.edittargetDuration, 'string');

set(handles.editWRVstartValues, 'string', '-20')

set(handles.edittargetDuration, 'string', '1')

MAPplot=get(handles.editMAPplot, 'string');

set(handles.editMAPplot, 'string', '1')

drawnow

run (handles)

set(handles.editMAPplot, 'string', MAPplot)

set(handles.editWRVstartValues, 'string', x)

set(handles.edittargetDuration, 'string', y)

% ------------------------------------------aReadAndCheckParameterBoxes

function errorMsg=aReadAndCheckParameterBoxes(handles)

global experiment  stimulusParameters betweenRuns  statsModel

global variableNames  LevittControl paramChanges betweenRunsVariables

% When the program sets the parameters all should be well

% But when the user changes them...

errorMsg='';

% name

experiment.name=get(handles.editName,'string');

% ear/models

option=get(handles.popupmenuEar,'value');

strings=get(handles.popupmenuEar,'string');

experiment.ear=strings{option};

switch experiment.ear

    case { 'MAPmodel', 'MAPmodelMultiCh', ...

            'MAPmodelSingleCh', 'MAPmodelListen'}

        % MAPmodel writes forced parameter settings to the screen

        %  so that they can be read from there

        % {'randomize within blocks', 'fixed sequence',...

        %  'randomize across blocks'}

        set(handles.popupmenuRandomize,'value',2)       % fixed sequence

        set(handles.editstimulusDelay,'string','0.01')  % no stimulus delay

        stimulusParameters.includeCue=0;                % no cue for MAP

end

% find tone type

option=get(handles.popupmenuTargetType,'value');

strings=get(handles.popupmenuTargetType,'string');

stimulusParameters.targetType=strings{option};

% find masker type

option=get(handles.popupmenuMaskerType,'value');

strings=get(handles.popupmenuMaskerType,'string');

stimulusParameters.maskerType=strings{option};

% find background type and level

option=get(handles.popupmenuBackgroundType,'value');

strings=get(handles.popupmenuBackgroundType,'string');

stimulusParameters.backgroundTypeValue=option;

stimulusParameters.backgroundLevel=...

    str2num(get(handles.editBackgroundLevel,'string'));

stimulusParameters.backgroundType=strings{option};

% Read all stimulus parameter boxes

for i=1:length(variableNames)

    cmd=['stimulusParameters.' variableNames{i} ...

        '= str2num(get(handles.edit' variableNames{i} ',''string'' ));'];

    eval(cmd);

end

% for multiple levels

stimulusParameters.targetLevels=stimulusParameters.targetLevel;

% check that all required values are in the edit boxes

for i=1:length(variableNames)-3 % do not include 'level limits'

    eval([ 'variableValues=stimulusParameters.' variableNames{i} ';'])

    if isempty(variableValues), errorMsg=[ variableNames{i} ...

            ' is an empty box']; return, end

end

chosenOption=get(handles.popupmenuVaryParameter1,'value');

betweenRuns.variableName1=betweenRunsVariables{chosenOption};

eval(['betweenRuns.variableList1 = stimulusParameters.' ...

    betweenRuns.variableName1 ';']);

chosenOption=get(handles.popupmenuVaryParameter2,'value');

betweenRuns.variableName2=variableNames{chosenOption};

eval(['betweenRuns.variableList2 = stimulusParameters.' ...

    betweenRuns.variableName2 ';']);

% Check that variable parameters 1 & 2 have different names

if strcmp(betweenRuns.variableName1,betweenRuns.variableName2) ...

        && ~strcmp(betweenRuns.variableName1,'none')

    errorMsg='variable parameters have the same name';

    return

end

% calibration

%  this is used to *reduce* the output signal from what it otherwise

%  would be

% signal values are between 1 - 2^23

%  these are interpreted as microPascals between -29 dB and 128 dB SPL

% calibrationdB adjusts these values to compensate for equipment

%  characteristics

%  this will change the range. e.g. a 7 dB calibration will yield

%   a range of -36 to 121 dB SPL

% Calibration is not used when modelling. Values are treated as dB SPL

stimulusParameters.calibrationdB=...

    str2num(get(handles.editcalibrationdB,'string'));

% check on cue

cueSetUp(handles)

stimulusParameters.WRVinitialStep=stimulusParameters.WRVsteps(1);

stimulusParameters.WRVsmallStep=stimulusParameters.WRVsteps(2);

% jitter start value set after reading in new parameters

switch experiment.paradigm

    case 'discomfort'

        stimulusParameters.jitterStartdB=0;

    otherwise

        stimulusParameters.jitterStartdB=abs(stimulusParameters.WRVsteps(1));

end

% stats model mean and slope

statsModelParameters= str2num(get(handles.editStatsModel,'string'));

switch experiment.ear

    case {'statsModelLogistic'}

        statsModel.logisticMean=statsModelParameters(1) ;

        statsModel.logisticSlope=statsModelParameters(2);

        statsModel.rareEvenGain=NaN ;

        statsModel.rareEventVmin=NaN;

    case 'statsModelRareEvent'

        statsModel.logisticMean=NaN ;

        statsModel.logisticSlope=NaN;

        statsModel.rareEvenGain=statsModelParameters(2) ;

        statsModel.rareEventVmin=statsModelParameters(1);

end

% MAP plotting

experiment.MAPplot=str2num(get(handles.editMAPplot,'string'));

% if ~isequal(experiment.MAPplot, 0),

%     % any other character will do it

%     experiment.MAPplot=1;

% end

% save data

experiment.saveData=str2num(get(handles.editSaveData,'string'));

if ~isequal(experiment.saveData, 0),

    % any other character will do it

    experiment.saveData=1;

end

% print tracks

experiment.printTracks=str2num(get(handles.editPrintTracks,'string'));

if ~isequal(experiment.printTracks, 0),

    % any other character will do it

    experiment.printTracks=1;

end

% catch trial rate

% (1)= start rate,   (2)= base rate,   (3)= time constant (trials)

stimulusParameters.catchTrialRates=...

    str2num(get(handles.editCatchTrialRate,'string'));

if stimulusParameters.catchTrialRates(1) ...

        < stimulusParameters.catchTrialRates(2)

    errorMsg=...

        'catch trial base rates must be less than catch trial start rate';

    return,

end

% force the decay rate for catchTrialRate

%  to avoid having to explain it to the user

stimulusParameters.catchTrialRates(3)=2;

% sample rate

% The sample rate is set in the paradigm file.

% Normally this is set in the startUp paradigm file and then left

% When the model is used, the multiThreshold sample rate

% overrules anything in the model

stimulusParameters.sampleRate=...

    str2num(get(handles.textsampleRate,'string'));

% music level

stimulusParameters.musicLeveldB=...

    str2num(get(handles.editMusicLevel,'string'));

%  set message box font size

experiment.msgFontSize=str2num(get(handles.editMsgFont,'string'));

experiment.subjGUIfontSize=str2num(get(handles.editSubjectFont,'string'));

% stop criteria

experiment.singleIntervalMaxTrials=...

    str2num(get(handles.editstopCriteriaBox,'string'));

experiment.maxTrials=experiment.singleIntervalMaxTrials(1);

% set up 2IFC is required (? better in atrhesholdAssessmentMethod)

switch experiment.threshEstMethod

    case {'2I2AFC++', '2A2AIFC+++'}

        experiment.peaksUsed=experiment.singleIntervalMaxTrials(2);

        experiment.PeakTroughCriterionSD=...

            experiment.singleIntervalMaxTrials(3);

        experiment.trialsToBeUsed= 5;

        LevittControl.maxTrials=experiment.singleIntervalMaxTrials(1);

        % start value for step until reduced

        LevittControl.startLevelStep= stimulusParameters.WRVsteps(1);

        % reduced step size

        LevittControl.steadyLevittStep= stimulusParameters.WRVsteps(2);

        LevittControl.TurnsToSmallSteps= 2;

        LevittControl.useLastNturns= 2*experiment.peaksUsed;

        LevittControl.minReversals= ...

            LevittControl.TurnsToSmallSteps+2*experiment.peaksUsed;

        LevittControl.targetsdPT = experiment.PeakTroughCriterionSD;

        LevittControl.maxLevittValue= stimulusParameters.WRVlimits(2);

        Levitt2;

end

% What kind of randomisation is required?

idx=get(handles.popupmenuRandomize,'value');

s=get(handles.popupmenuRandomize,'string');

betweenRuns.randomizeSequence=s{idx};

% Make small adjustments to variable values

%  to keep values unique against later sorting

x=betweenRuns.variableList1;

[y idx]= sort(x);

for i=1:length(y)-1, if y(i+1)==y(i); y(i+1)=y(i)*1.001; end, end

x(idx)=y;

betweenRuns.variableList1=x;

x=betweenRuns.variableList2;

[y idx]= sort(x);

for i=1:length(y)-1, if y(i+1)==y(i); y(i+1)=y(i)*1.001; end, end

x(idx)=y;

betweenRuns.variableList2=x;

% Checks: after reading and setting parameters  ------------------------------------------

% check for finger trouble (more checks possible

if stimulusParameters.maskerDuration>10

    errorMsg='maskerDuration is too long'; return, end

if stimulusParameters.gapDuration>10

    errorMsg='gapDuration is too long'; return, end

if stimulusParameters.targetDuration>10

    errorMsg='targetDuration is too long'; return, end

% rare event slope check

if experiment.psyFunSlope<0 ...

        && strcmp(experiment.functionEstMethod,'rareEvent')

    errorMsg='cannot use rareEvent option for negative psychometr slopes';

    return

end

% check ramps

if stimulusParameters.rampDuration*2> stimulusParameters.targetDuration

    errorMsg=' ramp duration is too long for the target';

    return

end

if max(stimulusParameters.maskerDuration)>0 ...

        && max(stimulusParameters.rampDuration...

        *2> stimulusParameters.maskerDuration)

    errorMsg=' ramp duration is too long for the masker';

    return

end

% Check WRVstartValues for length and compatibility with randomization

% only one start value supplied so all start values are the same

if length(stimulusParameters.WRVstartValues)==1

    stimulusParameters.WRVstartValues= ...

        repmat(stimulusParameters.WRVstartValues, 1, ...

        length(betweenRuns.variableList1)...

        *length(betweenRuns.variableList2));

elseif ~isequal(length(stimulusParameters.WRVstartValues), ...

        length(betweenRuns.variableList1)...

        *length(betweenRuns.variableList2))

    % otherwise we need one value for each combination of var1/ var2

    errorMsg='WRVstartValues not the same length as number of runs';

    return

elseif strcmp(betweenRuns.randomizeSequence, 'randomize within blocks')...

        && length(stimulusParameters.WRVstartValues)>1

    errorMsg='multiple WRVstartValues inappropriate';

    return

end

% set the target level in advance for every wnticipated run

switch experiment.paradigm

    case {'trainingIFMC', 'TMC','TMC_16ms', 'TMC - ELP', ...

            'IFMC', 'IFMC_16ms'}

        % For TMC and IFMC multiple target levels can be set

        if length(stimulusParameters.targetLevel)==1

            betweenRuns.targetLevels=...

                repmat(stimulusParameters.targetLevel, 1, ...

                length(betweenRuns.variableList1)...

                *length(betweenRuns.variableList2));

        elseif length(stimulusParameters.targetLevel)==...

                length(betweenRuns.variableList2)

            % only one level specified, so use this throughout

            x=stimulusParameters.targetLevel;

            x=repmat(x,length(betweenRuns.variableList1),1);

            x=reshape(x,1,length(betweenRuns.variableList1)*length(betweenRuns.variableList2));

            betweenRuns.targetLevels=x;

        else

            errorMsg='targetLevels not the same length as number of targetFrequencies';

        end

end

switch experiment.paradigm

    case  {'gapDetection', 'frequencyDiscrimination'}

        if ~isequal(stimulusParameters.targetDuration, ...

                stimulusParameters.maskerDuration)

            addToMsg(...

                'Warning: masker and target duration not the same.',1,1)

        end

        if ~isequal(stimulusParameters.maskerLevel, ...

                stimulusParameters.targetLevel)

            addToMsg(['Warning: masker and target level different.'...

                'They will be changed to be equal',1,1]);

        end

end

% identify model parameter changes if any

paramChanges=get(handles.editparamChanges,'string');

if ~strcmp(paramChanges, ';'), paramChanges=[paramChanges ';']; end

try

eval(paramChanges);

catch

    error('Problems with suggested parameter changes')

end

% -------------------------------------------- aSetSampleRate

function aSetSampleRate(sampleRate, handles)

global  stimulusParameters

stimulusParameters.sampleRate=sampleRate;

set(handles.textsampleRate,'string',num2str(stimulusParameters.sampleRate))

% -------------------------------------------- popupmenuEar_Callback

function popupmenuEar_Callback(hObject, eventdata, handles)

global experiment

option=get(handles.popupmenuEar,'value');

options=get(handles.popupmenuEar,'string');			% left/right/model

experiment.ear=options{option};

set(handles.editparamChanges,'visible','off')

switch experiment.ear

    case 'statsModelLogistic'

        set(handles.editStatsModel,'string', '15, 0.5')

    case 'statsModelRareEvent'

        set(handles.editStatsModel,'string', '20 1')

    case {'MAPmodelListen',  'MAPmodelMultiCh', 'MAPmodelSingleCh'}

        set(handles.editparamChanges,'visible','on')

end

earSetUp(handles)

% -------------------------------------------- earSetUp

function	earSetUp(handles)

global experiment stimulusParameters

% option=get(handles.popupmenuEar,'value');

% options=get(handles.popupmenuEar,'string');			% left/right/model

% experiment.ear=options{option};

switch experiment.ear

    case {'statsModelLogistic'}