FastLZeroSpikeInference (version 2018.12.10)

simulate_ar1: Simulate fluorescence trace based on simple AR(1) generative model

Description

Simulate fluorescence trace based on simple AR(1) generative model

Usage

simulate_ar1(n, gam, poisMean, sd, seed)

Arguments

n

number of timesteps

gam

AR(1) decay rate

poisMean

mean for Poisson distributed spikes

sd

standard deviation

seed

random seed

Value

spikes, fluorescence, and calcium concentration

Details

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)

Examples

# NOT RUN {
sim <- simulate_ar1(n = 500, gam = 0.998, poisMean = 0.009, sd = 0.05, seed = 1)
plot(sim)
# }