lagplot {sgeostat} | R Documentation |
Create a spatially lagged scatter plot, e.g. plot z(s) versus z(s+h), where h is a lag in a pair object.
lagplot(point.obj, pair.obj, a1, a2, lag=1, std=FALSE, query.a, xlim, ylim)
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 |
NULL
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.
http://www.gis.iastate.edu/SGeoStat/homepage.html
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)