> surved.households <- c(39, 179, 107, 56, 32, 22, 11, 15, 38)
> percentage.household<-surved.households*100/sum(surved.households)
> m2<-cbind(surved.households,percentage.household)
> dimnames(m2)<- list(row.names,c("surved.households","percentage.household"))
> income <- c(56126.36, 290124.32, 281048.35, 225586.54, 162442.64, 115243.11, 64389.11, 136743.99, 784834.07)
> average.per.capita.income<-income/surved.households
> percent.share.aggregate.income<-income*100/sum(income)
> cumulative.percent.share.aggregate.income <-cumsum(percent.share.aggregate.income)
> a<-cumulative.percent.share.aggregate.income
> successive.cumulative.percent.share.aggregate.income <- c(a[1], a[1]+a[2], a[2]+a[3], a[3]+a[4], a[4]+a[5], a[5]+a[6], a[6]+a[7], a[7]+a[8], a[8]+a[9])
> m3 <-cbind(m2,income,average.per.capita.income,percent.share.aggregate.income,cumulative.percent.share.aggregate.income,successive.cumulative.percent.share.aggregate.income)
> lorenz <- (1 - (sum(percentage.household*successive.cumulative.percent.share.aggregate.income))/10000)
> lorenz <- (1 - sum( m3[,2] * m3[,7] )/10000) # Using the matrix column




Note: These commands of R Console was documented in 9:53 AM 12/13/2003 by Mohammad Ehsanul Karim, Institute of Statistical Research and Training, Dhaka - 1000, Bangladesh.