CORRECTIONS

**Geisler, W. S., Najemnik, J., & Ing, A. D. (2009). Optimal stimulus encoders for natural tasks.**

*CORRECTIONS TO:**Journal of Vision, 9*(13):17, 1–16, http://journalofvision.org/9/13/17/, doi:10.1167/9.13.17.

There is a mathematical error in our publication entitled “Optimal stimulus encoders for natural tasks” (

*Journal of Vision, 9*(13):17, 1–16). The error is in text Equation 10 and 1, where Equation 10 is derived. The error has no affect on the first example application concerning “image patch identification,” and only a minor effect on the second example application concerning “foreground identification.” Nonetheless, under some circumstances, the mathematical error might have more substantial consequences. Here we provide replacements for text on p. 5 concerning Equation 10 and for 1.**Replacement text on p. 5:**

Under the above assumptions and expanding where

*p*(*k*∣**r**_{ q}(*k, l*)) using Bayes rule (see 1) we have:$p ( k | r q ( k , l ) )= 1 Z \u2211 j = 1 n k ( \u220f t = 1 q \sigma t ( k , j ) ) \u2212 1exp [ \u2212 1 2 \u2211 t = 1 q [ r t ( k , l ) \u2212 r t ( k , j ) ] 2 \sigma t ( k , j ) 2 ]$

(10)

*n*_{ k}is the number of training samples from category*k,*and*Z*is a normalization factor. In keeping with the approximation in Equation 4, the logarithm of this formula gives the average relative entropy when the stimulus is**s**(*k, l*). Thus, Equations 5– 10 provide a closed-form expression for the average relative entropy of the posterior probability distribution (that the ideal observer computes) for arbitrary samples from the joint probability distribution of environmental categories and associated stimuli,*p*_{0}(*k, l*).To estimate the optimal linear receptive fields we use a ‘greedy’ procedure. In other words, neurons are added to the population one at a time, with each neuron's receptive field being selected to produce the biggest decrease in decoding error. Specifically, we proceed sequentially by first finding the encoding function

*r*_{1}(*k, l*) that minimizes$ D \u2015$

_{1}(see Equation 5); then we substitute the estimated*r*_{1}(*k, l*) into Equation 10 and find the encoding function*r*_{2}(*k, l*) that minimizes$ D \u2015$

_{2}; then we substitute the estimated*r*_{1}(*k, l*) and*r*_{2}(*k, l*) into Equation 10, and find the encoding function*r*_{3}(*k, l*) that minimizes$ D \u2015$

_{3}, and so on.**Replacement for Appendix A:**

Appendix A

Here we derive formulas for the posterior probability distribution that is computed by the ideal Bayesian observer when receiving a population response

**R**_{ q}(*k, l*) to a presentation of stimulus**s**(*k, l*). (Keep in mind that the ideal observer does not know that the stimulus is**s**(*k, l*), but does know the mean response of each neuron in the population to each stimulus in the training set.) According to Bayes' rule:$p ( x | R q ( k , l ) )= p ( R q ( k , l ) | x ) p ( x ) \u2211 i = 1 m p ( R q ( k , l ) | i ) p ( i )$

To derive text Equation 10 we expand the above equation using the definition of conditional probability:

$p ( x | R q ( k , l ) )= p ( x ) \u2211 j = 1 n x p ( R q ( k , l ) | x , j ) p ( j | x ) \u2211 i = 1 m p ( i ) \u2211 j = 1 n i p ( R q ( k , l ) | i , j ) p ( j | i )$

Assuming that the samples are representative of the natural world, the prior probability of a category is the fraction of training samples from the category, and the prior probability of a particular sample from a category is the inverse of the number of training samples within the category:

*p*(*i*) =*n*_{ i}/*n, p*(*j*∣*i*) = 1/*n*_{ i}. Thus,$p ( x | R q ( k , l ) )= \u2211 j = 1 n x \u220f t = 1 q p ( R t ( k , l ) | x , j ) \u2211 i = 1 m \u2211 j = 1 n i \u220f t = 1 q p ( R t ( k , l ) | i , j )$

(A1)

By substitution of Equation 8 we have:

$p ( x | R q ( k , l ) )= \u2211 j = 1 n x ( \u220f t = 1 q \sigma t ( x , j ) ) \u2212 1 exp [ \u2212 1 2 \u2211 t = 1 q [ R t ( k , l ) \u2212 r t ( x , j ) ] 2 \sigma t ( x , j ) 2 ] \u2211 i = 1 m \u2211 j = 1 n i ( \u220f t = 1 q \sigma t ( i , j ) ) \u2212 1 exp [ \u2212 1 2 \u2211 t = 1 q [ R t ( k , l ) \u2212 r t ( i , j ) ] 2 \sigma t ( i , j ) 2 ]$

(A2)

Equation 10 follows by substitution from Equation 4. Note that

*Z*in Equation 10 corresponds to the denominator of Equation A2. In other words, the complete equation is$p ( k | r q ( k , l ) )= \u2211 j = 1 n k ( \u220f t = 1 q \sigma t ( k , j ) ) \u2212 1 exp [ \u2212 1 2 \u2211 t = 1 q [ r t ( k , l ) \u2212 r t ( k , j ) ] 2 \sigma t ( k , j ) 2 ] \u2211 i = 1 m \u2211 j = 1 n i ( \u220f t = 1 q \sigma t ( i , j ) ) \u2212 1 exp [ \u2212 1 2 \u2211 t = 1 q [ r t ( k , l ) \u2212 r t ( i , j ) ] 2 \sigma t ( i , j ) 2 ]$

(A3)

Citation

Geisler W. S.
Najemnik J.
Ing A. D.
(2010). Corrections to: Optimal stimulus encoders for natural tasks. Journal of Vision, 10, (2):27, 2, http://journalofvision.org/10/2/27/, doi:10.1167/10.2.27. [CrossRef]