 |
||||
|
||||
|
||||
|
||||
|
||||
The Journal of Neuroscience, June 1, 2003, 23(11):4677-4688
Previous Article | Next Article
Frequency-Specific Interaural Level Difference Tuning Predicts Spatial Response Patterns of Space-Specific Neurons in the Barn Owl Inferior Colliculus
Michael L. Spezio andTerry T. Takahashi Institute of Neuroscience, University of Oregon, Eugene, Oregon 97403
Abstract
Top
Abstract
Introduction
Space-specific neurons in the barn owl's inferior colliculus have spatial receptive fields (RFs) because of sensitivity to interaural time difference and frequency-specific interaural level difference (ILD). These neurons are assumed to be tuned to the frequency-specific ILDs occurring at their spatial RFs, but attempts to assess this tuning with traditional narrowband stimuli have had limited success. Indeed, tuning assessed in this manner, when processed via a linear model of spectral integration, typically explains only approximately half the variance in spatial response patterns. Here we report our findings that frequency-specific ILD tuning of space-specific neurons, when assessed from responses to broadband stimuli, predicted nearly 75% of the variance in spatial responses, using a linear model of spectral integration (p < 0.0001; n = 97 neurons). Furthermore, when we tested neurons using only those frequencies we found to be spatially relevant, we saw that their responses were similar to those elicited by broadband stimuli. When we used frequencies not identified as spatially relevant, such similarity was lacking. Furthermore, spectral components that elicited high firing rates when presented as narrowband stimuli were found in several cases to be irrelevant for or detrimental to the definition of spatial RFs. Thus, neurons achieved sharp spatial tuning by selecting for ILDs of a subset of spectral components in noise, some of which were not identified using narrowband stimuli.
Key words: sound localization; auditory space; reverse correlation; spatial hearing; receptive field; virtual auditory space; binaural
Introduction
Sound localization cues, such as the average binaural sound level (ABL) or the interaural time difference (ITD) and frequency-specific interaural level difference (ILD), vary by frequency because of the directional filtering of the head and external ears. Individual frequency components can only be localized to wide regions of space. To localize sounds precisely, the auditory system must integrate these cues across frequency.Studies of the neural mechanisms of sound localization have often assumed that neuronal tuning to spatial cues could be measured accurately by presenting narrowband stimuli (Kuwada and Yin, 1983
; Yin and Kuwada, 1983
The space-specific neurons of the auditory space map in the barn owl's inferior colliculus (IC) are a case in point. These neurons integrate ITD and ILD across frequency to compute their spatial receptive fields (RFs) (Knudsen and Konishi, 1978a
One possible explanation for this modest predictivity is that responses of space-specific neurons to tones are not good estimates of the frequency-specific ILD tuning that is operative in forming responses to broadband stimuli. Indeed, space-specific neurons respond poorly to tone pips. To achieve reliable firing, Euston and Takahashi (2002
In the present study, we measured ILD-alone responses using broadband noise and identified the frequency-specific ILD tuning directly from these responses. This tuning yielded better predictions of the ILD-alone response pattern of the neuron than did the tuning measured with tones. Moreover, presentation of only those frequencies identified from the response to noise produced ILD-alone responses that resembled those obtained with broadband noise. This suggests that our approach effectively identified the spatially relevant spectral components in the noise, and that these relevant components were not in fact as reliably identified using narrowband stimuli.
Materials and Methods
Subjects. Eight adult barn owls (Tyto alba) weighing 390–500 gm were used in this study. All procedures were approved by the Institutional Animal Care and Use Committee at the University of Oregon.Stimuli in virtual auditory space. All stimuli were presented in virtual auditory space (VAS). The stimulus waveforms were filtered with each bird's own head-related transfer functions (HRTFs) measured and processed as described by Keller et al. (1998
Terminology. Fully-cued refers to stimuli that simulate free-field conditions in all respects and to the neuronal responses to such stimuli. In fully-cued stimuli, spatial cues, including ITD, ILD, and ABL, vary in the frequency-specific manner characteristic of each location in the frontal hemisphere.
ILD-alone refers to stimuli in which only ILD was available as a spatial cue and to the neuronal responses to such stimuli. The ABL and ITD were held constant across space.
ILD-frequency tuning refers to the selectivity of a neuron for combinations of frequency and ILD.
Neurophysiology and stimuli. Surgery and neurophysiological recordings in the IC proceeded according to the methods of Keller and Takahashi (1996
) were inserted into the IC at known stereotaxic coordinates. Responses were recorded from isolated units in the IC, typically between 14 and 17 mm below the telencephalic surface. Once isolated, neurons were selected for further analysis on the basis of quantitative criteria for spatial specificity. All of the neurons included in this study exhibited a clear single response peak in the frontal hemisphere that was spatially restricted in both azimuth and elevation. Units that did not meet space-specific requirements were not included for further analysis. That such neurons were indeed in the IC was verified in one owl by inserting an electrode coated with the fluorescent marker DiI (Molecular Probes, Eugene, OR; owl 897) and recovering the electrode track histologically (Euston, 2001
Stimuli consisted of fixed-amplitude, random-phase broadband noise (2–11 kHz) and one-third octave narrowband noise, tone pips, and gamma tones. Noise stimuli had a uniformly distributed amplitude spectrum that deviated by <0.5 dB from the mean amplitude between specified frequencies and rolled off at 0.5 dB/Hz outside this range. Gamma tone filters were fourth-order and 100 Hz wide at half-height. All stimuli were synthesized de novo for each stimulus presentation and used in 3–10 repetitions per trial type.
On isolation, the responses of a unit to ABL, ITD, and ILD were initially characterized using broadband noise. A rate–ITD curve was measured at 5 µsec intervals using broadband noise set to the initial ILD preference. A rate–ILD curve was then determined at 5 dB intervals using broadband noise set to the ITD of the peak in the rate–ITD curve. A rate–ABL curve subsequently was measured at 5 dB intervals using broadband noise set to the ITD and ILD indicated by the previous two tests. All other tests, including those with tones and gamma tones, were done at 10–20 dB above the threshold of the neuron, as determined in the ABL test [30–60 dB sound pressure level (SPL), peak A-weighted]. In all but a few cases, tones, gamma tones, or both were presented at the same overall level as noise stimuli. All stimuli had 5 msec on–off ramps and were 100 msec in duration, with 200 or 300 msec interstimulus intervals. Spontaneous activity was determined during a 100 msec period without stimulation before the onset of each stimulus.
After the initial characterization, the spatial RF of the neuron was measured using fully cued stimuli at 359 locations in the frontal hemisphere. An ILD-alone response was measured using individualized ILD-alone filters, sampling the same locations that were sampled with fully cued stimuli. ITD cues were removed in the frequency domain first. We replaced measured phases for the left-and right-ear filters at all locations with one linear set of phases. This set ITD = 0 for each location. The optimal ITD of a neuron was then imposed in the time domain during unit recording via delay impulse responses. To correct for ABL, it was necessary to perform calculations in the frequency domain (i.e., on the HRTFs). This is because ABL varies in a frequency-dependent manner. Therefore, for each frequency value (e.g., 2–11 kHz, 1024 points in all), we calculated the average ABL across all locations. For each frequency value along the spectrum of each location, we imposed this average ABL so that the spectra maintained the original ILD value at that frequency. Thus, we modified the left and right HRTF spectra for each location so that (1) at each frequency value, the average of the left and right levels was the average ABL across all locations for that frequency; and (2) at each frequency value, subtracting the left level from the right level yielded the correct ILD value for that location and that frequency.
The ILD–frequency response of a neuron was recorded using tone pips, narrow gamma tones, or one-third octave noise at 4–5 dB ILD intervals between ±30 dB in ILD and at 450 Hz frequency intervals between 2 and 11 kHz. Finally, additional ILD-alone tests were performed for validation of the spatially relevant ILD cues (see below).
Identification of frequency-specific ILD tuning from responses to noise. Our approach is based on the idea that if the activity of a neuron is influenced by a particular frequency band in a broadband ILD-alone stimulus, the firing rate at each location would be related to the frequency-specific ILD value that occurs at that location. More precisely, the firing rate at a given location should be proportional to the difference between the ILD for that frequency at that location and the ILD for that frequency at the location that elicited the maximal response, i.e., the "optimal" location of the neuron.
The approach is illustrated in Figure 1. The diamond in Figure 1 A represents the frontal hemisphere in VAS, and the colors represent the firing evoked (red, high; blue, low) by ILD-alone stimuli in one cell. The optimal location of this neuron was at -25° elevation and 15° azimuth. The ILD spectrum giving rise to the activity at this location was defined as the optimal ILD spectrum and is shown as the solid gray line in Figure 1 B. For each of 15 one-sixth octave frequency bands (2–11 kHz), we computed the root mean square difference between the ILD spectrum at the optimal location and that of each of the remaining spatial locations. We appended a negative sign to this difference, if, overall, the ILD spectrum for the location in question was less than that of the optimal location and a positive sign otherwise. The resulting signed quantity is the "ILD distance" (dILD). Had we used the ordinary mean difference instead of dILD, we would have underestimated the true difference in ILD spectra, because the negative and positive differences within a frequency band would have canceled out. In Figure 1 B, the ILD spectra of the optimal location and locations that elicited low (Fig 1 A, white circle, B, blue dotted line) and high (Fig. 1A, black circle, B, magenta dashed line) firing rates are superimposed. The thickened segments of each ILD spectra in Figure 1 B indicate 3 of the 15 one-sixth octave bands (f4, f12, f14) for which we calculated the dILD.
View larger version (58K):
[in this window]
[in a new window]
Figure 1. Identification of frequency-specific ILD tuning from ILD-alone response surfaces measured using broadband noise. ILD-alone responses using broadband stimuli were measured (A) and decomposed by frequency band according to dILD, calculated from the location-dependent ILD spectra. B, Location-dependent ILD spectra are shown for 3 of the 359 locations from which responses were obtained: gray, spectrum for optimal location; magenta, spectrum for another location that elicits a high spike rate; blue, spectrum for a location that elicits a low spike rate. For each frequency band, the firing rate at a given spatial location was plotted against the dILD occurring at that location. Such plots are shown for three representative one-sixth octave bands in C–E. Local regression (C–E, gray lines) yielded a curve that estimated the relationship between firing rate and ILD distance for each frequency band. F, Surface showing the local regression functions for all frequency bands. G, Principal frequency components of the surface shown in F. For clarity, only three principal components are shown. The proportion of variance that each principal component explains is shown in the inset. H–J, Plot shown in F weighted by the three highest principal components. K–M, Principal components regression. The dILD–frequency surfaces weighted by the principal components were transformed into spatial plots via linear spectral integration, and multiple regression was used to compute the weight ( i) for each principal component. This is shown by the weighted summation of K–M, compared against the measured ILD-alone response surface, shown in A and reproduced in N. The principal component matrix (represented in G) was then multiplied by the vector of
Thus, for each frequency band, we computed a set of dILD values that were associated with the locations on the measured ILD-alone surface, resulting in a frequency decomposition in which frequencies were assumed to be independent. At low frequencies, at which the head and facial ruff imposed little sound shadowing, the range of ILDs was limited, whereas at higher frequencies, the range was wider (Keller et al., 1998
If a neuron is responsive to a particular frequency band, its location-dependent firing rate should be correlated with the dILD values of that band. We therefore plotted, for each frequency, the firing rate measured at the various locations in the ILD-alone surface against the associated dILD values. Figure 1C–E shows the plots for bands f4, f12, and f14. The gray lines are the estimated activity versus dILD curves obtained from zero-order local regression, using 15% of the data around each point (i.e.,
= 0.15; Loader, 1999
The contribution of each frequency band (i.e., the frequency band weighting) was determined using principal-component regression (MATLAB princomp and regress functions) as illustrated in Figure 1 F–O. First, the principal frequency components were calculated using the covariance matrix of the activity versus dILD curves (Fig. 1 F). Figure 1G shows the first three principal frequency components, which accounted for 41, 20, and 11% of the variance in the covariance matrix, respectively (inset). All components that uniquely accounted for at least 1% of the total variance in the covariance matrix of the broadband dILD–frequency surface were kept for further processing. For the neuron whose results are shown in Figure 1, the first six principal components met this criterion (although only the first three are shown for clarity). The frequency weights thus derived are closely related to the goodness of fit between spike rate and dILD (Fig. 1C–E). Weighting with the goodness of fit itself, although possible, poses the problem of multicollinearity, in which dILD tuning curves of neighboring frequency bands are correlated. Principal-component regression circumvents this problem.
The unweighted dILD–frequency response surface (Fig. 1 F) was then multiplied by each of the retained principal components to yield a set of component-weighted response surfaces (Fig. 1 H–J) as follows:
where ActivitydILD,frequencyuntransformed is an element in a matrix whose columns are the frequency-specific activity–dILD curves, and PCdILD,frequency is an element in a matrix of identical columns whose rows are the frequency weights of a given principal component. The three principal component-transformed surfaces for neuron N10 are shown in Figure 1 H–J.
(1) Each response surface in Figure 1 H–J was transformed into space using the method described by Euston and Takahashi (2002
Finally, multiple regression analysis was performed:
in which the dependent variable was the measured ILD-alone response surface (Fig. 1 N) and the independent variables were the principal-component estimates of the ILD-alone response surface (Fig. 1 K–M). Only those components that yielded regression coefficients (
(2) Each frequency-specific activity–dILD curve (Fig. 1 F) was multiplied by the appropriate frequency band weight to yield the final dILD–frequency response surface (Fig. 1 P). This surface was transformed into an ILD-alone response surface (Fig. 1Q), which, in essence, was a prediction of the ILD-alone surface measured with noise under the assumption that the space-specific neuron integrated the contribution of each frequency band independently and linearly (Euston and Takahashi, 2002
(3) In Figure 1, the ILD–frequency surface was obtained using the responses from all the locations of the measured ILD-alone response. The ILD–frequency surface may be thought of as a filter that transforms responses in the frequency domain into the spatial domain. The standard method of evaluating the quality of such a filter is to build it with a random subset of the measured locations and to test its output on the locations that were not selected for filter estimation. We took this approach to assess how well the modeled responses predicted the measured data. The square of the Pearson product–moment correlation coefficient (r) between the measured ILD-alone surface and that predicted using the method described above were computed. Mean r 2 values were calculated for 20 repetitions. For each repetition, we randomly selected n% of the locations and tested the predictivity on the remaining (1 n)%. It is important to note that this random procedure is a more robust test of the identified ILD–frequency tuning than, say, always selecting the same locations for parameterization and testing.
Neurophysiological validations. We evaluated the extent to which the frequency bands identified from the response of a neuron to noise were necessary and sufficient to account for the ILD-alone response surface measured with the broadband noise. We therefore constructed "minimal" stimuli consisting only of those spectral components identified from the noise response as being spatially relevant. These bands were defined as those that showed at least half the maximum frequency weight in the ILD–frequency response surface. Complementary stimuli, termed "negative" stimuli, consisted of only those frequency bands between 2 and 11 kHz that were not identified as spatially relevant. Minimal and negative stimuli were used to record ILD-alone response surfaces that were compared with the ILD-alone response surfaced obtained with the full-spectrum stimuli. ILD-alone response surfaces using the minimal and negative stimuli were measured using the same overall ABL level used to record responses to broadband stimuli. When comparing two measured ILD-alone response surfaces, we included all 359 measured locations.
We also compared the effectiveness with which the ILD tuning identified from the tonal stimuli and noise could account for the measured ILD-alone response surface. To this end, we constructed "tone-indicated" stimuli, comprising only those frequencies that elicited more than half-maximal firing in the ILD–frequency response surfaces assessed with the tone pips and gamma tones. Tone-indicated stimuli were analogous to the minimal stimuli just described, except that their spectral components were identified using the more standard method of presenting tone pips and gamma tones of differing center frequencies and ILDs. ILD-alone surfaces were measured using the tone-indicated stimuli, and r2 values were computed to compare them with those generated by broadband noise and by the minimal stimuli.
Comparison of frequency-specific ILD tuning assessed with tones and noise. We compared the ILD–frequency tuning estimated from responses to noise and tones. When using tones, we measured responses at a broad range of ILDs (±30 dB) and frequencies (2–11 kHz), and some ILD values were outside the ILD range encountered naturally. By contrast, when inferring the tuning from responses to noise, response estimates were necessarily restricted to the ethological range. To make the comparisons, therefore, we only considered frequency-specific ILD values common to the two surfaces. The ILD–frequency response surface estimated from the noise response was linearly interpolated in both the frequency and ILD dimensions to match the shared frequency and ILD values. The Pearson r between the two ILD-frequency surfaces was then computed.
We also compared the frequency tuning curves identified with noises and tones in terms of their shape, bandwidth, and peak location. To assess the similarity between simple tuning curves, we first generated rate versus frequency tuning curves by selecting, for each frequency, the maximal firing rate across ILD values. This was done for both the tone-and noise-based ILD–frequency surfaces. The correlation between these curves was then computed. To compare bandwidths, we measured the widths of each curve at half-height, after the curves were linearly interpolated to 1000 samples. To compare the peak positions of the interpolated curves, we calculated a "peak shift ratio," which is the absolute value of the difference in peak frequencies divided by the difference in peak bandwidths:
(4)
Results
Spatial properties of space-specific neurons tested
We analyzed 135 single units, but only 103 were categorized as space-specific neurons and included in the data set analyzed here. Twenty-nine of these were also included in a previous study (Euston and Takahashi, 2002ILD–frequency tuning identified from responses to broadband stimuli
We sought to identify the ILD selectivity at each frequency that could best explain the features of the ILD-alone response surface of a neuron measured with broadband noise. We first determined the fraction of locations of the ILD-alone response surface that was required for robust identification of spatially relevant frequency-specific ILD tuning. Figure 2A plots the fractional variance (r2) between measured and predicted ILD-alone responses against the percentage of loci used to derive the predicted values. The reported r2 values do not include the loci used to derive the prediction. When 30% of the data were used, the prediction accounted for 72% of the variance in the remaining locations. Incorporating up to 90% of the data yielded only a 5% increase.
View larger version (18K):
[in this window]
[in a new window]
Figure 2. Predictive capacity of the frequency-specific ILD tuning identified by responses to noise. All reported r 2 values are calculated using only those locations not included in the parameterization data set. A, Variance in the measured ILD-alone response surface accounted for by the frequency-specific ILD tuning extracted from noise responses plotted against the percent of the total data set (359 locations) used to calculate that tuning. Each point is the median r 2 value across all 102 neurons, calculated from mean values for each neuron obtained using 10 randomly selected fractional data sets. Error bars estimate the variance in the median values and are the SDs around the mean r 2 values obtained across the 10 randomly selected data sets for the 102 neurons. B, Mean ± SD of 20 randomly selected fractional data sets, using 30% of the total number of data points to predict the remaining 70%, plotted versus neuron number. C, Distribution of the mean r 2 values obtained using 30% of the data across 102 neurons.
Figure 2B assesses the robustness of our method of identifying frequency-specific ILD tuning for 102 cells. Each point represents the r2 value of a neuron averaged over the 20 iterations, randomly selecting 30% of the loci each time. The error bars represent the SDs. The SDs are quite restricted, suggesting that the tuning inferred from any set of 108 locations will yield consistent values of r2.
Figure 2C shows the distribution of r2 values for our population of neurons. The value of r2 was ≥0.75 (r ≥ 0.87) for 39% of the neurons analyzed here, and 69% of the neurons had r2 values ≥0.64 (r ≥ 0.8).
Comparison of the predictive capacities of ILD tuning identified from responses to broadband and tonal stimuli
We compared the degree to which the frequency-specific ILD tuning derived from responses to tones and noise could predict the measure ILD-alone response in 97 neurons.For each neuron, we used the tuning estimated from 30% of the responses to broadband ILD-alone stimuli to predict the ILD-alone responses at the remaining 70% of locations. This yielded a median r2 of 0.72 ± 0.14. Sixty-seven percent of these 97 neurons had r2 values ≥0.64. The same integration algorithm using the tuning obtained from responses to tonal stimuli produced a median r2 of 0.51 ± 0.19, consistent with the results of Euston and Takahashi (2002
View larger version (22K):
[in this window]
[in a new window]
Figure 3. Proportion of variance, r 2, in the measured ILD-alone response surfaces explained by frequency-specific ILD tuning inferred from broadband responses compared with r 2 values obtained from tuning measured with tones. The bootstrapping procedure used to infer the tuning from broadband responses was based on 30% of the locations picked at random, with replacement, and averaged over 20 repetitions.
A more detailed comparison of the predictive capacities of ILD–frequency tuning identified from responses to broadband and tonal stimuli is shown in Figure 4 for two representative neurons. The top row depicts results from neuron N14, showing a high level of agreement between the measured ILD-alone response surface and the prediction using noise-indicated ILD tuning (r2 = 0.86). Peaks in the ILD–frequency surfaces obtained from responses to noises and tones matched well (Fig. 4, top, first column). However, the tonal ILD–frequency surface indicated that the neuron was driven by a narrower range of frequencies (6–7 kHz) than was shown by the ILD–frequency tuning identified from responses to noise (5–8 kHz). As a result, there was poor overall agreement (r2 = 0.33) between measured ILD-alone responses and predicted responses using the ILD–frequency tuning based on responses to tones.
View larger version (49K):
[in this window]
[in a new window]
Figure 4. Predictive capacities of frequency-specific ILD tuning obtained from noise and tonal responses of two representative neurons. First column, Frequency-specific ILD tuning obtained from responses to broadband noise (top) and the tuning measured using tone or gamma tone pips (bottom). The neuron number is shown above the ILD–frequency tuning surface determined from responses to noise. Second column, Measured ILD-alone response surfaces to broadband noise. Third column, ILD-alone response surface predicted from the frequency-specific ILD tuning identified from responses to noise, shown in the first column. The reported r 2 above each surface is the mean variance accounted for in the 70% of locations that were not used to infer the frequency-specific ILD tuning. Fourth column, ILD-alone response surface predicted from the tonal response shown in the first column. The r 2 above each of these surfaces is the variance accounted for using the tone-indicated inputs.
Figure 4, bottom row, shows results from neuron N43, for which the ILD–frequency tuning obtained using noise stimuli accounted for 69% of the variance. In the case of N43, both the tonal and broadband ILD–frequency tuning indicated that this neuron relied on a broad frequency range for its spatial specificity. However, the ILD–frequency tuning identified using noise was shifted to higher frequencies and was broader in the ILD dimension than was the tuning estimated with tones. Because the ILD range of the tone-indicated tuning was so restricted, the corresponding predicted ILD-alone surface was actually more restricted than the measured surface (Fig. 4, row 3, compare second, fourth columns). The two estimates predicted more activity in the right hemisphere than seen in the measured ILD-alone surface, yet the noise-indicated ILD–frequency tuning was better at predicting the spatial response pattern of the neuron.
For both neurons N14 and N43, the ILD–frequency surfaces estimated by the two methods showed clear differences in identified spectral components. This suggests that spectral components that elicited maximal spike rates under narrowband conditions were not necessarily the most important in defining the spatial response of a neuron to complex stimuli.
Comparison of ILD–frequency tuning surfaces identified from responses to broadband and tonal stimuli
The ILD–frequency response surfaces that were obtained from responses to broadband and tonal stimuli differed in several key aspects. We grouped these differences into four general classes. For many neurons, examples of more than one class were observed. Thus, it was often impossible to classify neurons as belonging to only one class.In the first class, responses to noise identified a wider frequency bandwidth in the ILD–frequency surface than did responses to tonal stimuli. Overall, the median peak bandwidth of noise-identified tuning was 3410 Hz, whereas the median peak bandwidth of tone-indicated tuning was 2007 Hz, yielding a median difference of 1058 ± 1863 Hz (p < 0.0001, Wilcoxon signed rank test; n = 97). A total of 62 of 97 neurons showed noise-identified tuning that was at least 500 Hz wider than their tonal tuning curves. However, only 37 of 97 neurons showed peak bandwidths in their noise-identified tuning curves that were at least 1 SD wider than the peak bandwidths in their tone-indicated tuning curves. Examples of these, N13 and N14, are shown in Figure 5, A and B, respectively.
View larger version (27K):
[in this window]
[in a new window]
Figure 5. A–D, Comparisons of the frequency-specific ILD tuning obtained with tones (Ton; bottom), one-third octave noises (1/3; middle), or broadband noise (BB; top). The r 2 values quantifying the similarity of the measured ILD alone surfaces and those predicted using the noise and tonal responses are shown above each set of plots.
Some neurons (24 of 97) showed a narrower frequency range in their noise-identified tuning curves compared with their tone-indicated curves (median difference, 510 ± 1167 Hz; p < 0.0001, Wilcoxon signed rank test). This was somewhat surprising, given that space-specific neurons are generally characterized as broadband (Knudsen and Konishi, 1978a
A third class of difference was a noticeable frequency shift in the peaks of the ILD–frequency tuning surfaces. Of our 97 neurons, 56 exhibited peak frequency shifts that were at least as great as 50% of the difference in peak bandwidths in their tuning curves (peak shift ratio, ≥0.5; see Materials and Methods, Eq. 4). Of these, 30 shifted to a higher frequency in the noise-identified tuning curve, and 26 shifted to a lower frequency. One such case is seen with N10 in Figure 5C (peak shift ratio, 1.03).
Finally, for many neurons (68 of 97), tonal stimuli elicited stronger activity at frequencies of >9.5 kHz than did broadband noise (p < 0.0001, Wilcoxon signed rank test; n = 68). Only 21 of 97 neurons showed stronger high-frequency activity with noise. In neuron N8 (Fig. 5D), for example, tonal stimuli at
10 kHz elicited strong responses, producing a "high-frequency tail" not indicated by responses to noise. Note that in Figure 5, C and D, the high-frequency tails obtained using tonal stimuli were markedly reduced when one-third octave noises were used instead of tonal stimuli.
Despite these four classes of difference between the tonal and broadband ILD–frequency surfaces, the surfaces did correspond, often overlapping considerably in frequency even when they differed in ILD. The median correlation between the tonal frequency tuning curves and the frequency tuning obtained from responses to broadband stimuli was 0.59 ± 0.38. The median correlation between the two types of ILD–frequency surfaces was 0.52 ± 0.24. The distribution of correlation coefficients between surfaces is shown in Figure 6A. Of the 97 neurons for which both tonal- and noise-identified ILD–frequency surfaces were available, only 28% showed correlations between ILD–frequency surfaces of >0.60.
View larger version (19K):
[in this window]
[in a new window]
Figure 6. Comparisons of the frequency-specific ILD tuning obtained from responses to tones and broadband stimuli. A, Distribution of correlation coefficients between the ILD–frequency surfaces obtained by the two methods. B, Scatter plot showing that the predictive capacity of tone-indicated ILD–frequency tuning is better when it is more similar to the tuning obtained with broadband stimuli. C, Scatter plot showing that there is no relationship between the predictive capacity of the noise-identified ILD-frequency tuning and its similarity to the tone-indicated tuning.
As another test of whether the estimate of ILD tuning identified using the noise method more accurately estimated the response of a neuron as stimulus bandwidth increased, we also measured frequency-specific ILD tuning from 13 neurons using one-third octave noise. For most of the neurons tested in this manner (n = 13), the response pattern using one-third octave noise resembled the tuning identified from responses to broadband noise more closely than did the tuning identified using tonal stimuli. Correlations between tuning curves identified using broadband noise and those obtained with tones or one-third octave noise increased from a median r = 0.55 (tones) to a median r = 0.83 (one-third octave noise; p = 0.0012, Wilcoxon signed rank test). Correlations between ILD–frequency tuning surfaces increased from a median r = 0.47 (broadband noise vs tones) to a median r = 0.67 (broadband noise vs one-third octave noise; p < 0.0001, Wilcoxon signed rank test).
An important question was whether the ability of tonal ILD–frequency tuning surfaces to predict the measured ILD-alone response surfaces was related to their similarity to the tuning identified with noise. If the ILD–frequency tuning identified from responses to noise accurately reflected the tuning contributing to the measured ILD-alone responses of a neuron, then the closer a tonal ILD–frequency surface came to the tuning obtained from broadband stimuli, the better the tonal result should have predicted the measured ILD-alone surface. Thus, we expected a correlation between the quality of the tonal ILD-alone prediction and the agreement between the ILD–frequency tuning surfaces identified using tonal and broadband stimuli. However, we did not expect an association between the predictivity of ILD–frequency tuning surfaces obtained using noise and how closely these surfaces resembled the tonal ILD–frequency tuning surfaces.
This was, in fact, observed (Fig. 6B,C). In Figure 6B, the r2 value of the ILD-alone prediction obtained using tonal tuning is plotted for each neuron (n = 97) against the correlation between the tonal and noise-identified ILD–frequency surfaces. There was a clear correlation (r = 0.34; p = 0.0006), suggesting that the more the tonal ILD–frequency surface resembled the tuning obtained from responses to noise, the better the tonal tuning was at predicting the broadband ILD-alone response surface. There was no correlation seen between the predictive capacity of the noise-identified ILD–frequency tuning and the agreement between the ILD–frequency surfaces (Fig. 6C; r = 0.12; p > 0.1).
Direct neurophysiological assessment of frequency tuning obtained from responses to broadband noise
Tests with minimal, negative, and tone-indicated stimuli were used to assess the degree to which the spatially relevant frequency bands identified for a neuron could elicit ILD-alone responses similar to those obtained with broadband stimuli.ILD-alone response surfaces using minimal stimuli were recorded from 13 neurons. The median r2 value, calculated using all 359 locations, between the measured broadband and measured minimal ILD-alone response surfaces was 0.50 ± 0.15. This indicates good agreement with responses to broadband stimuli, keeping in mind that the trial-to-trial variance in firing can diminish the apparent agreement between the two measured ILD-alone surfaces.
Of these 13 neurons, 10 also were tested using negative stimuli, which yielded a median r2 value of 0.086 ± 0.20. Direct comparison of the r2 values obtained using the minimal and negative stimuli revealed that the former showed significantly higher predictive capacity than the latter (p = 0.002, Wilcoxon signed rank test). Thus, those frequency bands identified as spatially relevant were necessary and sufficient to produce ILD-alone responses similar to those obtained with broadband stimuli, whereas the remaining frequency bands, when presented alone, did not elicit such similar responses.
For 9 of 13 neurons, we compared the effects of minimal and tone-indicated stimuli. In these 9 cases, the ILD-frequency tuning surfaces obtained with tonal and broadband stimuli overlapped considerably, with a correlation coefficient of 0.75, a value significantly higher than the median value of 0.59 across all 97 neurons. Because of this significant overlap in frequencies and the subsequent similarity in the minimal and tone-indicated narrowband stimuli, we did not expect the minimal stimuli to be significantly better than tonal stimuli in eliciting ILD-alone surfaces similar to broadband ILD-alone surfaces. In fact, the median r2 value for the minimal stimuli was 0.50, whereas that for the tonal stimuli was 0.42, but this difference did not reach statistical significance (p > 0.1, Wilcoxon signed rank test). However, the ILD-frequency tuning identified from responses to broadband noise yielded better predictions of the measured minimal ILD-alone response surfaces than did the tone-indicated tuning. The median r2 value between minimal ILD-alone surface and that predicted using the ILD–frequency tuning identified from responses to broadband noise was 0.66, whereas that between the measured minimal and tonally predicted surfaces was 0.52 (p = 0.033; Wilcoxon signed rank test). Thus, even when there was significant overlap between tone- and noise-identified frequency bands, the noise-identified ILD–frequency tuning was better at predicting the ILD-alone response surface measured with minimal stimuli than was the tone-indicated ILD–frequency tuning. There was no difference between noise- and tone-indicated tuning in predicting measured tone-indicated ILD-alone response surfaces (r2 = 0.54 for each; p > 0.1). This crossvalidation test showed that the ILD–frequency tuning surfaces identified from responses to broadband stimuli yielded the tuning that most strongly determined the spatial response of a neuron.
In three of the nine neurons for which ILD-alone responses were obtained with both minimal and tone-indicated stimuli, the ILD–frequency tuning differed considerably, each showing a correlation coefficient ≤0.59 (N69, N72, and N135). In these neurons, the noise-identified frequency bands were better than the tone-indicated bands at eliciting ILD-alone response surfaces similar to those measured with broadband stimuli.
Results from one of these neurons, N69, are shown in Figure 7. Here, the ILD–frequency tuning surfaces obtained from responses to noise (Fig. 7A, top) and tonal (Fig. 7A, bottom) stimuli exhibited a marked lack of overlap (r = 0.59), attributable to a peak shift along the frequency axis. The ILD-alone response surface obtained using the minimal stimulus showed a better overall match to the broadband ILD-alone surface (Fig. 7C; r2 = 0.31) compared with the ILD-alone response surface obtained using the tone-indicated stimulus (Fig. 7D; r2 = 0.09).
View larger version (38K):
[in this window]
[in a new window]
Figure 7. Contribution of frequency bands identified by tonal and broadband stimuli, shown for neuron N69.A, Tuning for ILD and frequency obtained with broadband stimuli (BB; top) and tonal stimuli (Ton; bottom). B, ILD-alone response surface measured with broadband stimuli. C, ILD-alone surface measured using the minimal stimulus (see Materials and Methods for details). D, ILD-alone surface measured using a tone-indicated stimulus (see Materials and Methods for details). E, Firing elicited at each location by the minimal stimulus plotted against the firing at each corresponding location in the ILD alone surface obtained with the broadband stimulus (B). F, Firing elicited at each location by the tone-indicated stimulus plotted against the firing at each corresponding location in the ILD-alone surface obtained with the broadband stimulus (B). The variances, r2, obtained from the scatter plots in E and F are shown above the ILD-alone surfaces in C and D, respectively.
Similar results were obtained for the other two neurons. Both N135 and N72 exhibited significant disagreement between tone-indicated and noise-identified ILD–frequency tuning surfaces (r =-0.41 and 0.31, respectively) although nonetheless sharing a major frequency band. In neuron N135, the frequency bands identified with noise elicited an ILD-alone response surface more similar to that observed with broadband stimuli (r2 = 0.51) compared with those identified with tonal stimuli (r2 = 0.39). In neuron N72, the minimal stimulus also gave a response pattern (r2 = 0.29) more similar than did the tone-indicated stimulus (r2 = 0.23).
Discussion
The spatial RFs of auditory neurons are determined by the cues and frequencies to which neurons are tuned and the distribution of the cues in space. We showed that a linear model for the integration of ILD across frequency works well for most neurons but only if the frequency-specific ILD tuning is correctly estimated. We also showed that this tuning can be reliably identified from the neuronal responses to noise. Our approach was robust for noises 10–20 dB above threshold, in that only 30% of the measured responses picked at random were sufficient to explain much of the ILD-alone response surface of a neuron.We also presented data suggesting that narrowband stimuli, such as tones, are not optimal in assessing the frequency-specific ILD tuning of space-specific neurons in the barn owl (Euston and Takahashi, 2002
Saturation effects occur because of processing within the effective bandwidth of auditory nerve fibers (Greenwood, 1991
Comparisons with earlier studies
A number of recent studies have manipulated HRTFs to study the individual contributions of sound localization cues. Delgutte et al. (1999In the owl's space map, Euston and Takahashi (2002
The spatial transformation used here and by Euston and Takahashi (2002
Analogies to other input–output models of neural coding
The current method of extracting frequency-specific spatial information from responses to broadband stimuli is reminiscent of reverse correlation, in which spectral properties of auditory neurons are estimated from responses to broadband stimuli on the basis of the timing of individual spikes (deBoer and Kuyper, 1968Schnupp et al. (2001
Implications for developmental plasticity of neural mechanisms for sound localization
Studies of developmental plasticity in the neural coding of sound localization rely on altering frequency-dependent cues during an animal's critical developmental period and assessing the subsequent effects on neural activity (Szczepaniak and Moller, 1996Therefore, it is still possible that responses to broadband noise from space-specific neurons of device-reared owls would yield a correspondence between the frequency-specific ILD effects of the acoustic device and spatially relevant ILD–frequency response surfaces. This could help advance studies of the mechanisms of plasticity in ILD tuning responsible for neural encoding of auditory space.
Footnotes
Received Oct. 10, 2002; revised Feb. 25, 2003; accepted Mar. 18, 2003.This work was supported by National Institute of General Medical Sciences GrantT32-GM07257 (M.L.S.), National Institute on Deafness and Other Communication Disorders Grant DC-03925, and National Science Foundation Grant LIS CMS9720334. We thank Dr. Kip Keller for assistance in measuring head-related transfer functions, Dr. David Euston for providing some of the data and for programming assistance, Dr. Klaus Hartung for assistance with virtual auditory space techniques, Dr. Kip Keller, Dr. David Euston, Dr. Matthew Spitzer, Dr. Avinash Bala, Dr. Richard Marrocco, Dr. Patrick Phillips, Dr. Helen Neville, Dr. Janis Weeks, and Dr. Ben Arthur for helpful discussions, and anonymous reviewers for beneficial comments.
Correspondence should be addressed to Michael L. Spezio, Center for Neuroscience, 1544 Newton Court, University of California, Davis, CA 95616. E-mail: mlspezio{at}ucdavis.edu.
Copyright © 2003 Society for Neuroscience 0270-6474/03/234677-12$15.00/0
References
Arthur B (2001) Neural computations leading to space-specific auditory responses in the barn owl. PhD dissertation, California Institute of Technology.
Brainard MS, Knudsen EI, Esterly SD (1992) Neural derivation of sound source location: resolution of spatial ambiguities in binaural cues. J Acoust Soc Am 91: 1015–1027.[Web of Science][Medline]
Calhoun BM, Miller RL, Wong JC, Young ED (1998) Rate encoding of stimulus spectra by auditory nerve fibers. In: Psychophysical and physiological advances in hearing (Palmer AR, Rees A, Summerfield Q, Meddis R, eds), pp 170–177. London: Whurr.
Carney LH, Friedman M (1998) Spatiotemporal tuning of low-frequency cells in the anteroventral cochlear nucleus. J Neurosci 18: 1096–1104.
[Abstract/Free Full Text] deBoer E, Kuyper P (1968) Triggered correlation. IEEE Trans Biomed Eng 15: 169–179.[Medline]
deCharms RC, Blake DT, Merzenich MM (1998) Optimizing sound features for cortical neurons. Science 280: 1439–1443.
[Abstract/Free Full Text] Delgutte B, Joris PX, Litovsky RY, Yin TC (1999) Receptive fields and binaural interactions for virtual-space stimuli in the cat inferior colliculus. J Neurophysiol 81: 2833–2851.
[Abstract/Free Full Text] Eggermont JJ, Aertsen AM, Johannesma PIM (1983a) Quantitative characterization procedure for auditory neurons based on the spectro-temporal receptive field. Hear Res 10: 167–190.[Web of Science][Medline]
Eggermont JJ, Aertsen AM, Johannesma PIM (1983b) Prediction of the responses of auditory neurons in the midbrain of the grass frog based on the spectro-temporal receptive field. Hear Res 10: 191–202.[Web of Science][Medline]
Eggermont JJ, Johannesma PIM, Aertsen AM (1983c) Reverse-correlation methods in auditory research. Quant Rev Biophys 16: 341–414.
Euston D (2001) From spectrum to space: the integration of frequency specific intensity cues to produce auditory spatial receptive fields in the barn owl inferior colliculus. PhD dissertation, University of Oregon.
Euston D, Takahashi TT (2002) From spectrum to space: the contribution of level difference cues to spatial receptive fields in the barn owl inferior colliculus. J Neurosci 22: 284–293.
[Abstract/Free Full Text] Fitzpatrick DC, Kuwada S (2001) Tuning to interaural time differences across frequency. J Neurosci 21: 4844–4851.
[Abstract/Free Full Text] Fuzessery ZM, Wenstrup JJ, Pollak GD (1985) A representation of horizontal sound location in the inferior colliculus of the mustache bat (Pteronotus p. parnellii). Hear Res 20: 85–89.[Web of Science][Medline]
Gao E, Suga N (2000) Experience-dependent plasticity in the auditory cortex and the inferior colliculus of bats: role of the corticofugal system. Proc Natl Acad Sci USA 97: 8081–8086.
[Abstract/Free Full Text] Gold JI, Knudsen EI (1999) Hearing impairment induces frequency-specific adjustments in auditory spatial tuning in the optic tectum of young owls. J Neurophysiol 82: 2197–2209.
[Abstract/Free Full Text] Gold JI, Knudsen EI (2000a) Abnormal auditory experience induces frequency-specific adjustments in unit tuning for binaural localization cues in the optic tectum of juvenile owls. J Neurosci 20: 862–877.
[Abstract/Free Full Text] Gold JI, Knudsen EI (2000b) A site of auditory experience-dependent plasticity in the neural representation of auditory space in the barn owl's inferior colliculus. J Neurosci 20: 3469–3486.
[Abstract/Free Full Text] Gold JI, Knudsen EI (2001) Adaptive adjustment of connectivity in the inferior colliculus revealed by focal pharmacological inactivation. J Neurophysiol 85: 1575–1584.
[Abstract/Free Full Text] Greenwood DD (1991) Critical bandwidth and consonance: their operational definitions in relation to cochlear nonlinearity and combination tones. Hear Res 54: 209–246.[Medline]
Irvine DR, Gago G (1990) Binaural interaction in high-frequency neurons in inferior colliculus of the cat: effects of variations in sound pressure level on sensitivity to interaural intensity differences. J Neurophysiol 63: 570–591.
[Abstract/Free Full Text] Jiang D, McAlpine D, Palmer AR (1997) Responses of neurons in the inferior colliculus to binaural masking level difference stimuli measured by rate-versus-level functions. J Neurophysiol 77: 3085–3106.
[Abstract/Free Full Text] Keller CH, Takahashi TT (1996) Binaural cross-correlation predicts the responses of neurons in the owl's auditory space map under conditions simulating summing localization. J Neurosci 16: 4300–4309.
[Abstract/Free Full Text] Keller CH, Takahashi TT (2000) Representation of temporal features of complex sounds by the discharge patterns of neurons in the owl's inferior colliculus. J Neurophysiol 84: 2638–2650.
[Abstract/Free Full Text] Keller CH, Hartung K, Takahashi TT (1998) Head-related transfer functions of the barn owl: measurement and neural responses. Hear Res 118: 13–34.[Web of Science][Medline]
Knipschild M, Dorrscheidt GJ, Rubsamen R (1992) Setting complex tasks to single units in the avian auditory forebrain. I: processing of complex artificial stimuli. Hear Res 57: 216–230.[Web of Science][Medline]
Knudsen EI (1982) Auditory and visual maps of space in the optic tectum of the owl. J Neurosci 2: 1177–1194.[Abstract]
Knudsen EI (1999) Mechanisms of experience-dependent plasticity in the auditory localization pathway of the barn owl. J Comp Physiol [A] 185: 305–321.[Web of Science][Medline]
Knudsen EI, Konishi M (1978a) A neural map of auditory space in the owl. Science 200: 795–797.
[Abstract/Free Full Text] Knudsen EI, Konishi M (1978b) Center-surround organization of auditory receptive fields in the owl. Science 202: 778–780.
[Abstract/Free Full Text] Knudsen EI, Konishi M (1978c) Space and frequency are represented separately in auditory midbrain of the owl. J Neurophysiol 41: 870–884.
[Abstract/Free Full Text] Knudsen EI, Zheng W, DeBello WM (2000) Traces of learning in the auditory localization pathway. Proc Natl Acad Sci USA 97: 11815–11820.
[Abstract/Free Full Text] Köppl C (1997) Frequency tuning and spontaneous activity in the auditory nerve and cochlear nucleus magnocellularis of the barn owl Tyto alba. J Neurophysiol 77: 364–377.
[Abstract/Free Full Text] Köppl C, Yates G (1999) Coding of sound pressure level in the barn owl's auditory nerve. J Neurosci 19: 9674–9686.
[Abstract/Free Full Text] Kuwada S, Yin TC (1983) Binaural interaction in low-frequency neurons in inferior colliculus of the cat. I. Effects of long interaural delays, intensity, and repetition rate on interaural delay function. J Neurophysiol 50: 981–999.
[Abstract/Free Full Text] Loader C (1999) Local regression and likelihood. New York: Springer.
Marmarelis PZ, Naka K (1972) White-noise analysis of a neuron chain: an application of the Wiener theory. Science 175: 1276–1278.
[Abstract/Free Full Text] Marmarelis PZ, Naka KI (1973a) Nonlinear analysis and synthesis of receptive-field responses in the catfish retina. I. Horizontal cell leads to ganglion cell chain. J Neurophysiol 36: 605–618.
[Free Full Text] Marmarelis PZ, Naka KI (1973b) Nonlinear analysis and synthesis of receptive-field responses in the catfish retina. II. One-input white-noise analysis. J Neurophysiol 36: 619–633.
[Free Full Text] Marmarelis PZ, Naka KI (1973c) Nonlinear analysis and synthesis of receptive-field responses in the catfish retina. III. Two-input white-noise analysis. J Neurophysiol 36: 634–648.
[Free Full Text] Moiseff A, Konishi M (1981) Neuronal and behavioral sensitivity to binaural time differences in the owl. J Neurosci 1: 40–48.[Abstract]
Pena JL, Konishi M (2001) Auditory spatial receptive fields created by multiplication. Science 292: 249–252.
[Abstract/Free Full Text] Schafer M, Rubsamen R, Dorrscheidt GJ, Knipschild M (1992) Setting complex tasks to single units in the avian auditory forebrain. II: do we really need natural stimuli to describe neuronal response characteristics? Hear Res 57: 231–244.[Web of Science][Medline]
Schnupp JW, Mrsic-Flogel TD, King AJ (2001) Linear processing of spatial cues in primary auditory cortex. Nature 414: 200–204.[Medline]
Shinn-Cunningham B (2001) Models of plasticity in spatial auditory processing. Audiol Neurootol 6: 187–191.[Medline]
Szczepaniak WS, Moller AR (1996) Evidence of neuronal plasticity within the inferior colliculus after noise exposure: a study of evoked potentials in the rat. Electroencephalogr Clin Neurophysiol 100: 158–164.[Medline]
Takahashi T, Konishi M (1986) Selectivity for interaural time difference in the owl's midbrain. J Neurosci 6: 3413–3422.[Abstract]
Takahashi TT (1989) The neural coding of auditory space. J Exp Biol 146: 307–322.
[Abstract/Free Full Text] Theunissen FE, Sen K, Doupe AJ (2000) Spectral-temporal receptive fields of nonlinear auditory neurons obtained using natural sounds. J Neurosci 20: 2315–2331.
[Abstract/Free Full Text] Theunissen FE, David SV, Singh NC, Hsu A, Vinje WE, Gallant JL (2001) Estimating spatio-temporal receptive fields of auditory and visual neurons from their responses to natural stimuli. Network 12: 289–316.[Web of Science][Medline]
Tollin DJ, Yin TC (2002a) The coding of spatial location by single units in the lateral superior olive of the cat. I. Spatial receptive fields in azimuth. J Neurosci 22: 1454–1467.
[Abstract/Free Full Text] Tollin DJ, Yin TC (2002b) The coding of spatial location by single units in the lateral superior olive of the cat. II. The determinants of spatial receptive fields in azimuth. J Neurosci 22: 1468–1479.
[Abstract/Free Full Text] Wagner H, Takahashi T, Konishi M (1987) Representation of interaural time difference in the central nucleus of the barn owl's inferior colliculus. J Neurosci 7: 3105–3116.[Abstract]
Wang J, Salvi RJ, Powers N (1996) Plasticity of response properties of inferior colliculus neurons following acute cochlear damage. J Neurophysiol 75: 171–183.
[Abstract/Free Full Text] Wickesberg RE, Dickson JW, Gibson MM, Geisler CD (1984) Wiener kernel analysis of responses from anteroventral cochlear nucleus neurons. Hear Res 14: 155–174.[Web of Science][Medline]
Xu L, Furukawa S, Middlebrooks JC (1999) Auditory cortical responses in the cat to sounds that produce spatial illusions. Nature 399: 688–691.[Medline]
Yamada WM, Lewis ER (1999) Predicting the temporal responses of nonphase-locking bullfrog auditory units to complex acoustic waveforms. Hear Res 130: 155–170.[Web of Science][Medline]
Yin TC, Kuwada S (1983) Binaural interaction in low-frequency neurons in inferior colliculus of the cat. II. Effects of changing rate and direction of interaural phase. J Neurophysiol 50: 1000–1019.
[Abstract/Free Full Text] Yin TC, Kuwada S, Sujaku Y (1984) Interaural time sensitivity of highfrequency neurons in the inferior colliculus. J Acoust Soc Am 76: 1401–1410.[Web of Science][Medline]
Young ED, Calhoun BM, Yu JJ, Nelken I (1998) Representation of complex stimuli in the peripheral auditory system. In: International Conference on Acoustics/Acoustical Society of America, pp 1373–1376.
This article has been cited by other articles:
B. J. Fischer, G. B. Christianson, and J. L. Pena
Cross-Correlation in the Auditory Coincidence Detectors of Owls
J. Neurosci., August 6, 2008; 28(32): 8107 - 8115.
[Abstract] [Full Text] [PDF]
C. H. Keller and T. T. Takahashi
Localization and Identification of Concurrent Sounds in the Owl's Auditory Space Map
J. Neurosci., November 9, 2005; 25(45): 10446 - 10461.
[Abstract] [Full Text] [PDF]
M. W. Spitzer, A. D. S. Bala, and T. T. Takahashi
A Neuronal Correlate of the Precedence Effect Is Associated With Spatial Selectivity in the Barn Owl's Auditory Midbrain
J Neurophysiol, October 1, 2004; 92(4): 2051 - 2070.
[Abstract] [Full Text] [PDF]
This Article Abstract 
Full Text (PDF) Submit an eLetter Alert me when this article is cited Alert me when eLetters are posted Alert me if a correction is posted Citation Map Services Email this article to a friend Similar articles in this journal Similar articles in ISI Web of Science Similar articles in PubMed Alert me to new issues of the journal Download to citation manager 
Citing Articles Citing Articles via HighWire Citing Articles via ISI Web of Science (5) Citing Articles via Google Scholar Google Scholar Articles by Spezio, M. L. Articles by Takahashi, T. T. Search for Related Content PubMed PubMed Citation Articles by Spezio, M. L. Articles by Takahashi, T. T.
|
|
|
|
 |
|