`method'

Values for method that make gstat deviate from the default action are:

gs

override the default kriging method to get Gaussian simulations (section 2.6)

is

override the default kriging method to get indicator simulations (section 2.6)

semivariogram

sample semivariance on output (see also fit)

covariogram

sample covariogram on output

trend

spatial trend estimation ( $x_0 \hat \beta$), using generalised least squares if variograms are specified, or else using ordinary or weighted least squares (see Appendix A.2)

map

report the value of $m$ mask maps at the prediction location of data() as the output values, provided that $n$ (dummy) input variable are defined, with $n \ge m/2$ (example 9)

distance

report distance to nearest observation as the predicted value, and the distance to the most distant observation (in the neighbourhood selection, if defined) as prediction variance

nr

report number of observations in neighbourhood as predicted value, and, if omax was set, the number of non-empty octants (or quadrants) as prediction variance

div

diversity (the number of distinct values) in a local neighbourhood is reported as the predicted value, the range of the values (largest value $-$ smallest value) is reported as prediction variance

med

local median or quantile estimation, see quantile

point-in-polygon

Provided a set of (closed) polygons is given with the edges command, the (first) polygon in which the prediction location is located is given in the output. Predicted value carries the file number of the polygon, prediction variance the polygon number in the file. See also Appendix B