We propose a method that would allow for a rigorous statistical analysis of neural responses to natural stimuli, which are non–Gaussian and exhibit strong correlations. We have in mind a model in which neurons are selective for a small number of stimulus dimensions out of the high dimensional stimulus space, but within this subspace the responses can be arbitrarily nonlinear. Existing analysis ...