我正在使用Adehabitat
软件包中的compana
功能来确定是否从可用栖息地中选择了使用的栖息地。
我在可用矩阵中使用的栖息地类别之一具有非常低的值(即0.000015000
(或零值。因此,当我用更大的值(0.001
(替换R代码中的任何零值时,会影响我的结果。当我用其他值(即0.0001
和0.00000001
(替换较大的值时,为了确认结果的一致性,我的结果不一致。
我应该使用适当的小价值数字吗?
我尝试使用多个不同的小值来替换矩阵中的零值。小于0.001
的数量较小(与9种不同的栖息地类别一起工作(,而是不断增加较小的值会改变其他栖息地类别的值。
bars<-compana(used, avail, test = c("randomisation"),
rnv = 0.01, nrep = 1000, alpha = 0.1)
从 0.01
更改 rnv
值 -> 0.0001
改变结果,并继续使用较小的数字
我希望与大多数或所有rnv少量值看到一致的结果。
我建议阅读Aebischer等人的论文。al。(C 构成分析的栖息地使用中的栖息地数据,《生态学》,第74卷,第5期(1993年7月(,第1313-1325页(在。读和解释运行栖息地的组成分析的许多方面很容易。
在第1320-1321页上,作者说以下内容:
[...] 0%表示用途太低而无法记录,所以应该是被一个明显小于最小非零值的值所取代:少级数可能适合大多数情况。
因此,我会检查数据集中最小的非零值,并选择一个级级的值以交换零。