当存在整数溢出时如何重新采样基尼系数?

时间:2015-02-26 20:21:51

标签: r integer-overflow resampling

我想计算基尼系数的置信区间,我试过 为此目的应用下面显示的boot()函数(as suggested here):

library(reldist) # just for the gini() function
library(boot) # for the boot() function
x <- c(1,2,2,3,4,99)
gini(x)
y <- boot(x, gini, 500)
quantile(y$t, probs=c(0.025, 0.975))

在这个例子中,它完全正常。但是当我尝试将它应用于我的实际数据时,我会收到一个警告(&#34;总和(权重):整数溢出 - 使用总和(as.numeric(。))&#34;)并且函数没有&#39;工作。而不是分位数我收到一个错误: &#34;在quantile.default中出错(y2 $ t,probs = c(0.025,0.975)):如果&#39; na.rm&#39;,则不允许缺少值和NaN。是假的&#34;

我不完全确定整数溢出警告和此错误是否相关,但我想是的。我读了几篇关于整数溢出的帖子(for example this one)并下载了Rmpfr包(as suggested here)。如果没有非常清楚地了解发生了什么,我试图输入x作为mpfr对象,但它没有用。此外,我的变量是数字,因此只需使用as.numeric转换(as suggested here)也不是解决方案。

是否有其他方法可以计算gini的误差范围,或解决此问题?

0 个答案:

没有答案