MAP

% expGUI_MT = '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':

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

% 2. Measurements using the MAP model as the subject

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

setLocationOfGUIs(handles)

function setLocationOfGUIs(handles)

scrnsize=get(0,'screensize');

checkForPreviousGUI=[];

% if isstruct(checkForPreviousGUI)...

% addpath (['..' filesep 'modules'], ['..' filesep 'utilities'], ...

%     ['..' filesep 'parameterStore'],  ['..' filesep 'wavFileStore'],...

%     ['..' filesep 'testPrograms'])

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

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

paradigm=paradigmNames{chosenOption};

experiment.paradigm=paradigm;

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

% 'handles' are only occasionally used

addpath ('paradigms')

eval (cmd);

% establish popup menus on the basis of the paradigm file

set(handles.popupmenuPhase,'string', stimulusParameters.maskerPhase)

if stimulusParameters.includeCue

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

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

function aResetPopupMenus(handles)

global targetTypes maskerTypes experiment backgroundTypes

switch experiment.threshEstMethod

    case {'MaxLikelihood','oneIntervalUpDown'}

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

            num2str(experiment.singleIntervalMaxTrials))

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

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

            num2str(experiment.stopCriteria2IFC))

for i=1:length(variableNames)

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

        variableParameter1ID=i;

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

        variableParameter2ID=i;

    end

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

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

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

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

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

        set(handles.editcalibrationdB, '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')

        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')

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

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

        set(handles.editcalibrationdB, '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')

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

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

        end

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

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

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

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

function pushbuttonRun_Callback(hObject, eventdata, handles)

global experiment betweenRuns paradigmNames errormsg

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

experiment.singleShot=0;

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

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

        run (handles)

            disp(errormsg)

            optionNo=strmatch('profile',paradigmNames);

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

            aParadigmSelection(handles)

            return

        end

        global resultsTable

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

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

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

        run (handles)

            disp(errormsg)

            optionNo=strmatch('profile',paradigmNames);

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

            aParadigmSelection(handles)

            return

        end

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

        % use these threshold for TMC

        aParadigmSelection(handles)

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

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

        run (handles)

            disp(errormsg)

            optionNo=strmatch('profile',paradigmNames);

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

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

        gaps=resultsTable(2:end,1);

        BFs=resultsTable(1, 2:end);

        % use these threshold for IFMC

        optionNo=strmatch('IFMC',paradigmNames);

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

        aParadigmSelection(handles)

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

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

        run (handles)

            disp(errormsg)

            optionNo=strmatch('profile',paradigmNames);

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

            experiment.stop=-0;

            aParadigmSelection(handles)

            return

    otherwise

        run (handles)

        experiment.stop=0;

tic

expGUIhandles=handles;

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

% message box white (removes any previous error message)

set(handles.textMSG,...

% --- Executes on button press in pushbuttonSingleShot.

function pushbuttonSingleShot_Callback(hObject, eventdata, handles)

experiment.singleShot=1;

% special test for spontaneous activity

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

function errorMsg=aReadAndCheckParameterBoxes(handles)

global experiment  stimulusParameters betweenRuns  statsModel

% When the program sets the parameters all should be well

% But when the user changes them...

switch experiment.ear

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

        %  so that they can be read from there

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

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

end

% calibration

%  would be

% signal values are between 1 - 2^23

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

%  characteristics

%   a range of -36 to 121 dB SPL

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

stimulusParameters.calibrationdB=...

        % start value for step until reduced

        LevittControl.startLevelStep= stimulusParameters.WRVsteps(1);

        % reduced step size

        LevittControl.TurnsToSmallSteps= 2;

        LevittControl.useLastNturns= 2*experiment.peaksUsed;

        LevittControl.minReversals= ...

% Check WRVstartValues for length and compatibility with randomization

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

    stimulusParameters.WRVstartValues= ...

        repmat(stimulusParameters.WRVstartValues, 1, ...

        length(betweenRuns.variableList1)...

                stimulusParameters.maskerDuration)

            addToMsg(...

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

        if ~isequal(stimulusParameters.maskerLevel, ...

                stimulusParameters.targetLevel)

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

        end

end

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

function aSetSampleRate(sampleRate, handles)

global  stimulusParameters

    case {'statsModelLogistic', 'statsModelRareEvent'}

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

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

        % default psychometric bin size and logistic slopes

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

        set(handles.editName,'string', 'statsModel')

        %         experiment.headphonesUsed=0;

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

        stimulusParameters.includeCue=0;						 % no cue

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

        set(handles.editCatchTrialRate,'string','0 0  2 ');%no catch trials

        set(handles.editName,'string', 'Normal')			% force name

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

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

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

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

        set(handles.editSubjectFont,'string', '10');

        experiment.MacGThreshold=0; % num MacG spikes to exceed threshold

            case 'training'

                experiment.possLogSlopes=0.5;

        end

    case 'oneIntervalUpDown'

        experiment.functionEstMethod='logisticLS';

        set(handles.textstopCriteria,'string', 'stop criteria \ maxTrials')

            otherwise

                experiment.allowCatchTrials= 1;

        end

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

        LevittControl.rule='++'; %  e.g. '++' or '+++'

        experiment.singleIntervalMaxTrials=experiment.stopCriteria2IFC;

% NB responsibility for this is now transferred to the paradigm file

switch experiment.threshEstMethod

    % only one value required for level change

        stimulusParameters.subjectText=...

            'did the tone occur in window 1 or 2?';

    case {'MaxLikelihood',  'oneIntervalUpDown'};

drawnow

function pushbuttonOME_Callback(hObject, eventdata, handles)

aReadAndCheckParameterBoxes(handles);

function pushbuttonBM_Callback(hObject, eventdata, handles)

aReadAndCheckParameterBoxes(handles);

relativeFrequencies=[0.25    .5   .75  1  1.25 1.5    2];

testBM(stimulusParameters.targetFrequency, ...

function pushbuttonAN_Callback(hObject, eventdata, handles)

aReadAndCheckParameterBoxes(handles);

% now carry out tests

showPSTHs=0;

targetFrequency=stimulusParameters.targetFrequency(1);

BFlist=targetFrequency;

function pushbuttonPhLk_Callback(hObject, eventdata, handles)

aReadAndCheckParameterBoxes(handles);

function pushbuttonSYN_Callback(hObject, eventdata, handles)

aReadAndCheckParameterBoxes(handles);

function pushbuttonFM_Callback(hObject, eventdata, handles)

aReadAndCheckParameterBoxes(handles);

showPSTHs=1;

function popupmenuPhase_Callback(hObject, eventdata, handles)

global stimulusParameters

function pushbuttonParams_Callback(hObject, eventdata, handles)

global experiment stimulusParameters

aReadAndCheckParameterBoxes(handles);

showParams=1; BFlist=-1;

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

paramFunctionName=['method=MAPparams' experiment.name ...

    '(BFlist, stimulusParameters.sampleRate, showParams,paramChanges);'];

% --- Executes on button press in pushbuttonRP.

function pushbuttonRP_Callback(hObject, eventdata, handles)

aReadAndCheckParameterBoxes(handles);

% now carry out test

% function handles % ??

end

function editparamChanges_Callback(hObject, eventdata, handles)

errormsg='';

% interrupt by 'stop' button

% -----------------------------------------choose catch trials at random

% catch trials are for subject threshold measurements only

global outerMiddleEarParams DRNLParams AN_IHCsynapseParams

savePath=path;

modelResponse=[];

MacGregorResponse=[];

        MAPparamsName, AN_spikesOrProbability);

if showPlotsAndDetails

    options.showModelOutput=1;

    options.printFiringRates=1;

    options.showACF=0;

    options.showEfferent=1;

end

% No response,  probably caused by hitting 'stop' button

paradigm_training(handles) % default

stimulusParameters.WRVname='maskerLevel';

stimulusParameters.WRVsteps= [-10 -2];

stimulusParameters.WRVlimits=[-30 110];

paradigm_training(handles) % default

stimulusParameters.WRVname='maskerLevel';

stimulusParameters.WRVsteps= [-10 -4];

stimulusParameters.WRVlimits=[-30 110];

betweenRuns.variableName1='targetFrequency';

betweenRuns.variableList1=[250 500 1000 2000 4000 8000 ];

betweenRuns.variableName2='targetDuration';

betweenRuns.variableList2=0.016;

betweenRuns.randomizeSequence=1; % 'random sequence'

stimulusParameters.targetType='tone';

stimulusParameters.targetPhase='sin';

stimulusParameters.targetDuration=0.016;

stimulusParameters.targetLevel=stimulusParameters.WRVstartValues(1);

    load profile

end

];

% absolute thresholds

figure(90), clf

subplot(2,1,2)

semilogx(BFs,longTone,'ko-','lineWidth',2); hold on

semilogx(BFs,shortTone,'bo-','lineWidth',2); hold on

ylim([0 100])

% TMC

for BFno=1:length(BFs)

    subplot(2,6,BFno)

    plot(gaps,TMC(:,BFno)-longTone(BFno),'r','lineWidth',3), hold on

    xlim([0.01 0.1])

    if BFno==1

        ylabel('masker dB SL')

        xlabel('gap')

    end

end

% IFMCs

for BFno=1:length(BFs)

    freq=offBFs'*BFs(BFno);

    subplot(2,1,2)

    semilogx(freq,IFMCs(:,BFno),'r','lineWidth',3), hold on

    ylim([0 100])

    xlim([100 12000])

    grid on

xlabel('frequency (Hz)')

ylabel('masker dB / probe dB')

set(gca,'XTick',BFs)

global experiment stimulusParameters betweenRuns withinRuns statsModel   audio

global LevittControl expGUIhandles

global inputStimulusParams OMEParams DRNLParams

global IHC_VResp_VivoParams IHCpreSynapseParams  AN_IHCsynapseParams

    printReportGuide.showTracks=experiment.printTracks;

    printReportGuide.fileName=[];

    if experiment.saveData

    else

        % save this data (just in case)

        saveFileName=['savedData/mostRecentResults'];

msg=printTabTable(resultsTable,  headers);

addToMsg(msg,0)

fprintf('\n')

% sort tracks into the same order

betweenRuns.levelTracks=betweenRuns.levelTracks(idx1);

end

fprintf('\nparadigm:\t%s\n ', experiment.paradigm)

% ------------------------------------------------------- sortTablesForPrinting

function table= sortTablesForPrinting(idx1,idx2, var1values,var2values, x)

    %  without waiting for button press

    startNewRun(handles)

    set(expGUIhandles.textsampleRate,'string',...

        num2str(stimulusParameters.sampleRate))

        AN_IHCsynapseParams.plotSynapseContents=1;

    else

        AN_IHCsynapseParams.plotSynapseContents=0;

function [modelResponse, MacGregorResponse]=MAPmodel

global experiment stimulusParameters audio withinRuns

savePath=path;

addpath(['..' filesep 'MAP'], ['..' filesep 'utilities'])

% ---------------------------------------------- run Model

MAPparamsName=experiment.name;

showPlotsAndDetails=experiment.MAPplot;

AN_spikesOrProbability='spikes';

AN_spikesOrProbability='probability';

% [response, method]=MAPsequenceSeg(audio, method, 1:8);

if sum(strcmp(experiment.ear,{'MAPmodelMultiCh', 'MAPmodelListen'}))

    % use BFlist specified in MAPparams file

    BFlist=stimulusParameters.targetFrequency;

end

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

MAP1_14(audio, stimulusParameters.sampleRate, BFlist,...

    MAPparamsName, AN_spikesOrProbability, paramChanges);