Adaptively optimizing experiments can significantly reduce the number of trials needed to characterize neural responses using parametric statistical models. However, the potential for these methods has been limited to date by severe computational challenges: choosing the stimulus which will provide the most information about the (typically high-dimensional) model parameters requires evaluating a high-dimensional integration and optimization in near-real time. Here we present a fast algorithm for choosing the optimal (most informative) stimulus based on a Fisher approximation of the Shannon information and specialized numerical linear algebra techniques. This algorithm requires only low-rank matrix manipulations and a one-dimensional linesearch to choose the stimulus and is therefore efficient even for high-dimensional stimulus and parameter spaces; for example, we require just 15 milliseconds on a desktop computer to optimize a 100-dimensional stimulus. Our algorithm therefore makes real-time adaptive experimental design feasible. Simulation results show that model parameters can be estimated much more efficiently using these adaptive techniques than by using random (nonadaptive) stimuli. Finally, we generalize the algorithm to efficiently handle both fast adaptation due to spike-history effects and slow, non-systematic drifts in the model parameters.
[1]
Eric R. Ziegel,et al.
Generalized Linear Models
,
2002,
Technometrics.
[2]
L. Cosmides.
From : The Cognitive Neurosciences
,
1995
.
[3]
Liam Paninski,et al.
Asymptotic Theory of Information-Theoretic Experimental Design
,
2005,
Neural Computation.
[4]
M. V. Rossum,et al.
In Neural Computation
,
2022
.
[5]
E. Bizzi,et al.
The Cognitive Neurosciences
,
1996
.
[6]
S. Eisenstat,et al.
A Stable and Efficient Algorithm for the Rank-One Modification of the Symmetric Eigenproblem
,
1994,
SIAM J. Matrix Anal. Appl..
[7]
D. Wilkin,et al.
Neuron
,
2001,
Brain Research.