我正在对R中的一些昆虫计数数据进行零膨胀负二项式GLM。我的问题是如何让R将我的物种数据读作一个堆积的列以便保持零膨胀。如果我小计并将其导入R作为标题为Abundance的单行，我会松开零并且模型不起作用。我已经尝试过：
自己堆叠数据（有80列* 47行）所以手动堆叠后有3760行，你可以想象使用 pscl zeroinfl()命令时R的速度有多慢（我的电脑需要20分钟！它仍然有效）
下一个问题涉及空间相关性。某些采样器从同一介质中采样，以致违反独立性。我可以将介质作为模型的一个因素吗？
答案 0 :(得分：0)
使用PSCL需要2060分钟3760行？我的上帝，我有战斗30.000行:)这就是为什么我的pscl计算没有完成...
然而，我接着使用GLMM，包括嵌套随机效应（lme / gamm）和负二项分布，将theta设置为较低值，以便分布处理零膨胀。我认为这取决于零的程度。在我的情况下，它是44％，残差看起来相当不错。