THE AUDITORY MODELING TOOLBOX

Applies to version: 1.4.0

View the code

Go to function

BRUCE2018 - Auditory-nerve filterbank (improved synapse)

Usage:

[output] = bruce2018(stim,fsstim, fc);
[output,info] = bruce2018(stim, fsstim, fc, varargin);

Input parameters:

stim Pressure waveform of stimulus (timeseries)
fsstim Sampling frequency of stimulus
fc Frequency vector containing the CFs. Use logspace(log10(flow), log10(fhigh),numCF) to replicate the results from Bruce et al. (2018).
varargin various flags and key-value pairs.

Output parameters:

output

a struct containing the various modelstage outputs The struct output contains:

  • cohcs: actually used COHCs in the simulations
  • cihcs : actually used CIHCs in the simulations
  • sponts : actually used spontanous rates in the simulations
  • tabss : actually used absolute timings in the simulations
  • trels : actually used relative timings in the simulations
  • neurogram_ft : (fine-timing) neurogram calculated from the synapse output [time CFs]
  • t_ft : time axis of neurogram_ft (s)
  • neurogram_mr : average of spike count in each PSTH bin [time CF]
  • t_mr : time axis of neurogram_mr (s)
  • neurogram_Sout : synapse output [time CF]
  • t_Sout : time axis of neurogram_Sout (s)
  • fc : actually used CFs (Hz)
  • psth_ft : fine-timing PSTH [time CF]
  • meanrate : average spiking rate (spikes/s) [time CF]
  • varrate : variance of the spiking rate (spikes/s) [time CF]

Description:

bruce2018(...) returns modeled responses of multiple AN fibers tuned to various characteristic frequencies characterstic.

Please cite the references below if you use this model.

This function takes the following optional key/value pairs:

'ag_fs',ag_fs Frequencies at which the audiogram should be evaluated. Used only in fitaudiogram.
'ag_db',ag_db Hearing loss (dB) at frequencies ag_fs. Used only in fitaudiogram.
'cohcs',cohcs OHC scaling factors: 1 denotes normal OHC function (default); 0 denotes complete OHC dysfunction. Can be a vector of the size of fc or a scalar. Not used in fitaudiogram.
'cihcs',cihcs IHC scaling factors: 1 denotes normal IHC function (default); 0 denotes complete IHC dysfunction. Can be a vector of the size of fc or a scalar. Not used in fitaudiogram.
'numL',nl number of nerve fibres with low SR. Used only in autoSR or specificSRautoTiming.
'numM',nm number of nerve fibres with medium SR. Used only in autoSR or specificSRautoTiming.
'numH',nh number of nerve fibres with high SR. Used only in autoSR or specificSRautoTiming.
'lossL',lls loss of low-SR fibres ranging from 0 (no fibres) to 1 (healthy, all fibres). Used only in autoSR or specificSRautoTiming.
'lossM',lms loss of medium-SR fibres ranging from 0 (no fibres) to 1 (healthy, all fibres). Used only in autoSR or specificSRautoTiming.
'lossH',lhs loss of high-SR fibres ranging from 0 (no fibres) to 1 (healthy, all fibres). Used only in autoSR or specificSRautoTiming.
'SRL',SRL SR of the low-SR fibres (spikes/s). Used only in specificSRautoTiming.
'SRM',SRM SR of the medium-SR fibres (spikes/s). Used only in specificSRautoTiming.
'SRH',SRH SR of the high-SR fibres (spikes/s). Used only in specificSRautoTiming.
'numsponts',n Overall numbers of fibers. Used only in specificSR.
'spont',spont SR (spikes/s). Can be scalar or size of fc. Used only in specificSR.
'tabs',tabs Absolute timings (s). Can be scalar or size of fc. Used only in specificSR.
'trel',trel Relative timings (s). Can be scalar or size of fc. Used only in specificSR.
'psthbinwidth_mr',psthbw
 mean-rate binwidth (s).
'windur_ft',winft fine-timing neurogram window length.
'windur_mr',winmr mean-rate neurogram window length.
'nrep',nrep Number of repetitions for the mean rate, rate variance & psth calculation. Default is 1.
'reptime',rt length of one repetition of the stimuli with pause. Default is 1.2 * stimulus duration.
'fsmod',fsmod Model sampling rate. It is possible to run the model at a range of fsmod between 100 kHz and 500 kHz. Default value is 200 kHz for cats and 100 kHz for humans.

bruce2018 accepts the following flags:

'fitaudiogram' Calculate the hearing-loss factors cihcs and cohcs from the frequencies ag_fs and threshold shifts ag_dbloss (dB) by using bruce2018_fitaudiogram. The default parameters reflect a healthy cochlea.
'no_fitaudiogram' Default. Use directly provided cihcs and cohcs either as scalars or size of fc. The default parameters reflect a healthy cochlea.
'human' Default. Use model parameters for humans.
'cat' Use model parameters for cats.
'fixedFGn' Default. Fractional Gaussian noise will be the same in every simulation.
'varFGn' Fractional Gaussian noise will be different in every simulation.
'approxPL' Default. Use approxiate implementation of the power-law functions.
'actualPL' Use actual implementation of the power-law functions.
'outputPerSynapse' Output the synapse output of each individual nerve fibre. This can considerably slow down the calculations.
'outputPerCF' Default. Output the average results over all synapses. This mode is faster than outputPerSynapse.
'autoSR' Generate the parameters for the AN population wil be generated by the function bruce2018_generateanpopulation based on the number of low, medium, and high SR fibres numL, numM, numH. Default.
'specificSRautoTiming'
 Generate the timing parameters for the AN population by the function bruce2018_generateanpopulation based on the number of low, medium, and high SR fibres numL, numM, numH. The spontanous rates are provided in SRL, SRM, and SRH.
'specificSR' Do not generate AN parameters. The overall number of fibres numsponts, spontanouse rate spont, absolute timing tabs, and relative timing trel must be provided.

If 'outputPerSynapse' is specified, psth_ft, meanrate, and varrate have the dimensions [time CF Syn] (with Syn as the number of synapses) and output additionally contains:

'synout' the output of a synapse [time CF Syn]
'info' a struct containing the parameter settings applied

References:

I. C. Bruce, Y. Erfani, and M. S. R. Zilany. A phenomenological model of the synapse between the inner hair cell and auditory nerve: Implications of limited neurotransmitter release sites. Hearing Research, 360:40--54, 2018. [ http ]