To obtain a cell's contrast response functions, mean spike rate was calculated for each stimulus type and each contrast. Mean spike rate for the size III stimulus was computed as the average spike rate across all 30 trials at that contrast. We counted spikes for a 40-ms observation window beginning after a latency of 15 ms after stimulus onset (for on-cells) or offset (for off-cells). The average number of spikes over the 40-ms window was expressed as firing rate (impulses/s). Spike firing rates did not increase when the window was made longer than 40 ms, consistent with a psychophysical critical duration of no more than 40 ms for luminance increments.
35 Mean spike rate for the frequency-doubling stimulus was computed as the amplitude of the fundamental component from Fourier analysis of averaged response histograms. We only considered first harmonic amplitude, which is dominant compared with second or higher harmonics for contrast levels near threshold. A model for central detection of sinusoidal stimuli using a window of fixed time width
36 yielded results similar to those presented here.
For each cell and stimulus, contrast response functions (i.e., mean spike rate, number of impulses per second minus baseline firing versus contrast), were fit with a Michaelis-Menten function:
where
R is spike rate at contrast
c minus baseline firing rate,
G is contrast gain (impulses per second per percentage contrast), and
R max is maximum spike rate minus baseline firing rate. At low contrasts, the spike rate
R increases linearly with contrast
c with a slope of
G, at high contrasts
R asymptotes to
R max, and when contrast equals the
semisaturation contrast (
R max/
G),
R reaches half of
R max. Some of the P cells showed only small amounts of saturation, which means that
R max is poorly constrained by the data and can become unreasonably high (over 1000 spikes/s). Therefore, a maximum of 500 spikes/s was set for
R max; for some poorly responsive P cells, this maximum was reached because the semisaturation contrast was higher than the maximum stimulus. In these cases, fits were repeated with
R max = 300 spikes/s, and in all cases
G did not change by more than 0.2 log unit.
For each cell and stimulus, the contrast response function was fitted by minimizing the χ2 value using the Levenberg-Marquardt algorithm, with graphing and data analysis software (Igor Pro, version 6.1.2; Wavemetrics, Portland, OR). The fitted contrast gains (in log units) were then used for statistical analysis. For descriptive analysis, the means and associated standard deviations are reported. Given that some cells were measured with both stimulus types (n = 25; 19 M cells and 6 P cells) while other cells were tested with only one stimulus type, we used a mixed linear model to account for potential correlation among repeated measurements on the same cells. The mixed linear model tested the effects of stimulus type, cell type, and their interaction on log contrast gain. A significant interaction between cell type and stimulus type would indicate a difference between size III and frequency-doubling stimuli in differentiating M cell versus P cell responses. We first conducted mixed linear modeling with all cell data and then with data only from those cells that were measured with both stimulus types. The differences in contrast gains and associated standard errors (SE) estimated from the mixed linear model are reported. The statistical analyses were conducted using statistical software (SAS 9.1.3; SAS Institute, Cary, NC).