pro (version 0.1.1)

pro: Fit the PRO model

Description

Fit the PRO model to data. Reference: X Luo, S Gee, V Sohal, D Small (In Press). A Point-process Response Model for Optogenetics Experiments on Neural Circuits. _Statistics in Medicine_.

Usage

pro(spike, flash, ...)

Arguments

spike

A binary vector represents spiking (1) or no spiking (0).

flash

A binary vector of the same length of spike, 1 for flashing and 0 for non-flashing.

...

Additional parameters, see model.pro.

Value

a glm object of the fitted PRO coefficients.

Examples

n <- 500
set.seed(100)
re <- sim.lif(n, rbinom(n, 1, 0.14), 7, 3)
fit.pro <- pro(re$sbin, re$I)
summary(fit.pro)