ジニ係数と・・・
R-bloggersより
Gini index and Lorenz curve with R
http://www.r-bloggers.com/gini-index-and-lorenz-curve-with-r/
lthough I did not explain it during my lectures, calculating a Gini index or displaying the Lorenz curve can be done very easily with R.
だそうです。
とはいえ、王様の場合
xking <- c(0,0,0,0,50)
ineq(xking)
plot(Lc(xking))
これは1になるはずだが、帰ってきた答えは0.8・・・・
データ数が多いと近似的に1になるが
コードをみると
> Gini
function (x, corr = FALSE)
{
x <- as.numeric(x)
n <- length(x)
x <- sort(x)
G <- sum(x * 1:n)
G <- 2 * G/(n * sum(x))
G <- G - 1 - (1/n)
if (corr)
G * n/(n - 1)
else G
}
・・・・
大青木先生のとこでは
http://aoki2.si.gunma-u.ac.jp/R/Gini-index.html
同じ・・・
やはり絶対王制下の王様独占ケースでも0.8なのか