��Ū

�������Ķ��������������˷�����׻����롣

����ˡ

Gini.index(dat, main="Lorenz curve", xlab="", ylab="")

����

dat     �ǡ����٥��ȥ�
main    �ޤΥ����ȥ롡��ά���ϲ���񤫤ʤ�
xlab    x ����̾������ά���ϲ���񤫤ʤ�
ylab    y ����̾������ά���ˤϲ���񤫤ʤ�

������

���󥹥ȡ���ϡ��ʲ��� 1 �Ԥ򥳥ԡ�����R ���󥽡���˥ڡ����Ȥ���
source("http://aoki2.si.gunma-u.ac.jp/R/src/Gini_index.R", encoding="euc-jp")

# �������Ķ��������������˷�����׻�����
Gini.index <- function(      y,                      # �ǡ����٥��ȥ�
                        main="",                # �ޤΥ����ȥ�ʾ�ά���ϲ���񤫤ʤ���
                        xlab="",                # x ����̾���ʾ�ά���ϲ���񤫤ʤ���
                        ylab="")                # y ����̾���ʾ�ά���ϲ���񤫤ʤ���
{
        stopifnot(y >= 0)                       # ����ǡ����Ǥʤ���Фʤ�ʤ�
        n <- length(y)                               # �ǡ����θĿ�
        y <- sort(y)                         # ����������¤٤�
        y <- cumsum(y)                               # �����ٿ���Ȥ�
        y <- c(0, y/y[n])                    # ���������ٿ�����Ƭ�� 0 ��ä����
        x <- seq(0, 1, length=n+1)           # 0 �� 1 �����ֳ֤˶��ڤä��٥��ȥ����
        old <- par(xaxs="i", yaxs="i")
        plot(x, y, type="l", col="blue",        # ������֤ȥ������Ķ���
                main=main, xlab=xlab, ylab=ylab)
        abline(0, 1)                            # �г����ʸ������̤롤���� 1 ��ľ���ˤ�����
        par(old)
        return(2*sum(x-y)/n)                    # ���˷���
}


�����㡦���Ϸ����

> x <- c(0, 3, 3, 5, 5, 5, 5, 7, 7, 10)
> Gini.index(x)
[1] 0.276
> x <- c(0, 1, 2, 3, 5, 5, 7, 8, 9, 10)
> Gini.index(x)
[1] 0.376
> x <- c(3, 4, 4, 5, 5, 5, 5, 6, 6, 7)
> Gini.index(x)
[1] 0.12