LZeroSpikeInference (version 1.0.1)

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

Description

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

Usage

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

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)

See Also

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

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

Simulation: simulateAR1, plot.simdata.

Examples

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

Run the code above in your browser using DataCamp Workspace