lagplot {sgeostat}

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

NULL

Note

When query.a is entered, the user will be prompted to identify points on the display device. Because each point in the plot represents a pair of locations, the user must identify each point twice, once for the "from" point and once for the "to" point. Querying is ended by pressing the middle mouse button on the mouse while the cursor is in the display window.

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)