bruce2018
Auditory-nerve filterbank with improved synapse model

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. Size: (time x CFs).
t_ft: Time axis (in s) of neurogram_ft.
neurogram_mr: Average of spike count in each PSTH bin. Size: (time x CF).
t_mr: Time axis (in s) of neurogram_mr.
neurogram_Sout: Synapse output. Size: (time x CF).
t_Sout: Time axis (in s) of neurogram_Sout.
fc: Actual CFs (in Hz) used in the simulations.
psth_ft: Fine-timing PSTH. Size: (time x CF).
meanrate: Average spiking rate (in spikes/s). Size: (time x CF).
varrate: Variance of the spiking rate (in spikes/s). Size: (time x CF).
synout: Output of a synapse. Size: (time x CF x 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).

bruce2018(...) also accepts the following flags:

`'no_fitaudiogram'`	Use directly provided cihcs and cohcs either as scalars or size of fc. The default parameters reflect a healthy cochlea. Default.
`'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. Opposite of `'no_fitaudiogram'`.
`'human'`	Use model parameters for humans. Default.
`'cat'`	Use model parameters for cats. Opposite of `'human'`.
`'fixedFGn'`	Fractional Gaussian noise will be the same in every simulation. Default.
`'varFGn'`	Fractional Gaussian noise will be different in every simulation. Opposite of `'fixedFGn'`.
`'approxPL'`	Use approxiate implementation of the power-law functions. Default.
`'actualPL'`	Use actual implementation of the power-law functions. Opposite of `'approxPL'`.
`'outputPerCF'`	Output the average results over all synapses. This mode is faster than `'outputPerSynapse'`. Default.
`'outputPerSynapse'`	Output the synapse output of each individual nerve fibre. This can considerably slow down the calculations. Opposite of `'outputPerCF'`.
`'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 spontaneous rates given in SRL, SRM, and SRH, respectively. Can not be used with `'autoSR'` or `'specificSR'`.
`'specificSR'`	Do not generate AN parameters. The overall number of fibres numsponts, spontanouse rate spont, absolute timing tabs, and relative timing trel must be specifically provided. Can not be used with `'autoSR'` or `'specificSRautoTiming'`.
`'progress'`	Display the calculation progress. Default.
`'silent'`	No progress display. Opposite of `'progress'`.

Note: If 'outputPerSynapse' is specified, psth_ft, meanrate, and varrate have the size (time x CF x syn), with syn being the number of synapses.

This function takes the following optional key-value pairs:

`'ag_fs',ag_fs`	Frequencies (in Hz) at which the audiogram should be evaluated. Used only when `'fitaudiogram'` is specified.
`'ag_db',ag_db`	Hearing loss (in dB) at frequencies ag_fs. Used only in when `'fitaudiogram'` is specified.
`'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. Used only in when `'no_fitaudiogram'` is specified.
`'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. Used only in when `'no_fitaudiogram'` is specified.
`'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 (in spikes/s) of the low-SR fibres. Used only when `'specificSRautoTiming'` specified.
`'SRM',SRM`	SR (in spikes/s) of the medium-SR fibres. Used only when `'specificSRautoTiming'` specified.
`'SRH',SRH`	SR (in spikes/s) of the high-SR fibres. 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 (in samples).
`'windur_mr',winmr`	Mean-rate neurogram window length (in samples).
`'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.

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 ]