ginimap {GeoXp} | R Documentation |
The function ginimap()
computes a Lorentz curve from var
and calculates the Gini Index associated to var
.
ginimap(long,lat,var,listvar=NULL, listnomvar=NULL,carte=NULL,criteria=NULL, label="", cex.lab=1, pch=5, col="green", xlab="", ylab="", axes=FALSE, lablong="", lablat="")
long |
a vector x of size n |
lat |
a vector y of size n |
var |
a vector of numeric values of size n |
listvar |
matrix of variables which permit to plot bubbles on map or add a graphic using the tk window |
listnomvar |
a list with names of variables listvar |
carte |
matrix with 2 columns for drawing spatial polygonal contours : x and y coordinates of the vertices of the polygon |
criteria |
a vector of size n of boolean which permit to represent preselected sites with a cross, using the tcltk window |
label |
a list of character of size n with name of site. Names are printed on map after a selection |
cex.lab |
character size of label |
pch |
16 by default, symbol for selected points |
col |
"blue" by default, color of bars on the histogram |
xlab |
a title for the graphic x-axis |
ylab |
a title for the graphic y-axis |
axes |
a boolean with TRUE for drawing axes on the map |
lablong |
name of the x-axis that will be printed on the map |
lablat |
name of the y-axis that will be printed on the map |
Users have the possibility to choose a threshold by mouse clicking on the Lorentz curve or by specifying it in the menu. The corresponding
pair (F,G) and the value of `var' are then printed on the graph and the sites
with a value of var
lower or equal to the threshold are then selected on the map.
The Gini Index is given in the tcltk window (see function gini.r for the formula used to calculate it).
Thomas-Agnan C., Aragon Y., Ruiz-Gazen A., Laurent T., Robidou L.
Aragon Yves, Perrin Olivier, Ruiz-Gazen Anne, Thomas-Agnan Christine (2009), Statistique et Econométrie pour données géoréférencées : modèles et études de cas
# données eire data(eire) eire.contours<-polylist2list(eire.polys.utm) obs<-ginimap(eire.coords.utm$V1,eire.coords.utm$V2,eire.df$INCOME, listvar=eire.df,listnomvar=names(eire.df),carte=eire.contours,col="orange", label=as.character(1:26),cex.lab=0.7,xlab="Income",ylab="Total Income")