使用零一膨胀 beta 回归 (zoib) 来拟合具有交互作用和 R 中单个随机项的概率响应数据

Question

我正在尝试拟合一个

zoib

模型，以占用概率作为响应变量，以对抗家庭收入与分类变量、饮食（例如食虫动物、食粮动物）之间的相互作用，并使用物种作为随机术语。这是为了探讨一个地区的物种占有概率与该地区家庭收入的关系是否因物种觅食行会/饮食类别而异。我使用

zoib

是因为我有百分比数据，并且我的数据是一和零膨胀的。

我对使用

zoib

或 beta 回归没有信心，并已联系

zoib

包的作者，为我提供拟合模型的正确方法。

这是我正在使用的数据结构：

'data.frame':   242 obs. of  6 variables:
 $ pentad     : chr  "2350_2940" "2350_2945" "2350_3015" "2355_2920"
 $ species    : chr  "Common Ostrich" "Common Ostrich" "Common Ostrich"
 $ occupancy  : num  1 0.1717 0.2012 0.0445 1 ...
 $ Diet       : chr  "HerbiGrani" "HerbiGrani" "HerbiGrani" "HerbiGrani"
 $ urban.cover: num  39.1 31.1 25.4 21.7 32.3 ...
 $ meanInc    : num  73470 33968 81768 124308 277171

这是我提供的代码：

m1 <- zoib(occupancy ~ meanInc*Diet|1|meanInc*Diet|1, random = 1,    
       EUID = OPdat$species,  zero.inflation = TRUE,  
       one.inflation = TRUE, data = OPdat,   n.iter=50, 
       n.thin=20, n.burn=50)

这是我得到的错误：

joint.1z01(y, n, q, xmu.1, p.xmu, xsum.1, p.xsum, x0.1, p.x0, 中的错误： 未找到对象“x1.1””

zoib

期刊对我来说很难理解。谁能帮助我理解这个错误？物体 x1.1 是什么？任何其他适合此模型的建议/包也将不胜感激。

Answer 1

来自

?zoib::zoib

的“详细信息”部分：

在公式右侧，最多可以包含 5 个部分 顺序如下： xb | xd| x0 | x1 | z，其中 xb 代表 链接函数中的固定效应协变量/因子 beta 分布的均值，xd 代表固定效应 两者之和的链接函数中的协变量/因子 beta 分布的形状参数，x0 表示 Pr(y=0) 链接函数中的固定效应协变量/因子， x1 表示链接中的固定效应协变量/因子 Pr(y=1|y>0) 的函数，z 表示随机效应 协变量/因素。

相关：https://stats.stackexchange.com/questions/494369/how-to-fit-a-zoib-regression-model-with-random-effects-from-zoib-package

我猜测您在 RHS 末尾还需要另一个

|1

，即

meanInc*Diet|1|meanInc*Diet|1|1

，以便您指定：

• 条件均值模型：

  meanInc*Diet

• 色散模型：

  1

  （微不足道/常数）
• 零通胀模型：

  meanInc*Diet

• 单膨胀模型：

  1

  （微不足道/恒定）
• 随机效应模型：

  1

  （变化截距）[这是我添加的部分]

我相信

brms

套装也适合ZOIB/Z01B型号。

Answer 2

您是否曾经让这个软件包为您工作？我正在运行相同的模型，并且我也在努力运行它。如果您有任何建议，我们将不胜感激。