rmeddis@38: function testLibermanMOC_DPOAE rmeddis@38: rmeddis@38: % compares MOC response to LIberman's 1996 data for DPOAE reduction with rmeddis@38: % contralateral tone stimulation. rmeddis@38: % This program is used mainly as a check on the time constants involved. rmeddis@38: % rmeddis@38: % NB very different time constants are required for 'spikes' and rmeddis@38: % 'probability' rmeddis@38: rmeddis@38: % global dt ANdt savedBFlist saveAN_spikesOrProbability saveMAPparamsName... rmeddis@38: % savedInputSignal OMEextEarPressure TMoutput OMEoutput ARattenuation ... rmeddis@38: % DRNLoutput IHC_cilia_output IHCrestingCiliaCond IHCrestingV... rmeddis@38: % IHCoutput ANprobRateOutput ANoutput savePavailable ANtauCas ... rmeddis@38: % CNtauGk CNoutput ICoutput ICmembraneOutput ICfiberTypeRates ... rmeddis@38: % MOCattenuation rmeddis@38: global dt dtSpikes saveAN_spikesOrProbability ANprobRateOutput ICoutput rmeddis@38: rmeddis@38: global DRNLParams rmeddis@38: rmeddis@38: LibermanData=[ rmeddis@38: 2 0.2; rmeddis@38: 2.1 0.19;2.2 0.18;2.3 0.18;2.4 0.16;2.5 0.15;2.6 0.15;2.7 0.15; rmeddis@38: 2.8 0.12;2.9 0.12;3 0.1;3.1 0.1;3.2 0.05;3.3 0.05;3.4 0;3.5 -0.1; rmeddis@38: 3.6 -0.4;3.7 -1.2;3.8 -1.6;3.9 -1.8;4 -1.85;4.1 -2;4.2 -2.05; rmeddis@38: 4.3 -2.05;4.4 -2.15;4.5 -2.2;4.6 -2.15;4.7 -2.1;4.8 -2.15;4.9 -2.2; rmeddis@38: 5 -2.1;5.1 -2.1;5.2 -2.25;5.3 -2.1;5.4 -2.15;5.5 -2.1;5.6 -2.15; rmeddis@38: 5.7 -2.1;5.8 -2.2;5.9 -2.05;6 -2.15;6.1 -2.05;6.2 -2;6.3 -2.2;6.4 -2.1; rmeddis@38: 6.5 -2.05;6.6 -2.05;6.7 -2.05;6.8 -2.2;6.9 -2.1;7 -2.05;7.1 -2.05;7.2 -0.7; rmeddis@38: 7.3 -0.1;7.4 0;7.5 0.1;7.6 0.2;7.7 0.35;7.8 0.2;7.9 0.15;8 0.2;8.1 0.15;8.2 0.15; rmeddis@38: 8.3 0.15;8.4 0.12;8.5 0.1;8.6 0.09;8.7 0.08;8.8 0.07;8.9 0.06;9 0.05; rmeddis@38: ]; rmeddis@38: rmeddis@38: % Backus2006Data: time bilateral contralateral ipsilateral rmeddis@38: % all % max microPascals rmeddis@38: Backus2006Data=[ rmeddis@38: 100 20 15 10; rmeddis@38: 200 36 24 14; rmeddis@38: 300 44 30 18; rmeddis@38: 400 46 32 20; rmeddis@38: 500 48 34 22; rmeddis@38: 1000 50 36 24; rmeddis@38: 1500 52 37 25; rmeddis@38: 2000 54 38 27 rmeddis@38: ]; rmeddis@38: rmeddis@38: steadyMinimum=mean(LibermanData(LibermanData(:,1)>4 & LibermanData(:,1)<7,2)); rmeddis@38: rmeddis@38: restorePath=path; rmeddis@38: addpath (['..' filesep 'MAP'], ['..' filesep 'wavFileStore'], ... rmeddis@38: ['..' filesep 'utilities']) rmeddis@38: rmeddis@38: %% #1 parameter file name rmeddis@38: MAPparamsName='Normal'; rmeddis@38: rmeddis@38: rmeddis@38: %% #2 probability (fast) or spikes (slow) representation rmeddis@38: AN_spikesOrProbability='spikes'; rmeddis@38: % or rmeddis@38: AN_spikesOrProbability='probability'; rmeddis@38: rmeddis@38: rmeddis@38: %% #3 pure tone, harmonic sequence or speech file input rmeddis@38: signalType= 'tones'; rmeddis@38: sampleRate= 50000; rmeddis@38: rampDuration=.005; % raised cosine ramp (seconds) rmeddis@38: toneFrequency= 1000; % or a pure tone (Hz) rmeddis@38: duration=3.6; % Lieberman test rmeddis@38: beginSilence=1; % 1 for Lieberman rmeddis@38: endSilence=1; % 1 for Lieberman rmeddis@38: rmeddis@38: %% #4 rms level rmeddis@38: % signal details rmeddis@38: leveldBSPL= 80; % dB SPL (80 for Lieberman) rmeddis@38: rmeddis@38: rmeddis@38: %% #5 number of channels in the model rmeddis@38: rmeddis@38: numChannels=1; rmeddis@38: BFlist=toneFrequency; rmeddis@38: rmeddis@38: rmeddis@38: %% #6 change model parameters rmeddis@38: paramChanges={}; rmeddis@38: rmeddis@38: %% delare 'showMap' options to control graphical output rmeddis@38: showMapOptions.printModelParameters=1; % prints all parameters rmeddis@38: showMapOptions.showModelOutput=1; % plot of all stages rmeddis@38: showMapOptions.printFiringRates=1; % prints stage activity levels rmeddis@38: showMapOptions.showACF=0; % shows SACF (probability only) rmeddis@38: showMapOptions.showEfferent=1; % tracks of AR and MOC rmeddis@38: showMapOptions.surfAN=1; % 2D plot of HSR response rmeddis@38: showMapOptions.surfSpikes=0; % 2D plot of spikes histogram rmeddis@38: showMapOptions.ICrates=0; % IC rates by CNtauGk rmeddis@38: rmeddis@38: %% Generate stimuli rmeddis@38: rmeddis@38: % Create pure tone stimulus rmeddis@38: dt=1/sampleRate; % seconds rmeddis@38: time=dt: dt: duration; rmeddis@38: inputSignal=sum(sin(2*pi*toneFrequency'*time), 1); rmeddis@38: amp=10^(leveldBSPL/20)*28e-6; % converts to Pascals (peak) rmeddis@38: inputSignal=amp*inputSignal; rmeddis@38: % apply ramps rmeddis@38: % catch rampTime error rmeddis@38: if rampDuration>0.5*duration, rampDuration=duration/2; end rmeddis@38: rampTime=dt:dt:rampDuration; rmeddis@38: ramp=[0.5*(1+cos(2*pi*rampTime/(2*rampDuration)+pi)) ... rmeddis@38: ones(1,length(time)-length(rampTime))]; rmeddis@38: inputSignal=inputSignal.*ramp; rmeddis@38: ramp=fliplr(ramp); rmeddis@38: inputSignal=inputSignal.*ramp; rmeddis@38: % add silence rmeddis@38: intialSilence= zeros(1,round(beginSilence/dt)); rmeddis@38: finalSilence= zeros(1,round(endSilence/dt)); rmeddis@38: inputSignal= [intialSilence inputSignal finalSilence]; rmeddis@38: rmeddis@38: %% run the model rmeddis@38: tic rmeddis@38: fprintf('\n') rmeddis@38: disp(['Signal duration= ' num2str(length(inputSignal)/sampleRate)]) rmeddis@38: disp([num2str(numChannels) ' channel model: ' AN_spikesOrProbability]) rmeddis@38: disp('Computing ...') rmeddis@38: rmeddis@38: MAP1_14(inputSignal, sampleRate, BFlist, ... rmeddis@38: MAPparamsName, AN_spikesOrProbability, paramChanges); rmeddis@38: rmeddis@38: rmeddis@38: %% the model run is now complete. Now display the results rmeddis@38: UTIL_showMAP(showMapOptions) rmeddis@38: rmeddis@38: if strcmp(signalType,'tones') rmeddis@38: disp(['duration=' num2str(duration)]) rmeddis@38: disp(['level=' num2str(leveldBSPL)]) rmeddis@38: disp(['toneFrequency=' num2str(toneFrequency)]) rmeddis@38: disp(['attenuation factor =' ... rmeddis@38: num2str(DRNLParams.rateToAttenuationFactor, '%5.3f') ]) rmeddis@38: disp(['attenuation factor (probability)=' ... rmeddis@38: num2str(DRNLParams.rateToAttenuationFactorProb, '%5.3f') ]) rmeddis@38: disp(AN_spikesOrProbability) rmeddis@38: end rmeddis@38: disp(paramChanges) rmeddis@38: rmeddis@38: rmeddis@38: rmeddis@38: %% superimpose Lieberman (1996) data rmeddis@38: rmeddis@38: global MOCattenuation rmeddis@38: MOCdB=20*log10(MOCattenuation); rmeddis@38: MOCtime=dt:dt:dt*length(MOCdB); rmeddis@38: rmeddis@38: % scale up DPOAE results to the running maximum MOC dB rmeddis@38: steadyMOCminimum=mean(MOCdB(MOCtime>2 & MOCtime<4.5)); rmeddis@38: scalar=steadyMOCminimum/steadyMinimum; rmeddis@38: rmeddis@38: figure(90), clf rmeddis@38: plot(MOCtime,MOCdB), hold on rmeddis@38: plot(LibermanData(:,1)-2.5,scalar*LibermanData(:,2),'r:','linewidth',4) rmeddis@38: legend({'MAP', 'DPOAE'},'location', 'east') rmeddis@38: title('Compare Liberman (1996) DPOAE data with MAP MOC') rmeddis@38: xlabel('time (s)'), ylabel('MOC attenuation/ DPOAE reduction') rmeddis@38: if strcmp(saveAN_spikesOrProbability,'probability') rmeddis@38: text(0,2,['MOCtau= ' num2str(DRNLParams.MOCtauProb)]) rmeddis@38: else rmeddis@38: text(0,2,['MOCtau= ' num2str(DRNLParams.MOCtau)]) rmeddis@38: end rmeddis@38: set(gcf, 'name', 'Liberman compare') rmeddis@38: rmeddis@38: PSTHbinwidth=0.001; rmeddis@38: rmeddis@38: % show the source of the MOC activity rmeddis@38: figure(89) rmeddis@38: if strcmp(saveAN_spikesOrProbability,'probability') rmeddis@38: % brainstem activity rmeddis@38: PSTH=UTIL_PSTHmaker... rmeddis@38: (ANprobRateOutput(2,:), dt, PSTHbinwidth)*dt/PSTHbinwidth; rmeddis@38: else rmeddis@38: % AN probability rmeddis@38: PSTH=UTIL_PSTHmaker(ICoutput(2,:), dtSpikes, PSTHbinwidth)*dt/PSTHbinwidth; rmeddis@38: end rmeddis@38: rmeddis@38: time=PSTHbinwidth:PSTHbinwidth:PSTHbinwidth*length(PSTH); rmeddis@38: plot(time, PSTH) rmeddis@38: set(gcf,'name', 'Lieberman') rmeddis@38: title(saveAN_spikesOrProbability) rmeddis@38: rmeddis@38: toc rmeddis@38: path(restorePath) rmeddis@38: rmeddis@38: % figure(88), plot(MOCattenuation)