mvariocloudmap {GeoXp}R Documentation

Interactive multivariate variocloud and map

Description

The function mvariocloudmap() draws a scatterplot of pairwise Mahalanobis distances and spatial distances with a map. It is a multivariate version of the variocloud. The number of couples of sites plotted can be reduced by considering couples above a quantile regression curve.

Usage

mvariocloudmap(long, lat, object, listvar=NULL, listnomvar=NULL,quantiles=NULL, criteria=NULL,
carte = NULL, label = "",cex.lab=1,pch = 16,col="blue", xlab = "Pairwise spatial distances",
ylab="Pairwise Mahalanobis distances", lablong = "", lablat = "", axes=FALSE)

Arguments

long a vector x of size n
lat a vector y of size n
object a spatial weight matrix of class nb
listvar matrix of variables which permit to add graphics such as histogram, etc. and plot bubbles on map using the tlclk window
listnomvar names of variables from listvar
quantiles the value of alpha for representing alpha-quantile regression
criteria a vector of size n of boolean which permit to represent preselected sites with a cross, using the tcltk window
carte matrix with 2 columns for drawing spatial polygonal contours : x and y coordinates of the vertices of the polygon
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 points on the two graphics
xlab a title for the graphic x-axis
ylab a title for the graphic y-axis
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
axes TRUE for drawing axes on the map

Details

The pairwise Mahalanobis distances are calculated using the Minimum Covariance Determinant (MCD) estimator associated with 75% of observations (function covMcd in the robustbase package). Users have the possibility to select some couples of sites on the scatterplot that are also highlightened on the map. Selection of observations on the map is also possible and leads to the selection of all the couples which contain the selected observations on the scatterplot.

Value

A matrix of boolean of size n x n

Author(s)

Fizmoser P., Thomas-Agnan C., Ruiz-Gazen A., Laurent T.

References

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

See Also

misolationmap

Examples

data(meuse)
data(meuse.riv)

dat.meuse<-log(1+meuse[,3:7])

# matrice des plus proches voisins
coords <- coordinates(cbind(meuse$x,meuse$y))
col.knn <- knearneigh(coords, k=7)
nb.meuse <- knn2nb(col.knn)

obs<-mvariocloudmap(meuse$x,meuse$y,nb.meuse,dat.meuse,names(dat.meuse),quantiles=0.95,col='violet',pch=7,
carte=meuse.riv[-c(1:20,73:98,156:176),],label=as.character(1:155),cex.lab=0.7)

[Package GeoXp version 1.4 Index]