lines.krige.bayes: Add a (Bayesian) Estimated Variogram to a Plot
Description
Adds a variogram model to a plot with an empirical variogram.
The estimate is a chosen summary (mean, mode or mean) of the
posterior distribution returned by the function krige.bayes.
Usage
## S3 method for class 'krige.bayes':
lines(x, max.dist,
summary.posterior = c("mode", "median", "mean"), ...)
Arguments
x
an object of the class krige.bayes, typically output
of the function krige.bayes.
max.dist
maximum distance for the x-axis.
summary.posterior
specify which summary of the posterior
distribution should be used as parameter estimation. Options are
"mean", "median" or
"mode". See DETAILS below.
A line with the estimated variogram plot is added to the plot in the
current graphics device. No values are returned.
Details
The functions krige.bayes produces samples of the
posterior distribution of the parameters $(\sigma^2, \phi,
\tau^{2}_{rel})$.
For the option mode, the mode of $(\phi,
\tau^{2}_{rel})$ is used with the conditional
(on thses two values) mode of $\sigma^2$.
For the options mean and median these summaries are
computed from the samples of the posterior.
References
Further information about geoR can be found at:
http://www.maths.lancs.ac.uk/~ribeiro/geoR.