R中的零膨胀分布

Question

除了“iZid”包外，R中有哪些不同的包包含了模拟零膨胀分布的函数，涉及到诸如泊松、负二项、康-泊松、泊松逆高斯、泊松-林德利等常用的离散模型。

Answer 1

如果您有一个从非零充气分布产生随机偏差的现有函数，您可以编写一个包装器(或装饰器)来创建一个零充气偏差模拟器。我在这里所做的唯一假设是，原始函数的第一个参数称为n，并指定要选择的随机偏差的数目。

例如，如果我们想要扩展rbinom以返回零膨胀的二项式偏差.

ziversion <- function(rfun) {
   f <- function(n, ..., zi) {
      x <- rfun(n, ...)
      x <- ifelse(runif(n) < zi, 0, x)
      return(x)
   }
   return(f)
}
rzibinom <- ziversion(rbinom)
set.seed(101)
rzibinom(10, size = 10, prob = 0.2, zi = 0.5)
##  [1] 1 0 3 2 0 1 2 0 0 0

zi是零通胀概率.只要付出一点努力，代码就能变得更有效率.

Answer 2

看一看分配的CRAN任务视图。这是对帮助您使用发行版的R包的精心设计。您可以在页面中搜索“膨胀”，以快速找到相关部分。