calcBurstSzFn: Calculate burst size function

Description Usage Arguments Details Value Examples

Description

Calculates burst size function Phi(tau) according to Eq.4 of Dobrzynski & Bruggeman (2009) Elongation dynamics shape bursty transcription and translation. PNAS, 106(8) Formula: Phi (tau) = total no. of intervals / number of intervals longer than tau

Usage

1
calcBurstSzFn(in.v, in.nthr = 100, in.sc = 1/2)

Arguments

in.v

Vector with intervals

in.nthr

Total number of thresholds

in.sc

Scaling factor for threshold interval. Provide as 1/x to obtain the maximum threshold size = x.

Details

The input should contain a single-column data with intervals between events such as firings of a neuron, or production of a protein. The function searches for the number of events longer than a threshold. The procedure is repeated for a series of threshold lengths. The output data frame contains two columns with the number of events for each of the thresholds.

Value

Data frame with threshold lengths and the value of burst size function.

Examples

1
2
3
4
# generate a series of exponentially distributed intervals
df = rexp(1000, 1)
phi=calcBurstSzFn(v$V1, 100, 1/2)
plot(phi$t, phi$b, type = 'l')

dmattek/tca-package documentation built on Dec. 10, 2019, 4:10 p.m.