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

In FastLZeroSpikeInference: Fast Nonconvex Deconvolution of Calcium Imaging Data

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)

Example output

FastLZeroSpikeInference documentation built on May 2, 2019, 4:02 p.m.