Gammatone stimuli.

<p><b>A</b>) Human cochleotopic map of center frequencies from 100 to 11,234 Hz. <b>B</b>) Gammatone signal at 1 kHz center frequency. <b>C</b>) Gammatone filter bank with 50 channels.</p

Discovery-based selection of lambda (<i>λ</i>).

<p>Discovery-selection by permutation of the response <i>r</i>(<i>t</i>) to identify the minimum lambda necessary to reduce the degrees of freedom to zero. Empirical distribution shows the histogram of minimum <i>λ</i> over 200 repeated permutations and serves as the null distribution. The red arrow denotes the median of the distribution that we designate as the optimal <i>λ</i>.</p

STRF estimated from STA and from GLM.

<p><b>(A</b>) STRF from STA. Stimulus center frequencies ranged from 97 to 11,234 Hz. Spike-times were binned at 1 ms resolution. (<b>B</b>) Evolution of Spike-count STRF from GLM as a function of increasing <i>λ</i> for L1 norm LASSO: (<b>a</b>) very low values of <i>λ</i> lead to noisy estimates, (<b>c</b>) at very high values of <i>λ</i> all covariates are zero valued, (<b>d</b>) at optimal value of <i>λ</i> chosen from discovery-based selection. (<b>C</b>) Evolution of Spike-count STRF from GLM as a function of increasing <i>λ</i> for L1/L2 norm group LASSO: (<b>e</b>) very low, (<b>g</b>) very high, and (<b>h</b>) optimal values of <i>λ</i>. Optimization selects or removes, simultaneously, all the covariates forming a group. Groups are composed of 4x4, adjacent and non-overlapping covariates. <b>(D)</b> Predicted discharge <i>λ</i>_<i>CIF</i>(<i>t</i>|<i>H</i>_<i>t</i>) from representative segment of gammatone stimulus, with and without sparse-group regularization. Group-sparse regularized GLM consistently improved the prediction of validation data over non-regularized GLM prediction of expected spike-counts. Correlation coefficients are 0.133 with regularization (red) and 0.066 without regularization (gray). Neural responses from S178, electrode contact #4.</p

Spike-Count History contribution to spike-count and High-<i>γ</i> band power STRFs.

<p>Neural responses from S178, electrode contact #11. Anatomical location shown in <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0137915#pone.0137915.g001" target="_blank">Fig 1</a>. (<b>A</b>) Spike-count STRF estimated with L1/L2 sparsity-inducing norm and group structure GLM with Poisson distribution to link responses to predictors, without and with 350 ms spike-count history. (<b>B</b>) Magnitude of GLM spike-count history coefficients decreases with increasing history (i.e. time elapsed since current spike-count prediction). Shading represents 95% central range of null distribution estimated from permuted random shuffling of responses. (<b>C</b>) Cumulative distribution of p-values testing the contribution of spike-count history to current spike-count activity driven by gammatone stimuli. All p-values were adjusted for false discovery rate. (<b>D</b>) High-<i>γ</i> power STRF estimated with L1/L2 sparsity-inducing norm and group structure GLM with Poisson distribution to link responses to predictors, without and with 350 ms spike-count history. (<b>E</b>) Magnitude of GLM coefficients decreases with increasing spike-count history (i.e. time elapsed since current High-<i>γ</i> band power prediction). (<b>F</b>) Cumulative distribution of p-values testing the contribution of spike-count history to current High-<i>γ</i> power activity driven by gammatone stimuli. All p-values were adjusted for false discovery rate.</p

Locations of electrode recording sites within the superior temporal plane.

<p><b>(A</b>) MRI lateral-view rendering of a typical human left hemisphere. The Sylvian fissure is not visible from the cortical surface. The superior temporal plane was revealed along a section oriented at an oblique horizontal plane (solid red line with razor blade inset). (<b>B</b>) MRI rendering of superior temporal plane viewed from superior aspect. Light blue shading denotes the location of the obliquely oriented Heschl’s gyrus. The estimated locations of four recording sites selected from three different subjects (S140, S151, and S178) were projected to the surface of this illustrative brain and marked with filled red circles. MRI cross-sectional images containing the recording sites were obtained from sections oriented at an oblique frontal plane (solid green lines with razor blade inset), approximately perpendicular to the long axis of Heschl’s gyrus. (<b>C</b>) Line drawings of MRI cross sections show the position of the recording sites within the grey matter of Heschl’s gyrus for individual subjects.</p

Best frequency and bandwidth estimated for sites within Heschl’s gyrus.

<p><b>(Left column</b>) Fourteen recording locations in each subject were projected to the surface and marked with open circles. Solid red line marks the transition from core to non-core fields estimated with click train stimuli. (<b>Middle column</b>) Open circles mark a single or multiple best-frequency (BF) value for each location estimated from spike-count STRF. Bandwidth (BW) is depicted by a solid black line centered on the best-frequency for each location. (<b>Right column</b>) BF and BW mapped using High-<i>γ</i> power STRF.</p

Spike-count and High-<i>γ</i> power STRFs derived with sparse GLM models.

<p>Neural responses from S178 electrode contact #4 used in left column, responses from S140 electrode contact #6 used in middle column and from S151 electrode contact #14 used in right column (See also <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0137915#pone.0137915.g001" target="_blank">Fig 1</a>). Optimal <i>λ</i> values are shown on insets. (<b>A</b>) Spike-count STRF using L1 sparsity-inducing norm. (<b>B</b>) Spike-count STRF using L1/L2 norm regularization that exploits group structure when covariates are partitioned into neighborhoods, or groups. In this case, optimization selects or removes all the variables forming a group. Groups are composed of 4x4, adjacent and non-overlapping covariates. (<b>C</b>) High-<i>γ</i> (70 to 150 Hz) band power STRF from L1/L2 norm regularization that exploits group structure. Groups are composed of 4x4, adjacent and non-overlapping covariates.</p