Simulate from the zero-inflated continuous parameter binomial — rzicpb • underdisp

Simulate from the zero-inflated continuous parameter binomial

Usage

rzicpb(n, lambda, alpha, pi, support = c("ceiling", "floor"))

Arguments

n: Number of units.
lambda: Intensity rate (the CPB's lambda); scalar or length-n.
alpha: CPB shape parameter in (0, 1).
pi: Structural-zero probability; scalar or length-n.
support: The CPB's support, "ceiling" (King's CPB, the default) or "floor"; see cpb().

Value

An integer vector: a structural zero with probability pi, otherwise an (untruncated) CPB draw, which may itself be zero.

Examples

set.seed(1)
y <- rzicpb(1000, lambda = 3, alpha = 0.5, pi = 0.3)
mean(y == 0)
#> [1] 0.31