simulate_ar1 {FastLZeroSpikeInference} | R Documentation |
Simulate fluorescence trace based on simple AR(1) generative model
simulate_ar1(n, gam, poisMean, sd, seed)
n |
number of timesteps |
gam |
AR(1) decay rate |
poisMean |
mean for Poisson distributed spikes |
sd |
standard deviation |
seed |
random seed |
Simulate fluorescence trace based on simple AR(1) generative model
y_t = c_t + eps, eps ~ N(0, sd)
c_t = gam * c_t-1 + s_t
s_t ~ Pois(poisMean)
spikes, fluorescence, and calcium concentration
sim <- simulate_ar1(n = 500, gam = 0.998, poisMean = 0.009, sd = 0.05, seed = 1) plot(sim)