BRUCE2018 - Auditory-nerve filterbank (improved synapse)

Usage

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

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.

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.

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.

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]

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: