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 with TRUE on specific sites (these for non interactive selection) |
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, choice of representation of 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 (2008), ``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")