MAP

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

% Begin initialization code - DO NOT EDIT

gui_Singleton = 1;

% Then wait for user action

global stimulusParameters experiment betweenRuns

global targetTypes maskerTypes backgroundTypes

% Specify order of fields in main structures

% identify as empty values or empty strings only

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

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

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

% Variable variables

%  (names of variable that can changed between runs)

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

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

% populate the 'between runs variable' menus

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

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

% Establish available paradigms by scanning paradigms folder

paradigmNames= what('paradigms');

paradigmNames= paradigmNames.m; % select m files only

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

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

end

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

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

experiment.stopCriteria2IFC=[75 3 5];

% ** editBoxes that are only set by hand

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

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

% Catch Trial Rate

% calibration

stimulusParameters.restoreCalibration=7;

% saveData

experiment.saveData=1;

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

% printTracks

experiment.printTracks=0;

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

    (stimulusParameters.backgroundLevel))

% 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);

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

function aResetPopupMenus(handles)

global   stimulusParameters betweenRuns variableNames

switch experiment.threshEstMethod

    case {'MaxLikelihood','oneIntervalUpDown'}

%set variables popupmenus as specified in betweenRuns

variableParameter1ID=0; variableParameter2ID=0;

        variableParameter1ID=i;

    end

        variableParameter2ID=i;

    end

end

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

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

switch experiment.ear

    case {'statsModelLogistic', 'statsModelRareEvent'}

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

global experiment expGUIhandles stimulusParameters

tic

expGUIhandles=handles;

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

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

% --- Executes on button press in pushbuttonSingleShot.

function pushbuttonSingleShot_Callback(hObject, eventdata, handles)

experiment.singleShot=1;

% special test for spontaneous activity

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

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

% ------------------------------------------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...

end

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

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

    betweenRuns.variableName1 ';']);

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

    return,

end

% sample rate

% The sample rate is set in the paradigm file.

% identify model parameter changes if any

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

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

eval(paramChanges);

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

function aSetSampleRate(sampleRate, handles)

        stimulusParameters.includeCue=0;						 % no cue

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

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

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

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

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

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

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

end

aResetPopupMenus(handles)

% ------------------------------------------------------------- cueSetUp

function cueSetUp(handles)

end

% -------------------------------------------- popupmenuThreshEst_Callback

end

% Handles everything concerned with the stimulus presentation

%  called from startNewRun in subjGUI

experiment.status='presentingStimulus';

errormsg='';

            case {'training','absThreshold', 'absThreshold_8',  ...

                    'TENtest', 'threshold_duration','discomfort',...

                    'overShoot','overShootB','overShootMB1', ...

                cueTargetLevel=targetLevel+cueTestDifference;

            case {'forwardMasking','forwardMaskingD','trainingIFMC', ...

                % cue masker is weaker to make target more audible

                cueMaskerLevel=maskerLevel-cueTestDifference;

        end

end

% ----------------------------- catch trial

if withinRuns.catchTrial

    targetLevel=-100;	% no target

            return

        end

        % legacy programming

        if gapDuration<0

            errormsg=['gapDuration (' num2str(gapDuration) ...

                ') is less than zero  ***'];

globalStimParams.nSignalPoints=...

    round(globalStimParams.overallDuration*globalStimParams.FS);

%           ----------------------------------------------cue stimulus

% cue masker

componentNo=1;

componentNo=2;

precedingSilence=precedingSilence + maskerDuration+cueGapDuration;

stimComponents(targetEar,componentNo).type=targetType;

stimComponents(targetEar,componentNo).toneDuration=targetDuration;

stimComponents(targetEar,componentNo).frequencies=cueTargetFrequency;

stimComponents(targetEar,componentNo).amplitudesdB=cueTargetLevel;

componentNo=2;

targetDelay=precedingSilence+ maskerDuration+ gapDuration;

stimComponents(targetEar,componentNo).type=targetType;

stimComponents(targetEar,componentNo).toneDuration=targetDuration;

stimComponents(targetEar,componentNo).frequencies=targetFrequency;

stimComponents(targetEar,componentNo).amplitudesdB=targetLevel;

        %     case 'SRT'

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

    otherwise

end

% ------------------------------------------------------------- play!

AN_spikesOrProbability='spikes';

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

    MAP1_14(audio, 1/method.dt, method.nonlinCF,...

        MAPparamsName, AN_spikesOrProbability);

    options.showModelOutput=1;

    options.printFiringRates=1;

    options.showACF=0;

    UTIL_showMAP(options)

end

% ---------------------------------------------------------- end model run

time=dt:dt:dt*length(MacGregorResponse);

% group delay on unit response

function paradigm_IFMC(handles)

global stimulusParameters experiment betweenRuns

stimulusParameters.WRVname='maskerLevel';

stimulusParameters.WRVstartValues=10;

function paradigm_TENtest(handles)

global stimulusParameters experiment betweenRuns

stimulusParameters.WRVname='targetLevel';

stimulusParameters.WRVstartValues=40;

%   single interval up/down no cue

stimulusParameters.instructions{1}=[{'YES if you hear the added click'}, { }, { 'NO if not (or you are uncertain'}];

%   single interval up/down with cue

    {'The clicks must be **clearly distinct** to count'}];

function paradigm_TMC(handles)

global stimulusParameters experiment betweenRuns

stimulusParameters.WRVname='maskerLevel';

stimulusParameters.WRVstartValues=50;

function paradigm_absThreshold(handles)

global stimulusParameters experiment betweenRuns

betweenRuns.variableName1='targetFrequency';

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

function paradigm_absThreshold16(handles)

global stimulusParameters experiment betweenRuns

betweenRuns.variableName1='targetFrequency';

betweenRuns.variableList1=1000;

function paradigm_discomfort(handles)

global stimulusParameters experiment betweenRuns

stimulusParameters.WRVname='targetLevel';

stimulusParameters.WRVstartValues=75 ;

function paradigm_forwardMasking(handles)

global stimulusParameters experiment betweenRuns

stimulusParameters.WRVname='targetLevel';

stimulusParameters.WRVstartValues=50;

function paradigm_overShoot(handles)

global stimulusParameters  betweenRuns experiment

stimulusParameters.WRVname='targetLevel';

stimulusParameters.WRVstartValues=50;

function paradigm_profile(handles)

global stimulusParameters experiment betweenRuns

betweenRuns.variableName1='targetFrequency';

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

function paradigm_psychometric(handles)

global stimulusParameters experiment betweenRuns

experiment.printTracks=1;

experiment.maxTrials=30;

function paradigm_threshold_duration(handles)

global stimulusParameters experiment betweenRuns

stimulusParameters.WRVname='targetLevel';

stimulusParameters.WRVstartValues=40;

function paradigm_training(handles)

global stimulusParameters experiment betweenRuns

% stimulusParameters.subjectSampleRate=128000; % compatible with file input

%  assessment method

global stimulusParameters experiment betweenRuns

stimulusParameters.WRVname='maskerLevel';

stimulusParameters.WRVstartValues=-10;

betweenRuns.variableList1=[ 0.5       .9 .7   1.3 ];

betweenRuns.variableName2='targetFrequency';

% keep old list of target frequencies

experiment.maskerInUse=1;

stimulusParameters.maskerType='tone';

%% plot profile

if nargin<1

    bgName = '';

end

    freq=foreground.MaskerRatio'*foreground.IFMCFreq(BFno);

    subplot(2,1,2)

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

    xlim([100 12000])

%     grid on

end

xlabel('frequency (Hz)')

ylabel('masker dB / probe dB')

set(gca,'XTick',foreground.IFMCFreq)

if ~isempty(bgName)

    for BFno=1:length(background.IFMCFreq)

printReportGuide.showPsychometric=0;

printReportGuide.HorizontalTracks=1;

if nargin==0

    % print new report

    printReportGuide.showTracks=experiment.printTracks;

    disp('Psychometric function')

    fprintf(' level  \tfreq\tprob\n')

    fprintf('%6.0f\t%6.2f\t%6.0f\n',[levelsBinVector; binFrequencies; psy])

    fprintf('\n')

end

        otherwise

            title([stimulusParameters.WRVname ' = ' ...

                num2str(withinRuns.variableValue, '%5.1f') ...

    end

end

xlim([0 experiment.maxTrials+withinRuns.beginningOfPhase2]);

ylim(stimulusParameters.WRVlimits)

grid on

% Panel 4: Summary of threshold estimates (not used here)

end

xlimRM([ min(betweenRuns.variableList1) max(betweenRuns.variableList1) ])

ylim(stimulusParameters.WRVlimits)

xlabel(betweenRuns.variableName1)

set(gca,'ytick', [0 20 40 60 80 100])

try

    % update experimenter GUI

    addToMsg('Experiment completed.',1)

    printReport

    experiment.status='endOfExperiment';

    return

global experiment stimulusParameters audio withinRuns

% global outerMiddleEarParams DRNLParams AN_IHCsynapseParams

global paramChanges

savePath=path;

    options.showModelOutput=1;

    options.printFiringRates=1;

    options.showACF=0;

    options.surfProbability=0;

    showMapOptions.surfSpikes=0;

    UTIL_showMAP(options)

if strcmp(AN_spikesOrProbability,'spikes')

    MacGregorResponse= sum(ICoutput,1);                 % use IC

    time=dt:dt:dt*length(MacGregorResponse);

else

    % for one channel, ANprobResponse=ANprobRateOutput

windowOnsetDelay= 0.004;

windowOffsetDelay= 0.022; % long ringing time

switch experiment.threshEstMethod

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

        idx= time>stimulusParameters.testTargetBegins+windowOnsetDelay ...

            num2str([withinRuns.variableValue nSpikesTrueWindow ...

            nSpikesFalseWindow  nSpikesDuringTarget], '%4.0f') ] )

    otherwise

        idx=find(time>stimulusParameters.testTargetBegins +windowOnsetDelay...

            & time<stimulusParameters.testTargetEnds+windowOffsetDelay);

        if strcmp(AN_spikesOrProbability,'spikes')

end

try

    serobj.Baudrate = 9600;           		% Set the baud rate at the specific value

    set(serobj, 'Parity', 'none') ;     	% Set parity as none

    set(serobj, 'Databits', 8) ;          	% set the number of data bits