sgeostat (version 1.0-27)

lagplot: Lag Scatter Plot

Description

Create a spatially lagged scatter plot, e.g. plot z(s) versus z(s+h), where h is a lag in a pair object.

Usage

lagplot(point.obj, pair.obj, a1, a2, lag=1, std=FALSE, query.a, xlim, ylim)

Arguments

point.obj

a point object generated by point()

pair.obj

a pair object generated by pair()

a1

a variable to plot

a2

an optional variable name, if entered the plot will be created between a1 and a2

lag

the lag to plot

std

a logical variable indicating whether the data should be standardized to their means and standard deviations before plotting

query.a

an optional variable name, if entered, the value of the variable will be displayed on the graphics device for points identified by the user.

xlim

a vector of length 2 indicating the x limits of the graphics page

ylim

a vector of length 2 indicating the y limits of the graphics page

Value

References

http://www.gis.iastate.edu/SGeoStat/homepage.html

See Also

point, pair

Examples

opar <- par(ask = interactive() && .Device == "X11")
lagplot(maas.point,maas.pair,'zinc')
# with identifying pairs:
lagplot(maas.point,maas.pair,'zinc',lag=2,query.a='zinc')
par(opar)