balci2019 {YEAB}R Documentation

Peak individual trial analysis using moving average

Description

Peak individual trial analysis using moving average

Usage

balci2019(tasa_norm, bins)

Arguments

tasa_norm

numeric, normalized response rate

bins

numeric

Details

Based on Balci et al 2010

Value

a list with params: a numeric vector with start, stop, spread and argmax (the bin at which response rate is max) mov_av: the moving average

Examples

data("r_times")
# binarize r_times to create response rate at 2 sec bins
bins <- get_bins(r_times, 0, 180, 2)
bin_res <- 6
tasa <- f_table(bins, 0, 180, bin_res)
tasa_norm <- tasa$prop / max(tasa$prop)
bins <- tasa$bins
balci_ind <- balci2019(tasa_norm, bins)

plot(bins, tasa_norm, xlab = "6 sec bins", )
lines(bins, balci_ind$mov_av, col = "blue", lwd = 2)
abline(v = balci_ind$params[c(1, 2, 4)], lwd = c(1, 1, 2), col = c(1, 1, "red4"))

[Package YEAB version 1.0.6 Index]