The wrapper function polyCub can be used to call specific cubature
methods via its method argument. It calls the polyCub.SV
function by default, which implements general-purpose product Gauss cubature.
The desired cubature function should usually be called directly.
polyCub(polyregion, f, method = c("SV", "midpoint", "iso", "exact.Gauss"),
..., plot = FALSE)The approximated integral of f over polyregion.
a polygonal domain.
The following classes are supported:
"owin" from package spatstat.geom,
"gpc.poly" from gpclib,
"SpatialPolygons",
"Polygons",
and "Polygon" from package sp, as well as
"(MULTI)POLYGON" from package sf.
(For these classes, polyCub knows how to get an xylist.)
a two-dimensional real-valued function to be integrated over
polyregion. As its first argument it must take a coordinate matrix,
i.e., a numeric matrix with two columns, and it must return a numeric vector
of length the number of coordinates.
For the "exact.Gauss" method,
f is ignored since it is specific to the bivariate normal density.
choose one of the implemented cubature methods (partial
argument matching is applied), see help("polyCub-package")
for an overview. Defaults to using product Gauss cubature
implemented in polyCub.SV.
arguments of f or of the specific method.
logical indicating if an illustrative plot of the numerical integration should be produced.
Details and examples in the vignette("polyCub")
and on the method-specific help pages.
Other polyCub-methods:
polyCub.SV(),
polyCub.exact.Gauss(),
polyCub.iso(),
polyCub.midpoint()