Define the intervals associated with <code>y = j</code> based on the flooring function.
The function returns <code>-Inf</code> for <code>j = 0</code> (or smaller) and <code>Inf</code> for
any <code>j &gt;= y_max + 1</code>, where <code>y_max</code> is a known upper bound on the data <code>y</code>
(if specified).

For Bayesian and classical inference and prediction with count-valued data,
Simultaneous Transformation and Rounding (STAR) Models provide a flexible, interpretable,
and easy-to-use approach. STAR models the observed count data using a rounded
continuous data model and incorporates a transformation for greater flexibility.
Implicitly, STAR formalizes the commonly-applied yet incoherent procedure of
(i) transforming count-valued data and subsequently
(ii) modeling the transformed data using Gaussian models.
STAR is well-defined for count-valued data, which is reflected in predictive accuracy,
and is designed to account for zero-inflation, bounded or censored data, and over- or underdispersion.
Importantly, STAR is easy to combine with existing MCMC or point estimation
methods for continuous data, which allows seamless adaptation of continuous data
models (such as linear regressions, additive models, BART, random forests,
and gradient boosting machines) for count-valued data. The package also includes several
methods for modeling count time series data, namely via warped Dynamic Linear Models.
For more details and background on these methodologies, see the works of
Kowal and Canale (2020) <doi:10.1214/20-EJS1707>,
Kowal and Wu (2022) <doi:10.1111/biom.13617>,
King and Kowal (2022) <arXiv:2110.14790>, and
Kowal and Wu (2023) <arXiv:2110.12316>.

Brian King

countSTAR

Flexible Modeling of Count Data

Dan Kowal

a_j function

<dl><dt>j</dt>
<dd>the integer-valued input(s)</dd>
<dt>y_max</dt>
<dd>a fixed and known upper bound for all observations; default is <code>Inf</code></dd></dl>

Arguments

Inverse rounding function — a_j

<dl>

<dt>j</dt>
<dd>the integer-valued input(s)</dd>


<dt>y_max</dt>
<dd>a fixed and known upper bound for all observations; default is <code>Inf</code></dd>

</dl>

a_j: Inverse rounding function

Description

Usage

Value

Arguments

Examples