dbledensitymap {GeoXp} | R Documentation |
The function dbledensitymap
plots two kernel density estimates from var1
and var2
,
computed with `bkde.R', and a map with sites of coordinates (long
,lat
). Each site is associated to a value of var1
and var2
and
the two windows are interactive.
dbledensitymap(long, lat, var1, var2, kernel='triweight',listvar=NULL, listnomvar=NULL, carte=NULL, criteria=NULL, label="", cex.lab=1, pch=16, col=c("grey","blue"), xlab=c("",""), ylab=c("density","density"),axes=FALSE, lablong="", lablat="")
long |
a vector x of size n |
lat |
a vector y of size n |
var1 |
a vector of numeric values of size n |
var2 |
a vector of numeric values of size n |
kernel |
Smoothing kernel (see help(bkde) for list of options) |
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 |
vector of character of size n with names of sites |
cex.lab |
character size of label on map |
pch |
16 by default, symbol for selected points |
col |
c("grey","blue") by default, color of the two density curves |
xlab |
a list of title for the two x-axis graphics |
ylab |
a list of title for the two y-axis graphics |
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 |
The user can choose an interval on the density curve by mouse clicking on the lower and upper boundaries of the interval or by giving directly these values. The selected sites are then represented on the map in red. A selection by `points' or `polygon' on the map results in the drawing of the density of the corresponding sub-distribution on the density plot. Finally, the user can modify the bandwith parameter with a cursor in the tk window (parameter alpha). alpha is the smoothing parameter for the kernel smooth : it represents the mean percentage of sample points involved in the local averaging (example : alpha=20 means that on average, n x 0.2 points are in any interval of length 2h where h is the usual bandwidth).
A vector of boolean of size n. TRUE if the site was in the last selection.
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
Venables, W. N. and Ripley, B. D. (2002) Modern Applied Statistics with S. New York: Springer.
Wand M.P. et Jones M.C. (1995), Kernel Smoothing, Chapman & Hall.
histomap
, histobarmap
, scattermap
, densitymap
# data auckland data(auckland) contours.auckland<-polylist2list(auckpolys) dbledensitymap(auckland$Easting,auckland$Northing,var1=auckland$Deaths.1977.85, var2=auckland$Under.5.1981,carte=contours.auckland, listvar=cbind(auckland,as.numeric(auckland$Deaths.1977.85> mean(auckland$Deaths.1977.85))), listnomvar=c(names(auckland),"sup to mean"), xlab=c("Deaths.1977.85","Under.5.1981"), criteria=(auckland$Deaths.1977.85>mean(auckland$Deaths.1977.85))) # data columbus x <- readShapePoly(system.file("shapes/columbus.shp", package="maptools")[1]) colombus.contours<-spdf2list(x)$poly coords <- coordinates(x) colombus<-x@data obs<-dbledensitymap(coords[,1],coords[,2],colombus$HOVAL, colombus$CRIME,label=colombus$NEIGNO,carte=colombus.contours, cex.lab=0.7,listvar=colombus,listnomvar=names(colombus), col=c("purple","pink"),xlab=c("hoval","crime"),pch=14) # data eire data(eire) eire.contours<-polylist2list(eire.polys.utm) dbledensitymap(eire.coords.utm$V1,eire.coords.utm$V2,eire.df$A,eire.df$towns, xlab=c("Taux d'individus au groupe sanguin A", "Surface urbaine"),carte=eire.contours,kernel="normal", listvar=eire.df,listnomvar=names(eire.df))