1 / 55

What is the neural code?

What is the neural code?. Reading out the neural code. Alan Litke, UCSD. Reading out the neural code. Reading out the neural code. Kwan, HFSP 2010. What is the neural code?. M. Berry. What is the neural code?. Single neurons. Populations. Encoding and decoding.

ranit
Download Presentation

What is the neural code?

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. What is the neural code?

  2. Reading out the neural code Alan Litke, UCSD

  3. Reading out the neural code

  4. Reading out the neural code Kwan, HFSP 2010

  5. What is the neural code? M. Berry

  6. What is the neural code? Single neurons Populations

  7. Encoding and decoding • Encoding: how does a stimulus cause a pattern of responses? • what are the responses and what are their characteristics? • neural models: • -- from stimulus to response • -- descriptive mechanistic models • Decoding: what do these responses tell us about the stimulus? • Implies some kind of decoding algorithm • How to evaluate how good our algorithm is?

  8. Basic coding model: linear filtering Spatial filter: r=f(x,y) I(x0-x,y0-y)dxdy

  9. Basic coding model: temporal filtering Linear filter: r(t) = f(t-t) s(t) dt …shortcomings?

  10. Nonlinearity: tuning curves Nonlinear function: r= g(s) Gaussian tuning curve of a cortical (V1) neuron

  11. Nonlinearity: tuning curves Hand reaching direction Nonlinear function: r= g(s) Cosine tuning curve of a motor cortical neuron

  12. Nonlinearity: tuning curves Hand reaching direction Nonlinear function: r= g(s) Sigmoidaltuning curve of a V1 cortical neuron

  13. Nonlinearity: tuning curves Quian Quiroga, Reddy, Kreiman, Koch and Fried, Nature (2005)

  14. What is s?

  15. Nonlinearity: tuning curves Quian Quiroga, Reddy, Kreiman, Koch and Fried, Nature (2005)

  16. What is s?

  17. Coding models encoding P(response | stimulus) P(stimulus | response) decoding • What is response? • What is stimulus? • What is the function P?

  18. Basic coding model: temporal filtering Linear filter: r(t) = f(t-t) s(t) dt

  19. Next most basic coding model s*f1 Linear filter & nonlinearity: r(t) = g( f(t-t) s(t) dt)

  20. How to find the components of this model s*f1

  21. Reverse correlation: the spike-triggered average Spike-conditional ensemble

  22. The spike-triggered average

  23. Dimensionality reduction More generally, one can conceive of the action of the neuron or neural system as selecting a low dimensional subset of its inputs. Start with a very high dimensional description (eg. an image or a time-varying waveform) and pick out a small set of relevant dimensions. P(response | stimulus)  P(response | s1, s2, .., sn) s2 s3 s(t) s1 s3 sn s4 s5 s.. s.. s.. s2 s1 S(t) = (S1,S2,S3,…,Sn)

  24. Linear filtering Linear filtering = convolution = projection s2 s3 s1 Stimulus feature is a vector in a high-dimensional stimulus space

  25. Determining linear features from white noise Gaussian prior stimulus distribution Spike-conditional distribution Spike-triggered average

  26. How to find the components of this model s*f1

  27. Determining the nonlinear input/output function The input/output function is: which can be derived from data using Bayes’ rule:

  28. Nonlinear input/output function Tuning curve: P(spike|s) = P(s|spike) P(spike) / P(s) Tuning curve: P(spike|s) = P(s|spike) P(spike) / P(s) P(s) P(s|spike) P(s) P(s|spike) s s

  29. Next most basic coding model s*f1 Linear filter & nonlinearity: r(t) = g( f(t-t) s(t) dt) …shortcomings?

  30. Less basic coding models Linear filters & nonlinearity: r(t) = g(f1*s, f2*s, …, fn*s)

  31. covariance Determining linear features from white noise Gaussian prior stimulus distribution Spike-conditional distribution Spike-triggered average

  32. Identifying multiple features The covariance matrix is Stimulus prior Spike-triggered stimulus correlation Spike-triggered average • Properties: • The number of eigenvalues significantly different from zero • is the number of relevant stimulus features • The corresponding eigenvectors are the relevant features • (or span the relevant subspace) Bialek et al., 1988; Brenner et al., 2000; Bialek and de Ruyter, 2005

  33. A toy example: a filter-and-fire model Let’s develop some intuition for how this works: a filter-and-fire threshold-crossing model with AHP Keat, Reinagel, Reid and Meister, Predicting every spike. Neuron (2001) • Spiking is controlled by a single filter, f(t) • Spikes happen generally on an upward threshold crossing of • the filtered stimulus •  expect 2 relevant features, the filter f(t) and its time derivative f’(t)

  34. STA f(t) f'(t) Covariance analysis of a filter-and-fire model

  35. Let’s try it Example: rat somatosensory (barrel) cortex Ras Petersen and Mathew Diamond, SISSA Record from single units in barrel cortex

  36. Normalisedvelocity Pre-spike time (ms) White noise analysis in barrel cortex Spike-triggered average:

  37. White noise analysis in barrel cortex Is the neuron simply not very responsive to a white noise stimulus?

  38. Covariance matrices from barrel cortical neurons Prior Spike- triggered Difference

  39. 0.4 0.3 0.2 0.1 Velocity 0 -0.1 -0.2 -0.3 -0.4 150 100 50 0 Pre-spike time (ms) Eigenspectrum from barrel cortical neurons Eigenspectrum Leading modes

  40. Input/output relations from barrel cortical neurons Input/output relations wrt first two filters, alone: and in quadrature:

  41. 0.4 0.3 0.2 Velocity (arbitrary units) 0.1 0 -0.1 -0.2 -0.3 150 100 50 0 Pre-spike time (ms) Less significant eigenmodes from barrel cortical neurons How about the other modes? Next pair with +veeigenvalues Pair with -veeigenvalues

  42. Negative eigenmode pair Input/output relations for negative pair Firing rate decreases with increasing projection: suppressive modes

  43. Salamander retinal ganglion cells perform a variety of computations Michael Berry

  44. Almost filter-and-fire-like

  45. Not a threshold-crossing neuron

  46. Complex cell like

  47. Bimodal: two separate features are encoded as either/or

  48. When have you done a good job? • When the tuning curve over your variable is interesting. • How to quantify interesting?

  49. When have you done a good job? Tuning curve: P(spike|s) = P(s|spike) P(spike) / P(s) Tuning curve: P(spike|s) = P(s|spike) P(spike) / P(s) Boring: spikes unrelated to stimulus Interesting: spikes are selective P(s) P(s|spike) P(s) P(s|spike) s s Goodness measure: DKL(P(s|spike) | P(s))

  50. Maximally informative dimensions Sharpee, Rust and Bialek, Neural Computation, 2004 Choose filter in order to maximize DKL between spike-conditional and prior distributions Equivalent to maximizing mutual information between stimulus and spike P(s) P(s|spike) Does not depend on white noise inputs Likely to be most appropriate for deriving models from natural stimuli s = I*f

More Related