lines.variomodel.krige.bayes {geoR} | R Documentation |
Adds a Bayesian Estimate of the Variogram to a Plot
Description
Adds a Bayesian estimate of the variogram model to a plot typically 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.variomodel(x, summary.posterior, max.dist, uvec,
posterior = c("variogram", "parameters"), ...)
Arguments
x |
an object of the class |
summary.posterior |
specify which summary of the posterior
distribution should be used as the parameter estimate.
Options are |
max.dist |
numerical, the maximum distance for the x-axis. |
uvec |
a numerical vector with support points to compute the
variogram values. Only used if |
posterior |
indicates whether the the variogram line is based on the posterior of the variogram function (default) or the posterior of the model parameters. |
... |
Details
The function krige.bayes
returns samples from the
posterior distribution of the parameters .
This function allows for two basic options to draw a line with a summary of the variogram function.
- "variogram":
for each sample of the parameters the variogram function is computed at the support points defined in the argument
uvec
. Then a function provided by the user in the argumentsummary.posterior
is used to compute a summary of the values obtained at each support point.- "parameters":
in this case summaries of the posterior distribution of the model parameters as “plugged-in” in the variogram function. One of the options
"mode"
(default) ,"median"
or"mean"
can be provided in the argumentsummary.posterior
. The optionmode
, uses the mode ofand the mode of of
conditional on the modes of the former parameters. For the options
mean
andmedian
these summaries are computed from the samples of the posterior.
Value
A line with the estimated variogram plot is added to the plot in the current graphics device. No values are returned.
Author(s)
Paulo J. Ribeiro Jr. paulojus@leg.ufpr.br,
Peter J. Diggle p.diggle@lancaster.ac.uk.
References
Further information on the package geoR can be found at:
http://www.leg.ufpr.br/geoR/.
See Also
lines.variomodel
, krige.bayes
and lines
.
Examples
#See examples in the documentation of the function krige.bayes().