Old name: smithbrown2003-example1
At:
https://github.com/sohale/smithbrown2003-example1
The actual project development was done later on here
Statistical modelling of Spike Trains as Point Processes
Note Λ(t) = ∫ λ(t) dt is the times-rescaling function. x[k] is the input (i.e. the "state"). In physiological terms, the input signal is the EPSP/IPSP.
As spike rate estimator: The derived statistical model (with parameters) can also be seen as a rigorous method of calculating the "spiking rates" based on observed spike trains. The estimated moodel generates λ(t), which is an estimation of the probability of spikes (also can be seen as spiking rate). It's a quantity that is difficult to estimate, or would be incorrect if estimated in the naïve way simply by averaging (over bins) or smoothing of spike trains.
The estimations are based on EM (Expectation Maximisation) (the procedure devised by Ghahramani & Hinton 1996) apllied to Point Processes.
x(t+δt) = A x(t) + αI + ε
λ = exp(βx + μ)
It is comprised of a State Space model ( x = state) and a Point Process ( λ = CIF )
┌───┐ ┌───┐ ┌──────────────────────────────..─┐
│ x │ → │ λ │ → │ │ │║ │║│ │ │ │║║ │ .. │
└───┘ └───┘ └──────────────────────────────..─┘
A State-State model is a Dynamical System, i.e. a system with a feelback loop.
┌───┐ ┌───┐ ┌──────────────────────────────··─┐
│ x │ ───→ │ λ │───→ │ │ │║ │║│ │ │ │║║ │ ·· │
└╥──┘══>╗ └───┘ └──────────────────────────────··─┘
║ ║*A
║ ╟─<─ α * I ←─ I
╚═+══<═╝
+ε
I is the input (driving stimulus) and ε is additive noise, both are injected into the loop.
┌───┐ ┌───┐ ┌──────────────────────────────···
│ x │ ───────→ │ λ │───→ │ │ │║ │║│ │ │ │║║ │ ···
└╥──┘══>╗ └───┘ └──────────────────────────────···
║ ║
I ─→*α──→╢ ╟*A
╚═+══<═╝
+ε
LaTex:
$x_{t+δt} = A^{δt} \ x_t + αI + ε$ $λ = \exp(βx + μ)$
Coloured:
`
color test . color test .
Misc:
$\log λ = βx + μ$ - λ ∝
$\Pr$ (spike) / dt
In fact,
- Pr(spike) ∝ ∫ λ dt
- Pr(spike|no spike before)=Pr(ISI) ∝ ∫ λ dt
.
Using yuml diagrams
Didn't work: ![d]("http://yuml.me/diagram/scruffy/class/[State]uses -.->[Spikes {bg:orange}]")