Abstract
In cochlear implant (CI) recipients, measurement of the electrically evoked compound action potential (ECAP) is used in clinical routine to prove the electrode-nerve interface and to support fitting, especially in very young children. To record an ECAP amplitude growth function (AGF), electrical stimuli are presented at increasing intensities up to the maximum acceptable loudness (MAL). However, a continuously monotonically rising sound can quickly be perceived as unpleasant and thus impair the success rate of ECAP threshold determination. The present study investigates whether perceptual and objective parameters, which are involved in ECAP measurements, depend on the type of stimulus presentation. In 27 subjects with a CI (28 implants), stimuli of five different AGF patterns were presented ten times in random order at a medial electrode contact, resulting in a total of 50 AGFs per implant. One pattern (Pclin) in which the stimuli increase strictly monotonically is already implemented in the clinical software (MAESTRO). The present study also employed rhythmic patterns and patterns with pauses. The behavioral threshold (THR), the maximum comfortable level (MCL), the MAL, and objective ECAP threshold were analyzed. ECAP threshold values did not change significantly when different stimulus paradigms were applied. Rhythmic AGF patterns were perceived as less unpleasant and enabled a higher MAL compared to Pclin. This in turn led to improved ECAP threshold determination, with the success rate rising to 86.7%, compared to 79.1% for Pclin.
1 Introduction
A cochlear implant (CI) is a medical device that restores hearing to individuals with severe to profound sensorineural hearing loss. The CI electrode array is surgically implanted in the scala tympani of the cochlea. The patient usually perceives an auditory sensation when neural tissue is activated by sufficient stimulus intensity. However, measuring the evoked action potentials requires a response from multiple nerve fibers. This response is referred to as an electrically evoked compound action potential (ECAP), and the patient often perceives a clearly audible to loud sound.
ECAP measurement in CI users has become a valuable objective method (He et al., 2017). Among many other applications, it enables intra- and postoperative verification of auditory nerve functionality, can help in detecting electrode migration (Gärtner et al., 2014) and device failure (Gärtner et al., 2016; Gärtner et al., 2021a), and supports speech processor fitting (Gärtner et al., 2022). CI manufacturers have integrated ECAP measurement capabilities into their systems. No additional equipment is needed beyond what is already required for speech processor fitting, making it easier to use in daily clinical practice.
There are several ways to estimate the ECAP threshold. Visual inspection of ECAP recordings at different stimulus levels by an experienced clinician seems to be the best solution. However, this is time-consuming and results vary depending on the examiner. Automated detection would therefore be preferable. In fact, many methods are already available (Botros et al., 2007). They all have in common that the stimulus strength increases over time during the ECAP recording. If the associated increase in loudness becomes too uncomfortable for the patient, ECAP threshold determination could fail.
The abbreviation TECAP is used differently in various publications. In the following, TECAP stands for ECAP threshold. In this study, we investigate the influence of different stimulus patterns on behavioral and objective parameters in CI users for the first time. It is known that in normal hearing subjects, rising sounds are perceived as more salient than receding sounds. This perceptual bias for looming sounds (Neuhoff, 1998) is attributed to the biological significance of survival: a continuously rising sound may be perceived as an impending threat. The brain may have evolved a bias or increased sensitivity to rising sounds in order to gain an advantage in escaping from predators. Therefore, the continuously monotonically rising tone during the recording of the ECAP amplitude growth function (AGF) can quickly be perceived as uncomfortable and thus lead to the termination of the measurement. However, the algorithm used for TECAP determination is based on sufficiently strong stimulus responses (see Materials and methods section). The results of this study show how these can be achieved without causing discomfort to the CI user.
2 Materials and methods
2.1 Ethics statement
The Ethics Committee of the Hannover Medical School, Germany, where the data were collected using the research software, approved the study (No. 6586). All participants gave written informed consent before any study-specific procedures were performed. All participants also gave written informed consent for the publication of the respective data.
2.2 Participants
This study comprises 28 cochlear implants (ID01 to ID28) in 27 participants (ID12 and ID17 belong to the same patient). Demographic data are presented in Table 1. Further data on etiology and audiological background are listed in Supplementary Table 1. All measurements were conducted between August 2018 and August 2019. Only adults were included in this study. Another inclusion criterion was that a CI with an i100 platform manufactured by MED-EL Medical Electronics (Innsbruck, Austria) had to be in use for compatibility with the ECAP measurement hardware and software. Users were asked to agree to participate in the study during their regular follow-up visit. There was no additional criterion for selection of participants. No additional visit was required for measurement. Participants were not compensated for their participation, and no travel expenses were paid.
TABLE 1
| Implant ID | Gender | Age at implantation (years) | Duration of hearing with CI at time of measurement (years) | Implant type | Electrode type | Implanted side |
|---|---|---|---|---|---|---|
| ID01 | M | 68.7 | 5.3 | Synchrony | Flex 28 | R |
| ID02 | M | 38.7 | 0.7 | Synchrony | Flex 28 | L |
| ID03 | M | 40.7 | 5.1 | Concerto | Flex 28 | L |
| ID04 | F | 49.0 | 0.4 | Synchrony | Flex 28 | R |
| ID05 | F | 46.3 | 3.2 | Synchrony | Flex 20 | R |
| ID06 | M | 72.4 | 0.7 | Synchrony | Flex 24 | R |
| ID07 | M | 52.8 | 0.4 | Synchrony | Flex 28 | R |
| ID08 | F | 55.5 | 1.1 | Synchrony | Flex 24 | R |
| ID09 | F | 44.5 | 3.1 | Synchrony | Flex 28 | R |
| ID10 | M | 48.8 | 4.0 | Synchrony | Flex 28 | R |
| ID11 | F | 56.0 | 3.0 | Synchrony | Flex 20 | L |
| ID12 | M | 63.2 | 1.2 | Synchrony | Flex 28 | R |
| ID13 | M | 35.5 | 8.5 | Sonata | Standard | L |
| ID14 | F | 63.7 | 6.1 | Concerto | Standard | L |
| ID15 | F | 69.6 | 0.8 | Synchrony | Flex 28 | R |
| ID16 | F | 49.9 | 0.1 | Synchrony | Flex 28 | L |
| ID17 | M | 64.3 | 0.1 | Synchrony | Flex 28 | L |
| ID18 | F | 27.2 | 0.1 | Synchrony | Flex 28 | R |
| ID19 | F | 66.4 | 2.9 | Synchrony | Flex 28 | L |
| ID20 | F | 19.5 | 13.8 | Pulsar | Standard | L |
| ID21 | M | 38.8 | 7.3 | Concerto | Standard | L |
| ID22 | F | 21.1 | 0.8 | Synchrony | Flex 28 | R |
| ID23 | F | 50.5 | 2.7 | Synchrony | Flex 16 | L |
| ID24 | M | 45.9 | 2.3 | Synchrony | Flex 28 | L |
| ID25 | M | 63.3 | 1.4 | Synchrony | Flex 24 | L |
| ID26 | M | 47.1 | 1.1 | Synchrony | Flex 24 | R |
| ID27 | F | 58.2 | 1.1 | Synchrony | Flex 28 | R |
| ID28 | M | 58.0 | 0.6 | Synchrony | Flex 28 | L |
Participant demographics.
M, male; F, female; R, right; L, left. Further audiological background including etiology, onset and duration of hearing loss, pre-implant hearing experience is provided by Supplementary Table 1. The length of the electrode array is as follows: Flex 16: 16 mm, Flex 20: 24 mm, Flex 28: 28 mm, Standard: 31.5 mm.
2.3 Stimuli and ECAP-Recordings
All measurements were conducted using the MAX USB interface box (MED-EL Company, Innsbruck, Austria) with the measurement coil placed over the implant. This system was driven by a personal computer running a research software. Bi-phasic, charge-balanced stimuli were presented to the medial electrode contact E06. Pulse duration was 40 μs (exceptions: 50 μs for ID20 and 60 μs for ID25), interphase gap 2.1 μs, and repetition rate 60 pulses per second. The measurement window captured responses up to 1.7 ms after stimulus onset. To reduce stimulus artifact, pulses were presented in pairs with alternating polarity and the response was subsequently averaged. A zero-amplitude template subtraction algorithm was applied to all recordings to compensate for the implant signature artifact. Four adjacent electrode contacts (E04, E05, E07, and E08) were used to record the response. Therefore, stimuli at a given stimulus strength were presented as alternating pairs four times in a row (0.13 s duration per level). The stimulus level was increased in 174 steps from 0 to 875 cu (where 1 current unit, cu, is equal to 1 μA), corresponding to a maximum charge of 35 nC.
For the first 11 implants, ID01–ID11, measurements were stored with reduced resolution only, and for the remaining 17 implants, ID12–ID28, the storage procedure was revised by making use of an adaptive Sigma-Delta modulation (Zierhofer, 2008) and ECAP measurements were stored identically to the clinical software MAESTRO, version 7 and higher (MED-EL Company, Innsbruck, Austria). For subsequent analysis, ECAP parameterization based on ECAP amplitudes showed no differences, nevertheless the signal-to-noise ratio of single traces was affected by the reduced resolution. This was evaluated for the 17 implants, ID12–ID28, by artificially reducing the data quality and comparing the results: correlation analysis using linear mixed models (LMM) and subject as random effect revealed a correlation of r = 0.984 between TECAPs calculated on the basis of the two storage qualities (see Supplementary Table 2).
The hardware setup and data processing are similar to those of the clinical setup and the clinical software, MAESTRO version 7 and higher (MED-EL Company, Innsbruck, Austria). The custom-made research software used enabled different ways to raise the stimulus level compared to the clinical software.
Five different patterns of raising up the stimulus level were under investigation (Figure 1): the Fine-grain paradigm (Gärtner et al., 2018), as implemented in the clinical software MAESTRO (version 7 and higher), presents stimuli with a continuously increasing level (pattern Pclin) (Figure 1A). As a variation, pauses of 1 s duration were added (pattern P1) (Figure 1B). In addition, three non-monotonic patterns were introduced. Pattern P2 (Figure 1C) had a light and P3 (Figure 1D) a stronger rhythm. To achieve this, pulses were rearranged within a set of 24 stimulus levels (i.e., a duration of 3.2 s). In the last paradigm, stimulus strength was periodically reduced with increasing duration (pattern P4) (Figure 1E), i.e., longer pauses (stimulus reduction) for stronger stimuli. A similar constant growth in stimulus intensity was present throughout the patterns in order to minimize bias for algorithms included in ECAP parameterization.
FIGURE 1

Different stimulus patterns for ECAP AGF recordings. (A) Strictly monotonic increasing, as is the current standard in clinical routine (Pclin). (B) Monotonically increasing interspersed with pauses (P1). (C) Stimulus oscillates with a less pronounced rhythm (P2). (D) Stimulus oscillates in a more pronounced rhythm (P3). (E) Stimulus oscillates with pauses whose duration depends on the stimulus strength (P4). A sound sample for each pattern is available in Supplementary material.
In Figure 1, the black dots and lines refer to the stimulus strength E in nC. The red lines refer to the theoretically perceived loudness L and were calculated using following equation posed by Fu and Shannon (1998):
where q = 2.72. k is a constant factor. As loudness integration window of 200 μs was used.
An AGF was presented 10 times with each of the 5 paradigms (Pclin, P1, P2, P3, P4). The order was pseudorandomized. This resulted in a total of 50 AGFs per implant ID. Sound samples of these five pattern are available as Supplementary material. The study was single-blinded. Therefore, only the conductor of the measurements knew which pattern was currently being tested.
2.4 ECAP parameterization
For all AGF patterns post-processing including filtering of ECAP traces, determination of N- and P-peaks of the ECAP response, parameterization of ECAP AGF with a sigmoid function, as well as selection of the most representative recording electrode was done using the AutoART algorithms (Strahl et al., 2018), as they are implemented in the clinical Software MAESTRO (version 7 and higher).
After recordings of a single run out of 50 for each implant, ECAP recordings for each stimulation pattern were sorted according to stimulus levels, separately for each recording electrode. The four resulting data sets were processed independently to derive four ECAP parameterizations. For each data set, ECAP amplitude at a given stimulus level was calculated using a moving average of ± 3 levels (7 values in total). ECAP amplitude [y(x)] as a function of stimulus level (x) was parameterized by fitting of a sigmoidal function
using the Levenberg-Marquardt algorithm (Marquardt, 1963). Here, y0 (noise level, i.e. the spontaneous activity of the neuron without electrical stimulation, in μV), c (the maximum observable ECAP amplitude at a given recording electrode, in μV), x0 (stimulus strength, in nC, corresponding to the inflection point of a symmetric sigmoidal function), and k (representing the rate of change per nC within the dynamic range of the modeled neuronal population) are the parameters to be fitted (van de Heyning et al., 2016), which were then used to estimate TECAP and the slope of the ECAP AGF at the steepest point of the sigmoid. The latter is defined by the value of the first order derivative at x0:
The stimulation intensity at which the slope intersects with y0 was used to denote TECAP:
From the four parameterizations derived in this way (one for each recording contact), the AutoART algorithms were also used to decide whether a TECAP value could be estimated (i.e., whether the quality of the approximation to a sigmoidal function was sufficient) and which recording electrode provided the best result. Therefore, the final data set for each stimulation electrode contains one recording electrode. If no TECAP could be reliably determined, the remaining objective parameters (y0, x0, c, k, ECAP slope) were also not included in the statistical evaluations.
2.5 Behavioral loudness perception tests
The ECAP AGF measurements were accompanied by behavioral perception tests. Subjects were instructed to press the space bar of the personal computer three times during the presentation of each AGF. The first response indicates the threshold level (THR), i.e., when sound becomes audible. The second response indicates the most comfortable level (MCL), i.e., when sound reaches a pleasantly loud perception. If sound becomes uncomfortably loud, a third response indicates the maximum acceptable loudness (MAL). Due to the limited maximum presentation level of 35 nC, reaching the MAL was not necessarily expected. At the end of the experiment, participants were asked to rate which pattern sounded most pleasant to them. Therefore, the ranking was given as a subjective response, where the subjects could replay representative 5 s for each paradigm.
2.6 Variability and intrinsic parameters of test methods
Loudness levels, such as threshold, are best identified using stepwise approaches that bracket the level of interest (Carhart and Jerger, 1959). Repeated short stimuli avoid effects such as loudness masking, and stepwise approaches facilitate the identification of different loudness levels. The present experiment, however, aims at identifying and quantifying these influencing factors such as looming (continuous, not stepwise increase of stimulus levels) and loudness masking. Ten repetitions for each paradigm were chosen as an a priori compromise between accuracy and duration of the experiment.
To estimate the variability of the determined parameters (THR, MCL, MAL, TECAP, ECAP slope, y0), the maximum, minimum, and mean values were first determined for each parameter in each of the 10 runs per paradigm. The variability is then calculated as the median across all 28 implants of the difference between the maximum and minimum values, divided by the mean value [(max-min)/mean].
2.7 Statistical analysis
Python programming software (version 3.6, 2017), was used to program the research tool and store databases using besides built-in functionality packages numpy (version 1.13), PyQt5 (version 5.10.1, 2018), sqlite3 (version 3.23, 2018), and the “Research Interface Box 2” dynamic link library RIB2.dll (version 1.71, provided by the University of Innsbruck) to stimulate with the CI via MAX Interface Box (Litovsky et al., 2017).
For statistical analyses, python (version 3.7.3, 2019) with updated packages of the research tool and in addition matplotlib (version 3.1.0, 2019) and scipy (version 1.7.1, 2021) were used. Via jupyter notebook (version 6.4.8, 2021) and the package rpy2 (version 3.4.5, 2021) the “R programming language and statistical environment” was used to perform statistical analyses (R Core Team, 2021)(R Core Team, 2021), version 4.1.1, with packages “lmerTest” (Kuznetsova et al., 2017), “lme4” (Bates et al., 2015), “MuMin” (Burnham and Anderson, 2004), and “performance” (Nakagawa et al., 2017), function “lmer” and “r2” for the LMM, “multcomp” (Hothorn et al., 2008), and functions “glht” and “squaredGLMM” for the post-hoc analysis.
For subsequent statistical comparisons in this study, LMMs using fixed and random factors were used. The details of the used factors are given in the appropriate section.
3 Results
3.1 General
The average duration of a complete study session was 43 min 5 s per subject (22 min 44 s—69 min 9 s). Occasionally, the subjects took a break, which averaged 2 min 16 s (8 s—7 min 24 s).
In 8 out of 1,400 cases (0.57%), no perceptual thresholds (THR, MCL, and MAL) could be determined because no data was saved due to an error. Similarly, in 71 out of 1,400 cases (5.07%), no objective measures could be determined (i.e., values for y0, c, k, x0, TECAP, and ECAP slope), again due to an error resulting in no data being saved.
The level of MCL in particular, but also THR and MAL, is chosen individually by each subject, so that comparisons of absolute levels between subjects are not valid. Here we are comparing between different stimulus patterns per subject.
In order to comply with the loudness safety guidelines of the ethics committee, a stimulation limit of 35 nC was set. A ceiling effect was considered to have occurred if a perceptual threshold was above 34 nC. Several subjects reached this limit (see Table 2), and it must be taken into account in the following statistics.
TABLE 2
| Parameter | Pclin | P1 | P2 | P3 | P4 |
|---|---|---|---|---|---|
| THR variability | 19.2% | 17.7% | 14.3% | 16.1% | 21.1% |
| THR ceiling frequency | 0 | 0 | 0 | 0 | 0 |
| MCL variability | 25.5% | 25.2% | 25.7% | 24.4% | 22.8% |
| MCL ceiling frequency | 2 (m: 2, f: 0) |
2 (m: 2, f: 0) |
2 (m: 2, f: 0) |
3 (m: 2, f: 1) |
4 (m: 3, f: 1) |
| MAL variability | 29.4% | 26.8% | 26.8% | 29.1% | 18.5% |
| MAL ceiling frequency | 13 (m: 8, f: 5) |
14 (m: 9, f: 5) |
21 (m: 11, f: 10) |
22 (m: 11, f: 11) |
20 (m: 11, f: 9) |
Variability of behavioral limits and frequency of ceiling effects for each of the five stimulus patterns.
The variability of threshold (THR), maximum comfort level (MCL), and maximum acceptable loudness (MAL) is expressed as median across implants of (max-min)/mean for each of the five stimulus paradigms (Pclin, P1, P2, P3, P4). If a behavioral limit exceeded a limit of 34 nC in more than 3 out of 10 runs, a ceiling effect was considered. The frequency of such ceiling effects is specified according to gender (m, male; f, female).
Besides averaged values, 10 repetitions as a test-retest scenario provide insights. Over the duration of the experiment, y0 showed a high variation (see Table 3). The algorithms work as expected in terms of TECAP accuracy as well as for ECAP slope accuracy, with the latter exhibiting greater variation. This is due to the fact that the selected parameter set is chosen among several recording electrodes. While the TECAP is less affected by the recording electrodes, the ECAP amplitude, and, consequently the ECAP slope change with distance, i.e., with the specific location of the stimulation and recording electrodes. For the evaluation of differences between the paradigms, the factor recording electrode is taken into account.
TABLE 3
| Parameter | Pclin | P1 | P2 | P3 | P4 |
|---|---|---|---|---|---|
| TECAP variability | 16.8% | 14.0% | 14.7% | 15.9% | 13.0% |
| ECAP slope variability | 45.8% | 41.1% | 40.1% | 39.7% | 33.7% |
| y0 variability | 169.0% | 123.7% | 149.3% | 196.0% | 179.5% |
| Success rate | 79.1% | 77.9% | 82.2% | 84.4% | 86.7% |
| Effective number of runs | 259 | 265 | 270 | 265 | 270 |
Variability of ECAP parametrization and success rate in TECAP determination for each of the five stimulus patterns.
The variability of TECAP, ECAP slope, and y0 is expressed as median across implants of (max-min)/mean for each of the five stimulus paradigms (Pclin, P1, P2, P3, P4). In very few cases, a series of measured values could not be saved, meaning that the effective number of runs was < 280.
Overall, the ECAP parameterization in this study was found to be robust in a test-retest scenario with a TECAP variability of 16.8% for the clinically implemented pattern Pclin.
3.2 ECAP measurements
The TECAP determination success rate was defined as the number of TECAP values determined divided by the number of runs (usually 10). For safety reasons, any unexpected event was handled in abortion of the current run and starting with the next one. In those instances, the attainment of 10 runs was impeded by data loss that occurred with the research software. The success rate was 79% for Pclin, and for the new AGF stimulus patterns 78% (P1), 82% (P2), 84% (P3), and 87% (P4). More details are shown in Table 3.
The influence of stimulus pattern on estimated TECAP is shown in Figure 2. Figures 2A–D show the TECAP of the Pclin pattern on the x-axis, while the TECAP of the P1, P2, P3, and P4 patterns are shown on the y-axis. Figure 2E shows the difference in TECAP of the 4 newly designed paradigms compared to the classic strict Fine-grain pattern Pclin. Correlation analysis revealed similar TECAP for all paradigms with Pearson’s correlation coefficients (vs. Pclin) = 0.899 (P1), 0.881 (P2), 0.885 (P3), and 0.908 (P4). Pairwise comparisons did not reveal any statistical difference of TECAP between Pclin and any of the other stimulus patterns.
FIGURE 2

Effect of stimulus pattern on TECAP. (A–D) The x-axis shows TECAP values (nC) determined by using the strict stimulus pattern Pclin. The y-axis shows TECAP values (nC) determined by using the newly introduced stimulus patterns P1, P2, P3, P4. Individual measurements (up to 10 per paradigm) are shown as colored dots (28 colors for 28 implants) without borders. The corresponding mean values per implant are shown as colored dots with black borders if a TECAP value could be determined in more than 3 cases, otherwise as stars with black borders. Lilac lines are shown to guide the eye. (E) Overview of data already shown in (A–D) as difference to pattern Pclin.
3.2.1 Linear mixed model
We performed LMM analysis as explained in Statistical Analysis. Investigations of ECAP parameter (TECAP, ECAP slope, y0) were performed using the model
Estimates of the effects and significance are shown in Table 4.
TABLE 4
| Parameter | P1 | P2 | P3 | P4 |
|---|---|---|---|---|
| TECAP (nC) | +0.0 | –0.1 | –0.2 | –0.1 |
| ECAP slope (μV/nC) | –0.7 | +0.5 | +1.0 | +0.5 |
| y0 (μV) | –1.3 | –1.0 | +1.7 | +3.2 |
Linear mixed model analysis of objective ECAP parameters TECAP, ECAP slope, and y0.
“ECAP parameter∼ paradigm + (1| subject) + (1| recording_electrode)” revealed no significant differences (p < 0.05) for tested parameters, except for: TECAP, ECAP slope for “P3” resp. y0 for “P4.”
In contrast to single pairwise comparisons, TECAP and ECAP slope showed a significant difference for paradigm P3 (p < 0.05) and y0 for paradigm P4 (p = < 0.05).
3.3 Behavioral loudness perception tests
In very few of the 1,400 runs were the perception levels reported by the test subjects clearly incorrect. This concerned three runs in which THR was reported as zero (ID11, one run with P2; ID23, one run each with P1 and P3). Three other runs in which THR was below 3 nC could be discussed, as this value appears too low (ID 23, one run each with P1 and Pclin; ID26, one run with P2). In one of these subjects (ID23), MCL and MAL were reported with values below 3 nC, which was highly unlikely (one run each with P1, P3, and Pclin). In one case (ID13, one run with pattern P3), the feedback at THR was apparently forgotten, as THR, MCL, and MAL were identical (34.77 nC) and THR was significantly lower in the other runs. Nevertheless, the data from each test subject was included in the evaluation as provided by the test subject without correction.
The results of the behavioral tests are shown in Figure 3. In Figures 3A–D, values of stimulus charge for THR, MCL, and MAL values are indicated by different symbols. The size of the symbol reflects the amount of data with the same value. The x-axis shows the stimulus charge of the Pclin paradigm. The y-axis shows the stimulus charge of the newly introduced AGF paradigms. The relative difference compared to the Pclin pattern is shown in Figure 3E. The differences become more pronounced with increasing stimulus strength. As expected, our motivated AGF patterns had an effect on loudness perception. All data are available in Supplementary Table 3.
FIGURE 3

Effect of stimulus pattern on loudness perception. (A–D) The x-axis shows the stimulus charge (nC) of behavioral limits [threshold (THR), maximum comfort level (MCL), and maximum acceptable loudness (MAL)] for the strict Pclin pattern. The y-axis shows the stimulus charge (nC) of the corresponding limits for the newly introduced AGF paradigms. The lilac line indicates the same stimulus strengths for orientation. (E) Overview of data already shown in (A–D) as difference to pattern Pclin.
Results of subjective preferences are shown in Figure 4. Paradigms with rhythm P2 and P3 were rated highest. Interestingly, the bilateral patient with implant ID12 and ID17 ranked both ears different. Since he was the only candidate who participated in the study with both of his CIs, one can only speculate about the causes. ID12 was measured first, followed by ID17 2 days later. Different preferences could have been due to daily moods. The different periods of use/acclimatization could also have had an influence. ID17 was implanted 1 year after ID 12 (see Table 1). The preference did not correlate with TECAP determination success rate.
FIGURE 4

Subjective preferences regarding the various stimulus paradigms. The graph shows how often each pattern was rated as the most preferred (1) to the least preferred (5).
3.3.1 Linear mixed model
We performed LMM analysis as explained in Statistical Analysis. Investigations of behavioral parameter (THR, MCL, and MAL) were performed using the model
once with all data, and once using only data below ceiling limit. Estimates of the effects and significance are shown in Table 5.
TABLE 5
| Parameter | P1 | P2 | P3 | P4 |
|---|---|---|---|---|
| THR (nC) | –0.2 (–0.2) | +1.1 (+1.1) | +1.3 (+1.3) | +1.7 (+1.7) |
| MCL (nC) | –0.0 (+0.2) | +1.0 (+1.1) | +2.1 (+2.2) | +4.0 (+3.9) |
| MAL (nC) | +0.1 (0.2) | +1.8 (+1.8) | +1.4 (+1.9) | +3.7 (+2.4) |
Linear mixed model analysis of behavioral parameters THR, MCL, and MAL.
“Behavioral parameter∼ paradigm + (1| subject)” revealed significant differences (p < 0.001) for all three parameters and paradigms P2, P3, and P4, but for no parameter Pclin vs. P1. Limiting data with respect to ceiling showed similar results. Data without ceiling effect in parentheses.
The effect of gender was examined using the LMM behavioral
and revealed a significant effect (p < 0.001) for all behavioral parameters of 1.6 nC for MAL (2.7 nC including ceiled data), 2.2 nC for MCL (3.5 nC including ceiled data), and 1.1 nC (1.1 nC including ceiled data) for THR. Pairwise tests between female and male groups separated by paradigms did not reveal statistical significance.
4 Discussion
The algorithm implemented in the MAESTRO clinical software (version 7 and higher) for determining TECAP uses an approximation of the recorded AGF data to a sigmoidal function. The accuracy increases with increasing stimulus intensity. However, the patient then also perceives signals that are continuously getting louder. In acoustic hearing, a perceptual bias for rising sound is already known (Neuhoff, 1998). Sounds that increase in intensity are perceived as louder and change in loudness more than sounds that decrease in intensity. The effect is more pronounced for tones than for noise (Neuhoff, 2001).
If a perceptual bias for rising tones is also present in CI users, measuring the ECAP AGF with strict up-ramp stimulus intensities (Pclin) may not be the best solution. The novelty of this study lies in the introduction of alternative AGF patterns by adding pauses and rhythm in an attempt to avoid continuous auditory looming, which should ultimately make the measurement more comfortable. The focus of this study was not to find out, whether a perceptual bias is also present in CI users, but to improve the ECAP measurement procedure and TECAP determination. Loudness integration of strict up-ramp stimuli may be reduced by masking effects: soft pulses following strong pulses cannot be heard and are perceived as pauses. Rearranging soft pulses after louder ones thus creates a rhythm and is hypothesized to be more comfortable to the user.
One limitation of this study is that only four patterns were tested in addition to the classic strict Fine-grain AGF pattern Pclin. There is still room for further studies addressing different AGF patterns.
The rationale for the four patterns presented in this study was to reduce a potential perceptual bias for looming sounds and to make overall impression more comfortable for the subjects. Paradigm P1 could not reduce the perceived loudness as much as the other test paradigms (difference at MCL < 0.5 nC), so either the pauses used were too few to inhibit auditory looming, respectively change loudness perception, or the looming effect is not dominant below MAL levels in paradigm Pclin. Paradigm P4, many pauses with stimulus-dependent duration, showed a strong effect (differences P4 – Pclin ∼ 4 nC, which is ∼20%) (see Figure 3E). We conclude that the method used to avoid looming, respectively, loudness masking, has a great potential. Although P4 had the largest loudness reduction, users did not rate it as pleasant (see Figure 4). Pattern P2 was rated as moderately pleasant and pattern P3 as most pleasant, suggesting that users recognize the intended rhythm and prefer stimulation with rhythmic patterns. Since the difference in loudness integration is only minor, auditory looming is more inhibited by the stronger rhythm, i.e., deviation from monotonic increment is greater. Our results show increased MCLs for both paradigms, but the difference P3 – Pclin = 2.1 nC is greater than P2 – Pclin = 1.0 nC. For MAL, ceiling effects bias the results: both paradigms increase MAL to a large extend. From this perspective, we conclude that loudness integration is a stronger effect than auditory looming.
Neuhoff and Heckel (2004) reported that looming has gender-specific effects. They found that as intensity increased, females perceived more loudness change than males, which is consistent with a more pronounced awareness of danger. In our study, we also find gender effects, but we were not able to quantify differences in this effect for individual paradigms with our data and experimental design.
The overall loudness of ECAP recordings can be influenced by different ways of estimating TECAP. Approximating a sigmoid function may require a higher stimulus level. However, there are other ways of estimating TECAP, for example by comparing signals in different time windows (Hoth et al., 2018), which has recently been made applicable to Fine-grain recordings (Gärtner et al., 2021b). We see effects on ECAP recordings for paradigm P4 on y0, which we connect to a larger number of recordings for subthreshold stimulation recorded at different times during the AGF. We also link these additional recordings to the improved test-retest accuracy for pattern P4 (TECAP variability: median across implants of (max-min)/mean < 13%), which was best among all paradigms. The differences for pattern P3 for TECAP and ECAP slope cannot be directly explained by this additional source of noise. In the future, other rhythmic patterns could be examined to investigate possible systematic effects, such as the frequency of changes between strong and soft pulses.
The Fine-grain method (Pclin) has already led to a 5% higher success rate in determining TECAP compared to the former standard ART method (Garrido et al., 2018). The present study shows that TECAP estimation can be further improved by skillfully re-sorting the stimuli, i.e., by introducing rhythm. The use of different non-monotonous AFG patterns may contribute to further improvement.
5 Conclusion
In CI recipients, a monotonically increasing stimulus intensity during ECAP recordings can be perceived as unpleasant. Our study shows that optimized loudness integration increases overall achievable levels, but not comfort. Rhythmic stimuli that deviate more from the monotonic increment show higher attainable levels as well as greater pleasantness. This suggests that in addition to loudness integration, other effects such as the hypothesized looming effect play a role in ECAP recordings. Loudness integration has been observed to have a strong influence on perceived loudness even at very low stimulation rates, such as those used for ECAP recordings. ECAP measurements may provide more comfort when rhythm is introduced without altering the TECAP measurement results.
Sounds that continuously increase in intensity, such as those perceived during ECAP AGF measurement, may compromise successful TECAP determination. A better alternative is to use rhythmic stimuli, which also sound more comfortable to most subjects. The underlying effect may be attributed to a perceptual bias for rising sounds and/or loudness integration. TECAP values did not change significantly while applying different stimulus paradigms.
Statements
Data availability statement
The original contributions presented in this study are included in this article/Supplementary material, further inquiries can be directed to the corresponding author.
Ethics statement
The studies involving humans were approved by Ethics Committee of the Hannover Medical School, Germany. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
LG: Conceptualization, Data curation, Formal analysis, Methodology, Validation, Writing – original draft, Writing – review & editing. KS: Conceptualization, Data curation, Formal analysis, Methodology, Validation, Visualization, Writing – review & editing. TB: Data curation, Investigation, Writing – review & editing. WL: Investigation, Writing – review & editing. TL: Funding acquisition, Supervision, Writing – review & editing. AB: Conceptualization, Funding acquisition, Methodology, Supervision, Writing – review & editing.
Funding
The author(s) declared that financial support was received for this work and/or its publication. This study was funded and supported by MED-EL Medical Electronics (Innsbruck, Austria). The funder provided support in the form of salaries for authors KS, TB, WL. MED-EL did not play a role in the decision to publish this work.
Conflict of interest
KS and TB were employed by MED-EL. WL was an intern at MED-EL when the study was conducted.
The remaining author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declared that generative AI was not used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnins.2025.1711456/full#supplementary-material
References
1
Bates D. Mächler M. Bolker B. Walker S. (2015). Fitting linear mixed-effects models using lme4.J. Stat. Softw.671–48. 10.18637/jss.v067.i01
2
Botros A. van Dijk B. Killian M. (2007). AutoNR: an automated system that measures ECAP thresholds with the Nucleus Freedom cochlear implant via machine intelligence.Artif. Intell. Med.4015–28. 10.1016/j.artmed.2006.06.003
3
Burnham K. P. Anderson D. R. (2004). in Model Selection and Multimodel Inference, edsBurnhamK. P.AndersonD. R. (New York, NY: Springer), 10.1007/b97636
4
Carhart R. Jerger J. F. (1959). Preferred method for clinical determination of pure-tone thresholds.J. Speech Hear. Disord.24330–345. 10.1044/jshd.2404.330
5
Fu Q. J. Shannon R. V. (1998). Effects of amplitude nonlinearity on phoneme recognition by cochlear implant users and normal-hearing listeners.J. Acoust. Soc. Am.1042570–2577. 10.1121/1.423912
6
Garrido L. C. Schwarz K. Lauss K. Vernetta C. P. Kontides A. Gomez M. D. et al (2018). Comparison of a traditional and novel evoked compound action potentials recording approach and evoked auditory brainstem responses in pediatric cochlear implants users.J. Int. Adv. Otol.14353–358. 10.5152/iao.2018.5591
7
Gärtner L. Bräcker T. Kals M. Penninger R. T. Billinger-Finke M. Lenarz T. et al (2022). ARTFit-A quick and reliable tool for performing initial fittings in users of MED-EL cochlear implants.Life12:269. 10.3390/life12020269
8
Gärtner L. Büchner A. Illg A. Lenarz T. (2021a). Hidden electrode failure in a cochlear implant user.Laryngoscope131E1275–E1278. 10.1002/lary.29023
9
Gärtner L. Lenarz T. Büchner A. (2018). Fine-grain recordings of the electrically evoked compound action potential amplitude growth function in cochlear implant recipients.Biomed. Eng. Online17:140. 10.1186/s12938-018-0588-z
10
Gärtner L. Lesinski-Schiedat A. Büchner A. Lenarz T. (2016). [Headache and degraded speech recognition in a patient with unilateral cochlear implant].Laryngorhinootologie95559–560. 10.1055/s-0042-105763
11
Gärtner L. Spitzer P. Lauss K. Takanen M. Lenarz T. Hoth S. (2021b). Optimized SNR-based ECAP threshold determination is comparable to the judgement of human evaluators.PLoS One16:e0259347. 10.1371/journal.pone.0259347
12
Gärtner L. Würfel W. Büchner A. Lenarz T. (2014). [Unexpected severe degradation of speech recognition in a patient with cochlear implant].Laryngorhinootologie93398–400. 10.1055/s-0034-1372590
13
He S. Teagle H. F. B. Buchman C. A. (2017). The electrically evoked compound action potential: from laboratory to clinic.Front. Neurosci.11:339. 10.3389/fnins.2017.00339
14
Hoth S. Spitzer P. Praetorius M. (2018). A new approach for the determination of ECAP thresholds.Cochlear Implants Int.19104–114. 10.1080/14670100.2017.1402472
15
Hothorn T. Bretz F. Westfall P. (2008). Simultaneous inference in general parametric models.Biom J.50346–363. 10.1002/bimj.200810425
16
Kuznetsova A. Brockhoff P. B. Christensen R. H. B. (2017). lmerTest package: tests in linear mixed effects models.J. Stat. Softw.821–26. 10.18637/jss.v082.i13
17
Litovsky R. Y. Goupell M. J. Kan A. Landsberger D. M. (2017). Use of research interfaces for psychophysical studies with cochlear-implant users.Trends Hear.21:2331216517736464. 10.1177/2331216517736464
18
Marquardt D. W. (1963). An algorithm for least-squares estimation of nonlinear parameters.J. Soc. Indust. Appl. Math.11431–441. 10.1137/0111030
19
Nakagawa S. Johnson P. C. D. Schielzeth H. (2017). The coefficient of determination R2 and intra-class correlation coefficient from generalized linear mixed-effects models revisited and expanded.J. R. Soc. Interface14:20170213. 10.1098/rsif.2017.0213
20
Neuhoff J. G. (1998). Perceptual bias for rising tones.Nature395123–124. 10.1038/25862
21
Neuhoff J. G. (2001). An adaptive bias in the perception of looming auditory motion.Ecol. Psychol.1387–110. 10.1207/S15326969ECO1302_2
22
Neuhoff J. G. Heckel T. (2004). “Sex differences in perceiving auditory “looming” produced by acoustic intensity change,” in Proceedings of the 10th International Conference on Auditory Display (ICAD2004), (Sydney).
23
Strahl S. Dierker A. Spitzer P. Schwarz K. (2018). “AutoART-A system for automatic determination of eCAP thresholds,” in 21st Annual Meeting of the German Society of Audiology, (Halle).
24
van de Heyning P. Arauz S. L. Atlas M. Baumgartner W. D. Caversaccio M. Chester-Browne R. et al (2016). Electrically evoked compound action potentials are different depending on the site of cochlear stimulation.Cochlear Implants Int.17251–262. 10.1080/14670100.2016.1240427
25
Zierhofer C. M. (2008). Adaptive delta-sigma modulation for enhanced input dynamic range.Eur. J. Adv. Signal Process.2008:439203. 10.1155/2008/439203
Summary
Keywords
auditory looming, cochlear implant, electrically evoked compound action potential, fine-grain, loudness perception
Citation
Gärtner L, Schwarz K, Bräcker T, Lamping W, Lenarz T and Büchner A (2026) Rhythm increases perceptual tolerance when measuring electrically evoked compound action potentials in cochlear implant recipients. Front. Neurosci. 19:1711456. doi: 10.3389/fnins.2025.1711456
Received
23 September 2025
Revised
29 November 2025
Accepted
30 November 2025
Published
05 January 2026
Volume
19 - 2025
Edited by
Bruna Mussoi, University of Tennessee Health Science Center (UTHSC), United States
Reviewed by
Colleen Psarros, Macquarie University, Australia
Florian Schertenleib, University of Oldenburg, Germany
Updates
Copyright
© 2026 Gärtner, Schwarz, Bräcker, Lamping, Lenarz and Büchner.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Lutz Gärtner, Gaertner.Lutz@mh-hannover.de
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.