| lines.variomodel.krige.bayes {geoR} | R Documentation |
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.
## S3 method for class 'krige.bayes':
lines.variomodel(x, summary.posterior, max.dist, uvec,
posterior = c("variogram", "parameters"), ...)
x |
an object of the class krige.bayes, typically an output
of the function krige.bayes. |
summary.posterior |
specify which summary of the posterior
distribution should be used as the parameter estimate.
Options are "mean", "median" or
"mode". See DETAILS below. |
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 = "variogram".
Defaults to seq(0, max.dist, length = 51).
|
posterior |
indicates whether the the variogram line is based on the posterior of the variogram function (default) or the posterior of the model parameters. |
... |
arguments passed to the functions lines or curve. |
The function krige.bayes returns samples from the
posterior distribution of the parameters (sigma^2, phi, tausq.rel).
This function allows for two basic options to draw a line with a summary of the variogram function.
uvec. Then a function provided by the user in the
argument summary.posterior is used to compute a summary of
the values obtained at each support point.}
item[2.]{"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 argument summary.posterior.
The option mode, uses the mode of (phi, tau^2.rel) and the mode of
of sigma^2 conditional on the modes of the former parameters.
For the options mean and median these summaries are
computed from the samples of the posterior.
}
A line with the estimated variogram plot is added to the plot in the current graphics device. No values are returned.
Paulo J. Ribeiro Jr. paulojus@leg.ufpr.br,
Peter J. Diggle p.diggle@lancaster.ac.uk.
Further information on the package geoR can be found at:
http://www.leg.ufpr.br/geoR.
lines.variomodel, krige.bayes and lines.
#See examples in the documentation of the function krige.bayes().