ジニ係数と・・・

 

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なのか