Learn R Programming

Rdistance (version 4.1.1)

integrateOneStepNumeric: Numeric Integration of One-step Function

Description

Compute integral of the one-step distance function using numeric integration. This function is only called for oneStep functions that contain expansion factors.

Usage

integrateOneStepNumeric(
  object,
  newdata = NULL,
  w.lo = NULL,
  w.hi = NULL,
  Units = NULL,
  expansions = NULL,
  series = NULL,
  isPoints = NULL
)

Value

A vector of areas under distance functions. If object is a distance function and newdata is specified, the returned vector's length is nrow(newdata). If object is a distance function and newdata is NULL, returned vector's length is length(distances(object)). If object is a matrix, return's length is nrow(object).

Arguments

object

Either an Rdistance fitted distance function (an object that inherits from class "dfunc"; usually produced by a call to dfuncEstim), or a matrix of canonical distance function parameters (e.g., matrix(fit$par,1)). If a matrix, each row corresponds to a distance function and each column is a parameter. If object is a matrix, it should not have measurement units. Only quantities derived from function parameters (e.g., ESW) have units. Rdistance function parameters themselves never have units.

newdata

A data frame containing new values for the distance function covariates. If NULL and object is a fitted distance function, the observed covariates stored in object are used (behavior similar to predict.lm). Argument newdata is ignored if object is a matrix.

w.lo

Minimum sighting distance or left-truncation value if object is a matrix. Ignored if object is a fitted distance function. Must have physical measurement units.

w.hi

Maximum sighting distance or right-truncation value if object is a matrix. Ignored if object is a fitted distance function. Must have physical measurement units.

Units

Physical units of sighting distances if object is a matrix. Sighting distance units can differ from units of w.lo or w.hi. Ignored if object is a fitted distance function.

expansions

A scalar specifying the number of terms in series to compute. Depending on the series, this could be 0 through 5. The default of 0 equates to no expansion terms of any type. No expansion terms are allowed (i.e., expansions is forced to 0) if covariates are present in the detection function (i.e., right-hand side of formula includes something other than 1).

series

If expansions > 0, this string specifies the type of expansion to use. Valid values at present are 'simple', 'hermite', and 'cosine'.

isPoints

Boolean. TRUE if integration is for point surveys. FALSE for line-transect surveys. Line-transect surveys integrate under the distance function, g(x), while point surveys integrate under the distance function times distances, xg(x).

Details

The oneStep.like function has an extremely large discontinuity at Theta. Accurate numeric integration requires inserting Theta and Theta+ (a value just larger than Theta) into the series of points being evaluated. Because this creates un-equal intervals, the Trapazoid rule must be used. Rdistance's Simpson's rule routine (integrateNumeric) will not work for oneStep likelihoods that have expansions.

See Also

integrateNumeric; integrateOneStepLines; integrateOneStepPoints