GeoScatterplot function

h-scatterplot for space and space-time data.

The function produces h-scatterplots for the spatial, spatio-temporal and bivariate setting. UTF-8

GeoScatterplot(data, coordx, coordy=NULL, coordz=NULL, coordt=NULL, coordx_dyn=NULL,
           distance='Eucl', grid=FALSE, maxdist=NULL,neighb=NULL,
           times=NULL, numbins=4, radius=6371, bivariate=FALSE,...)

Arguments

• data: A $d$-dimensional vector (a single spatial realisation) or a ($n x d$)-matrix ($n$ iid spatial realisations) or a ($d x d$)-matrix (a single spatial realisation on regular grid) or an ($d x d x n$)-array ($n$ iid spatial realisations on regular grid) or a ($t x d$)-matrix (a single spatial-temporal realisation) or an ($t x d x n$)-array ($n$ iid spatial-temporal realisations) or or an ($d x d x t$)-array (a single spatial-temporal realisation on regular grid) or an ($d x d x t x n$)-array ($n$ iid spatial-temporal realisations on regular grid). See GeoFit for details.

• coordx: A numeric ($d x 2$)-matrix or ($d x 3$)-matrix Coordinates on a sphere for a fixed radius radius

  are passed in lon/lat format expressed in decimal degrees.

• coordy: A numeric vector giving 1-dimension of spatial coordinates; Optional argument, the default is NULL.

• coordz: A numeric vector giving 1-dimension of spatial coordinates; Optional argument, the default is NULL.

• coordt: A numeric vector assigning 1-dimension of temporal coordinates. Optional argument, the default is NULL then a spatial random field is expected.

• coordx_dyn: A list of $m$ numeric ($d x 2$)-matrices containing dynamical (in time) spatial coordinates. Optional argument, the default is NULL

• distance: String; the name of the spatial distance. The default is Eucl, the euclidean distance. See the Section Details of GeoFit.

• grid: Logical; if FALSE (the default) the data are interpreted as spatial or spatial-temporal realisations on a set of non-equispaced spatial sites.

• maxdist: A numeric value denoting the spatial maximum distance, see the Section Details .

• neighb: Numeric; an optional positive integer indicating the order of neighborhood. See the Section Details for more information.

• times: A numeric vector denoting the temporal instants involved Details .

• numbins: A numeric value denoting the numbers of bins, see the Section Details .

• radius: Numeric; a value indicating the radius of the sphere when using the great circle distance. Default value is the radius of the earth in Km (i.e. 6371)

• bivariate: Logical; if FALSE (the default) the data are interpreted as univariate spatial or spatial-temporal realisations. Otherwise they are intrepreted as a a realization from a bivariate field.

• ...: Optional parameters passed to the plot function.

Returns

Produces a plot. No values are returned.

Details

h-scatterplot is the plot of the pair values that are neighborhood of a certain order or with distances belonging to a certain interval. In the first case a (vector of) neighborhood must be specified. In the second case a maximum distance (maxdist) and a number of lag-bins (numbins) must be specified. The method based on neighborhoods is recommended in particular for large datasets.

Author(s)

Moreno Bevilacqua, moreno.bevilacqua89@gmail.com ,https://sites.google.com/view/moreno-bevilacqua/home, Víctor Morales Oñate, victor.morales@uv.cl , https://sites.google.com/site/moralesonatevictor/, Christian", Caamaño-Carrillo, chcaaman@ubiobio.cl ,https://www.researchgate.net/profile/Christian-Caamano

Examples

library(GeoModels)
set.seed(514)

NN = 600
coords = cbind(runif(NN),runif(NN))

param = list(mean=0,sill=1,nugget=0,power2=4,scale=0.4,smooth=0)

corrmodel = "GenWend"; model = "Gaussian"

data = GeoSim(coordx = coords,corrmodel = corrmodel,
                  model = model,param = param)$data

# h-scatterplots for given a vector of neighborhoods
GeoScatterplot(data,coords,neighb=c(2,4))