Lorenz curves and Gini index
Wednesday, May 26th, 2010La Lorenz curve y el Gini Index se usan mucho como indicadores de la concentración de la riqueza (miden el grado de concentración de la renta en un país). En la literatura pesquera ha sido utilizada unas pocas veces para ilustrar la concentración de la abundancia. Nosotros queremos utilizarlos para analizar la concentración del esfuerzo de pesca, lo que en principio no se ha hecho antes.
#PARA HACER CURVAS DE LORENZ E INDICES GINI: #1. Hace falta tener cargado el package "ineq" library(ineq) Lorenz <- read.table(file = "C:\\R.gps\\Lorenz.txt",header = TRUE) names(Lorenz) summary(Lorenz) #Loop que hizo Inés para hacer las curvas y estimar los indices de las 31 patches de erizo #y guardar los primeros en un pdf (sin el loop tendría que hacer un código para cada uno de los 31 patches): giniIndex <- NULL pdf("GraficosLorentz.pdf", width=11.69, height=8.26) op <- par(mfrow=c(3,7)) for(i in 1:31){ Mtemp <- subset(Lorenz,Mancha==as.character(i), select=4)#Columna4=n.celda plot(Lc(Mtemp[,1]), main=paste("Mancha", i, sep=" ")) giniIndex[i] <- ineq(Mtemp[,1]) } par(op) dev.off() #Para guardar los índices: write(giniIndex,file="clipboard",ncolumns=1) #o write.csv(giniIndex,file="ginis.csv") #Resultó mejor la 2º pq da los datos en una unica columna con el nº #de patch asociado ##FIGURA DE LOS 3 PATCHES MÁS VISITADOS (7,15Y22): M7<-subset(Lorenz,Mancha=="7") M15<-subset(Lorenz,Mancha=="15") M22<-subset(Lorenz,Mancha=="22") Lc.7<-Lc(M7$n.celda) Lc.15<-Lc(M15$n.celda) Lc.22<-Lc(M22$n.celda) plot(Lc.7,ylab="",xlab="",main="") lines(Lc.15,lty="dashed") lines(Lc.22,lty="dotted")