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: