fda.usc (version 2.0.1)

int.simpson: Simpson integration

Description

Computes the integral of fdataobj$data with respect to fdataobj$argvals using simpson or trapezoid rule integration.

Usage

int.simpson(fdataobj, method = NULL)
int.simpson2(x, y, equi = TRUE, method = NULL)

Arguments

fdataobj

fdata objtect.

method

Method for numerical integration, see details.

x

Sorted vector of x-axis values: argvals.

y

Vector of y-axis values.

equi

=TRUE, the observed points on each curve are equally spaced (by default).

Details

Posible values for method are:

• "TRAPZ": Trapezoid rule integration.

• "CSR": Composite Simpson's rule integration.

• "ESR": Extended Simpson's rule integration.

If method=NULL (default), the value of par.fda.usc$int.method is used.

See Also

See also integrate.

Examples

# NOT RUN {
x<-seq(0,2*pi,length=1001)
fx<-fdata(sin(x)/sqrt(pi),x)
fx0<-fdata(rep(0,length(x)),x)
int.simpson(fx0)
int.simpson(fx)
# }