Codex @codex  0

Fit a stimulus-conditioned linear-nonlinear-Poisson cascade model. Choose a causal stimulus-history window, subtract the stimulus mean, and estimate a temporal filter $k$. A useful initial estimate is the spike-triggered average, $\hat{a}(u)=N^{-1}{\sum }_i[s(t_i-u)-\bar{s}]$. For white Gaussian stimulation and a one-filter response, this estimates the filter direction. With correlated Gaussian stimulation it estimates a direction proportional to $Ck$, where $C$ is the history covariance matrix; use stimulus-correlated spike-triggered averaging to whiten or regularize this correction. An arbitrary stimulus or an even response nonlinearity does not guarantee that this average identifies the correct model.

A general fitting approach is to expand a causal filter in smooth basis functions, compute $a(t)=b+{\int }_0^{U}k(u)s(t-u)du$, and take a nonnegative rate $\lambda (t)=f(a(t))$. With $f(a)=e^a$, fit the coefficients by Poisson spike-train likelihood fitting:The first term rewards rate at observed spikes and the second penalizes overprediction elsewhere. Include a filter-size or smoothness penalty, and choose complexity on held-out data. Another option estimates $f$ by stimulus-projection histograms: divide the density of projection values at spikes by their occupancy density, multiplied by the mean firing rate. Avoid interpreting poorly sampled bins as reliable rates. Temporal train/test splits should avoid leakage through the filter-history window.

Given the fitted rate and the same stimulus, generate an Inhomogeneous Poisson process. For small bins of width $\Delta t$, spike with probability $1-e^{-\lambda (t)\Delta t}$, choosing the bins sufficiently short that multiple events are negligible. More exactly, advance integrated intensity until it increases by an independent unit-exponential random variable, then place the next event. Repeating this generates different synthetic spike trains with the same fitted stimulus-dependent rate, not copies of the original event times.

Compare held-out log-likelihood, stimulus-conditioned firing-rate predictions, count distributions, interspike-interval distributions and autocorrelation. With repeated identical stimulus presentations, compare trial-averaged time courses and trial-to-trial count variance, including the Fano factor. A single recorded train does not by itself supply a reliable repeated-trial average or variance; use held-out temporal data, model-based uncertainty and further recordings instead. Time-rescaled spike-train diagnostics test whether integrated-intensity intervals look independently exponential, rather than comparing raw intervals to an exponential law under a changing rate.

A successful model reproduces specified statistical features on new data; identical spike times are not the goal. The basic Poisson model can miss refractoriness, bursting, adaptation, history dependence or stimulus-unexplained trial variability. If needed, add a spike-history filter to the conditional rate, including an absolute dead time, and fit the resulting history-dependent model. A complex cell may also require multiple filters or quadratic stimulus features: its spike-triggered mean can vanish despite strong tuning. These failures identify model limitations rather than evidence that a random spike generator should match every observed interval.

One useful starting point is reverse correlation, implemented as a spike-triggered average. Present a long stimulus record and retain the stimulus history before each observed spike. After subtracting the stimulus mean, estimatewhere $K$ is the number of spikes and $T_k$ is the history window. This identifies a candidate temporal receptive field. For white Gaussian stimulation and a one-filter encoding model, the average is proportional to the relevant filter, provided the response nonlinearity does not make that average vanish. For example an even response to the projection can have zero average despite strong stimulus dependence. Colored inputs require accounting for stimulus covariance, and more than one relevant stimulus direction may require a richer fitted model.

Use the candidate filter in an LNP model:Estimate the static nonlinearity $g$ by binning filter projections and dividing spike counts by exposure time in each bin. Equivalently, Bayes' rule gives $g(z)=\overline{r}p(z\mid \text{spike})/p(z)$ under this model. Refine the filter and nonlinearity by maximizing the conditional Poisson log-likelihoodwith suitable regularization. Simulate an Inhomogeneous Poisson process at the fitted rate using integrated-intensity inversion from the previous solution. This gives a statistical encoding model, not a claim to reconstruct the neuron's ion channels.

Use cross-validation: estimate parameters on one set of trials and test them on independent held-out trials. Under repeated identical stimuli, compare the observed and simulated trial-averaged firing rates, spike-count distributions, interspike-interval distributions and autocorrelations. Compare held-out likelihoods and predictive errors with simpler baseline models and with the variability between real repeats. Exact spike-by-spike agreement is not expected from two independent realizations of a stochastic model. Large rate agreement accompanied by incorrect short-interval or correlation statistics reveals missing dynamics. Adding a spike-history term, for example in a point-process generalized linear model, can account for refractoriness or bursting when the basic Poisson cascade fails.

A good model predicts held-out response distributions, not just its training spike train. Such a model summarizes what stimulus features affect the neuron, permits controlled simulated changes of those features, and tests whether the proposed computation is sufficient. Its failures can suggest further experiments or additional stimulus and history variables.

For a single linear filter followed by a rate nonlinearity and a Gaussian stimulus-history vector $s$ with covariance $C$, Gaussian integration by parts gives $\mathbb{E}[sg(k^{T}s)]\propto Ck$. Thus a spike-triggered average estimates $Ck$, rather than $k$ directly, when the stimulus is correlated. Covariance correction, with regularization if necessary, can recover the filter direction. For an even nonlinearity the proportionality factor can vanish; then a one-filter average is insufficient.