Generating distributions

The ARPobservation package generates behavior streams that follow an alternating renewal process with specified generating distributions. The package provides two approaches to generating the initial interim time and initial event duration, which we explain further below. Subsequent event durations A₂, A₃, A₄, ... are generated independently, following a specified probability distribution with mean μ and cumulative distribution function (cdf) F(t; μ). Subsequent interim times B₁, B₂, B₃, ... are generated independently, following a specified probability distribution with mean λ and cdf G(t; λ).

The package currently includes functions for exponential distributions, gamma distributions, mixtures of two gamma distributions, Weibull distributions, uniform distributions, and constant values. Each distribution is implemented as an object of class eq_dist, which provides functions for generating random deviates from the specified distribution and the corresponding equilibrium distribution. For distributions involving more than a single parameter, all parameters except for the mean must be specified.

Initial conditions

ARPobservation provides two approaches to generating the initial interim time and initial event duration. The first approach involves the following steps:

1. Generate a random number X from a Bernoulli distribution with user-specified probability p₀.
2. If X = 0, then generate B₀ from the same distribution as subsequent interim times, i.e., from G(t; λ). If X = 1, then set B₀ = 0.
3. Generate A₁ from the same distribution as subsequent event durations, i.e., from F(t; μ).

If p₀ is not specified by the user, the default value of p₀ = 0 is used, so that behavior streams always begin with an interim time. This approach produces behavior streams that are initially out of equilibrium.

The other approach uses initial conditions chosen so that the resulting process is in equilibrium. This involves the following steps:

1. Generate a random number X from a Bernoulli distribution with probability μ/(μ + λ).
2. If X = 0, then generate B₀ from the distribution with cdf $$ \tilde{G}(t; \lambda) = \frac{1}{\lambda} \int_0^t \left[1 - G(x; \lambda)\right] dx. $$ If X = 1, then set B₀ = 0.
3. If X = 0, then generate A₁ from F(t; μ). If X = 1, then generate A₁ from the distribution with cdf $$ \tilde{F}(t; \mu) = \frac{1}{\mu} \int_0^t \left[1 - F(x; \mu) \right] dx. $$

Event counting

Event counting produces a measurement Y^E equal to the number of events that begin during the observation session. Let J denote the number of last behavioral event that begins during the observation session, which can be calculated by finding the integer that satisfies the inequalities $$ \sum_{j=0}^{J-1} \left(A_j + B_j \right) \leq L < \sum_{j=0}^{J} \left(A_j + B_j \right), $$ where we define A₀ = 0 for notational convenience. It follows that Y^E = J.

Continuous recording

Continuous recording produces a measurement Y^C equal to the proportion of the observation session during which the behavior occurs. In order to calculate this quantity from the behavior stream, we must account for the possibility that the last event beginning during the observation session may have a duration that extends beyond when the session ends. The measurement based on continuous recording can be calculated as $$ Y^C = \begin{cases} \frac{1}{L} \sum_{j=1}^J A_j & \text{if}\quad \sum_{j=1}^{J} \left(B_{j-1} + A_j\right) \leq L \\ 1 - \frac{1}{L} \sum_{j=0}^{J-1} B_j & \text{if}\quad \sum_{j=1}^{J} \left(B_{j-1} + A_j\right) > L \end{cases} $$

Momentary time sampling

In momentary time sampling, an observer divides the observation session into K intervals of equal length and notes whether the behavior is present or absent at the very end of each interval. The summary measurement Y^M then corresponds to the proportion of moments during which the behavior is observed. Let X_k = 1 if the behavior is occurring at the end of interval k for k = 1, ..., K. The value of X_k can be calculated from the behavior stream as follows. Let I(X) denote the indicator function, equal to one if condition X is true and zero otherwise. Let m_k be the number of the last event that ends before the k^th interval ends, defined formally as $$ m_k = \sum_{i=1}^J I\left[\sum_{j=1}^i \left(B_{j-1} + A_j\right) < k L \right] $$ for k = 1, ..., K. If interim time B_mk concludes before the end of interval k (or equivalently, if event A_mk + 1 begins before the end of interval k), then X_k = 1; formally, $$ X_k = I\left[\sum_{j=0}^{m_k} \left(A_j + B_j\right) < k L \right] $$ for k = 1, ..., K. The summary measurement is then calculated as $\displaystyle{Y^M = \sum_{k=1}^K X_k / K}$.

Partial interval recording

Like momentary time sampling, partial interval recording is also based on a set of K intervals of equal length, but a different rule is used to score each interval. In partial interval recording, the observer counts a behavior as present if it occurs at any point during the first c time units of the interval, where c ≤ L/K; the remaining L/K − c time units are used to record notes or rest. Let U_k = 1 if the behavior occurs at any point during the k^th interval, U_k = 0 otherwise. The k^th interval will be equal to one if and only if interim time m_k − 1 ends during the active part of the interval. Noting that interim time m_k − 1 ends at time $\sum_{j=0}^{m_{k-1}} \left(A_j + B_j\right)$ and that the active part of the k^th interval ends at time (k − 1)L + c, it can be seen that $$ U_k = I \left[\sum_{j=0}^{m_{k-1}} \left(A_j + B_j\right) < (k-1)L + c \right], $$ for k = 1, ..., K. The summary measurement Y^P is then calculated as the proportion of intervals during which the behavior is observed at any point: $\displaystyle{Y^P = \sum_{k=1}^K U_k / K}$.

Whole interval recording

Whole interval recording is similar to partial interval recording but uses yet a different rule to score each interval. Specifically, the observer counts a behavior as present only if it occurs for all c time units at the beginning of the interval. Let W_k = 1 if the behavior occurs for the duration, with W_k = 0 otherwise. Let n_k be the number of the last event that begins before the k^th interval begins, defined formally as $$ n_k = \sum_{i=1}^J I\left[\sum_{j=0}^i \left(A_j + B_j\right) < (k - 1) L \right] $$ for k = 1, ..., K. It follows that W_k will be equal to one if and only if event n_k ends after the active part of interval k: $$ W_k = I \left[\sum_{j=1}^{n_k} \left(B_{j-1} + A_j\right) \geq (k - 1) L + c \right], $$ for k = 1, ..., K. The summary measurement Y^W is then calculated as the proportion of intervals during which the behavior is observed at any point: $\displaystyle{Y^W = \sum_{k=1}^K W_k / K}$.