Simulate fluorescence trace based on a simple AR-1 generative model
simulate_ar1.RdSimulate fluorescence trace based on a simple AR-1 generative model
simulate_ar1(n, gam, poisMean, sd, seed, c0 = 0)
Arguments
| n | Numeric; length of the time series |
|---|---|
| gam | Numeric; AR-1 decay rate |
| poisMean | Numeric; mean for Poisson distributed spikes |
| sd | Numeric; standard deviation |
| seed | Numeric; random seed |
| c0 | Numeric; initial calcium concentration, default to 0 |
Value
spikes-- A list of timesteps at which a spike occursfl-- The noisy fluorescence \(y_t\)conc-- The true calcium concentration \(c_t\)
Details
Simulate fluorescence trace based on a simple AR-1 generative model: $$y_t = c_t + \epsilon_t, \epsilon_t \sim N(0, \sigma^2),$$ $$c_t = gam * c_{t-1} + s_t,$$ $$s_t \sim Poisson(poisMean).$$
