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

Description Usage Arguments Details Value See Also Examples

Description

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

Usage

1
simulateAR1(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

See Also

Estimate spikes: estimateSpikes, print.estimatedSpikes, plot.estimatedSpikes.

Cross validation: cv.estimateSpikes, print.cvSpike, plot.cvSpike.

Simulation: simulateAR1, plot.simdata.

Examples

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

jewellsean/LZeroSpikeInference documentation built on May 19, 2019, 7:16 a.m.