simulate_ar1 {FastLZeroSpikeInference}R Documentation

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

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)

Value

spikes, fluorescence, and calcium concentration

Examples

sim <- simulate_ar1(n = 500, gam = 0.998, poisMean = 0.009, sd = 0.05, seed = 1)
plot(sim)

[Package FastLZeroSpikeInference version 2018.12.10 Index]