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: