我试图拟合一个zoib
模型,其中占用概率作为响应变量,以对抗家庭收入与分类变量,饮食(例如,食虫动物,花岗岩动物)之间的相互作用,并使用物种作为随机项。本研究旨在探讨物种在一个地区的占有概率与该地区家庭收入的关系是否因物种觅食行会/饮食类别的不同而不同。我使用zoib
,因为我有百分比数据,因为我的数据是1和0膨胀。
我对使用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)
这是我得到的错误:
关节错误。1z01(y, n, q, x)1, p, xmu, xsum。1、p.xsum, x0.1, p.x0,对象'x1.1'未找到">
zoib
rjournal对我来说很难理解。有人能帮我理解这个错误吗?对象x1.1是什么?如有任何其他适合此模型的建议/包,也将不胜感激。
From "Details"?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
我猜你需要另一个|1
在你的RHS结束,即。meanInc*Diet|1|meanInc*Diet|1|1
,以便您指定:
- 条件平均模型:
meanInc*Diet
- 色散模型:
1
(trivial/constant) - 零膨胀模型:
meanInc*Diet
- 单膨胀模型:
1
(trivial/constant) - 随机效应模型:
1
(可变截距)[这是我添加的部分]
我相信brms
包也将适合ZOIB/Z01B型号。