感谢分享。ineq中的Gini函数也可以算,不过在算相对基尼系数时需要调整。我这个无需调整,算是一个小改进吧。
另外还推荐一个更针对区域经济研究的REAT,不过他那个gini.conc参数设定实在太过复杂。有兴趣的去看看吧,还有很多有用的函数如算区域专业化程度的gini.spec等。
至于是自编还是用现成的函数,这个仁者见仁、智者见智吧。我这个函数的优点是参数设定简明,且有明确文献支持。他人函数你不知道他是基于什么公式来写的。
下面是各个函数的比较分析结果。
> x <- 1:5
> y <- seq(50,10,-10)
> z <- (x/sum(x)) / (y/(sum(y)))
>
> geo.gini(x)
[1] 0.2666667
> ineq::Gini(x)
[1] 0.2666667
> REAT::gini.conc(x, e_j = rep(100, 5))
[1] 0.2666667
>
> geo.gini(z)
[1] 0.5103448
> ineq::Gini(z)
[1] 0.5103448
> REAT::gini.conc(x, e_j = y)
[1] 0.5103448