GeoScatterplot: h-scatterplot for space and space-time data.

View source: R/GeoScatterplot.R

GeoScatterplotR Documentation

h-scatterplot for space and space-time data.

Description

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

Usage

GeoScatterplot(data, coordx, coordy=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 \times d)-matrix (n iid spatial realisations) or a (d \times d)-matrix (a single spatial realisation on regular grid) or an (d \times d \times n)-array (n iid spatial realisations on regular grid) or a (t \times d)-matrix (a single spatial-temporal realisation) or an (t \times d \times n )-array (n iid spatial-temporal realisations) or or an (d \times d \times t \times n )-array (a single spatial-temporal realisation on regular grid) or an (d \times d \times t \times n )-array (n iid spatial-temporal realisations on regular grid). See GeoFit for details.

coordx

A numeric (d \times 2)-matrix (where d is the number of spatial sites) assigning 2-dimensions of spatial coordinates or a numeric d-dimensional vector assigning 1-dimension of spatial coordinates. Coordinates on a sphere for a fixed radius radius are passed in lon/lat format expressed in decimal degrees.

coordy

A numeric vector assigning 1-dimension of spatial coordinates; coordy is interpreted only if coordx is a numeric vector or grid=TRUE otherwise it will be ignored. Optional argument, the default is NULL then coordx is expected to be numeric a (d \times 2)-matrix.

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_t \times 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.

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.

Value

Produces a plot. No values are returned.

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))


GeoModels documentation built on June 22, 2024, 10:10 a.m.