THE AUDITORY MODELING TOOLBOX

Applies to version: 1.5.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 (time series).
fsstim Sampling frequency of stimulus (in Hz).
fc Frequency vector containing the CFs (in Hz). Use logspace(log10(flow), log10(fhigh),numCF) to replicate the results from Bruce et al. (2018).
varargin Various optional flags and key-value pairs.

Output parameters:

output

A struct containing the various modelstage outputs. It contains the following fields:

cohcs: Actual COHCs used in the simulations.

cihcs: Actual CIHCs used in the simulations.

sponts: Actual spontanous rates (SRs) used in the simulations.

tabss: Actual absolute timings used in the simulations.

trels: Actual relative timings used in the simulations.

neurogram_ft: Fine-timing neurogram calculated from the synapse output [time CFs].

t_ft: Time axis of neurogram_ft (in s).

neurogram_mr: Average of spike count in each PSTH bin [time CF].

t_mr: Time axis of neurogram_mr (in s).

neurogram_Sout: Synapse output [time CF].

t_Sout: Time axis of neurogram_Sout (in s).

fc: Actual CFs used in the simulations (in Hz).

psth_ft: Fine-timing PSTH [time CF].

meanrate: Average spiking rate (in spikes/s) [time CF].

varrate: Variance of the spiking rate (in spikes/s) [time CF].

synout: Output of a synapse [time CF Syn]. Available only if 'outputPerSynapse' is specified.

info A struct containing the parameter settings applied. Available only if 'outputPerSynapse' is specified.

Description:

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

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 (in 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. cohcs 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. cihcs 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 when 'autoSR' or 'specificSRautoTiming' specified.
'numM',nm Number of nerve fibres with medium SR. Used only when 'autoSR' or 'specificSRautoTiming' specified.
'numH',nh Number of nerve fibres with high SR. Used only when 'autoSR' or 'specificSRautoTiming' specified.
'lossL',lls Loss of low-SR fibres ranging from 0 (no fibres) to 1 (healthy, all fibres). Used only when 'autoSR' or 'specificSRautoTiming' specified.
'lossM',lms Loss of medium-SR fibres ranging from 0 (no fibres) to 1 (healthy, all fibres). Used only when 'autoSR' or 'specificSRautoTiming' specified.
'lossH',lhs Loss of high-SR fibres ranging from 0 (no fibres) to 1 (healthy, all fibres). Used only when 'autoSR' or 'specificSRautoTiming' specified.
'SRL',SRL SR of the low-SR fibres (in spikes/s). Used only when 'specificSRautoTiming' specified.
'SRM',SRM SR of the medium-SR fibres (in spikes/s). Used only when 'specificSRautoTiming' specified.
'SRH',SRH SR of the high-SR fibres (in spikes/s). Used only when 'specificSRautoTiming' specified.
'numsponts',n Overall numbers of fibers. Used only when 'specificSR' specified.
'spont',spont SR (in spikes/s). Can be scalar or size of fc. Used only when 'specificSR' specified.
'tabs',tabs Absolute timings (in s). Can be scalar or size of fc. Used only when 'specificSR' specified.
'trel',trel Relative timings (in s). Can be scalar or size of fc. Used only when 'specificSR' specified.
'psthbinwidth_mr',psthbw
 Mean-rate binwidth (in 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, and PSTH calculation. Default is 1.
'reptime',rt Length of one repetition of the stimuli with pause. Default is 1.2 times the stimulus duration.
'fsmod',fsmod Model sampling rate (in Hz). The model was tested for fsmod between 100 kHz and 500 kHz. Default is 200 kHz for cats and 100 kHz for humans.

bruce2018 also accepts the following flags:

'fitaudiogram' Calculate the hearing-loss factors cihcs and cohcs from the frequencies ag_fs and threshold shifts ag_dbloss (in 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, respectively. 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, and the SRs SRL, SRM, and SRH, respectively.
'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.

Note: If 'outputPerSynapse' is specified, psth_ft, meanrate, and varrate have the dimensions [time CF Syn] (with Syn being the number of synapses).

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 ]