[output] = bruce2018(stim,fsstim, fc); [output,info] = bruce2018(stim, fsstim, fc, varargin);
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 | 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. |
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).
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 ]